The inherent fragility of collective proliferative control.

Tissues achieve and maintain their sizes through active feedback, whereby cells collectively regulate proliferation and differentiation so as to facilitate homeostasis and the ability to respond to disturbances. One of the best understood feedback mechanisms-renewal control-achieves remarkable feats of robustness in determining and maintaining desired sizes. Yet in a variety of biologically relevant situations, we show that stochastic effects should cause rare but catastrophic failures of renewal control. We define the circumstances under which this occurs and raise the possibility such events account for important non-genetic steps in the development of cancer. We further suggest that the spontaneous stochastic reversal of these events could explain cases of cancer normalization or dormancy following treatment. Indeed, we show that the kinetics of post-treatment recurrence for many cancers are often better fit by a model of stochastic re-emergence due to loss of collective proliferative control, than by deterministic models of cancer relapse.

Uncovering Minimal Pathways in Melanoma Initiation.

Cutaneous melanomas are clinically and histologically heterogeneous. Most display activating mutations in Braf or Nras and complete loss of function of one or more tumor suppressor genes. Mouse models that replicate such mutations produce fast-growing, pigmented tumors. However, mice that combine Braf activation with only heterozygous loss of Pten also produce tumors and, as we show here, in an Albino background this occurs even with Braf activation alone. Such tumors arise rarely, grow slowly, and express low levels of pigmentation genes. The timing of their appearance was consistent with a single step stochastic event, but no evidence could be found that it required de novo mutation, suggesting instead the involvement of an epigenetic transition. Single-cell transcriptomic analysis revealed such tumors to be heterogeneous, including a minor cell type we term LNM ( L ow-pigment, N eural- and extracellular M atrix-signature) that displays gene expression resembling "neural crest"-like cell subsets detected in the fast-growing tumors of more heavily-mutated mice, as well as in human biopsy and xenograft samples. We provide evidence that LNM cells pre-exist in normal skin, are expanded by Braf activation, can transition into malignant cells, and persist with malignant cells through multiple rounds of transplantation. We discuss the possibility that LNM cells not only serve as a pre-malignant state in the production of some melanomas, but also as an important intermediate in the development of drug resistance.

Predictive nonlinear modeling of malignant myelopoiesis and tyrosine kinase inhibitor therapy.

Chronic myeloid leukemia (CML) is a blood cancer characterized by dysregulated production of maturing myeloid cells driven by the product of the Philadelphia chromosome, the BCR-ABL1 tyrosine kinase. Tyrosine kinase inhibitors (TKIs) have proved effective in treating CML, but there is still a cohort of patients who do not respond to TKI therapy even in the absence of mutations in the BCR-ABL1 kinase domain that mediate drug resistance. To discover novel strategies to improve TKI therapy in CML, we developed a nonlinear mathematical model of CML hematopoiesis that incorporates feedback control and lineage branching. Cell-cell interactions were constrained using an automated model selection method together with previous observations and new in vivo data from a chimeric BCR-ABL1 transgenic mouse model of CML. The resulting quantitative model captures the dynamics of normal and CML cells at various stages of the disease and exhibits variable responses to TKI treatment, consistent with those of CML patients. The model predicts that an increase in the proportion of CML stem cells in the bone marrow would decrease the tendency of the disease to respond to TKI therapy, in concordance with clinical data and confirmed experimentally in mice. The model further suggests that, under our assumed similarities between normal and leukemic cells, a key predictor of refractory response to TKI treatment is an increased maximum probability of self-renewal of normal hematopoietic stem cells. We use these insights to develop a clinical prognostic criterion to predict the efficacy of TKI treatment and design strategies to improve treatment response. The model predicts that stimulating the differentiation of leukemic stem cells while applying TKI therapy can significantly improve treatment outcomes.

Modeling collective cell behavior in cancer: Perspectives from an interdisciplinary conversation.

Collective cell behavior contributes to all stages of cancer progression. Understanding how collective behavior emerges through cell-cell interactions and decision-making will advance our understanding of cancer biology and provide new therapeutic approaches. Here, we summarize an interdisciplinary discussion on multicellular behavior in cancer, draw lessons from other scientific disciplines, and identify future directions.

Leveraging gene correlations in single cell transcriptomic data.

