Search

Abstract Background The Epithelial-Mesenchymal Transition (EMT) endows epithelial-looking cells with enhanced migratory ability during embryonic development and tissue repair. EMT can also be co-opted by cancer cells to acquire metastatic potential and drug-resistance. Recent research has argued that epithelial (E) cells can undergo either a partial EMT to attain a hybrid epithelial/mesenchymal (E/M) phenotype that typically displays collective migration, or a complete EMT to adopt a mesenchymal (M) phenotype that shows individual migration. The core EMT regulatory network - miR-34/SNAIL/miR-200/ZEB1 - has been identified by various studies, but how this network regulates the transitions among the E, E/M, and M phenotypes remains controversial. Two major mathematical models – ternary chimera switch (TCS) and cascading bistable switches (CBS) - that both focus on the miR-34/SNAIL/miR-200/ZEB1 network, have been proposed to elucidate the EMT dynamics, but a detailed analysis of how well either or both of these two models can capture recent experimental observations about EMT dynamics remains to be done. Results Here, via an integrated experimental and theoretical approach, we first show that both these two models can be used to understand the two-step transition of EMT - E→E/M→M, the different responses of SNAIL and ZEB1 to exogenous TGF-β and the irreversibility of complete EMT. Next, we present new experimental results that tend to discriminate between these two models. We show that ZEB1 is present at intermediate levels in the hybrid E/M H1975 cells, and that in HMLE cells, overexpression of SNAIL is not sufficient to initiate EMT in the absence of ZEB1 and FOXC2. Conclusions These experimental results argue in favor of the TCS model proposing that miR-200/ZEB1 behaves as a three-way decision-making switch enabling transitions among the E, hybrid E/M and M phenotypes.

One of the most important roles of cells is performing their cellular tasks properly for survival. Cells usually achieve robust functionality, for example, cell-fate decision-making and signal transduction, through multiple layers of regulation involving many genes. Despite the combinatorial complexity of gene regulation, its quantitative behavior has been typically studied on the basis of experimentally verified core gene regulatory circuitry, composed of a small set of important elements. It is still unclear how such a core circuit operates in the presence of many other regulatory molecules and in a crowded and noisy cellular environment. Here we report a new computational method, named random circuit perturbation (RACIPE), for interrogating the robust dynamical behavior of a gene regulatory circuit even without accurate measurements of circuit kinetic parameters. RACIPE generates an ensemble of random kinetic models corresponding to a fixed circuit topology, and utilizes statistical tools to identify generic properties of the circuit. By applying RACIPE to simple toggle-switch-like motifs, we observed that the stable states of all models converge to experimentally observed gene state clusters even when the parameters are strongly perturbed. RACIPE was further applied to a proposed 22-gene network of the Epithelial-to-Mesenchymal Transition (EMT), from which we identified four experimentally observed gene states, including the states that are associated with two different types of hybrid Epithelial/Mesenchymal phenotypes. Our results suggest that dynamics of a gene circuit is mainly determined by its topology, not by detailed circuit parameters. Our work provides a theoretical foundation for circuit-based systems biology modeling. We anticipate RACIPE to be a powerful tool to predict and decode circuit design principles in an unbiased manner, and to quantitatively evaluate the robustness and heterogeneity of gene expression.

Neural networks, despite their highly interconnected nature, exhibit distinctly localized and gated activation. Modularity, a distinctive feature of neural networks, has been recently proposed as an important parameter determining the manner by which networks support activity propagation. Here we use an engineered biological model, consisting of engineered rat cortical neurons, to study the role of modular topology in gating the activity between cell populations. We show that pairs of connected modules support conditional propagation (transmitting stronger bursts with higher probability), long delays and propagation asymmetry. Moreover, large modular networks manifest diverse patterns of both local and global activation. Blocking inhibition decreased activity diversity and replaced it with highly consistent transmission patterns. By independently controlling modularity and disinhibition, experimentally and in a model, we pose that modular topology is an important parameter affecting activation localization and is instrumental for population-level gating by disinhibition.

Financial markets are partially composed of sectors dominated by external driving forces, such as commodity prices, infrastructure and other indices. We characterize the statistical properties of such sectors and present a novel model for the coupling of the stock prices and their dominating driving forces, inspired by mean reverting stochastic processes. Using the model we were able to explain the market sectors' long term behavior and estimate the coupling strength between stocks in financial markets and the sector specific driving forces. Notably, the analysis was successfully applied to the shipping market, in which the Baltic dry index (BDI), an assessment of the price of transporting the major raw materials by sea, influences the shipping financial market. We also present the analysis of other sectors-the gold mining market and the food production market, for which the model was also successfully applied. The model can serve as a general tool for characterizing the coupling between external forces and affected financial variables and therefore for estimating the risk in sectors and their vulnerability to external stress.

BACKGROUND:Regulation of hepatic glucose production has been a target for antidiabetic drug development, due to its major contribution to glucose homeostasis. Previous pre-clinical study demonstrated that peripheral electrical stimulation (PES) may stimulate glucose utilization and improve hepatic insulin sensitivity. The aim of the present study was to evaluate safety, tolerability, and the glucose-lowering effect of this approach in patients with type 2 diabetes (T2DM). METHODS:Twelve patients with T2DM were recruited for an open label, interventional, randomized trial. Eleven patients underwent, in a crossover design, an active, and a no-intervention control periods, separated with a two-week washout phase. During the active period, the patients received a daily lower extremity PES treatment (1.33Hz/16Hz burst mode), for 14 days. Study endpoints included changes in glucose levels, number of hypoglycemic episodes, and other potential side effects. Endpoints were analyzed based on continuous glucose meter readings, and laboratory evaluation. RESULTS:We found that during the active period, the most significant effect was on nocturnal glucose control (P < 0.0004), as well as on pre-meal mean glucose levels (P < 0.02). The mean daily glucose levels were also decreased although it did not reach clinical significance (P = 0.07). A reduction in serum cortisol (P < 0.01) but not in insulin was also detected after 2 weeks of treatment. No adverse events were recorded. CONCLUSIONS:These results indicate that repeated PES treatment, even for a very short duration, can improve blood glucose control, possibly by suppressing hepatic glucose production. This effect may be mediated via hypothalamic-pituitary-adrenal axis modulation. TRIAL REGISTRATION:ClinicalTrials.gov NCT02727790.

