Search

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., Competing Interests: The authors declare no conflict of interest.

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., Competing Interests: The authors declare no conflict of interest.

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.

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.

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.

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.

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.

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.

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.

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.

In the heterogeneous environment surrounding plant roots (the rhizosphere), microorganisms both compete and cooperate. Here, we show that two very different inhabitants of the rhizosphere, the nonmotile fungus Aspergillus fumigatus and the swarming bacterium Paenibacillus vortex, can facilitate each other's dispersal. A. fumigatus conidia (nonmotile asexual fungal spores) can be transported by P. vortex swarms over distances of at least 30 cm and at rates of up to 10.8 mm h(-1). Moreover, conidia can be rescued and transported by P. vortex from niches of adverse growth conditions. Potential benefit to the bacteria may be in crossing otherwise impenetrable barriers in the soil: fungal mycelia seem to act as bridges to allow P. vortex to cross air gaps in agar plates. Transport of conidia was inhibited by proteolytic treatment of conidia or the addition of purified P. vortex flagella, suggesting specific contacts between flagella and proteins on the conidial surface. Conidia were transported by P. vortex into locations where antibiotics inhibited bacteria growth, and therefore, growth and sporulation of A. fumigatus were not limited by bacterial competition. Conidia from other fungi, similar in size to those fungi from A. fumigatus, were not transported as efficiently by P. vortex. Conidia from a range of fungi were not transported by another closely related rhizosphere bacterium, Paenibacillus polymyxa, or the more distantly related Proteus mirabilis, despite both being efficient swarmers.

Sibling Paenibacillus dendritiformis bacterial colonies grown on low-nutrient agar medium mutually inhibit growth through secretion of a lethal factor. Analysis of secretions reveals the presence of subtilisin (a protease) and a 12 kDa protein, termed sibling lethal factor (Slf). Purified subtilisin promotes the growth and expansion of P. dendritiformis colonies, whereas Slf is lethal and lyses P. dendritiformis cells in culture. Slf is encoded by a gene belonging to a large family of bacterial genes of unknown function, and the gene is predicted to encode a protein of approximately 20 kDa, termed dendritiformis sibling bacteriocin. The 20 kDa recombinant protein was produced and found to be inactive, but exposure to subtilisin resulted in cleavage to the active, 12 kDa form. The experimental results, combined with mathematical modeling, show that subtilisin serves to regulate growth of the colony. Below a threshold concentration, subtilisin promotes colony growth and expansion. However, once it exceeds a threshold, as occurs at the interface between competing colonies, Slf is then secreted into the medium to rapidly reduce cell density by lysis of the bacterial cells. The presence of genes encoding homologs of dendritiformis sibling bacteriocin in other bacterial species suggests that this mechanism for self-regulation of colony growth might not be limited to P. dendritiformis.

Bacteria serve as the central arena for understanding how gene networks and proteins process information and control cellular behaviors. Recently, much effort has been devoted to the investigation of specific bacteria gene circuits as functioning modules. The next challenge is the integrative modeling of complex cellular networks composed of many such modules. A tractable integrative model of the sophisticated decision-making signal transduction system that determines the fate between sporulation and competence is presented. This model provides an understanding of how information is sensed and processed to reach an "informative" decision in the context of cell state and signals from other cells. The competence module (ComK dynamics) is modeled as a stochastic switch whose transition rate is controlled by a quorum-sensing unit. The sporulation module (Spo0A dynamics) is modeled as a timer whose clock rate is adjusted by a stress-sensing unit. The interplay between these modules is mediated via the Rap assessment system, which gates the sensing units, and the AbrB-Rok decision module, which creates an opportunity for competence within a specific window of the sporulation timer. The timer is regulated via a special repressilator-like inhibition of Spo0A* by Spo0E, which is itself inhibited by AbrB. For some stress and input signals, this repressilator can generate a frustration state with large variations (fluctuations or oscillations) in Spo0A* and AbrB concentrations, which might serve an important role in generating cell variability. This integrative framework is a starting point that can be extended to include transition into cannibalism and the role of colony organization.