BACKGROUND: Many approaches have been developed to overcome technical noise in single cell RNA-sequencing (scRNAseq). As researchers dig deeper into data-looking for rare cell types, subtleties of cell states, and details of gene regulatory networks-there is a growing need for algorithms with controllable accuracy and fewer ad hoc parameters and thresholds. Impeding this goal is the fact that an appropriate null distribution for scRNAseq cannot simply be extracted from data when ground truth about biological variation is unknown (i.e., usually). RESULTS: We approach this problem analytically, assuming that scRNAseq data reflect only cell heterogeneity (what we seek to characterize), transcriptional noise (temporal fluctuations randomly distributed across cells), and sampling error (i.e., Poisson noise). We analyze scRNAseq data without normalization-a step that skews distributions, particularly for sparse data-and calculate p-values associated with key statistics. We develop an improved method for selecting features for cell clustering and identifying gene-gene correlations, both positive and negative. Using simulated data, we show that this method, which we call BigSur (Basic Informatics and Gene Statistics from Unnormalized Reads), captures even weak yet significant correlation structures in scRNAseq data. Applying BigSur to data from a clonal human melanoma cell line, we identify thousands of correlations that, when clustered without supervision into gene communities, align with known cellular components and biological processes, and highlight potentially novel cell biological relationships. CONCLUSIONS: New insights into functionally relevant gene regulatory networks can be obtained using a statistically grounded approach to the identification of gene-gene correlations.

When oncogenes do not cause cancer.

Environmental cues, not oncogene-induced senescence, may stop melanocytes with an activating mutation in the BRAF gene from turning into melanoma.

Dynamics of nevus development implicate cell cooperation in the growth arrest of transformed melanocytes.

Mutational activation of the BRAF proto-oncogene in melanocytes reliably produces benign nevi (pigmented 'moles'), yet the same change is the most common driver mutation in melanoma. The reason nevi stop growing, and do not progress to melanoma, is widely attributed to a cell-autonomous process of 'oncogene-induced senescence'. Using a mouse model of Braf-driven nevus formation, analyzing both proliferative dynamics and single-cell gene expression, we found no evidence that nevus cells are senescent, either compared with other skin cells, or other melanocytes. We also found that nevus size distributions could not be fit by any simple cell-autonomous model of growth arrest, yet were easily fit by models based on collective cell behavior, for example in which arresting cells release an arrest-promoting factor. We suggest that nevus growth arrest is more likely related to the cell interactions that mediate size control in normal tissues, than to any cell-autonomous, 'oncogene-induced' program of senescence.

Melanocytes are pigment-producing cells found throughout the skin. Mutations that activate a gene called BRAF cause these cells to divide and produce melanocytic nevi, also known as "moles". These mutations are oncogenic, meaning they can cause cancer. Indeed, BRAF is the most commonly mutated gene in melanoma, a deadly skin cancer that arises from melanocytes. Yet, moles hardly ever progress to melanoma. A proposed explanation for this behavior is that, once activated, BRAF initiates a process called "oncogene-induced senescence" in each melanocyte. This process, likened to premature aging, is thought to be what causes cells in a mole to quit dividing. Although this hypothesis is widely accepted, it has proved difficult to test directly. To investigate this notion, Ruiz-Vega et al. studied mice with hundreds of moles created by the same BRAF mutation found in human moles. Analyzing the activity of genes in individual cells revealed that nevus melanocytes that have stopped growing are no more senescent than other skin cells, including non-mole melanocytes. Ruiz-Vega et al. then analyzed the sizes at which moles stopped growing, estimating the number of cells in each mole. The data were then compared with the results of a simulation and mathematical modeling. This revealed that any model based on the idea of cells independently shutting down after a number of random events could not reproduce the distribution of mole sizes that had been experimentally observed. On the other hand, models based on melanocytes acting collectively to shut down each other's growth fit the observed data much better. These findings suggest that moles do not stop growing as a direct result of the activation of BRAF, but because they sense and respond to their own overgrowth. The same kind of collective sensing is observed in normal tissues that maintain a constant size. Discovering that melanocytes do this not only sheds light on why moles stop growing, it could also help researchers devise new ways to prevent melanomas from forming.

Expanding signaling-molecule wavefront model of cell polarization in the Drosophila wing primordium.