The rapid increase of wealth inequality in the past few decades is one of the most disturbing social and economic issues of our time. Studying its origin and underlying mechanisms is essential for policy aiming to control and even reverse this trend. In that context, controlling the distribution of income, using income tax or other macroeconomic policy instruments, is generally perceived as effective for regulating the wealth distribution. We provide a theoretical tool, based on the realistic modeling of wealth inequality dynamics, to describe the effects of personal savings and income distribution on wealth inequality. Our theoretical approach incorporates coupled equations, solved using iterated maps to model the dynamics of wealth and income inequality. Notably, using the appropriate historical parameter values we were able to capture the historical dynamics of wealth inequality in the United States during the course of the 20th century. It is found that the effect of personal savings on wealth inequality is substantial, and its major decrease in the past 30 years can be associated with the current wealth inequality surge. In addition, the effect of increasing income tax, though naturally contributing to lowering income inequality, might contribute to a mild increase in wealth inequality and vice versa. Plausible changes in income tax are found to have an insignificant effect on wealth inequality, in practice. In addition, controlling the income inequality, by progressive taxation, for example, is found to have a very small effect on wealth inequality in the short run. The results imply, therefore, that controlling income inequality is an impractical tool for regulating wealth inequality.

A new paradigm has recently emerged in brain science whereby communications between glial cells and neuron-glia interactions should be considered together with neurons and their networks to understand higher brain functions. In particular, astrocytes, the main type of glial cells in the cortex, have been shown to communicate with neurons and with each other. They are thought to form a gap-junction-coupled syncytium supporting cell-cell communication via propagating Ca(2+) waves. An identified mode of propagation is based on cytoplasm-to-cytoplasm transport of inositol trisphosphate (IP(3)) through gap junctions that locally trigger Ca(2+) pulses via IP(3)-dependent Ca(2+)-induced Ca(2+) release. It is, however, currently unknown whether this intracellular route is able to support the propagation of long-distance regenerative Ca(2+) waves or is restricted to short-distance signaling. Furthermore, the influence of the intracellular signaling dynamics on intercellular propagation remains to be understood. In this work, we propose a model of the gap-junctional route for intercellular Ca(2+) wave propagation in astrocytes. Our model yields two major predictions. First, we show that long-distance regenerative signaling requires nonlinear coupling in the gap junctions. Second, we show that even with nonlinear gap junctions, long-distance regenerative signaling is favored when the internal Ca(2+) dynamics implements frequency modulation-encoding oscillations with pulsating dynamics, while amplitude modulation-encoding dynamics tends to restrict the propagation range. As a result, spatially heterogeneous molecular properties and/or weak couplings are shown to give rise to rich spatiotemporal dynamics that support complex propagation behaviors. These results shed new light on the mechanisms implicated in the propagation of Ca(2+) waves across astrocytes and the precise conditions under which glial cells may participate in information processing in the brain.

Acetylcholine (ACh) is a regulator of neural excitability and one of the neurochemical substrates of sleep. Amongst the cellular effects induced by cholinergic modulation are a reduction in spike-frequency adaptation (SFA) and a shift in the phase response curve (PRC). We demonstrate in a biophysical model how changes in neural excitability and network structure interact to create three distinct functional regimes: localized asynchronous, traveling asynchronous, and traveling synchronous. Our results qualitatively match those observed experimentally. Cortical activity during slow wave sleep (SWS) differs from that during REM sleep or waking states. During SWS there are traveling patterns of activity in the cortex; in other states stationary patterns occur. Our model is a network composed of Hodgkin-Huxley type neurons with a M-current regulated by ACh. Regulation of ACh level can account for dynamical changes between functional regimes. Reduction of the magnitude of this current recreates the reduction in SFA the shift from a type 2 to a type 1 PRC observed in the presence of ACh. When SFA is minimal (in waking or REM sleep state, high ACh) patterns of activity are localized and easily pinned by network inhomogeneities. When SFA is present (decreasing ACh), traveling waves of activity naturally arise. A further decrease in ACh leads to a high degree of synchrony within traveling waves. We also show that the level of ACh determines how sensitive network activity is to synaptic heterogeneity. These regimes may have a profound functional significance as stationary patterns may play a role in the proper encoding of external input as memory and traveling waves could lead to synaptic regularization, giving unique insights into the role and significance of ACh in determining patterns of cortical activity and functional differences arising from the patterns.

ABSTRACT Swarming bacteria are challenged by the need to invade hostile environments. Swarms of the flagellated bacterium Paenibacillus vortex can collectively transport other microorganisms. Here we show that P. vortex can invade toxic environments by carrying antibiotic-degrading bacteria; this transport is mediated by a specialized, phenotypic subpopulation utilizing a process not dependent on cargo motility. Swarms of beta-lactam antibiotic (BLA)-sensitive P. vortex used beta-lactamase-producing, resistant, cargo bacteria to detoxify BLAs in their path. In the presence of BLAs, both transporter and cargo bacteria gained from this temporary cooperation; there was a positive correlation between BLA resistance and dispersal. P. vortex transported only the most beneficial antibiotic-resistant cargo (including environmental and clinical isolates) in a sustained way. P. vortex displayed a bet-hedging strategy that promoted the colonization of nontoxic niches by P. vortex alone; when detoxifying cargo bacteria were not needed, they were lost. This work has relevance for the dispersal of antibiotic-resistant microorganisms and for strategies for asymmetric cooperation with agricultural and medical implications. IMPORTANCE Antibiotic resistance is a major health threat. We show a novel mechanism for the local spread of antibiotic resistance. This involves interactions between different bacteria: one species provides an enzyme that detoxifies the antibiotic (a sessile cargo bacterium carrying a resistance gene), while the other (Paenibacillus vortex) moves itself and transports the cargo. P. vortex used a bet-hedging strategy, colonizing new environments alone when the cargo added no benefit, but cooperating when the cargo was needed. This work is of interest in an evolutionary context and sheds light on fundamental questions, such as how environmental antibiotic resistance may lead to clinical resistance and also microbial social organization, as well as the costs, benefits, and risks of dispersal in the environment.

