Share, but unequally: a plausible mechanism for emergence and maintenance of intratumour heterogeneity

A circuit for secretion-coupled cellular autonomy in multicellular eukaryotic cells

Cancers represent complex autonomous systems, displaying self-sufficiency in growth signaling. Autonomous growth is fueled by a cancer cell's ability to "secrete-and-sense" growth factors (GFs): a poorly understood phenomenon. Using an integrated computational and experimental approach, here we dissect the impact of a feedback-coupled GTPase circuit within the secretory pathway that imparts secretion-coupled autonomy. The circuit is assembled when the Ras-superfamily monomeric GTPase Arf1, and the heterotrimeric GTPase Giαβγ and their corresponding GAPs and GEFs are coupled by GIV/Girdin, a protein that is known to fuel aggressive traits in diverse cancers. One forward and two key negative feedback loops within the circuit create closed-loop control, allow the two GTPases to coregulate each other, and convert the expected switch-like behavior of Arf1-dependent secretion into an unexpected dose-response alignment behavior of sensing and secretion. Such behavior translates into cell survival that is self-sustained by stimulus-proportionate secretion. Proteomic studies and protein-protein interaction network analyses pinpoint GFs (e.g., the epidermal GF) as key stimuli for such self-sustenance. Findings highlight how the enhanced coupling of two biological switches in cancer cells is critical for multiscale feedback control to achieve secretion-coupled autonomy of growth factors.

Scales of Cancer Evolution: Selfish Genome or Cooperating Cells?

Simple Summary

Cancer continuously evolves its ability to survive in time-varying microenvironment, which results, regarding the therapeutic context, in its therapeutic resistance. As it is accepted that the development of resistance is the direct consequence of intratumour heterogeneity, its evolutionary etiology is intensively studied. Models of carinogenesis are often assessed accordingly to how well they fit into the evolutionary scenario. In the paper, the relevant observations and concepts in cancer research, such as intratumour heterogeneity, cell plasticity, and Markov cell state dynamics, are reviewed and integrated into an evolutionary model. The possibility that the interaction between cancer cells can be interpreted as cooperation is proposed.

Abstract

The exploitation of the evolutionary modus operandi of cancer to steer its progression towards drug sensitive cancer cells is a challenging research topic. Integrating evolutionary principles into cancer therapy requires properly identified selection level, the relevant timescale, and the respective fitness of the principal selection unit on that timescale. Interpretation of some features of cancer progression, such as increased heterogeneity of isogenic cancer cells, is difficult from the most straightforward evolutionary view with the cancer cell as the principal selection unit. In the paper, the relation between the two levels of intratumour heterogeneity, genetic, due to genetic instability, and non-genetic, due to phenotypic plasticity, is reviewed and the evolutionary role of the latter is outlined. In analogy to the evolutionary optimization in a changing environment, the cell state dynamics in cancer clones are interpreted as the risk diversifying strategy bet hedging, optimizing the balance between the exploitation and exploration of the cell state space.

A mathematical model for phenotypic heterogeneity in breast cancer with implications for therapeutic strategies

Inevitably, almost all cancer patients develop resistance to targeted therapy. Intratumour heterogeneity is a major cause of drug resistance. Mathematical models that explain experiments quantitatively are useful in understanding the origin of intratumour heterogeneity, which then could be used to explore scenarios for efficacious therapy. Here, we develop a mathematical model to investigate intratumour heterogeneity in breast cancer by exploiting the observation that HER2+ and HER2- cells could divide symmetrically or asymmetrically. Our predictions for the evolution of cell fractions are in quantitative agreement with single-cell experiments. Remarkably, the colony size of HER2+ cells emerging from a single HER2- cell (or vice versa), which occurs in about four cell doublings, also agrees with experimental results, without tweaking any parameter in the model. The theory explains experimental data on the responses of breast tumours under different treatment protocols. We then used the model to predict that, not only the order of two drugs, but also the treatment period for each drug and the tumour cell plasticity could be manipulated to improve the treatment efficacy. Mathematical models, when integrated with data on patients, make possible exploration of a broad range of parameters readily, which might provide insights in devising effective therapies.

Intratumor Heterogeneity in Hepatocellular Carcinoma: Challenges and Opportunities