In developing tissues, cell polarization and proliferation are regulated by morphogens and signaling pathways. Cells throughout the Drosophila wing primordium typically show subcellular localization of the unconventional myosin Dachs on the distal side of cells (nearest the center of the disc). Dachs localization depends on the spatial distribution of bonds between the protocadherins Fat (Ft) and Dachsous (Ds), which form heterodimers between adjacent cells; and the Golgi kinase Four-jointed (Fj), which affects the binding affinities of Ft and Ds. The Fj concentration forms a linear gradient while the Ds concentration is roughly uniform throughout most of the wing pouch with a steep transition region that propagates from the center to the edge of the pouch during the third larval instar. Although the Fj gradient is an important cue for polarization, it is unclear how the polarization is affected by cell division and the expanding Ds transition region, both of which can alter the distribution of Ft-Ds heterodimers around the cell periphery. We have developed a computational model to address these questions. In our model, the binding affinity of Ft and Ds depends on phosphorylation by Fj. We assume that the asymmetry of the Ft-Ds bond distribution around the cell periphery defines the polarization, with greater asymmetry promoting cell proliferation. Our model predicts that this asymmetry is greatest in the radially-expanding transition region that leaves polarized cells in its wake. These cells naturally retain their bond distribution asymmetry after division by rapidly replenishing Ft-Ds bonds at new cell-cell interfaces. Thus we predict that the distal localization of Dachs in cells throughout the pouch requires the movement of the Ds transition region and the simple presence, rather than any specific spatial pattern, of Fj.

Organ-level quorum sensing directs regeneration in hair stem cell populations.

Coordinated organ behavior is crucial for an effective response to environmental stimuli. By studying regeneration of hair follicles in response to patterned hair plucking, we demonstrate that organ-level quorum sensing allows coordinated responses to skin injury. Plucking hair at different densities leads to a regeneration of up to five times more neighboring, unplucked resting hairs, indicating activation of a collective decision-making process. Through data modeling, the range of the quorum signal was estimated to be on the order of 1 mm, greater than expected for a diffusible molecular cue. Molecular and genetic analysis uncovered a two-step mechanism, where release of CCL2 from injured hairs leads to recruitment of TNF-α-secreting macrophages, which accumulate and signal to both plucked and unplucked follicles. By coupling immune response with regeneration, this mechanism allows skin to respond predictively to distress, disregarding mild injury, while meeting stronger injury with full-scale cooperative activation of stem cells.

The protective role of symmetric stem cell division on the accumulation of heritable damage.

Stem cell divisions are either asymmetric-in which one daughter cell remains a stem cell and one does not-or symmetric, in which both daughter cells adopt the same fate, either stem or non-stem. Recent studies show that in many tissues operating under homeostatic conditions stem cell division patterns are strongly biased toward the symmetric outcome, raising the question of whether symmetry confers some benefit. Here, we show that symmetry, via extinction of damaged stem-cell clones, reduces the lifetime risk of accumulating phenotypically silent heritable damage (mutations or aberrant epigenetic changes) in individual stem cells. This effect is greatest in rapidly cycling tissues subject to accelerating rates of damage accumulation over time, a scenario that describes the progression of many cancers. A decrease in the rate of cellular damage accumulation may be an important factor favoring symmetric patterns of stem cell division.

Predicting mechanism of biphasic growth factor action on tumor growth using a multi-species model with feedback control.

A large number of growth factors and drugs are known to act in a biphasic manner: at lower concentrations they cause increased division of target cells, whereas at higher concentrations the mitogenic effect is inhibited. Often, the molecular details of the mitogenic effect of the growth factor are known, whereas the inhibitory effect is not. Hepatoctyte Growth Factor, HGF, has recently been recognized as a strong mitogen that is present in the microenvironment of solid tumors. Recent evidence suggests that HGF acts in a biphasic manner on tumor growth. We build a multi-species model of HGF action on tumor cells using different hypotheses for high dose-HGF activation of a growth inhibitor and show that the shape of the dose-response curve is directly related to the mechanism of inhibitor activation. We thus hypothesize that the shape of a dose-response curve is informative of the molecular action of the growth factor on the growth inhibitor.

Dynamics and precision in retinoic acid morphogen gradients.

