Toward understanding cancer stem cell heterogeneity in the tumor microenvironment.

The epithelial-mesenchymal transition (EMT) and cancer stem cell (CSC) formation are two paramount processes driving tumor progression, therapy resistance, and cancer metastasis. Recent experiments show that cells with varying EMT and CSC phenotypes are spatially segregated in the primary tumor. The underlying mechanisms generating such spatiotemporal dynamics in the tumor microenvironment, however, remain largely unexplored. Here, we show through a mechanism-based dynamical model that the diffusion of EMT-inducing signals such as TGF-ß, together with noncell autonomous control of EMT and CSC decision making via the Notch signaling pathway, can explain experimentally observed disparate localization of subsets of CSCs with varying EMT phenotypes in the tumor. Our simulations show that the more mesenchymal CSCs lie at the invasive edge, while the hybrid epithelial/mesenchymal (E/M) CSCs reside in the tumor interior. Further, motivated by the role of Notch-Jagged signaling in mediating EMT and stemness, we investigated the microenvironmental factors that promote Notch-Jagged signaling. We show that many inflammatory cytokines such as IL-6 that can promote Notch-Jagged signaling can (i) stabilize a hybrid E/M phenotype, (ii) increase the likelihood of spatial proximity of hybrid E/M cells, and (iii) expand the fraction of CSCs. To validate the predicted connection between Notch-Jagged signaling and stemness, we knocked down JAG1 in hybrid E/M SUM149 human breast cancer cells in vitro. JAG1 knockdown significantly restricted tumor organoid formation, confirming the key role that Notch-Jagged signaling can play in tumor progression. Together, our integrated computational-experimental framework reveals the underlying principles of spatiotemporal dynamics of EMT and CSCs.

Astrocytes restore connectivity and synchronization in dysfunctional cerebellar networks.

Evidence suggests that astrocytes play key roles in structural and functional organization of neuronal circuits. To understand how astrocytes influence the physiopathology of cerebellar circuits, we cultured cells from cerebella of mice that lack the ATM gene. Mutations in ATM are causative of the human cerebellar degenerative disease ataxia-telangiectasia. Cerebellar cultures grown from Atm-/- mice had disrupted network synchronization, atrophied astrocytic arborizations, reduced autophagy levels, and higher numbers of synapses per neuron than wild-type cultures. Chimeric circuitries composed of wild-type astrocytes and Atm-/- neurons were indistinguishable from wild-type cultures. Adult cerebellar characterizations confirmed disrupted astrocyte morphology, increased GABAergic synaptic markers, and reduced autophagy in Atm-/- compared with wild-type mice. These results indicate that astrocytes can impact neuronal circuits at levels ranging from synaptic expression to global dynamics.

Immunoproteasome deficiency is a feature of non-small cell lung cancer with a mesenchymal phenotype and is associated with a poor outcome.

The immunoproteasome plays a key role in generation of HLA peptides for T cell-mediated immunity. Integrative genomic and proteomic analysis of non-small cell lung carcinoma (NSCLC) cell lines revealed significantly reduced expression of immunoproteasome components and their regulators associated with epithelial to mesenchymal transition. Low expression of immunoproteasome subunits in early stage NSCLC patients was associated with recurrence and metastasis. Depleted repertoire of HLA class I-bound peptides in mesenchymal cells deficient in immunoproteasome components was restored with either IFNγ or 5-aza-2'-deoxycytidine (5-aza-dC) treatment. Our findings point to a mechanism of immune evasion of cells with a mesenchymal phenotype and suggest a strategy to overcome immune evasion through induction of the immunoproteasome to increase the cellular repertoire of HLA class I-bound peptides.

Interrogating the topological robustness of gene regulatory circuits by randomization.

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.

Mechanically-driven phase separation in a growing bacterial colony.