The rapid increase of wealth inequality in the past few decades is a most disturbing social and economic issue of our time. In order to control, and even reverse that surge, its origin and underlying mechanisms should be revealed. One of the challenges in studying these mechanisms is to incorporate realistic individual dynamics in the population level in a self-consistent manner. Our theoretical approach meets the challenge by using interacting multi-agent master-equations to model the dynamics of wealth inequality. The model is solved using stochastic multi-agent iterated maps. Taking into account growth rate, return on capital, private savings and economic mobility, we were able to capture the historical dynamics of wealth inequality in the United States during the course of the 20th century. We show that the fraction of capital income in the national income and the fraction of private savings are the critical factors that govern the wealth inequality dynamics. In addition, we found that economic mobility plays a crucial role in wealth accumulation. Notably, we found that the major decrease in private savings since the 1980s could be associated primarily with the recent surge in wealth inequality and if nothing changes in this respect we predict further increase in wealth inequality in the future. However, the 2007-08 financial crisis brought an opportunity to restrain the wealth inequality surge by increasing private savings. If this trend continues, it may lead to prevention, and even reversing, of the ongoing inequality surge.

Fibromyalgia Syndrome (FMS) is a persistent and debilitating disorder estimated to impair the quality of life of 2-4% of the population, with 9:1 female-to-male incidence ratio. FMS is an important representative example of central nervous system sensitization and is associated with abnormal brain activity. Key symptoms include chronic widespread pain, allodynia and diffuse tenderness, along with fatigue and sleep disturbance. The syndrome is still elusive and refractory. The goal of this study was to evaluate the effect of hyperbaric oxygen therapy (HBOT) on symptoms and brain activity in FMS.A prospective, active control, crossover clinical trial. Patients were randomly assigned to treated and crossover groups: The treated group patients were evaluated at baseline and after HBOT. Patients in the crossover-control group were evaluated three times: baseline, after a control period of no treatment, and after HBOT. Evaluations consisted of physical examination, including tender point count and pain threshold, extensive evaluation of quality of life, and single photon emission computed tomography (SPECT) imaging for evaluation of brain activity. The HBOT protocol comprised 40 sessions, 5 days/week, 90 minutes, 100% oxygen at 2ATA. Sixty female patients were included, aged 21-67 years and diagnosed with FMS at least 2 years earlier. HBOT in both groups led to significant amelioration of all FMS symptoms, with significant improvement in life quality. Analysis of SPECT imaging revealed rectification of the abnormal brain activity: decrease of the hyperactivity mainly in the posterior region and elevation of the reduced activity mainly in frontal areas. No improvement in any of the parameters was observed following the control period.The study provides evidence that HBOT can improve the symptoms and life quality of FMS patients. Moreover, it shows that HBOT can induce neuroplasticity and significantly rectify abnormal brain activity in pain related areas of FMS patients.ClinicalTrials.gov NCT01827683.

Neuron-astrocyte communication is an important regulatory mechanism in various brain functions but its complexity and role are yet to be fully understood. In particular, the temporal pattern of astrocyte response to neuronal firing has not been fully characterized. Here, we used neuron-astrocyte cultures on multi-electrode arrays coupled to Ca2+ imaging and explored the range of neuronal stimulation frequencies while keeping constant the amount of stimulation. Our results reveal that astrocytes specifically respond to the frequency of neuronal stimulation by intracellular Ca2+ transients, with a clear onset of astrocytic activation at neuron firing rates around 3-5 Hz. The cell-to-cell heterogeneity of the astrocyte Ca2+ response was however large and increasing with stimulation frequency. Astrocytic activation by neurons was abolished with antagonists of type I metabotropic glutamate receptor, validating the glutamate-dependence of this neuron-to-astrocyte pathway. Using a realistic biophysical model of glutamate-based intracellular calcium signaling in astrocytes, we suggest that the stepwise response is due to the supralinear dynamics of intracellular IP3 and that the heterogeneity of the responses may be due to the heterogeneity of the astrocyte-to-astrocyte couplings via gap junction channels. Therefore our results present astrocyte intracellular Ca2+ activity as a nonlinear integrator of glutamate-dependent neuronal activity.

It has recently been discovered that single neuron stimulation can impact network dynamics in immature and adult neuronal circuits. Here we report a novel mechanism which can explain in neuronal circuits, at an early stage of development, the peculiar role played by a few specific neurons in promoting/arresting the population activity. For this purpose, we consider a standard neuronal network model, with short-term synaptic plasticity, whose population activity is characterized by bursting behavior. The addition of developmentally inspired constraints and correlations in the distribution of the neuronal connectivities and excitabilities leads to the emergence of functional hub neurons, whose stimulation/deletion is critical for the network activity. Functional hubs form a clique, where a precise sequential activation of the neurons is essential to ignite collective events without any need for a specific topological architecture. Unsupervised time-lagged firings of supra-threshold cells, in connection with coordinated entrainments of near-threshold neurons, are the key ingredients to orchestrate population activity.