Hepatocellular carcinoma (HCC) represents a leading cause of cancer-related death, but it remains difficult to treat. Intratumor genetic and phenotypic heterogeneity are inherent properties of breast, skin, lung, prostate, and brain tumors, and intratumor heterogeneity (ITH) helps define prognosis and therapeutic response in these cancers. Several recent studies estimate that ITH is inherent to HCC and attribute the clinical intractability of HCC to this heterogeneity. In this review, we examine the evidence for genomic, phenotypic, and tumor microenvironment ITH in HCC, with a focus on two of the top molecular drivers of HCC: β-catenin (CTNNB1) and Telomerase reverse transcriptase (TERT). We discuss the influence of ITH on HCC diagnosis, prognosis, and therapy, while highlighting the gaps in knowledge and possible future directions.

Self-generated persistent random forces drive phase separation in growing tumors

A single solid tumor, composed of nearly identical cells, exhibits heterogeneous dynamics. Dynamics of cells in the core is glass-like, whereas those in the periphery undergoes diffusive or super-diffusive behavior. Quantification of heterogeneity using the mean square displacement or the self-intermediate scattering function, which involves averaging over the cell population, hides the complexity of the collective movement. Using the t-distributed stochastic neighbor embedding (t-SNE), a popular unsupervised machine learning dimensionality reduction technique, we show that the phase space structure of an evolving colony of cells, driven by cell division and apoptosis, partitions into nearly disjoint sets composed principally of the core and periphery cells. The non-equilibrium phase separation is driven by the differences in the persistence of self-generated active forces induced by cell division. Extensive heterogeneity revealed by t-SNE paves the way toward understanding the origins of intratumor heterogeneity using experimental imaging data.

A Circuit for Secretion-coupled Cellular Autonomy in Multicellular Eukaryotes.. [DOI: 10.1101/2021.03.18.436048]

Cancers represent complex autonomous systems, displaying self-sufficiency in growth signaling. Autonomous growth is fueled by a cancer cell’s ability to ‘secrete-and-sense’ growth factors: a poorly understood phenomenon. Using an integrated systems and experimental approach, here we dissect the impact of a feedback-coupled GTPase circuit within the secretory pathway that imparts secretion-coupled autonomy. The circuit is assembled when the Ras-superfamily monomeric GTPase Arf1, and the heterotrimeric GTPase Giαβγ and their corresponding GAPs and GEFs are coupled by GIV/Girdin, a protein that is known to fuel aggressive traits in diverse cancers. One forward and two key negative feedback loops within the circuit create closed-loop control (CLC), allow the two GTPases to coregulate each other, and convert the expected switch-like behavior of Arf1-dependent secretion into an unexpected dose response alignment behavior of sensing and secretion. Such behavior translates into cell survival that is self-sustained by stimulus-proportionate secretion. Proteomic studies and protein-protein interaction network analyses pinpoint growth factors (e.g., the epidermal growth factor; EGF) as a key stimuli for such self-sustenance. Findings highlight how enhanced coupling of two biological switches in cancer cells is critical for multiscale feedback control to achieve secretion-coupled autonomy of growth factors.

SYNOPSIS IMAGE

STANDFIRST TEXT

This work defines the inner workings of a Golgi-localized molecular circuitry comprised of coupled GTPases, which empowers cells to achieve self-sufficiency in growth factor signaling by creating a secrete-and-sense autocrine loop.

HIGHLIGHTS/MAIN FINDINGS

Modeling and experimental approaches were used to dissect a coupled GTPase circuit.

Coupling enables closed loop feedback and mutual control of GTPases.

Coupling generates dose response alignment behavior of sensing and secretion of growth factors.

Coupling is critical for multiscale feedback control to achieve secretion-coupled autonomy.

Spatially heterogeneous dynamics of cells in a growing tumor spheroid: comparison between theory and experiments

Collective cell movement, characterized by multiple cells that are in contact for substantial periods of time and undergo correlated motion, plays a central role in cancer and embryogenesis. Recent imaging experiments have provided time-dependent traces of individual cells, thus providing an unprecedented picture of tumor spheroid growth. By using simulations of a minimal cell model, we analyze the experimental data that map the movement of cells in a fibrosarcoma tumor spheroid embedded in a collagen matrix. Both simulations and experiments show that cells in the core of the spheroid exhibit subdiffusive glassy dynamics (mean square displacement, Δ(t) ≈ tα with α < 1), whereas cells in the periphery exhibit superdiffusive motion, Δ(t) ≈ tα with α > 1. The motion of most of the cells near the periphery is highly persistent and correlated directional motion due to cell doubling and apoptosis rates, thus explaining the observed superdiffusive behavior. The α values for cells in the core and periphery, extracted from simulations and experiments, are in near quantitative agreement with each other, which is surprising given that no parameter in the model was used to fit the measurements. The qualitatively different dynamics of cells in the core and periphery is captured by the fourth order susceptibility, introduced to characterize metastable states in glass forming systems. Analyses of the velocity autocorrelation of individual cells show remarkable spatial heterogeneity with no two cells exhibiting similar behavior. The prediction that α should depend on the location of the cells in the tumor is amenable to experimental testing. The highly heterogeneous dynamics of cells in the tumor spheroid provides a plausible mechanism for the origin of intratumor heterogeneity.