Retinoic acid (RA) regulates many cellular behaviors during embryonic development and adult homeostasis. Like other morphogens, RA forms gradients through the use of localized sources and sinks, feedback, and interactions with other signals; this has been particularly well studied in the context of hindbrain segmentation in vertebrate embryos. Yet, as a small lipophilic molecule derived from a dietary source-vitamin A-RA differs markedly from better-studied polypeptide morphogens in its mechanisms of transport, signaling, and removal. Computational models suggest that the distinctive features of RA gradients make them particularly robust to large perturbations. Such features include combined positive and negative feedback effects via intracellular fatty acid binding proteins and RA-degrading enzymes. Here, we discuss how these features, together with feedback interactions among RA target genes, help enable RA to specify multiple, accurate pattern elements in the developing hindbrain, despite operating in an environment of high cellular and biochemical uncertainty and noise.

Activin and GDF11 collaborate in feedback control of neuroepithelial stem cell proliferation and fate.

Studies of the olfactory epithelium model system have demonstrated that production of neurons is regulated by negative feedback. Previously, we showed that a locally produced signal, the TGFß superfamily ligand GDF11, regulates the genesis of olfactory receptor neurons by inhibiting proliferation of the immediate neuronal precursors (INPs) that give rise to them. GDF11 is antagonized by follistatin (FST), which is also produced locally. Here, we show that Fst(-/-) mice exhibit dramatically decreased neurogenesis, a phenotype that can only be partially explained by increased GDF11 activity. Instead, a second FST-binding factor, activin ßB (ACTßB), inhibits neurogenesis by a distinct mechanism: whereas GDF11 inhibits expansion of INPs, ACTßB inhibits expansion of stem and early progenitor cells. We present data supporting the concept that these latter cells, previously considered two distinct types, constitute a dynamic stem/progenitor population in which individual cells alternate expression of Sox2 and/or Ascl1. In addition, we demonstrate that interplay between ACTßB and GDF11 determines whether stem/progenitor cells adopt a glial versus neuronal fate. Altogether, the data indicate that the transition between stem cells and committed progenitors is neither sharp nor irreversible and that GDF11, ACTßB and FST are crucial components of a circuit that controls both total cell number and the ratio of neuronal versus glial cells in this system. Thus, our findings demonstrate a close connection between the signals involved in the control of tissue size and those that regulate the proportions of different cell types.

Heparan sulfate acts as a bone morphogenetic protein coreceptor by facilitating ligand-induced receptor hetero-oligomerization.

Cell surface heparan sulfate (HS) not only binds several major classes of growth factors but also sometimes potentiates their activities--an effect usually termed "coreception." A view that coreception is due to the stabilization of growth factor-receptor interactions has emerged primarily from studies of the fibroblast growth factors (FGFs). Recent in vivo studies have strongly suggested that HS also plays an important role in regulating signaling by the bone morphogenetic proteins (BMPs). Here, we provide evidence that the mechanism of coreception for BMPs is markedly different from that established for FGFs. First, we demonstrate a direct, stimulatory role for cell surface HS in the immediate signaling activities of BMP2 and BMP4, and we provide evidence that HS-BMP interactions are required for this effect. Next, using several independent assays of ligand binding and receptor assembly, including coimmunoprecipitation, cross-linking, and fluorescence fluctuation microscopy, we show that HS does not affect BMP binding to type I receptor subunits but instead enhances the subsequent recruitment of type II receptor subunits to BMP-type I receptor complexes. This suggests a view of HS as a catalyst of the formation of signaling complexes, rather than as a stabilizer of growth factor binding.

The 'stem cell' concept: is it holding us back?

Developmental biology, regenerative medicine and cancer biology are increasingly occupied with the molecular characterization of stem cells. Yet recent work adds to a growing body of literature suggesting that 'stemness' cannot be reduced to the molecular features of cell types, and is instead an emergent property of cell lineages under feedback control.

Discovery of a novel molecule that regulates tumor growth and metastasis.

The heparan sulfate proteoglycan, Glypican-1 (GPC1), significantly impacts the growth of pancreatic cancer cells in vivo and markedly attenuates tumor angiogenesis and metastasis in athymic mice. Interestingly, both cancer cell-derived and host-derived GPC1 play an important role in tumor development and spread. These data suggest that GPC1 may be a valid therapeutic target for pancreatic cancer.

