Real-time cell cycle imaging during melanoma growth, invasion, and drug response

Patient-derived melanoma models

Melanoma is a very aggressive, rapidly metastasizing tumor that has been studied intensively in the past regarding the underlying genetic and molecular mechanisms. More recently developed treatment modalities have improved response rates and overall survival of patients. However, the majority of patients suffer from secondary treatment resistance, which requires in depth analyses of the underlying mechanisms. Here, melanoma models based on patients-derived material may play an important role. Consequently, a plethora of different experimental techniques have been developed in the past years. Among these are 3D and 4D culture techniques, organotypic skin reconstructs, melanoma-on-chip models and patient-derived xenografts, Every technique has its own strengths but also weaknesses regarding throughput, reproducibility, and reflection of the human situation. Here, we provide a comprehensive overview of currently used techniques and discuss their use in different experimental settings.

Characterization of two different alginate-based bioinks and the influence of melanoma growth within

Extrusion-based bioprinting is an established method in biofabrication. Suitable bioinks have fundamentally different compositions and characteristics, which should be examined, in order to find a perfect model system. Here, we investigate the effect of two alginate-based, yet unalike 3D-printed bioinks, pre-crosslinked alginate-dialdehyde gelatin (ADA-GEL) and a mixture of alginate, hyaluronic acid, and gelatin (Alg/HA/Gel), on the melanoma cell line Mel Im and vice versa in terms of stiffness, shrinkage, cellular behavior and colony formation over 15 days. Rheological stiffness measurements revealed two soft gels with similar storage moduli. The cells did not have a significant impact on the overall stiffness, whereas ADA-GEL (2.5/2.5%) was significantly stiffer than Alg/HA/Gel (0.5/0.1/3%). Regarding the shrinkage of printed constructs, cells had a significant influence, especially in ADA-GEL, which has covalent bonds between the oxidized alginate and gelatin. Multi-photon microscopy exhibited proliferation, cell spreading and migration in ADA-GEL with cell-cell and cell-matrix interaction, dissimilarly to Alg/HA/Gel, in which cells formed spherical, encapsulated colonies. Scanning electron microscopy and histology showed degradation and multi-layered growth on ADA-GEL and fewer examples of escaped cells on Alg/HA/Gel. Both gels serve as proliferation bioink for melanoma with more necrosis in deeper Alg/HA/Gel colonies and differences in spreading and matrix interaction. These findings show the importance of proper characterization of the bioinks for different applications.

Cell cycle perturbation uncouples mitotic progression and invasive behavior in a post-mitotic cell

The acquisition of the post-mitotic state is crucial for the execution of many terminally differentiated cell behaviors during organismal development. However, the mechanisms that maintain the post-mitotic state in this context remain poorly understood. To gain insight into these mechanisms, we used the genetically and visually accessible model of C. elegans anchor cell (AC) invasion into the vulval epithelium. The AC is a terminally differentiated uterine cell that normally exits the cell cycle and enters a post-mitotic state before initiating contact between the uterus and vulva through a cell invasion event. Here, we set out to identify the set of negative cell cycle regulators that maintain the AC in this post-mitotic, invasive state. Our findings revealed a critical role for CKI-1 (p21CIP1/p27KIP1) in redundantly maintaining the post-mitotic state of the AC, as loss of CKI-1 in combination with other negative cell cycle regulators-including CKI-2 (p21CIP1/p27KIP1), LIN-35 (pRb/p107/p130), FZR-1 (Cdh1/Hct1), and LIN-23 (β-TrCP)-resulted in proliferating ACs. Remarkably, time-lapse imaging revealed that these ACs retain their ability to invade. Upon examination of a node in the gene regulatory network controlling AC invasion, we determined that proliferating, invasive ACs do so by maintaining aspects of pro-invasive gene expression. We therefore report that the requirement for a post-mitotic state for invasive cell behavior can be bypassed following direct cell cycle perturbation.

Generation of 3D melanoma models using an assembloid-based approach

While 3D tumor models have greatly evolved over the past years, there is still a strong requirement for more biosimilar models which are capable of recapitulating cellular crosstalk within the tumor microenvironment while equally displaying representative levels of tumor aggressiveness and invasion. Herein, we disclose an assembloid melanoma model based on the fusion of individual stromal multicellular spheroids (MCSs). In contrast to more traditional tumor models, we show that it is possible to develop self-organizing, heterotypic melanoma models where tumor cells present stem-cell like features like up-regulated pluripotency master regulators SOX2, POU5F1 and NANOG. Additionally, these assembloids display high levels of invasiveness while embedded in 3D matrices as evidenced by stromal cell promotion of melanoma cell invasion via metalloproteinase production. Furthermore, sensitivity to anticancer drug doxorubicin was demonstrated for the melanoma assembloid model. These findings suggest that melanoma assembloids may play a significant role in the field of 3D cancer models as they more closely mimic the tumor microenvironment when compared to more traditional MCSs, opening the doors to a better understanding of the role of tumor microenvironment in supporting tumor progression. STATEMENT OF SIGNIFICANCE: The development of complex 3D tumor models that better recapitulate the tumor microenvironment is crucial for both an improved comprehension of intercellular crosstalk and for more efficient drug screening. We have herein developed a self-organizing heterotypic assembloid-based melanoma model capable of closely mimicking the tumor microenvironment. Key features recapitulated were the preservation of cancer cell stemness, sensitivity to anti-cancer agents and tumor cell invasion promoted by stromal cells. The approach of pre-establishing distinct stromal domains for subsequent combination into more complex tumor constructs provides a route for developing superior tumor models with a higher degree of similarity to native cancer tissues.

Dormancy of cutaneous melanoma

Many cancer-related deaths including melanoma result from metastases that develop months or years after the initial cancer therapy. Even the most effective drugs and immune therapies rarely eradicate all tumor cells. Instead, they strongly reduce cancer burden, permitting dormant cancer cells to persist in niches, where they establish a cellular homeostasis with their host without causing clinical symptoms. Dormant cancers respond poorly to most drugs and therapies since they do not proliferate and hide in niches. It therefore remains a major challenge to develop novel therapies for dormant cancers. In this review we focus on the mechanisms regulating the initiation of cutaneous melanoma dormancy as well as those which are involved in reawakening of dormant cutaneous melanoma cells. In recent years the role of neutrophils and niche components in reawakening of melanoma cells came into focus and indicate possible future therapeutic applications. Sophisticated in vitro and in vivo melanoma dormancy models are needed to make progress in this field and are discussed.

Cell cycle perturbation uncouples mitotic progression and invasive behavior in a post-mitotic cell

The acquisition of the post-mitotic state is crucial for the execution of many terminally differentiated cell behaviors during organismal development. However, the mechanisms that maintain the post-mitotic state in this context remain poorly understood. To gain insight into these mechanisms, we used the genetically and visually accessible model of C. elegans anchor cell (AC) invasion into the vulval epithelium. The AC is a terminally differentiated uterine cell that normally exits the cell cycle and enters a post-mitotic state, initiating contact between the uterus and vulva through a cell invasion event. Here, we set out to identify the set of negative cell cycle regulators that maintain the AC in this post-mitotic, invasive state. Our findings revealed a critical role for CKI-1 (p21CIP1/p27KIP1) in redundantly maintaining the post-mitotic state of the AC, as loss of CKI-1 in combination with other negative cell cycle regulators-including CKI-2 (p21CIP1/p27KIP1), LIN-35 (pRb/p107/p130), FZR-1 (Cdh1/Hct1), and LIN-23 (β-TrCP)-resulted in proliferating ACs. Remarkably, time-lapse imaging revealed that these ACs retain their ability to invade. Upon examination of a node in the gene regulatory network controlling AC invasion, we determined that proliferating, invasive ACs do so by maintaining aspects of pro-invasive gene expression. We therefore report that the requirement for a post-mitotic state for invasive cell behavior can be bypassed following direct cell cycle perturbation.