The immune system is essential to body defense and maintenance. Specific antibodies to foreign invaders function in body defense, and it has been suggested that autoantibodies binding to self molecules are important in body maintenance. Recently, the autoantibody repertoires in the bloods of healthy mothers and their newborns were studied using an antigen microarray containing hundreds of self molecules. It was found that the mothers expressed diverse repertoires for both IgG and IgM autoantibodies. Each newborn shares with its mother a similar repertoire of IgG antibodies, which cross the placental but its IgM repertoire is more similar to those of other newborns. Here, we took a system-level approach and analyzed the correlations between autoantibody reactivities of the previous data and extended the study to new data from newborns at birth and a week later, and from healthy young women. For the young women, we found modular organization of both IgG and IgM isotypes into antigen cliques-subgroups of highly correlated antigen reactivities. In contrast, the newborns were found to share a universal congenital IgM profile with no modular organization. Moreover, the IgG autoantibodies of the newborns manifested buds of the mothers' antigen cliques, but they were noticeably less structured. These findings suggest that the natural autoantibody repertoire of humans shows relatively little organization at birth, but, by young adulthood, it becomes sorted out into a modular organization of subgroups (cliques) of correlated antigens. These features revealed by antigen microarrays can be used to define personal states of autoantibody organizational motifs.

Bacteria can secrete a wide array of antibacterial compounds when competing with other bacteria for the same resources. Some of these compounds, such as bacteriocins, can affect bacteria of similar or closely related strains. In some cases, these secretions have been found to kill sibling cells that belong to the same colony. Here, we present experimental observations of competition between 2 sibling colonies of Paenibacillus dendritiformis grown on a low-nutrient agar gel. We find that neighboring colonies (growing from droplet inoculation) mutually inhibit growth through secretions that become lethal if the level exceeds a well-defined threshold. In contrast, within a single colony developing from a droplet inoculation, no growth inhibition is observed. However, growth inhibition and cell death are observed if material extracted from the agar between 2 growing colonies is introduced outside a growing single colony. To interpret the observations, we devised a simple mathematical model for the secretion of an antibacterial compound. Simulations of this model illustrate how secretions from neighboring colonies can be deadly, whereas secretions from a single colony growing from a droplet are not.

The role of stochasticity and noise in controlling genetic circuits is investigated in the context of transitions into and from competence in Bacillus subtilis. Recent experiments have demonstrated that bistability is not necessary for this function, but that the existence of one stable fixed point (vegetation) and an excitable unstable one (competence) is sufficient. Stochasticity therefore plays a crucial role in this excitation. Noise can be generated by discrete events such as RNA and protein synthesis and their degradation. We consider an alternative noise source connected with the protein binding/unbinding to the DNA. A theoretical model that includes this "nonadiabatic" mechanism appears to produce a better agreement with experiments than models where only the adiabatic limit is considered, suggesting that this nonconventional stochasticity source may be important for biological functions.

The mitotic spindle checkpoint halts the cell cycle until all chromosomes are attached to the mitotic spindles. Evidence suggests that the checkpoint prevents cell-cycle progression by inhibiting the activity of the APC-Cdc20 complex, but the precise mechanism underlying this inhibition is not yet known. Here, we use mathematical modeling to compare several mechanisms that could account for this inhibition. We describe the interplay between the capacities to strongly inhibit cell-cycle progression before spindle attachment on one hand and to rapidly resume cell-cycle progression once the last kinetochore is attached on the other hand. We find that inhibition that is restricted to the kinetochore region is not sufficient for supporting both requirements when realistic diffusion constants are considered. A mechanism that amplifies the checkpoint signal through autocatalyzed inhibition is also insufficient. In contrast, amplifying the signal through the release of a diffusible inhibitory complex can support reliable checkpoint function. Our results suggest that the design of the spindle checkpoint network is limited by physical constraints imposed by realistic diffusion constants and the relevant spatial and temporal dimensions where computation is performed.