Time scales and wave formation in non-linear spatial public goods games

The co-evolutionary dynamics of competing populations can be strongly affected by frequency-dependent selection and spatial population structure. As co-evolving populations grow into a spatial domain, their initial spatial arrangement and their growth rate differences are important factors that determine the long-term outcome. We here model producer and free-rider co-evolution in the context of a diffusive public good (PG) that is produced by the producers at a cost but evokes local concentration-dependent growth benefits to all. The benefit of the PG can be non-linearly dependent on public good concentration. We consider the spatial growth dynamics of producers and free-riders in one, two and three dimensions by modeling producer cell, free-rider cell and public good densities in space, driven by the processes of birth, death and diffusion (cell movement and public good distribution). Typically, one population goes extinct, but the time-scale of this process varies with initial conditions and the growth rate functions. We establish that spatial variation is transient regardless of dimensionality, and that structured initial conditions lead to increasing times to get close to an extinction state, called ε-extinction time. Further, we find that uncorrelated initial spatial structures do not influence this ε-extinction time in comparison to a corresponding well-mixed (non-spatial) system. In order to estimate the ε-extinction time of either free-riders or producers we derive a slow manifold solution. For invading populations, i.e. for populations that are initially highly segregated, we observe a traveling wave, whose speed can be calculated. Our results provide quantitative predictions for the transient spatial dynamics of cooperative traits under pressure of extinction.

Evolutionary public good (PG) games capture the essence of production of growth-beneficial factors that are vulnerable to exploitation by free-riders who do not carry the cost of production. PGs emerge in cellular populations, for example in growing bacteria and cancer cells. We study the eco-evolutionary dynamics of a PG in populations that grow in space. In our model, PG-producer cells and free-rider cells can grow according to their own birth and death rates. Co-evolution occurs due to public good-driven surplus in the intrinsic growth rates at a cost to producers. A net growth rate-benefit to free-riders leads to the well-known tragedy of the commons in which producers go extinct. What is often omitted from discussions is the time scale on which this extinction can occur, especially in spatial populations. Here, we derive analytical estimates of the ε-extinction time in different spatial settings. As we do not consider a stochastic process, the ε-extinction time captures the time needed to approach an extinction state. We identify spatial scenarios in which extinction takes long enough such that the tragedy of the commons never occurs within a meaningful lifetime of the system. Using numerical simulations we analyze the deviations from our analytical predictions.

Glioblastoma: Role of Mitochondria N-acetylserotonin/Melatonin Ratio in Mediating Effects of miR-451 and Aryl Hydrocarbon Receptor and in Coordinating Wider Biochemical Changes

A wide array of different factors and processes have been linked to the biochemical underpinnings of glioblastoma multiforme (GBM) and glioblastoma stem cells (GSC), with no clear framework in which these may be integrated. Consequently, treatment of GBM/GSC is generally regarded as very poor. This article provides a framework that is based on alterations in the regulation of the melatonergic pathways within mitochondria of GBM/GSC. It is proposed that the presence of high levels of mitochondria-synthesized melatonin is toxic to GBM/GSC, with a number of processes in GBM/GSC acting to limit melatonin’s synthesis in mitochondria. One such factor is the aryl hydrocarbon receptor, which increases cytochrome P450 (CYP)1b1 in mitochondria, leading to the ‘backward’ conversion of melatonin to N-acetylserotonin (NAS). N-acetylserotonin has some similar, but some important differential effects compared with melatonin, including its activation of the tyrosine receptor kinase B (TrkB) receptor. TrkB activation is important to GBM/GSC survival and proliferation. A plethora of significant, but previously disparate, data on GBM/GSC can then be integrated within this framework, including miR-451, AMP-activated protein kinase (AMPK)/mTOR, 14-3-3 proteins, sirtuins, tryptophan 2,3-dioxygenase, and the kynurenine pathways. Such a conceptualization provides a framework for the development of more effective treatment for this poorly managed condition.