Real-Time Cell Cycle Imaging in a 3D Cell Culture Model of Melanoma, Quantitative Analysis, Optical Clearing, and Mathematical Modeling

Aberrant cell cycle progression is a hallmark of solid tumors. Therefore, cell cycle analysis is an invaluable technique to study cancer cell biology. However, cell cycle progression has been most commonly assessed by methods that are limited to temporal snapshots or that lack spatial information. In this chapter, we describe a technique that allows spatiotemporal real-time tracking of cell cycle progression of individual cells in a multicellular context. The power of this system lies in the use of 3D melanoma spheroids generated from melanoma cells engineered with the fluorescent ubiquitination-based cell cycle indicator (FUCCI). This technique, combined with mathematical modeling, allows us to gain further and more detailed insight into several relevant aspects of solid cancer cell biology, such as tumor growth, proliferation, invasion, and drug sensitivity.

Formation and Growth of Co-Culture Tumour Spheroids: New Compartment-Based Mathematical Models and Experiments

Co-culture tumour spheroid experiments are routinely performed to investigate cancer progression and test anti-cancer therapies. Therefore, methods to quantitatively characterise and interpret co-culture spheroid growth are of great interest. However, co-culture spheroid growth is complex. Multiple biological processes occur on overlapping timescales and different cell types within the spheroid may have different characteristics, such as differing proliferation rates or responses to nutrient availability. At present there is no standard, widely-accepted mathematical model of such complex spatio-temporal growth processes. Typical approaches to analyse these experiments focus on the late-time temporal evolution of spheroid size and overlook early-time spheroid formation, spheroid structure and geometry. Here, using a range of ordinary differential equation-based mathematical models and parameter estimation, we interpret new co-culture experimental data. We provide new biological insights about spheroid formation, growth, and structure. As part of this analysis we connect Greenspan's seminal mathematical model to co-culture data for the first time. Furthermore, we generalise a class of compartment-based spheroid mathematical models that have previously been restricted to one population so they can be applied to multiple populations. As special cases of the general model, we explore multiple natural two population extensions to Greenspan's seminal model and reveal biological mechanisms that can describe the internal dynamics of growing co-culture spheroids and those that cannot. This mathematical and statistical modelling-based framework is well-suited to analyse spheroids grown with multiple different cell types and the new class of mathematical models provide opportunities for further mathematical and biological insights.

H3K4me3 remodeling induced acquired resistance through O-GlcNAc transferase

AIMS

Drivers of the drug tolerant proliferative persister (DTPP) state have not been well investigated. Histone H3 lysine-4 trimethylation (H3K4me3), an active histone mark, might enable slow cycling drug tolerant persisters (DTP) to regain proliferative capacity. This study aimed to determine H3K4me3 transcriptionally active sites identifying a key regulator of DTPPs.

METHODS

Deploying a model of adaptive cancer drug tolerance, H3K4me3 ChIP-Seq data of DTPPs guided identification of top transcription factor binding motifs. These suggested involvement of O-linked N-acetylglucosamine transferase (OGT), which was confirmed by metabolomics analysis and biochemical assays. OGT impact on DTPPs and adaptive resistance was explored in vitro and in vivo.

RESULTS

H3K4me3 remodeling was widespread in CPG island regions and DNA binding motifs associated with O-GlcNAc marked chromatin. Accordingly, we observed an upregulation of OGT, O-GlcNAc and its binding partner TET1 in chronically treated cancer cells. Inhibition of OGT led to loss of H3K4me3 and downregulation of genes contributing to drug resistance. Genetic ablation of OGT prevented acquired drug resistance in in vivo models. Upstream of OGT, we identified AMPK as an actionable target. AMPK activation by acetyl salicylic acid downregulated OGT with similar effects on delaying acquired resistance.

CONCLUSION

Our findings uncover a fundamental mechanism of adaptive drug resistance that governs cancer cell reprogramming towards acquired drug resistance, a process that can be exploited to improve response duration and patient outcomes.

Fluorescent biosensors illuminate the spatial regulation of cell signaling across scales

As cell signaling research has advanced, it has become clearer that signal transduction has complex spatiotemporal regulation that goes beyond foundational linear transduction models. Several technologies have enabled these discoveries, including fluorescent biosensors designed to report live biochemical signaling events. As genetically encoded and live-cell compatible tools, fluorescent biosensors are well suited to address diverse cell signaling questions across different spatial scales of regulation. In this review, methods of examining spatial signaling regulation and the design of fluorescent biosensors are introduced. Then, recent biosensor developments that illuminate the importance of spatial regulation in cell signaling are highlighted at several scales, including membranes and organelles, molecular assemblies, and cell/tissue heterogeneity. In closing, perspectives on how fluorescent biosensors will continue enhancing cell signaling research are discussed.

Applications and Advances of Multicellular Tumor Spheroids: Challenges in Their Development and Analysis

Biomedical research requires both in vitro and in vivo studies in order to explore disease processes or drug interactions. Foundational investigations have been performed at the cellular level using two-dimensional cultures as the gold-standard method since the early 20th century. However, three-dimensional (3D) cultures have emerged as a new tool for tissue modeling over the last few years, bridging the gap between in vitro and animal model studies. Cancer has been a worldwide challenge for the biomedical community due to its high morbidity and mortality rates. Various methods have been developed to produce multicellular tumor spheroids (MCTSs), including scaffold-free and scaffold-based structures, which usually depend on the demands of the cells used and the related biological question. MCTSs are increasingly utilized in studies involving cancer cell metabolism and cell cycle defects. These studies produce massive amounts of data, which demand elaborate and complex tools for thorough analysis. In this review, we discuss the advantages and disadvantages of several up-to-date methods used to construct MCTSs. In addition, we also present advanced methods for analyzing MCTS features. As MCTSs more closely mimic the in vivo tumor environment, compared to 2D monolayers, they can evolve to be an appealing model for in vitro tumor biology studies.

Advancements in melanoma cancer metastasis models

Metastatic melanoma is a complex and deadly disease. Due to its complexity, the development of novel therapeutic strategies to inhibit metastatic melanoma remains an outstanding challenge. Our ability to study metastasis is advanced with the development of in vitro and in vivo models that better mimic the different steps of the metastatic cascade beginning from primary tumor initiation to final metastatic seeding. In this review, we provide a comprehensive summary of in vitro models, in vivo models, and in silico platforms to study the individual steps of melanoma metastasis. Furthermore, we highlight the advantages and limitations of each model and discuss the challenges of how to improve current models to enhance translation for melanoma cancer patients and future therapies.

Growth and adaptation mechanisms of tumour spheroids with time-dependent oxygen availability