Secretion of extracellular polymeric substances (EPSs) by growing bacteria is an integral part of forming biofilm-like structures. In such dense systems, mechanical interactions among the structural components can be expected to significantly contribute to morphological properties. Here, we use a particle-based modeling approach to study the self-organization of nonmotile rod-shaped bacterial cells growing on a solid substrate in the presence of self-produced EPSs. In our simulation, all of the components interact mechanically via repulsive forces, occurring as the bacterial cells grow and divide (via consuming diffusing nutrient) and produce EPSs. Based on our simulation, we show that mechanical interactions control the collective behavior of the system. In particular, we find that the presence of nonadsorbing EPSs can lead to spontaneous aggregation of bacterial cells by a depletion attraction and thereby generates phase separated patterns in the nonequilibrium growing colony. Both repulsive interactions between cell and EPSs and the overall concentration of EPSs are important factors in the self-organization in a nonequilibrium growing colony. Furthermore, we investigate the interplay of mechanics with the nutrient diffusion and consumption by bacterial cells and observe that suppression of branch formation occurs due to EPSs compared with the case where no EPS is produced.

Jagged mediates differences in normal and tumor angiogenesis by affecting tip-stalk fate decision.

Angiogenesis is critical during development, wound repair, and cancer progression. During angiogenesis, some endothelial cells adopt a tip phenotype to lead the formation of new branching vessels; the trailing stalk cells proliferate to develop the vessel. Notch and VEGF signaling mediate the selection of these tip endothelial cells. However, how Jagged, a Notch ligand that is overexpressed in cancer, affects angiogenesis remains elusive. Here, by developing a theoretical framework for Notch-Delta-Jagged-VEGF signaling, we found that higher production levels of Jagged destabilizes the tip and stalk cell fates and can give rise to a hybrid tip/stalk phenotype that leads to poorly perfused and chaotic angiogenesis, which is a hallmark of cancer. Consistently, the signaling interactions that restrict Notch-Jagged signaling, such as Fringe, cis-inhibition, and increased production of Delta, stabilize tip and stalk fates and limit the existence of hybrid tip/stalk phenotype. Our results underline how overexpression of Jagged can transform physiological angiogenesis into pathological one.

Jagged-Delta asymmetry in Notch signaling can give rise to a Sender/Receiver hybrid phenotype.

Notch signaling pathway mediates cell-fate determination during embryonic development, wound healing, and tumorigenesis. This pathway is activated when the ligand Delta or the ligand Jagged of one cell interacts with the Notch receptor of its neighboring cell, releasing the Notch Intracellular Domain (NICD) that activates many downstream target genes. NICD affects ligand production asymmetrically--it represses Delta, but activates Jagged. Although the dynamical role of Notch-Jagged signaling remains elusive, it is widely recognized that Notch-Delta signaling behaves as an intercellular toggle switch, giving rise to two distinct fates that neighboring cells adopt--Sender (high ligand, low receptor) and Receiver (low ligand, high receptor). Here, we devise a specific theoretical framework that incorporates both Delta and Jagged in Notch signaling circuit to explore the functional role of Jagged in cell-fate determination. We find that the asymmetric effect of NICD renders the circuit to behave as a three-way switch, giving rise to an additional state--a hybrid Sender/Receiver (medium ligand, medium receptor). This phenotype allows neighboring cells to both send and receive signals, thereby attaining similar fates. We also show that due to the asymmetric effect of the glycosyltransferase Fringe, different outcomes are generated depending on which ligand is dominant: Delta-mediated signaling drives neighboring cells to have an opposite fate; Jagged-mediated signaling drives the cell to maintain a similar fate to that of its neighbor. We elucidate the role of Jagged in cell-fate determination and discuss its possible implications in understanding tumor-stroma cross-talk, which frequently entails Notch-Jagged communication.

Operating principles of tristable circuits regulating cellular differentiation.