The ability of cells to coordinately migrate in groups is crucial to enable them to travel long distances during embryonic development, wound healing and tumorigenesis, but the fundamental mechanisms underlying intercellular coordination during collective cell migration remain elusive despite considerable research efforts. A novel analytical framework is introduced here to explicitly detect and quantify cell clusters that move coordinately in a monolayer. The analysis combines and associates vast amount of spatiotemporal data across multiple experiments into transparent quantitative measures to report the emergence of new modes of organized behavior during collective migration of tumor and epithelial cells in wound healing assays. First, we discovered the emergence of a wave of coordinated migration propagating backward from the wound front, which reflects formation of clusters of coordinately migrating cells that are generated further away from the wound edge and disintegrate close to the advancing front. This wave emerges in both normal and tumor cells, and is amplified by Met activation with hepatocyte growth factor/scatter factor. Second, Met activation was found to induce coinciding waves of cellular acceleration and stretching, which in turn trigger the emergence of a backward propagating wave of directional migration with about an hour phase lag. Assessments of the relations between the waves revealed that amplified coordinated migration is associated with the emergence of directional migration. Taken together, our data and simplified modeling-based assessments suggest that increased velocity leads to enhanced coordination: higher motility arises due to acceleration and stretching that seems to increase directionality by temporarily diminishing the velocity components orthogonal to the direction defined by the monolayer geometry. Spatial and temporal accumulation of directionality thus defines coordination. The findings offer new insight and suggest a basic cellular mechanism for long-term cell guidance and intercellular communication during collective cell migration.

The principal interactions leading to the emergence of order in swarms of marching locust nymphs was studied both experimentally, using small groups of marching locusts in the lab, and using computer simulations. We utilized a custom tracking algorithm to reveal fundamental animal-animal interactions leading to collective motion. Uncovering this behavior introduced a new agent-based modeling approach in which pause-and-go motion is pivotal. The behavioral and modeling findings are largely based on motion-related visual sensory inputs obtained by the individual locust. Results suggest a generic principle, in which intermittent animal motion can be considered as a sequence of individual decisions as animals repeatedly reassess their situation and decide whether or not to swarm. This interpretation implies, among other things, some generic characteristics regarding the build-up and emergence of collective order in swarms: in particular, that order and disorder are generic meta-stable states of the system, suggesting that the emergence of order is kinetic and does not necessarily require external environmental changes. This work calls for further experimental as well as theoretical investigation of the neural mechanisms underlying locust coordinative behavior.

The characterization of asset price returns is an important subject in modern finance. Traditionally, the dynamics of stock returns are assumed to lack any temporal order. Here we present an analysis of the autocovariance of stock market indices and unravel temporal order in several major stock markets. We also demonstrate a fundamental difference between developed and emerging markets in the past decade - emerging markets are marked by positive order in contrast to developed markets whose dynamics are marked by weakly negative order. In addition, the reaction to financial crises was found to be reversed among developed and emerging markets, presenting large positive/negative autocovariance spikes following the onset of these crises. Notably, the Chinese market shows neutral or no order while being regarded as an emerging market. These findings show that despite the coupling between international markets and global trading, major differences exist between different markets, and demonstrate that the autocovariance of markets is correlated with their stability, as well as with their state of development.

Traumatic brain injury (TBI) is the leading cause of death and disability in the US. Approximately 70-90% of the TBI cases are classified as mild, and up to 25% of them will not recover and suffer chronic neurocognitive impairments. The main pathology in these cases involves diffuse brain injuries, which are hard to detect by anatomical imaging yet noticeable in metabolic imaging. The current study tested the effectiveness of Hyperbaric Oxygen Therapy (HBOT) in improving brain function and quality of life in mTBI patients suffering chronic neurocognitive impairments.The trial population included 56 mTBI patients 1-5 years after injury with prolonged post-concussion syndrome (PCS). The HBOT effect was evaluated by means of prospective, randomized, crossover controlled trial: the patients were randomly assigned to treated or crossover groups. Patients in the treated group were evaluated at baseline and following 40 HBOT sessions; patients in the crossover group were evaluated three times: at baseline, following a 2-month control period of no treatment, and following subsequent 2-months of 40 HBOT sessions. The HBOT protocol included 40 treatment sessions (5 days/week), 60 minutes each, with 100% oxygen at 1.5 ATA. "Mindstreams" was used for cognitive evaluations, quality of life (QOL) was evaluated by the EQ-5D, and changes in brain activity were assessed by SPECT imaging. Significant improvements were demonstrated in cognitive function and QOL in both groups following HBOT but no significant improvement was observed following the control period. SPECT imaging revealed elevated brain activity in good agreement with the cognitive improvements.HBOT can induce neuroplasticity leading to repair of chronically impaired brain functions and improved quality of life in mTBI patients with prolonged PCS at late chronic stage.ClinicalTrials.gov NCT00715052.

Recovery after stroke correlates with non-active (stunned) brain regions, which may persist for years. The current study aimed to evaluate whether increasing the level of dissolved oxygen by Hyperbaric Oxygen Therapy (HBOT) could activate neuroplasticity in patients with chronic neurologic deficiencies due to stroke.A prospective, randomized, controlled trial including 74 patients (15 were excluded). All participants suffered a stroke 6-36 months prior to inclusion and had at least one motor dysfunction. After inclusion, patients were randomly assigned to "treated" or "cross" groups. Brain activity was assessed by SPECT imaging; neurologic functions were evaluated by NIHSS, ADL, and life quality. Patients in the treated group were evaluated twice: at baseline and after 40 HBOT sessions. Patients in the cross group were evaluated three times: at baseline, after a 2-month control period of no treatment, and after subsequent 2-months of 40 HBOT sessions. HBOT protocol: Two months of 40 sessions (5 days/week), 90 minutes each, 100% oxygen at 2 ATA. We found that the neurological functions and life quality of all patients in both groups were significantly improved following the HBOT sessions while no improvement was found during the control period of the patients in the cross group. Results of SPECT imaging were well correlated with clinical improvement. Elevated brain activity was detected mostly in regions of live cells (as confirmed by CT) with low activity (based on SPECT) - regions of noticeable discrepancy between anatomy and physiology.The results indicate that HBOT can lead to significant neurological improvements in post stroke patients even at chronic late stages. The observed clinical improvements imply that neuroplasticity can still be activated long after damage onset in regions where there is a brain SPECT/CT (anatomy/physiology) mismatch.