Tumours are subject to external environmental variability. However, in vitro tumour spheroid experiments, used to understand cancer progression and develop cancer therapies, have been routinely performed for the past fifty years in constant external environments. Furthermore, spheroids are typically grown in ambient atmospheric oxygen (normoxia), whereas most in vivo tumours exist in hypoxic environments. Therefore, there are clear discrepancies between in vitro and in vivo conditions. We explore these discrepancies by combining tools from experimental biology, mathematical modelling, and statistical uncertainty quantification. Focusing on oxygen variability to develop our framework, we reveal key biological mechanisms governing tumour spheroid growth. Growing spheroids in time-dependent conditions, we identify and quantify novel biological adaptation mechanisms, including unexpected necrotic core removal, and transient reversal of the tumour spheroid growth phases.

Slow-cycling murine melanoma cells display plasticity and enhanced tumorigenicity in syngeneic transplantation assay

Slow-cycling cancer cells (SCC) contribute to the aggressiveness of many cancers, and their invasiveness and chemoresistance pose a great therapeutic challenge. However, in melanoma, their tumor-initiating abilities are not fully understood. In this study, we used the syngeneic transplantation assay to investigate the tumor-initiating properties of melanoma SCC in the physiologically relevant in vivo settings. For this we used B16-F10 murine melanoma cell line where we identified a small fraction of SCC. We found that, unlike human melanoma, the murine melanoma SCC were not marked by the high expression of the epigenetic enzyme Jarid1b. At the same time, their slow-cycling phenotype was a temporary state, similar to what was described in human melanoma. Progeny of SCC had slightly increased doxorubicin resistance and altered expression of melanogenesis genes, independent of the expression of cancer stem cell markers. Single-cell expansion of SCC revealed delayed growth and reduced clone formation when compared to non-SCC, which was further confirmed by an in vitro limiting dilution assay. Finally, syngeneic transplantation of 10 cells in vivo established that SCC were able to initiate growth in primary recipients and continue growth in the serial transplantation assay, suggesting their self-renewal nature. Together, our study highlights the high plasticity and tumorigenicity of murine melanoma SCC and suggests their role in melanoma aggressiveness.

From simplicity to complexity in current melanoma models

Despite the recent impressive clinical success of immunotherapy against melanoma, development of primary and adaptive resistance against immune checkpoint inhibitors remains a major issue in a large number of treated patients. This highlights the need for melanoma models that replicate the tumor's intricate dynamics in the tumor microenvironment (TME) and associated immune suppression to study possible resistance mechanisms in order to improve current and test novel therapeutics. While two-dimensional melanoma cell cultures have been widely used to perform functional genomics screens in a high-throughput fashion, they are not suitable to answer more complex scientific questions. Melanoma models have also been established in a variety of experimental (humanized) animals. However, due to differences in physiology, such models do not fully represent human melanoma development. Therefore, fully human three-dimensional in vitro models mimicking melanoma cell interactions with the TME are being developed to address this need for more physiologically relevant models. Such models include melanoma organoids, spheroids, and reconstructed human melanoma-in-skin cultures. Still, while major advances have been made to complement and replace animals, these in vitro systems have yet to fully recapitulate human tumor complexity. Lastly, technical advancements have been made in the organ-on-chip field to replicate functions and microstructures of in vivo human tissues and organs. This review summarizes advancements made in understanding and treating melanoma and specifically aims to discuss the progress made towards developing melanoma models, their applications, limitations, and the advances still needed to further facilitate the development of therapeutics.

Roles of IGFBP-3 in cell migration and growth in an endophytic tongue squamous cell carcinoma cell line

Insulin-like growth factor binding protein-3 (IGFBP-3) is a member of the IGFBP family that has high affinity for IGFs and functions as either an oncogene or tumor suppressor in various types of cancer. We previously found that IGFBP3 mRNA levels are higher in endophytic-type human tongue squamous cell carcinoma (TSCC) that is more invasive and more prone to metastasis than exophytic and superficial types. This finding prompted us to investigate the roles of IGFBP-3 in TSCC using SAS cells, which were originally derived from endophytic-type TSCC. Specifically, we used SAS cells that express a fluorescent ubiquitination-based cell-cycle indicator (Fucci). RNA-sequencing analysis indicated that IGFBP-3 is associated with cell migration and cell growth. In fact, IGFBP-3 knockdown downregulates cell migration and causes cells to arrest in G₁. This migratory potential appears to be cell cycle–independent. IGFBP-3 knockdown also reduced levels of secreted IGFBP-3; however, decreased migratory potential was not rescued by exogenous recombinant human IGFBP-3. Furthermore, ERK activity was downregulated by IGFBP-3 depletion, which suggests that MEK/ERK signaling may be involved in IGFBP-3-mediated cell migration. We therefore conclude that intracellular IGFBP-3 enhances cell migration independently of the cell cycle in TSCC with a higher metastatic potential.

Persister state-directed transitioning and vulnerability in melanoma

Melanoma is a highly plastic tumor characterized by dynamic interconversion of different cell identities depending on the biological context. Melanoma cells with high expression of the H3K4 demethylase KDM5B (JARID1B) rest in a slow-cycling, yet reversible persister state. Over time, KDM5Bhigh cells can promote rapid tumor repopulation with equilibrated KDM5B expression heterogeneity. The cellular identity of KDM5Bhigh persister cells has not been studied so far, missing an important cell state-directed treatment opportunity in melanoma. Here, we have established a doxycycline-titratable system for genetic induction of permanent intratumor expression of KDM5B and screened for chemical agents that phenocopy this effect. Transcriptional profiling and cell functional assays confirmed that the dihydropyridine 2-phenoxyethyl 4-(2-fluorophenyl)-2,7,7-trimethyl-5-oxo-1,4,5,6,7,8-hexa-hydro-quinoline-3-carboxylate (termed Cpd1) supports high KDM5B expression and directs melanoma cells towards differentiation along the melanocytic lineage and to cell cycle-arrest. The high KDM5B state additionally prevents cell proliferation through negative regulation of cytokinetic abscission. Moreover, treatment with Cpd1 promoted the expression of the melanocyte-specific tyrosinase gene specifically sensitizing melanoma cells for the tyrosinase-processed antifolate prodrug 3-O-(3,4,5-trimethoxybenzoyl)-(-)-epicatechin (TMECG). In summary, our study provides proof-of-concept for a dual hit strategy in melanoma, in which persister state-directed transitioning limits tumor plasticity and primes melanoma cells towards lineage-specific elimination.

A novel role for aspirin in enhancing the reprogramming function of miR-302/367 cluster and breast tumor suppression