Glypican-1 modulates the angiogenic and metastatic potential of human and mouse cancer cells.

Cells isolated from many types of human cancers express heparin-binding growth factors (HBGFs) that drive tumor growth, metastasis, and angiogenesis. The heparan sulfate proteoglycan glypican-1 (GPC1) is a coreceptor for HBGFs. Here we show that both cancer cell-derived and host-derived GPC1 are crucial for efficient growth, metastasis, and angiogenesis of human and mouse cancer cells. Thus downregulation of GPC1 in the human pancreatic cancer cell line PANC-1, using antisense approaches, resulted in prolonged doubling times and decreased anchorage-independent growth in vitro as well as attenuated tumor growth, angiogenesis, and metastasis when these cells were transplanted into athymic mice. Moreover, athymic mice that lacked GPC1 exhibited decreased tumor angiogenesis and metastasis following intrapancreatic implantation with either PANC-1 or T3M4 human pancreatic cancer cells and fewer pulmonary metastases following intravenous injection of murine B16-F10 melanoma cells. In addition, hepatic endothelial cells isolated from these mice exhibited an attenuated mitogenic response to VEGF-A. These data indicate that cancer cell- and host-derived GPC1 are crucial for full mitogenic, angiogenic, and metastatic potential of cancer cells. Thus targeting GPC1 might provide new avenues for cancer therapy and for the prevention of cancer metastasis.

Growth factor-induced shedding of syndecan-1 confers glypican-1 dependence on mitogenic responses of cancer cells.

The cell surface heparan sulfate proteoglycan (HSPG) glypican-1 is up-regulated by pancreatic and breast cancer cells, and its removal renders such cells insensitive to many growth factors. We sought to explain why the cell surface HSPG syndecan-1, which is also up-regulated by these cells and is a known growth factor coreceptor, does not compensate for glypican-1 loss. We show that the initial responses of these cells to the growth factor FGF2 are not glypican dependent, but they become so over time as FGF2 induces shedding of syndecan-1. Manipulations that retain syndecan-1 on the cell surface make long-term FGF2 responses glypican independent, whereas those that trigger syndecan-1 shedding make initial FGF2 responses glypican dependent. We further show that syndecan-1 shedding is mediated by matrix metalloproteinase-7 (MMP7), which, being anchored to cells by HSPGs, also causes its own release in a complex with syndecan-1 ectodomains. These results support a specific role for shed syndecan-1 or MMP7-syndecan-1 complexes in tumor progression and add to accumulating evidence that syndecans and glypicans have nonequivalent functions in vivo.

Widespread defects in the primary olfactory pathway caused by loss of Mash1 function.

MASH1, a basic helix-loop-helix transcription factor, is widely expressed by neuronal progenitors in the CNS and PNS, suggesting that it plays a role in the development of many neural regions. However, in mice lacking a functional Mash1 gene, major alterations have been reported in only a few neuronal populations; among these is a generalized loss of olfactory receptor neurons of the olfactory epithelium. Here, we use a transgenic reporter mouse line, in which the cell bodies and growing axons of subsets of central and peripheral neurons are marked by expression of a tau-lacZ reporter gene (the Tattler-4 allele), to look both more broadly and deeply at defects in the nervous system of Mash1-/- mice. In addition to the expected lack of olfactory receptor neurons in the main olfactory epithelium, developing Mash1-/-;Tattler-4+/- mice exhibited reductions in neuronal cell number in the vomeronasal organ and in the olfactory bulb; the morphology of the rostral migratory stream, which gives rise to olfactory bulb interneurons, was also abnormal. Further examination of cell proliferation, cell death, and cell type-specific markers in Mash1-/- animals uncovered parallels between the main olfactory epithelium and the vomeronasal organ in the regulation of sensory neuron development. Interestingly, this analysis also revealed that, in the olfactory epithelium of Mash1-/- animals, there is an overproduction of proliferating cells that co-express markers of both neuronal progenitors and supporting cells. This finding suggests that olfactory receptor neurons and olfactory epithelium supporting cells may share a common progenitor, and that expression of Mash1 may be an important factor in determining whether these progenitors ultimately generate neurons or glia.