Insulin resistance (IR) is one of the most widespread health problems in modern times. The gold standard for quantification of IR is the hyperinsulinemic-euglycemic glucose clamp technique. During the test, a regulated glucose infusion is delivered intravenously to maintain a constant blood glucose concentration. Current control algorithms for regulating this glucose infusion are based on feedback control. These models require frequent sampling of blood, and can only partly capture the complexity associated with regulation of glucose. Here we present an improved clamp control algorithm which is motivated by the stochastic nature of glucose kinetics, while using the minimal need in blood samples required for evaluation of IR. A glucose pump control algorithm, based on artificial neural networks model was developed. The system was trained with a data base collected from 62 rat model experiments, using a back-propagation Levenberg-Marquardt optimization. Genetic algorithm was used to optimize network topology and learning features. The predictive value of the proposed algorithm during the temporal period of interest was significantly improved relative to a feedback control applied at an equivalent low sampling interval. Robustness to noise analysis demonstrates the applicability of the algorithm in realistic situations.

Collective cell migration plays a major role in embryonic morphogenesis, tissue remodeling, wound repair and cancer invasion. Despite many decades of extensive investigations, only few analytical tools have been developed to enhance the biological understanding of this important phenomenon. Here we present a novel quantitative approach to analyze long term kinetics of bright field time-lapse wound healing. Fully-automated spatiotemporal measures and visualization of cells' motility and implicit morphology were proven to be sound, repetitive and highly informative compared to single-cell tracking analysis. We study cellular collective migration induced by tyrosine kinase-growth factor signaling (Met-Hepatocyte Growth Factor/Scatter Factor (HGF/SF)). Our quantitative approach is applied to demonstrate that collective migration of the adenocarcinoma cell lines is characterized by simple morpho-kinetics. HGF/SF induces complex morpho-kinetic coordinated collective migration: cells at the front move faster and are more spread than those further away from the wound edge. As the wound heals, distant cells gradually accelerate and enhance spread and elongation -resembling the epithelial to mesenchymal transition (EMT), and then the cells become more spread and maintain higher velocity than cells located closer to the wound. Finally, upon wound closure, front cells halt, shrink and round up (resembling mesenchymal to epithelial transition (MET) phenotype) while distant cells undergo the same process gradually. Met inhibition experiments further validate that Met signaling dramatically alters the morpho-kinetic dynamics of the healing wound. Machine-learning classification was applied to demonstrate the generalization of our findings, revealing even subtle changes in motility patterns induced by Met-inhibition. It is concluded that activation of Met-signaling induces an elaborated model in which cells lead a coordinated increased motility along with gradual differentiation-based collective cell motility dynamics. Our quantitative phenotypes may guide future investigation on the molecular and cellular mechanisms of tyrosine kinase-induced coordinate cell motility and morphogenesis in metastasis.

BACKGROUND: In the current era of strong worldwide market couplings the global financial village became highly prone to systemic collapses, events that can rapidly sweep throughout the entire village. METHODOLOGY/PRINCIPAL FINDINGS: We present a new methodology to assess and quantify inter-market relations. The approach is based on the correlations between the market index, the index volatility, the market Index Cohesive Force and the meta-correlations (correlations between the intra-correlations.) We investigated the relations between six important world markets--U.S., U.K., Germany, Japan, China and India--from January 2000 until December 2010. We found that while the developed "western" markets (U.S., U.K., Germany) are highly correlated, the interdependencies between these markets and the developing "eastern" markets (India and China) are volatile and with noticeable maxima at times of global world events. The Japanese market switches "identity"--it switches between periods of high meta-correlations with the "western" markets and periods when it behaves more similarly to the "eastern" markets. CONCLUSIONS/SIGNIFICANCE: The methodological framework presented here provides a way to quantify the evolvement of interdependencies in the global market, evaluate a world financial network and quantify changes in the world inter market relations. Such changes can be used as precursors to the agitation of the global financial village. Hence, the new approach can help to develop a sensitive "financial seismograph" to detect early signs of global financial crises so they can be treated before they develop into worldwide events.

Short-term presynaptic plasticity designates variations of the amplitude of synaptic information transfer whereby the amount of neurotransmitter released upon presynaptic stimulation changes over seconds as a function of the neuronal firing activity. While a consensus has emerged that the resulting decrease (depression) and/or increase (facilitation) of the synapse strength are crucial to neuronal computations, their modes of expression in vivo remain unclear. Recent experimental studies have reported that glial cells, particularly astrocytes in the hippocampus, are able to modulate short-term plasticity but the mechanism of such a modulation is poorly understood. Here, we investigate the characteristics of short-term plasticity modulation by astrocytes using a biophysically realistic computational model. Mean-field analysis of the model, supported by intensive numerical simulations, unravels that astrocytes may mediate counterintuitive effects. Depending on the expressed presynaptic signaling pathways, astrocytes may globally inhibit or potentiate the synapse: the amount of released neurotransmitter in the presence of the astrocyte is transiently smaller or larger than in its absence. But this global effect usually coexists with the opposite local effect on paired pulses: with release-decreasing astrocytes most paired pulses become facilitated, namely the amount of neurotransmitter released upon spike i+1 is larger than that at spike i, while paired-pulse depression becomes prominent under release-increasing astrocytes. Moreover, we show that the frequency of astrocytic intracellular Ca(2+) oscillations controls the effects of the astrocyte on short-term synaptic plasticity. Our model explains several experimental observations yet unsolved, and uncovers astrocytic gliotransmission as a possible transient switch between short-term paired-pulse depression and facilitation. This possibility has deep implications on the processing of neuronal spikes and resulting information transfer at synapses.

Collective navigation and swarming have been studied in animal groups, such as fish schools, bird flocks, bacteria, and slime molds. Computer modeling has shown that collective behavior of simple agents can result from simple interactions between the agents, which include short range repulsion, intermediate range alignment, and long range attraction. Here we study collective navigation of bacteria-inspired smart agents in complex terrains, with adaptive interactions that depend on performance. More specifically, each agent adjusts its interactions with the other agents according to its local environment--by decreasing the peers' influence while navigating in a beneficial direction, and increasing it otherwise. We show that inclusion of such performance dependent adaptable interactions significantly improves the collective swarming performance, leading to highly efficient navigation, especially in complex terrains. Notably, to afford such adaptable interactions, each modeled agent requires only simple computational capabilities with short-term memory, which can easily be implemented in simple swarming robots.