Recent studies have provided evidence for tumor suppressive function of the embryonic stem cell-specific miR-302/367 cluster through induction of a reprogramming process. Aspirin has been found to induce reprogramming factors of mesenchymal-to-epithelial transition in breast cancer cells. Therefore, we aimed to investigate whether overexpression of miR-302/367 cluster and aspirin treatment cooperate in the induction of reprogramming and tumor suppression in breast cancer cells. MDA-MB-231 and SK-BR-3 human breast cancer cell lines were transfected with a miR-302/367 expressing vector and treated with aspirin. The cells were evaluated for indices of apoptosis, proliferation, migration, and invasion. In both cell lines, treatment of miR-302/367-transfected cells with aspirin upregulated expression of some main pluripotency factors such as OCT4, SOX2, NANOG, and KLF4, and downregulated expression of some invasion and angiogenesis markers at gene and protein levels. Aspirin increased the apoptotic rate in both cell lines transfected with miR-302/367. Both miR-302/367 and aspirin upregulated the expression of FOXD3 protein which is a known inducer of OCT4 and NANOG. Our results demonstrate that aspirin can enhance miR-302/367-induced reprogramming of breast cancer cells possibly through upregulation of FOXD3 expression. This can further augment the reversal of epithelial-mesenchymal transition and inhibits migration, invasion, and angiogenic signaling in breast cancer cells reprogrammed by miR-302/367. Therefore, aspirin may serve as a useful adjuvant for reprogramming of cancer cells.

MITF deficiency accelerates GNAQ-driven uveal melanoma

Cutaneous melanoma (CM) and uveal melanoma (UM) both originate from the melanocytic lineage but are primarily driven by distinct oncogenic drivers, BRAF/NRAS or GNAQ/GNA11, respectively. The melanocytic master transcriptional regulator, MITF, is essential for both CM development and maintenance, but its role in UM is largely unexplored. Here, we use zebrafish models to dissect the key UM oncogenic signaling events and establish the role of MITF in UM tumors. Using a melanocytic lineage expression system, we showed that patient-derived mutations of GNAQ (GNAQQ209L) or its upstream CYSLTR2 receptor (CYSLTR2L129Q) both drive UM when combined with a cooperating mutation, tp53M214K/M214K. The tumor-initiating potential of the major GNAQ/11 effector pathways, YAP, and phospholipase C-β (PLCβ)–ERK was also investigated in this system and thus showed that while activated YAP (YAPAA) induced UM with high potency, the patient-derived PLCβ4 mutation (PLCB4D630Y) very rarely yielded UM tumors in the tp53M214K/M214K context. Remarkably, mitfa deficiency was profoundly UM promoting, dramatically accelerating the onset and progression of tumors induced by Tg(mitfa:GNAQQ209L);tp53M214K/M214K or Tg(mitfa:CYSLTR2L129Q);tp53M214K/M214K. Moreover, mitfa loss was sufficient to cooperate with GNAQQ209L to drive tp53–wild type UM development and allowed Tg(mitfa:PLCB4D630Y);tp53M214K/M214K melanocyte lineage cells to readily form tumors. Notably, all of the mitfa−/− UM tumors, including those arising in Tg(mitfa:PLCB4D630Y);tp53M214K/M214K;mitfa−/− zebrafish, displayed nuclear YAP while lacking hyperactive ERK indicative of PLCβ signaling. Collectively, these data show that YAP signaling is the major mediator of UM and that MITF acts as a bona fide tumor suppressor in UM in direct opposition to its essential role in CM.

A stochastic mathematical model of 4D tumour spheroids with real-time fluorescent cell cycle labelling

In vitro tumour spheroids have been used to study avascular tumour growth and drug design for over 50 years. Tumour spheroids exhibit heterogeneity within the growing population that is thought to be related to spatial and temporal differences in nutrient availability. The recent development of real-time fluorescent cell cycle imaging allows us to identify the position and cell cycle status of individual cells within the growing spheroid, giving rise to the notion of a four-dimensional (4D) tumour spheroid. We develop the first stochastic individual-based model (IBM) of a 4D tumour spheroid and show that IBM simulation data compares well with experimental data using a primary human melanoma cell line. The IBM provides quantitative information about nutrient availability within the spheroid, which is important because it is difficult to measure these data experimentally.

Real-time visualization and quantitation of cell death and cell cycle progression in 2D and 3D cultures utilizing genetically encoded probes

Cancer cells grown as 3D-structures are better models for mimicking in vivo conditions than the 2D-culture systems employable in drug discovery applications. Cell cycle and cell death are important determinants for preclinical drug screening and tumor growth studies in laboratory conditions. Though several 3D-models and live-cell compatible approaches are available, a method for simultaneous real-time detection of cell cycle and cell death is required. Here we demonstrate a high-throughput adaptable method using genetically encoded fluorescent probes for the real-time quantitative detection of cell death and cell cycle. The cell-cycle indicator cdt1-Kusabira orange (KO) is stably integrated into cancer cells and further transfected with the Fluorescence Resonance Energy Transfer-based ECFP-DEVD-EYFP caspase activation sensor. The nuclear cdt1-KO expression serves as the readout for cell-cycle, and caspase activation is visualized by ECFP/EYFP ratiometric imaging. The image-based platform allowed imaging of growing spheres for prolonged periods in 3D-culture with excellent single-cell resolution through confocal microscopy. High-throughput screening (HTS) adaptation was achieved by targeting the caspase-sensor at the nucleus, which enabled the quantitation of cell death in 3D-models. The HTS using limited compound libraries, identified two lead compounds that induced caspase-activation both in 2D and 3D-cultures. This is the first report of an approach for noninvasive stain-free quantitative imaging of cell death and cell cycle with potential drug discovery applications.

Designing and interpreting 4D tumour spheroid experiments

Tumour spheroid experiments are routinely used to study cancer progression and treatment. Various and inconsistent experimental designs are used, leading to challenges in interpretation and reproducibility. Using multiple experimental designs, live-dead cell staining, and real-time cell cycle imaging, we measure necrotic and proliferation-inhibited regions in over 1000 4D tumour spheroids (3D space plus cell cycle status). By intentionally varying the initial spheroid size and temporal sampling frequencies across multiple cell lines, we collect an abundance of measurements of internal spheroid structure. These data are difficult to compare and interpret. However, using an objective mathematical modelling framework and statistical identifiability analysis we quantitatively compare experimental designs and identify design choices that produce reliable biological insight. Measurements of internal spheroid structure provide the most insight, whereas varying initial spheroid size and temporal measurement frequency is less important. Our general framework applies to spheroids grown in different conditions and with different cell types.

Quantitative analysis of tumour spheroid structure

Tumour spheroids are common in vitro experimental models of avascular tumour growth. Compared with traditional two-dimensional culture, tumour spheroids more closely mimic the avascular tumour microenvironment where spatial differences in nutrient availability strongly influence growth. We show that spheroids initiated using significantly different numbers of cells grow to similar limiting sizes, suggesting that avascular tumours have a limiting structure; in agreement with untested predictions of classical mathematical models of tumour spheroids. We develop a novel mathematical and statistical framework to study the structure of tumour spheroids seeded from cells transduced with fluorescent cell cycle indicators, enabling us to discriminate between arrested and cycling cells and identify an arrested region. Our analysis shows that transient spheroid structure is independent of initial spheroid size, and the limiting structure can be independent of seeding density. Standard experimental protocols compare spheroid size as a function of time; however, our analysis suggests that comparing spheroid structure as a function of overall size produces results that are relatively insensitive to variability in spheroid size. Our experimental observations are made using two melanoma cell lines, but our modelling framework applies across a wide range of spheroid culture conditions and cell lines.

Fluorescence-Based Quantitative and Spatial Analysis of Tumour Spheroids: A Proposed Tool to Predict Patient-Specific Therapy Response