Many cell-fate decisions during embryonic development are governed by a motif comprised of two transcription factors (TFs) A and B that mutually inhibit each other and may self-activate. This motif, called as a self-activating toggle switch (SATS), can typically have three stable states (phenotypes)-two corresponding to differentiated cell fates, each of which has a much higher level of one TF than the other-[Formula: see text] or [Formula: see text]-and the third state corresponding to an 'undecided' stem-like state with similar levels of both A and B-[Formula: see text]. Furthermore, two or more SATSes can be coupled together in various topologies in different contexts, thereby affecting the coordination between multiple cellular decisions. However, two questions remain largely unanswered: (a) what governs the co-existence and relative stability of these three stable states? (b) What orchestrates the decision-making of coupled SATSes? Here, we first demonstrate that the co-existence and relative stability of the three stable states in an individual SATS can be governed by the relative strength of self-activation, external signals activating and/or inhibiting A and B, and mutual degradation between A and B. Simultaneously, we investigate the effects of these factors on the decision-making of two coupled SATSes. Our results offer novel understanding into the operating principles of individual and coupled tristable self-activating toggle switches (SATSes) regulating cellular differentiation and can yield insights into synthesizing three-way genetic circuits and understanding of cellular reprogramming.

Modularity Induced Gating and Delays in Neuronal Networks.

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.

Tumor-associated and disease-associated autoantibody repertoires in healthy colostrum and maternal and newborn cord sera.

In this work, we studied autoantibody repertoires and Ig isotypes in 71 mothers and their 104 healthy newborns (including twins and triplets delivered term or premature). Newborns receive maternal IgG Abs via the placenta before birth, but developing infants must produce their own IgM and IgA Abs. We used an Ag microarray analysis to detect binding to a selection of 295 self-Ags, compared with 27 standard foreign Ags. The magnitude of binding to specific self-Ags was found to be not less than that to the foreign Ags. As expected, each newborn shared with its mother a similar IgG repertoire-manifest as early as the 24th week of gestation. IgM and IgA autoantibody repertoires in cord sera were highly correlated among the newborns and differed from their mothers' repertoires; the latter differed in sera and milk. The autoantibodies bound to self-Ags known to be associated with tumors and to autoimmune diseases. Thus, autoantibody repertoires in healthy humans--the immunological homunculus--arise congenitally, differ in maternal milk and sera, and mark the potential of the immune system to attack tumors, beneficially, or healthy tissues, harmfully; regulation of the tissue site, the dynamics, and the response phenotype of homuncular autoimmunity very likely affects health.

Growth feedback as a basis for persister bistability.

A small fraction of cells in many bacterial populations, called persisters, are much less sensitive to antibiotic treatment than the majority. Persisters are in a dormant metabolic state, even while remaining genetically identical to the actively growing cells. Toxin and antitoxin modules in bacteria are believed to be one possible cause of persistence. A two-gene operon, HipBA, is one of many chromosomally encoded toxin and antitoxin modules in Escherichia coli and the HipA7 allelic variant was the first validated high-persistence mutant. Here, we present a stochastic model that can generate bistability of the HipBA system, via the reciprocal coupling of free HipA to the cellular growth rate. The actively growing state and the dormant state each correspond to a stable state of this model. Fluctuations enable transitions from one to the other. This model is fully in agreement with experimental data obtained with synthetic promoter constructs.

Modeling putative therapeutic implications of exosome exchange between tumor and immune cells.

Development of effective strategies to mobilize the immune system as a therapeutic modality in cancer necessitates a better understanding of the contribution of the tumor microenvironment to the complex interplay between cancer cells and the immune response. Recently, effort has been directed at unraveling the functional role of exosomes and their cargo of messengers in this interplay. Exosomes are small vesicles (30-200 nm) that mediate local and long-range communication through the horizontal transfer of information, such as combinations of proteins, mRNAs and microRNAs. Here, we develop a tractable theoretical framework to study the putative role of exosome-mediated cell-cell communication in the cancer-immunity interplay. We reduce the complex interplay into a generic model whose three components are cancer cells, dendritic cells (consisting of precursor, immature, and mature types), and killer cells (consisting of cytotoxic T cells, helper T cells, effector B cells, and natural killer cells). The framework also incorporates the effects of exosome exchange on enhancement/reduction of cell maturation, proliferation, apoptosis, immune recognition, and activation/inhibition. We reveal tristability-possible existence of three cancer states: a low cancer load with intermediate immune level state, an intermediate cancer load with high immune level state, and a high cancer load with low immune-level state, and establish the corresponding effective landscape for the cancer-immunity network. We illustrate how the framework can contribute to the design and assessments of combination therapies.