Many eukaryotic cells are able to crawl on surfaces and guide their motility based on environmental cues. These cues are interpreted by signaling systems which couple to cell mechanics; indeed membrane protrusions in crawling cells are often accompanied by activated membrane patches, which are localized areas of increased concentration of one or more signaling components. To determine how these patches are related to cell motion, we examine the spatial localization of RasGTP in chemotaxing Dictyostelium discoideum cells under conditions where the vertical extent of the cell was restricted. Quantitative analyses of the data reveal a high degree of spatial correlation between patches of activated Ras and membrane protrusions. Based on these findings, we formulate a model for amoeboid cell motion that consists of two coupled modules. The first module utilizes a recently developed two-component reaction diffusion model that generates transient and localized areas of elevated concentration of one of the components along the membrane. The activated patches determine the location of membrane protrusions (and overall cell motion) that are computed in the second module, which also takes into account the cortical tension and the availability of protrusion resources. We show that our model is able to produce realistic amoeboid-like motion and that our numerical results are consistent with experimentally observed pseudopod dynamics. Specifically, we show that the commonly observed splitting of pseudopods can result directly from the dynamics of the signaling patches.

MotivationNew antigen microarray technology enables parallel recording of antibody reactivities with hundreds of antigens. Such data affords system level analysis of the immune system's organization using methods and approaches from network theory. Here we measured the reactivity of 290 antigens (for both the IgG and IgM isotypes) of 10 healthy mothers and their term newborns. We constructed antigen correlation networks (or immune networks) whose nodes are the antigens and the edges are the antigen-antigen reactivity correlations, and we also computed their corresponding minimum spanning trees (MST)--maximal information reduced sub-graphs. We quantify the network organization (topology) in terms of the network theory divergence rate measure and rank the antigen importance in the full antigen correlation networks by the eigen-value centrality measure. This analysis makes possible the characterization and comparison of the IgG and IgM immune networks at birth (newborns) and adulthood (mothers) in terms of topology and node importance.ResultsComparison of the immune network topology at birth and adulthood revealed partial conservation of the IgG immune network topology, and significant reorganization of the IgM immune networks. Inspection of the antigen importance revealed some dominant (in terms of high centrality) antigens in the IgG and IgM networks at birth, which retain their importance at adulthood.

BACKGROUND: Semantic memory has generated much research. As such, the majority of investigations have focused on the English language, and much less on other languages, such as Hebrew. Furthermore, little research has been done on search processes within the semantic network, even though they are abundant within cognitive semantic phenomena. METHODOLOGY/PRINCIPAL FINDINGS: We examine a unique dataset of free association norms to a set of target words and make use of correlation and network theory methodologies to investigate the global and local features of the Hebrew lexicon. The global features of the lexicon are investigated through the use of association correlations--correlations between target words, based on their association responses similarity; the local features of the lexicon are investigated through the use of association dependencies--the influence words have in the network on other words. CONCLUSIONS/SIGNIFICANCE: Our investigation uncovered Small-World Network features of the Hebrew lexicon, specifically a high clustering coefficient and a scale-free distribution, and provides means to examine how words group together into semantically related 'free categories'. Our novel approach enables us to identify how words facilitate or inhibit the spread of activation within the network, and how these words influence each other. We discuss how these properties relate to classical research on spreading activation and suggest that these properties influence cognitive semantic search processes. A semantic search task, the Remote Association Test is discussed in light of our findings.

Living cells of many types need to move in response to external stimuli in order to accomplish their functional tasks; these tasks range from wound healing to immune response to fertilization. While the directional motion is typically dictated by an external signal, the actual motility is also restricted by physical constraints, such as the presence of other cells and the extracellular matrix. The ability to successfully navigate in the presence of obstacles is not only essential for organisms, but might prove relevant in the study of autonomous robotic motion.We study a computational model of amoeboid chemotactic navigation under differing conditions, from motion in an obstacle-free environment to navigation between obstacles and finally to moving in a maze. We use the maze as a simple stand-in for a motion task with severe constraints, as might be expected in dense extracellular matrix. Whereas agents using simple chemotaxis can successfully navigate around small obstacles, the presence of large barriers can often lead to agent trapping. We further show that employing a simple memory mechanism, namely secretion of a repulsive chemical by the agent, helps the agent escape from such trapping.Our main conclusion is that cells employing simple chemotactic strategies will often be unable to navigate through maze-like geometries, but a simple chemical marker mechanism (which we refer to as "self-assistance") significantly improves success rates. This realization provides important insights into mechanisms that might be employed by real cells migrating in complex environments as well as clues for the design of robotic navigation strategies. The results can be extended to more complicated multi-cellular systems and can be used in the study of mammalian cell migration and cancer metastasis.

BACKGROUND: The 2007-2009 financial crisis, and its fallout, has strongly emphasized the need to define new ways and measures to study and assess the stock market dynamics. METHODOLOGY/PRINCIPAL FINDINGS: The S&P500 dynamics during 4/1999-4/2010 is investigated in terms of the index cohesive force (ICF--the balance between the stock correlations and the partial correlations after subtraction of the index contribution), and the Eigenvalue entropy of the stock correlation matrices. We found a rapid market transition at the end of 2001 from a flexible state of low ICF into a stiff (nonflexible) state of high ICF that is prone to market systemic collapses. The stiff state is also marked by strong effect of the market index on the stock-stock correlations as well as bursts of high stock correlations reminiscence of epileptic brain activity. CONCLUSIONS/SIGNIFICANCE: The market dynamical states, stability and transition between economic states was studies using new quantitative measures. Doing so shed new light on the origin and nature of the current crisis. The new approach is likely to be applicable to other classes of complex systems from gene networks to the human brain.