Tumour spheroids are widely used to pre-clinically assess anti-cancer treatments. They are an excellent compromise between the lack of microenvironment encountered in adherent cell culture conditions and the great complexity of in vivo animal models. Spheroids recapitulate intra-tumour microenvironment-driven heterogeneity, a pivotal aspect for therapy outcome that is, however, often overlooked. Likely due to their ease, most assays measure overall spheroid size and/or cell death as a readout. However, as different tumour cell subpopulations may show a different biology and therapy response, it is paramount to obtain information from these distinct regions within the spheroid. We describe here a methodology to quantitatively and spatially assess fluorescence-based microscopy spheroid images by semi-automated software-based analysis. This provides a fast assay that accounts for spatial biological differences that are driven by the tumour microenvironment. We outline the methodology using detection of hypoxia, cell death and PBMC infiltration as examples, and we propose this procedure as an exploratory approach to assist therapy response prediction for personalised medicine.

Estimating parameters of a stochastic cell invasion model with fluorescent cell cycle labelling using approximate Bayesian computation

We develop a parameter estimation method based on approximate Bayesian computation (ABC) for a stochastic cell invasion model using fluorescent cell cycle labelling with proliferation, migration and crowding effects. Previously, inference has been performed on a deterministic version of the model fitted to cell density data, and not all parameters were identifiable. Considering the stochastic model allows us to harness more features of experimental data, including cell trajectories and cell count data, which we show overcomes the parameter identifiability problem. We demonstrate that, while difficult to collect, cell trajectory data can provide more information about the parameters of the cell invasion model. To handle the intractability of the likelihood function of the stochastic model, we use an efficient ABC algorithm based on sequential Monte Carlo. Rcpp and MATLAB implementations of the simulation model and ABC algorithm used in this study are available at https://github.com/michaelcarr-stats/FUCCI.

Rapid initiation of cell cycle reentry processes protects neurons from amyloid-β toxicity

Neurons are postmitotic cells. Reactivation of the cell cycle by neurons has been reported in Alzheimer's disease (AD) brains and models. This gave rise to the hypothesis that reentering the cell cycle renders neurons vulnerable and thus contributes to AD pathogenesis. Here, we use the fluorescent ubiquitination-based cell cycle indicator (FUCCI) technology to monitor the cell cycle in live neurons. We found transient, self-limited cell cycle reentry activity in naive neurons, suggesting that their postmitotic state is a dynamic process. Furthermore, we observed a diverse response to oligomeric amyloid-β (oAβ) challenge; neurons without cell cycle reentry activity would undergo cell death without activating the FUCCI reporter, while neurons undergoing cell cycle reentry activity at the time of the oAβ challenge could maintain and increase FUCCI reporter signal and evade cell death. Accordingly, we observed marked neuronal FUCCI positivity in the brains of human mutant Aβ precursor protein transgenic (APP23) mice together with increased neuronal expression of the endogenous cell cycle control protein geminin in the brains of 3-mo-old APP23 mice and human AD brains. Taken together, our data challenge the current view on cell cycle in neurons and AD, suggesting that pathways active during early cell cycle reentry in neurons protect from Aβ toxicity.

Endolysosomal Cation Channels and MITF in Melanocytes and Melanoma

Microphthalmia-associated transcription factor (MITF) is the principal transcription factor regulating pivotal processes in melanoma cell development, growth, survival, proliferation, differentiation and invasion. In recent years, convincing evidence has been provided attesting key roles of endolysosomal cation channels, specifically TPCs and TRPMLs, in cancer, including breast cancer, glioblastoma, bladder cancer, hepatocellular carcinoma and melanoma. In this review, we provide a gene expression profile of these channels in different types of cancers and decipher their roles, in particular the roles of two-pore channel 2 (TPC2) and TRPML1 in melanocytes and melanoma. We specifically discuss the signaling cascades regulating MITF and the relationship between endolysosomal cation channels, MAPK, canonical Wnt/GSK3 pathways and MITF.

Plasticity of Cancer Stem Cell: Origin and Role in Disease Progression and Therapy Resistance

In embryonic development and throughout life, there are some cells can exhibit phenotypic plasticity. Phenotypic plasticity is the ability of cells to differentiate into multiple lineages. In normal development, plasticity is highly regulated whereas cancer cells re-activate this dynamic ability for their own progression. The re-activation of these mechanisms enables cancer cells to acquire a cancer stem cell (CSC) phenotype- a subpopulation of cells with increased ability to survive in a hostile environment and resist therapeutic insults. There are several contributors fuel CSC plasticity in different stages of disease progression such as a complex network of tumour stroma, epidermal microenvironment and different sub-compartments within tumour. These factors play a key role in the transformation of tumour cells from a stable condition to a progressive state. In addition, flexibility in the metabolic state of CSCs helps in disease progression. Moreover, epigenetic changes such as chromatin, DNA methylation could stimulate the phenotypic change of CSCs. Development of resistance to therapy due to highly plastic behaviour of CSCs is a major cause of treatment failure in cancers. However, recent studies explored that plasticity can also expose the weaknesses in CSCs, thereby could be utilized for future therapeutic development. Therefore, in this review, we discuss how cancer cells acquire the plasticity, especially the role of the normal developmental process, tumour microenvironment, and epigenetic changes in the development of plasticity. We further highlight the therapeutic resistance property of CSCs attributed by plasticity. Also, outline some potential therapeutic options against plasticity of CSCs. Graphical Abstract .

Epithelial-to-mesenchymal-like transition events in melanoma

Epithelial-to-mesenchymal transition (EMT), a process through which epithelial tumor cells acquire mesenchymal phenotypic properties, contributes to both metastatic dissemination and therapy resistance in cancer. Accumulating evidence indicates that nonepithelial tumors, including melanoma, can also gain mesenchymal-like properties that increase their metastatic propensity and decrease their sensitivity to therapy. In this review, we discuss recent findings, illustrating the striking similarities-but also knowledge gaps-between the biology of mesenchymal-like state(s) in melanoma and mesenchymal state(s) from epithelial cancers. Based on this comparative analysis, we suggest hypothesis-driven experimental approaches to further deepen our understanding of the EMT-like process in melanoma and how such investigations may pave the way towards the identification of clinically relevant biomarkers for prognosis and new therapeutic strategies.

Superlinear growth reveals the Allee effect in tumors

Integrating experimental data into ecological models plays a central role in understanding biological mechanisms that drive tumor progression where such knowledge can be used to develop new therapeutic strategies. While the current studies emphasize the role of competition among tumor cells, they fail to explain recently observed superlinear growth dynamics across human tumors. Here we study tumor growth dynamics by developing a model that incorporates evolutionary dynamics inside tumors with tumor-microenvironment interactions. Our results reveal that tumor cells' ability to manipulate the environment and induce angiogenesis drives superlinear growth-a process compatible with the Allee effect. In light of this understanding, our model suggests that, for high-risk tumors that have a higher growth rate, suppressing angiogenesis can be the appropriate therapeutic intervention.

Mathematical Model of Tumour Spheroid Experiments with Real-Time Cell Cycle Imaging