Polarity mechanisms such as contact inhibition of locomotion regulate persistent rotational motion of mammalian cells on micropatterns.

Pairs of endothelial cells on adhesive micropatterns rotate persistently, but pairs of fibroblasts do not; coherent rotation is present in normal mammary acini and kidney cells but absent in cancerous cells. Why? To answer this question, we develop a computational model of pairs of mammalian cells on adhesive micropatterns using a phase field method and study the conditions under which persistent rotational motion (PRM) emerges. Our model couples the shape of the cell, the cell's internal chemical polarity, and interactions between cells such as volume exclusion and adhesion. We show that PRM can emerge from this minimal model and that the cell-cell interface may be influenced by the nucleus. We study the effect of various cell polarity mechanisms on rotational motion, including contact inhibition of locomotion, neighbor alignment, and velocity alignment, where cells align their polarity to their velocity. These polarity mechanisms strongly regulate PRM: Small differences in polarity mechanisms can create significant differences in collective rotation. We argue that the existence or absence of rotation under confinement may lead to insight into the cell's methods for coordinating collective cell motility.

Formation and Dynamics of Waves in a Cortical Model of Cholinergic Modulation.

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.

MicroRNA-based regulation of epithelial-hybrid-mesenchymal fate determination.

Forward and backward transitions between epithelial and mesenchymal phenotypes play crucial roles in embryonic development and tissue repair. Aberrantly regulated transitions are also a hallmark of cancer metastasis. The genetic network that regulates these transitions appears to allow for the existence of a hybrid phenotype (epithelial/mesenchymal). Hybrid cells are endowed with mixed epithelial and mesenchymal characteristics, enabling specialized capabilities such as collective cell migration. Cell-fate determination between the three phenotypes is in fact regulated by a circuit composed of two highly interconnected chimeric modules--the miR-34/SNAIL and the miR-200/ZEB mutual-inhibition feedback circuits. Here, we used detailed modeling of microRNA-based regulation to study this core unit. More specifically, we investigated the functions of the two isolated modules and subsequently of the combined unit when the two modules are integrated into the full regulatory circuit. We found that miR-200/ZEB forms a tristable circuit that acts as a ternary switch, driven by miR-34/SNAIL, that is a monostable module that acts as a noise-buffering integrator of internal and external signals. We propose to associate the three stable states--(1,0), (high miR-200)/(low ZEB); (0,1), (low miR-200)/(high ZEB); and (1/2,1/2), (medium miR-200)/(medium ZEB)--with the epithelial, mesenchymal, and hybrid phenotypes, respectively. Our (1/2,1/2) state hypothesis is consistent with recent experimental studies (e.g., ZEB expression measurements in collectively migrating cells) and explains the lack of observed mesenchymal-to-hybrid transitions in metastatic cells and in induced pluripotent stem cells. Testable predictions of dynamic gene expression during complete and partial transitions are presented.

Engineering Trojan-horse bacteria to fight cancer.

In this issue of Blood, Massa et al demonstrated that tumor-tropic bacteria engineered to express specific cancer-related antibodies can recognize and specifically destroy the malignant cells by infecting them with toxic cargo. The findings provide an essential step in establishing the novel approach of enlisting engineered bacteria that speak cancer's language to outsmart and defeat cancer

Antibiotic-induced anomalous statistics of collective bacterial swarming.

Under sublethal antibiotics concentrations, the statistics of collectively swarming Bacillus subtilis transitions from normal to anomalous, with a heavy-tailed speed distribution and a two-step temporal correlation of velocities. The transition is due to changes in the properties of the bacterial motion and the formation of a motility-defective subpopulation that self-segregates into regions. As a result, both the colonial expansion and the growth rate are not affected by antibiotics. This phenomenon suggests a new strategy bacteria employ to fight antibiotic stress.

Clique of functional hubs orchestrates population bursts in developmentally regulated neural networks.

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.

Propagating waves of directionality and coordination orchestrate collective cell migration.

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.

Glutamate mediated astrocytic filtering of neuronal activity.

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.