Confocal microscopy analysis of fluorescence and morphology is becoming the standard tool in cell biology and molecular imaging. Accurate quantification algorithms are required to enhance the understanding of different biological phenomena. We present a novel approach based on image-segmentation of multi-cellular regions in bright field images demonstrating enhanced quantitative analyses and better understanding of cell motility. We present MultiCellSeg, a segmentation algorithm to separate between multi-cellular and background regions for bright field images, which is based on classification of local patches within an image: a cascade of Support Vector Machines (SVMs) is applied using basic image features. Post processing includes additional classification and graph-cut segmentation to reclassify erroneous regions and refine the segmentation. This approach leads to a parameter-free and robust algorithm. Comparison to an alternative algorithm on wound healing assay images demonstrates its superiority. The proposed approach was used to evaluate common cell migration models such as wound healing and scatter assay. It was applied to quantify the acceleration effect of Hepatocyte growth factor/scatter factor (HGF/SF) on healing rate in a time lapse confocal microscopy wound healing assay and demonstrated that the healing rate is linear in both treated and untreated cells, and that HGF/SF accelerates the healing rate by approximately two-fold. A novel fully automated, accurate, zero-parameters method to classify and score scatter-assay images was developed and demonstrated that multi-cellular texture is an excellent descriptor to measure HGF/SF-induced cell scattering. We show that exploitation of textural information from differential interference contrast (DIC) images on the multi-cellular level can prove beneficial for the analyses of wound healing and scatter assays. The proposed approach is generic and can be used alone or alongside traditional fluorescence single-cell processing to perform objective, accurate quantitative analyses for various biological applications.

The spontaneous activity of engineered quadruple cultured neural networks (of four coupled sub-networks) exhibits a repertoire of different types of mutual synchronization events. Each event corresponds to a specific activity propagation mode defined by the order of activity propagation between the sub-networks. We statistically characterized the frequency of spontaneous appearance of the different types of activity propagation modes (APMs). The relative frequencies of the APMs were then examined for their power law properties. We found that the frequencies of appearance of the leading (most frequent) APMs have close to constant algebraic ratio reminiscent of Zipf's scaling of words. We show that the observations are consistent with a simplified wrestling model. This model represents an extension of the boxing arena model which was previously proposed to describe the ratio between the two activity modes in two coupled sub-networks. The additional new element in the wrestling model presented here is that the firing within each network is modeled by a time interval generator with similar intra-network Lévy distribution. We modeled the different burst-initiation-zones’ (BIZ) interaction by competition between the stochastic generators with Gaussian inter-network variability. Estimation of the model parameters revealed similarity across different cultures while the inter-burst-interval (IBI) of the cultures was similar across different APMs as numerical simulation of the model predicts.

What are the dominant stocks which drive the correlations present among stocks traded in a stock market? Can a correlation analysis provide an answer to this question? In the past, correlation based networks have been proposed as a tool to uncover the underlying backbone of the market. Correlation based networks represent the stocks and their relationships, which are then investigated using different network theory methodologies. Here we introduce a new concept to tackle the above question--the partial correlation network. Partial correlation is a measure of how the correlation between two variables, e.g., stock returns, is affected by a third variable. By using it we define a proxy of stock influence, which is then used to construct partial correlation networks. The empirical part of this study is performed on a specific financial system, namely the set of 300 highly capitalized stocks traded at the New York Stock Exchange, in the time period 2001-2003. By constructing the partial correlation network, unlike the case of standard correlation based networks, we find that stocks belonging to the financial sector and, in particular, to the investment services sub-sector, are the most influential stocks affecting the correlation profile of the system. Using a moving window analysis, we find that the strong influence of the financial stocks is conserved across time for the investigated trading period. Our findings shed a new light on the underlying mechanisms and driving forces controlling the correlation profile observed in a financial market.

BACKGROUND: Information processing in neuronal networks relies on the network's ability to generate temporal patterns of action potentials. Although the nature of neuronal network activity has been intensively investigated in the past several decades at the individual neuron level, the underlying principles of the collective network activity, such as the synchronization and coordination between neurons, are largely unknown. Here we focus on isolated neuronal clusters in culture and address the following simple, yet fundamental questions: What is the minimal number of cells needed to exhibit collective dynamics? What are the internal temporal characteristics of such dynamics and how do the temporal features of network activity alternate upon crossover from minimal networks to large networks? METHODOLOGY/PRINCIPAL FINDINGS: We used network engineering techniques to induce self-organization of cultured networks into neuronal clusters of different sizes. We found that small clusters made of as few as 40 cells already exhibit spontaneous collective events characterized by innate synchronous network oscillations in the range of 25 to 100 Hz. The oscillation frequency of each network appeared to be independent of cluster size. The duration and rate of the network events scale with cluster size but converge to that of large uniform networks. Finally, the investigation of two coupled clusters revealed clear activity propagation with master/slave asymmetry. CONCLUSIONS/SIGNIFICANCE: The nature of the activity patterns observed in small networks, namely the consistent emergence of similar activity across networks of different size and morphology, suggests that neuronal clusters self-regulate their activity to sustain network bursts with internal oscillatory features. We therefore suggest that clusters of as few as tens of cells can serve as a minimal but sufficient functional network, capable of sustaining oscillatory activity. Interestingly, the frequencies of these oscillations are similar those observed in vivo.

The Spindle Assembly Checkpoint (SAC) is an intracellular mechanism that ensures proper chromosome segregation. By inhibiting Cdc20, a co-factor of the Anaphase Promoting Complex (APC), the checkpoint arrests the cell cycle until all chromosomes are properly attached to the mitotic spindle. Inhibition of Cdc20 is mediated by a conserved network of interacting proteins. The individual functions of these proteins are well characterized, but understanding of their integrated function is still rudimentary. We here describe our attempts to reverse-engineer the SAC network based on gene deletion phenotypes. We begun by formulating a general model of the SAC which enables us to predict the rate of chromosomal missegregation for any putative set of interactions between the SAC proteins. Next the missegregation rates of seven yeast strains are measured in response to the deletion of one or two checkpoint proteins. Finally, we searched for the set of interactions that correctly predicted the observed missegregation rates of all deletion mutants. Remarkably, although based on only seven phenotypes, the consistent network we obtained successfully reproduces many of the known properties of the SAC. Further insights provided by our analysis are discussed.