Three-dimensional (3D) in vitro tumour spheroid experiments are an important tool for studying cancer progression and potential cancer drug therapies. Standard experiments involve growing and imaging spheroids to explore how different conditions lead to different rates of spheroid growth. These kinds of experiments, however, do not reveal any information about the spatial distribution of the cell cycle within the expanding spheroid. Since 2008, a new experimental technology called fluorescent ubiquitination-based cell cycle indicator (FUCCI) has enabled real-time in situ visualisation of the cell cycle progression. Observations of 3D tumour spheroids with FUCCI labelling reveal significant intratumoural structure, as the cell cycle status can vary with location. Although many mathematical models of tumour spheroid growth have been developed, none of the existing mathematical models are designed to interpret experimental observations with FUCCI labelling. In this work, we adapt the mathematical framework originally proposed by Ward and King (Math Med Biol 14:39-69, 1997. https://doi.org/10.1093/imammb/14.1.39 ) to produce a new mathematical model of FUCCI-labelled tumour spheroid growth. The mathematical model treats the spheroid as being composed of three subpopulations: (i) living cells in G1 phase that fluoresce red; (ii) living cells in S/G2/M phase that fluoresce green; and (iii) dead cells that are not fluorescent. We assume that the rates at which cells pass through different phases of the cell cycle, and the rate of cell death, depend upon the local oxygen concentration. Parameterising the new mathematical model using experimental measurements of cell cycle transition times, we show that the model can qualitatively capture important experimental observations that cannot be addressed using previous mathematical models. Further, we show that the mathematical model can be used to qualitatively mimic the action of anti-mitotic drugs applied to the spheroid. All software programs required to solve the nonlinear moving boundary problem associated with the new mathematical model are available on GitHub. at https://github.com/wang-jin-mathbio/Jin2021.

A novel mathematical model of heterogeneous cell proliferation

We present a novel mathematical model of heterogeneous cell proliferation where the total population consists of a subpopulation of slow-proliferating cells and a subpopulation of fast-proliferating cells. The model incorporates two cellular processes, asymmetric cell division and induced switching between proliferative states, which are important determinants for the heterogeneity of a cell population. As motivation for our model we provide experimental data that illustrate the induced-switching process. Our model consists of a system of two coupled delay differential equations with distributed time delays and the cell densities as functions of time. The distributed delays are bounded and allow for the choice of delay kernel. We analyse the model and prove the nonnegativity and boundedness of solutions, the existence and uniqueness of solutions, and the local stability characteristics of the equilibrium points. We find that the parameters for induced switching are bifurcation parameters and therefore determine the long-term behaviour of the model. Numerical simulations illustrate and support the theoretical findings, and demonstrate the primary importance of transient dynamics for understanding the evolution of many experimental cell populations.

Synchronized oscillations in growing cell populations are explained by demographic noise

Understanding synchrony in growing populations is important for applications as diverse as epidemiology and cancer treatment. Recent experiments employing fluorescent reporters in melanoma cell lines have uncovered growing subpopulations exhibiting sustained oscillations, with nearby cells appearing to synchronize their cycles. In this study, we demonstrate that the behavior observed is consistent with long-lasting transient phenomenon initiated and amplified by the finite-sample effects and demographic noise. We present a novel mathematical analysis of a multistage model of cell growth, which accurately reproduces the synchronized oscillations. As part of the analysis, we elucidate the transient and asymptotic phases of the dynamics and derive an analytical formula to quantify the effect of demographic noise in the appearance of the oscillations. The implications of these findings are broad, such as providing insight into experimental protocols that are used to study the growth of asynchronous populations and, in particular, those investigations relating to anticancer drug discovery.

To Divide or Invade: A Look Behind the Scenes of the Proliferation-Invasion Interplay in the Caenorhabditis elegans Anchor Cell

Cell invasion is defined by the capability of cells to migrate across compartment boundaries established by basement membranes (BMs). The development of complex organs involves regulated cell growth and regrouping of different cell types, which are enabled by controlled cell proliferation and cell invasion. Moreover, when a malignant tumor takes control over the body, cancer cells evolve to become invasive, allowing them to spread to distant sites and form metastases. At the core of the switch between proliferation and invasion are changes in cellular morphology driven by remodeling of the cytoskeleton. Proliferative cells utilize their actomyosin network to assemble a contractile ring during cytokinesis, while invasive cells form actin-rich protrusions, called invadopodia that allow them to breach the BMs. Studies of developmental cell invasion as well as of malignant tumors revealed that cell invasion and proliferation are two mutually exclusive states. In particular, anchor cell (AC) invasion during Caenorhabditis elegans larval development is an excellent model to study the transition from cell proliferation to cell invasion under physiological conditions. This mini-review discusses recent insights from the C. elegans AC invasion model into how G1 cell-cycle arrest is coordinated with the activation of the signaling networks required for BM breaching. Many regulators of the proliferation-invasion network are conserved between C. elegans and mammals. Therefore, the worm may provide important clues to better understand cell invasion and metastasis formation in humans.

Persistence as an Optimal Hedging Strategy

Bacteria invest in a slow-growing subpopulation, called persisters, to ensure survival in the face of uncertainty. This hedging strategy is remarkably similar to financial hedging, where diversifying an investment portfolio protects against economic uncertainty. We provide a new, to our knowledge, theoretical foundation for understanding cellular hedging by unifying the study of biological population dynamics and the mathematics of financial risk management through optimal control theory. Motivated by the widely accepted role of volatility in the emergence of persistence, we consider several models of environmental volatility described by continuous-time stochastic processes. This allows us to study an emergent cellular hedging strategy that maximizes the expected per capita growth rate of the population. Analytical and simulation results probe the optimal persister strategy, revealing results that are consistent with experimental observations and suggest new opportunities for experimental investigation and design. Overall, we provide a new, to our knowledge, way of conceptualizing and modeling cellular decision making in volatile environments by explicitly unifying theory from mathematical biology and finance.

Intravital imaging reveals cell cycle-dependent myogenic cell migration during muscle regeneration

During muscle regeneration, extracellular signal-regulated kinase (ERK) promotes both proliferation and migration. However, the relationship between proliferation and migration is poorly understood in this context. To elucidate this complex relationship on a physiological level, we established an intravital imaging system for measuring ERK activity, migration speed, and cell-cycle phases in mouse muscle satellite cell-derived myogenic cells. We found that in vivo, ERK is maximally activated in myogenic cells two days after injury, and this is then followed by increases in cell number and motility. With limited effects of ERK activity on migration on an acute timescale, we hypothesized that ERK increases migration speed in the later phase by promoting cell-cycle progression. Our cell-cycle analysis further revealed that in myogenic cells, ERK activity is critical for G1/S transition, and cells migrate more rapidly in S/G2 phase 3 days after injury. Finally, migration speed of myogenic cells was suppressed after CDK1/2-but not CDK1-inhibitor treatment, demonstrating a critical role of CDK2 in myogenic cell migration. Overall, our study demonstrates that in myogenic cells, the ERK-CDK2 axis promotes not only G1/S transition but also migration, thus providing a novel mechanism for efficient muscle regeneration.

FUCCI Real-Time Cell-Cycle Imaging as a Guide for Designing Improved Cancer Therapy: A Review of Innovative Strategies to Target Quiescent Chemo-Resistant Cancer Cells

Simple Summary

Chemotherapy of solid tumors has made very slow progress over many decades. A major problem has been that solid tumors very often contain non-dividing cells due to lack of oxygen deep in the tumor and these non-dividing cells resist most currently-used chemotherapy which usually only targets dividing cells. The present review demonstrates how a unique imaging system, FUCCI, which color codes cells depending on whether they are in a dividing or non-dividing phase, is being used to design very novel therapy that targets non-dividing cancer cells which can greatly improve the efficacy of cancer chemotherapy.

Abstract

Progress in chemotherapy of solid cancer has been tragically slow due, in large part, to the chemoresistance of quiescent cancer cells in tumors. The fluorescence ubiquitination cell-cycle indicator (FUCCI) was developed in 2008 by Miyawaki et al., which color-codes the phases of the cell cycle in real-time. FUCCI utilizes genes linked to different color fluorescent reporters that are only expressed in specific phases of the cell cycle and can, thereby, image the phases of the cell cycle in real-time. Intravital real-time FUCCI imaging within tumors has demonstrated that an established tumor comprises a majority of quiescent cancer cells and a minor population of cycling cancer cells located at the tumor surface or in proximity to tumor blood vessels. In contrast to most cycling cancer cells, quiescent cancer cells are resistant to cytotoxic chemotherapy, most of which target cells in S/G₂/M phases. The quiescent cancer cells can re-enter the cell cycle after surviving treatment, which suggests the reason why most cytotoxic chemotherapy is often ineffective for solid cancers. Thus, quiescent cancer cells are a major impediment to effective cancer therapy. FUCCI imaging can be used to effectively target quiescent cancer cells within tumors. For example, we review how FUCCI imaging can help to identify cell-cycle-specific therapeutics that comprise decoy of quiescent cancer cells from G₁ phase to cycling phases, trapping the cancer cells in S/G₂ phase where cancer cells are mostly sensitive to cytotoxic chemotherapy and eradicating the cancer cells with cytotoxic chemotherapy most active against S/G₂ phase cells. FUCCI can readily image cell-cycle dynamics at the single cell level in real-time in vitro and in vivo. Therefore, visualizing cell cycle dynamics within tumors with FUCCI can provide a guide for many strategies to improve cell-cycle targeting therapy for solid cancers.

Abrogation of RAB27A expression transiently affects melanoma cell proliferation

The role of the small GTPase RAB27A as an essential melanosome trafficking regulator in melanocytes is well-accepted. A decade ago, RAB27A was identified as a tumor dependency gene that promotes melanoma cell proliferation. RAB27A has since been linked to another propeller of cancer progression: exosome secretion. We have recently demonstrated that RAB27A is overexpressed in a subset of melanomas. High RAB27A gene and protein expression correlate with poor prognosis in melanoma patients. Mechanistic investigations revealed that the generation of pro-invasive exosomes was RAB27A-dependent and, therefore, silencing RAB27A reduced melanoma cell invasion in vitro and in vivo. However, previous studies have implicated RAB27A to be involved in both proliferation and invasion of melanoma cells. Employing four human cell lines, stratified by RAB27A expression, and one RAB27A-high mouse cell line, we demonstrate in this study that the effects of abrogating RAB27A expression on proliferation are only temporary, in contrast to our previously reported persistent effects on tumor invasion and metastasis. Therefore, we assist in the dissection of the short-term effects of RAB27A knockdown on melanoma cell proliferation versus long-term effects on melanoma invasion and metastasis. We believe that our findings provide novel insights into the effects of RAB27A blockade.

The EMT Transcription Factor ZEB2 Promotes Proliferation of Primary and Metastatic Melanoma While Suppressing an Invasive, Mesenchymal-Like Phenotype

Epithelial-to-mesenchymal transition (EMT)-inducing transcription factors (TF) are well known for their ability to induce mesenchymal states associated with increased migratory and invasive properties. Unexpectedly, nuclear expression of the EMT-TF ZEB2 in human primary melanoma has been shown to correlate with reduced invasion. We report here that ZEB2 is required for outgrowth for primary melanomas and metastases at secondary sites. Ablation of Zeb2 hampered outgrowth of primary melanomas in vivo, whereas ectopic expression enhanced proliferation and growth at both primary and secondary sites. Gain of Zeb2 expression in pulmonary-residing melanoma cells promoted the development of macroscopic lesions. In vivo fate mapping made clear that melanoma cells undergo a conversion in state where ZEB2 expression is replaced by ZEB1 expression associated with gain of an invasive phenotype. These findings suggest that reversible switching of the ZEB2/ZEB1 ratio enhances melanoma metastatic dissemination. SIGNIFICANCE: ZEB2 function exerts opposing behaviors in melanoma by promoting proliferation and expansion and conversely inhibiting invasiveness, which could be of future clinical relevance. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/14/2983/F1.large.jpg.

Practical parameter identifiability for spatio-temporal models of cell invasion

We examine the practical identifiability of parameters in a spatio-temporal reaction-diffusion model of a scratch assay. Experimental data involve fluorescent cell cycle labels, providing spatial information about cell position and temporal information about the cell cycle phase. Cell cycle labelling is incorporated into the reaction-diffusion model by treating the total population as two interacting subpopulations. Practical identifiability is examined using a Bayesian Markov chain Monte Carlo (MCMC) framework, confirming that the parameters are identifiable when we assume the diffusivities of the subpopulations are identical, but that the parameters are practically non-identifiable when we allow the diffusivities to be distinct. We also assess practical identifiability using a profile likelihood approach, providing similar results to MCMC with the advantage of being an order of magnitude faster to compute. Therefore, we suggest that the profile likelihood ought to be adopted as a screening tool to assess practical identifiability before MCMC computations are performed.

miR-16 enhances miR-302/367-induced reprogramming and tumor suppression in breast cancer cells

Overexpression of either miR-302 or miR-302/367 cluster induces reprogramming of cancer cells and exerts tumor-suppressive effects by induction of mesenchymal-to-epithelial transition, apoptosis and a less proliferative capacity. Several reports have described miR-16 as a tumor suppressor microRNA (miRNA). Here, we studied the impact of exogenous induction of miR-16 in MDA-MB-231 and SK-BR-3 breast cancer cells following overexpression of miR-302/367 cluster and investigated whether transfection of these cells by a mature miR-16 mimic could affect the reprogramming state of the cells and their tumorigenicity. miR-16 enhanced the expression levels of OCT4A, SOX2, and NANOG, generally known as transcription or pluripotency factors, and suppressed proliferation and invasiveness of these cells. Meanwhile, inhibition of miR-16 counteracted both the reprogramming effect and the antitumor function of miR-302/367 in the breast cancer cells. Current results indicate that miR-16 can work as an adjuvant to improve both cancer cell reprogramming and tumor-suppressive function of miR-302/367 cluster in MDA-MB-231 and SK-BR-3 cells, while its inhibition counteracts all of these effects. Combined application of miRNAs that share some common targets in cancer cell signaling pathways may provide new approaches for repression of multiple hallmarks of cancer.

Examining Go-or-Grow Using Fluorescent Cell-Cycle Indicators and Cell-Cycle-Inhibiting Drugs

The go-or-grow hypothesis states that adherent cells undergo reversible phenotype switching between migratory and proliferative states, with cells in the migratory state being more motile than cells in the proliferative state. Here, we examine go-or-grow in two-dimensional in vitro assays using melanoma cells with fluorescent cell-cycle indicators and cell-cycle-inhibiting drugs. We analyze the experimental data using single-cell tracking to calculate mean diffusivities and compare motility between cells in different cell-cycle phases and in cell-cycle arrest. Unequivocally, our analysis does not support the go-or-grow hypothesis. We present clear evidence that cell motility is independent of the cell-cycle phase and that nonproliferative arrested cells have the same motility as cycling cells.

Effect of heterogeneity and spatial correlations on the structure of a tumor invasion front in cellular environments

Analysis of invasion front has been widely used to decipher biological properties, as well as the growth dynamics of the corresponding populations. Likewise, the invasion front of tumors has been investigated, from which insights into the biological mechanisms of tumor growth have been gained. We develop a model to study how tumors' invasion front depends on the relevant properties of a cellular environment. To do so, we develop a model based on a nonlinear reaction-diffusion equation, the Fisher-Kolmogorov-Petrovsky-Piskunov equation, to model tumor growth. Our study aims to understand how heterogeneity in the cellular environment's stiffness, as well as spatial correlations in its morphology, the existence of both of which has been demonstrated by experiments, affects the properties of tumor invasion front. It is demonstrated that three important factors affect the properties of the front, namely the spatial distribution of the local diffusion coefficients, the spatial correlations between them, and the ratio of the cells' duplication rate and their average diffusion coefficient. Analyzing the scaling properties of tumor invasion front computed by solving the governing equation, we show that, contrary to several previous claims, the invasion front of tumors and cancerous cell colonies cannot be described by the well-known models of kinetic growth, such as the Kardar-Parisi-Zhang equation.

Melanoma plasticity and phenotypic diversity: therapeutic barriers and opportunities

An incomplete view of the mechanisms that drive metastasis, the primary cause of cancer-related death, has been a major barrier to development of effective therapeutics and prognostic diagnostics. Increasing evidence indicates that the interplay between microenvironment, genetic lesions, and cellular plasticity drives the metastatic cascade and resistance to therapies. Here, using melanoma as a model, we outline the diversity and trajectories of cell states during metastatic dissemination and therapy exposure, and highlight how understanding the magnitude and dynamics of nongenetic reprogramming in space and time at single-cell resolution can be exploited to develop therapeutic strategies that capitalize on nongenetic tumor evolution.

Mathematical models incorporating a multi-stage cell cycle replicate normally-hidden inherent synchronization in cell proliferation

We present a suite of experimental data showing that cell proliferation assays, prepared using standard methods thought to produce asynchronous cell populations, persistently exhibit inherent synchronization. Our experiments use fluorescent cell cycle indicators to reveal the normally hidden cell synchronization, by highlighting oscillatory subpopulations within the total cell population. These oscillatory subpopulations would never be observed without these cell cycle indicators. On the other hand, our experimental data show that the total cell population appears to grow exponentially, as in an asynchronous population. We reconcile these seemingly inconsistent observations by employing a multi-stage mathematical model of cell proliferation that can replicate the oscillatory subpopulations. Our study has important implications for understanding and improving experimental reproducibility. In particular, inherent synchronization may affect the experimental reproducibility of studies aiming to investigate cell cycle-dependent mechanisms, including changes in migration and drug response.

SMAD signaling promotes melanoma metastasis independently of phenotype switching

The development of metastatic melanoma is thought to require the dynamic shifting of neoplastic cells between proliferative and invasive phenotypes. Contrary to this conventional "phenotype switching" model, we now show that disease progression can involve malignant melanoma cells simultaneously displaying proliferative and invasive properties. Using a genetic mouse model of melanoma in combination with in vitro analyses of melanoma cell lines, we found that conditional deletion of the downstream signaling molecule Smad4, which abrogates all canonical TGF-β signaling, indeed inhibits both tumor growth and metastasis. Conditional deletion of the inhibitory signaling factor Smad7, however, generated cells that are both highly invasive and proliferative, indicating that invasiveness is compatible with a high proliferation rate. In fact, conditional Smad7 deletion led to sustained melanoma growth and at the same time promoted massive metastasis formation, a result consistent with data indicating that low SMAD7 levels in patient tumors are associated with a poor survival. Our findings reveal that modulation of SMAD7 levels can overcome the need for phenotype switching during tumor progression and may thus represent a novel therapeutic target in metastatic disease.

A Nonquiescent "Idling" Population State in Drug-Treated, BRAF-Mutated Melanoma

Targeted therapy is an effective standard of care in BRAF-mutated malignant melanoma. However, the duration of tumor remission varies unpredictably among patients, and relapse is almost inevitable. Here, we examine the responses of several BRAF-mutated melanoma cell lines (including isogenic subclones) to BRAF inhibitors. We observe complex response dynamics across cell lines, with short-term responses (<100 h) varying from cell line to cell line. In the long term, however, we observe equilibration of all drug-treated populations into a nonquiescent state characterized by a balanced rate of death and division, which we term the "idling" state, and to our knowledge, this state has not been previously reported. Using mathematical modeling, we propose that the observed population-level dynamics are the result of cells transitioning between basins of attraction within a drug-modified phenotypic landscape. Each basin is associated with a drug-induced proliferation rate, a recently introduced metric of an antiproliferative drug effect. The idling population state represents a new dynamic equilibrium in which cells are distributed across the landscape such that the population achieves zero net growth. By fitting our model to experimental drug-response data, we infer the phenotypic landscapes of all considered melanoma cell lines and provide a unifying view of how BRAF-mutated melanomas respond to BRAF inhibition. We hypothesize that the residual disease observed in patients after targeted therapy is composed of a significant number of idling cells. Thus, defining molecular determinants of the phenotypic landscape that idling populations occupy may lead to "targeted landscaping" therapies based on rational modification of the landscape to favor basins with greater drug susceptibility.

Magnolol induces cell death through PI3K/Akt-mediated epigenetic modifications boosting treatment of BRAF- and NRAS-mutant melanoma

Most BRAF‐mutant melanoma patients experience a fulminate relapse after several months of treatment with BRAF/MEK inhibitors. To improve therapeutic efficacy, natural plant‐derived compounds might be considered as potent additives. Here, we show that magnolol, a constituent of Magnolia officinalis, induced G1 arrest, apoptosis and cell death in BRAF‐ and NRAS‐mutant melanoma cells at low concentration, with no effect in BRAF‐ and NRAS wild‐type melanoma cells and human keratinocytes. This was confirmed in a 3D spheroid model. The apoptosis‐inducing effect of magnolol was completely rescued by activating Akt suggesting a mechanism relying primarily on Akt signaling. Magnolol significantly downregulated the PI3K/Akt pathway which led to a global decrease of the active histone mark H3K4me3. Alongside, the repressive histone mark H3K9me3 was increased as a response to DNA damage. Magnolol‐induced alterations of histone modifications are reversible upon activation of the Akt pathway. Magnolol‐induced a synergistic effect in combination with either BRAF/MEK inhibitors dabrafenib/trametinib or docetaxel at a lower concentration than usually applied in melanoma patients. Combination of magnolol with targeted therapy or chemotherapy also led to analogous effects on histone marks, which was rescued by Akt pathway activation. Our study revealed a novel epigenetic mechanism of magnolol‐induced cell death in melanoma. Magnolol might therefore be a clinically useful addition to BRAF/MEK inhibitors with enhanced efficacy delaying or preventing disease recurrence.