We present here assessment of the latent market information embedded in the raw, affinity (normalized), and partial correlations. We compared the Zipf plot, spectrum, and distribution of the eigenvalues for each matrix with the results of the corresponding random matrix. The analysis was performed on stocks belonging to the New York and Tel Aviv Stock Exchange, for the time period of January 2000 to March 2009. Our results show that in comparison to the raw correlations, the affinity matrices highlight the dominant factors of the system, and the partial correlation matrices contain more information. We propose that significant stock market information, which cannot be captured by the raw correlations, is embedded in the affinity and partial correlations. Our results further demonstrate the differences between NY and TA markets.

DNA chips allow simultaneous measurements of genome-wide response of thousands of genes, i.e. system level monitoring of the gene-network activity. Advanced analysis methods have been developed to extract meaningful information from the vast amount of raw gene-expression data obtained from the microarray measurements. These methods usually aimed to distinguish between groups of subjects (e.g., cancer patients vs. healthy subjects) or identifying marker genes that help to distinguish between those groups. We assumed that motifs related to the internal structure of operons and gene-networks regulation are also embedded in microarray and can be deciphered by using proper analysis.The analysis presented here is based on investigating the gene-gene correlations. We analyze a database of gene expression of Bacillus subtilis exposed to sub-lethal levels of 37 different antibiotics. Using unsupervised analysis (dendrogram) of the matrix of normalized gene-gene correlations, we identified the operons as they form distinct clusters of genes in the sorted correlation matrix. Applying dimension-reduction algorithm (Principal Component Analysis, PCA) to the matrices of normalized correlations reveals functional motifs. The genes are placed in a reduced 3-dimensional space of the three leading PCA eigen-vectors according to their corresponding eigen-values. We found that the organization of the genes in the reduced PCA space recovers motifs of the operon internal structure, such as the order of the genes along the genome, gene separation by non-coding segments, and translational start and end regions. In addition to the intra-operon structure, it is also possible to predict inter-operon relationships, operons sharing functional regulation factors, and more. In particular, we demonstrate the above in the context of the competence and sporulation pathways.We demonstrated that by analyzing gene-gene correlation from gene-expression data it is possible to identify operons and to predict unknown internal structure of operons and gene-networks regulation.

Notch pathway is an evolutionarily conserved cell–cell communication mechanism governing cell-fate during development and tumor progression. It is activated when Notch receptor of one cell binds to either of its ligand—Delta or Jagged—of another cell. Notch–Delta (ND) signaling forms a two-way switch, and two cells interacting via ND signaling adopt different fates—Sender (high ligand, low receptor) and Receiver (low ligand, high receptor). Notch–Delta–Jagged signaling (NDJ) behaves as a three-way switch and enables an additional fate—hybrid Sender/Receiver (S/R) (medium ligand, medium receptor). Here, by extending our framework of NDJ signaling for a two-cell system, we show that higher production rate of Jagged, but not that of Delta, expands the range of parameters for which both cells attain the hybrid S/R state. Conversely, glycosyltransferase Fringe and cis -inhibition reduces this range of conditions, and reduces the relative stability of the hybrid S/R state, thereby promoting cell-fate divergence and consequently lateral inhibition-based patterns. Lastly, soluble Jagged drives the cells to attain the hybrid S/R state, and soluble Delta drives them to be Receivers. We also discuss the critical role of hybrid S/R state in promoting cancer metastasis by enabling collective cell migration and expanding cancer stem cell (CSC) population.

Modelling the behaviour of stock markets has been of major interest in the past century. The market can be treated as a network of many investors reacting in accordance to their group behaviour, as manifested by the index and effected by the flow of external information into the system. Here we devise a model that encapsulates the behaviour of stock markets. The model consists of two terms, demonstrating quantitatively the effect of the individual tendency to follow the group and the effect of the individual reaction to the available information. Using the above factors we were able to explain several key features of the stock market: the high correlations between the individual stocks and the index; the Epps effect; the high fluctuating nature of the market, which is similar to real market behaviour. Furthermore, intricate long term phenomena are also described by this model, such as bursts of synchronized average correlation and the dominance of the index as demonstrated through partial correlation.

Bacterial swarming resulting in collective navigation over surfaces provides a valuable example of cooperative colonization of new territories. The social bacterium Paenibacillus vortex exhibits successful and diverse swarming strategies. When grown on hard agar surfaces with peptone , P. vortex develops complex colonies of vortices (rotating bacterial aggregates). In contrast, during growth on Mueller–Hinton broth gelled into a soft agar surface, a new strategy of multi-level organization is revealed: the colonies are organized into a special network of swarms (or ‘snakes’ of a fraction of millimeter in width) with intricate internal traffic. More specifically, cell movement is organized in two or three lanes of bacteria traveling between the back and the front of the swarm. This special form of cellular logistics suggests new methods in which bacteria can share resources and risk while searching for food or migrating into new territories. While the vortices-based organization on hard agar surfaces has been modeled before, here, we introduce a new multi-agent bacterial swarming model devised to capture the swarms-based organization on soft surfaces. We test two putative generic mechanisms that may underlie the observed swarming logistics: (i) chemo-activated taxis in response to chemical cues and (ii) special align-and-push interactions between the bacteria and the boundary of the layer of lubricant collectively generated by the swarming bacteria. Using realistic parameters, the model captures the observed phenomena with semi-quantitative agreement in terms of the velocity as well as the dynamics of the swarm and its envelope. This agreement implies that the bacteria interactions with the swarm boundary play a crucial role in mediating the interplay between the collective movement of the swarm and the internal traffic dynamics.