VEGF/VEGFR axis and its signaling in melanoma: Current knowledge toward therapeutic targeting agents and future perspectives.

Melanoma is responsible for most skin cancer-associated deaths globally. The progression of melanoma is influenced by a number of pathogenic processes. Understanding the VEGF/VEGFR axis, which includes VEGF-A, PlGF, VEGF-B, VEGF-C, and VEGF-D and their receptors, VEGFR-1, VEGFR-2, and VEGFR-3, is of great importance in melanoma due to its crucial role in angiogenesis. This axis generates multifactorial and complex cellular signaling, engaging the MAPK/ERK, PI3K/AKT, PKC, PLC-γ, and FAK signaling pathways. Melanoma cell growth and proliferation, migration and metastasis, survival, and acquired resistance to therapy are influenced by this axis. The VEGF/VEGFR axis was extensively examined for their potential as diagnostic/prognostic biomarkers in melanoma patients and results showed that VEGF overexpression can be associated with unfavorable prognosis, higher level of tumor invasion and poor response to therapy. MicroRNAs linking to the VEGF/VEGFR axis were identified and, in this review, divided into two categories according to their functions, some of them promote melanoma angiogenesis (promotive group) and some restrict melanoma angiogenesis (protective group). In addition, the approach of treating melanoma by targeting the VEGF/VEGFR axis has garnered significant interest among researchers. These agents can be divided into two main groups: anti-VEGF and VEGFR inhibitors. These therapeutic options may be a prominent step along with the modern targeting and immune therapies for better coverage of pathological processes leading to melanoma progression and therapy resistance.

Targeting 14-3-3θ-mediated TDP-43 pathology in amyotrophic lateral sclerosis and frontotemporal dementia mice.

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are characterized by cytoplasmic deposition of the nuclear TAR-binding protein 43 (TDP-43). Although cytoplasmic re-localization of TDP-43 is a key event in the pathogenesis of ALS/FTD, the underlying mechanisms remain unknown. Here, we identified a non-canonical interaction between 14-3-3θ and TDP-43, which regulates nuclear-cytoplasmic shuttling. Neuronal 14-3-3θ levels were increased in sporadic ALS and FTD with TDP-43 pathology. Pathogenic TDP-43 showed increased interaction with 14-3-3θ, resulting in cytoplasmic accumulation, insolubility, phosphorylation, and fragmentation of TDP-43, resembling pathological changes in disease. Harnessing this increased affinity of 14-3-3θ for pathogenic TDP-43, we devised a gene therapy vector targeting TDP-43 pathology, which mitigated functional deficits and neurodegeneration in different ALS/FTD mouse models expressing mutant or non-mutant TDP-43, including when already symptomatic at the time of treatment. Our study identified 14-3-3θ as a mediator of cytoplasmic TDP-43 localization with implications for ALS/FTD pathogenesis and therapy.

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.

The MC1R r allele does not increase melanoma risk in MITF E318K carriers.

BACKGROUND: Population-wide screening for melanoma is not cost-effective, but genetic characterization could facilitate risk stratification and targeted screening. Common Melanocortin-1 receptor (MC1R) red hair colour (RHC) variants and Microphthalmia-associated transcription factor (MITF) E318K separately confer moderate melanoma susceptibility, but their interactive effects are relatively unexplored. OBJECTIVES: To evaluate whether MC1R genotypes differentially affect melanoma risk in MITF E318K+ vs. E318K- individuals. MATERIALS AND METHODS: Melanoma status (affected or unaffected) and genotype data (MC1R and MITF E318K) were collated from research cohorts (five Australian and two European). In addition, RHC genotypes from E318K+ individuals with and without melanoma were extracted from databases (The Cancer Genome Atlas and Medical Genome Research Bank, respectively). χ2 and logistic regression were used to evaluate RHC allele and genotype frequencies within E318K+/- cohorts depending on melanoma status. Replication analysis was conducted on 200 000 general-population exomes (UK Biobank). RESULTS: The cohort comprised 1165 MITF E318K- and 322 E318K+ individuals. In E318K- cases MC1R R and r alleles increased melanoma risk relative to wild type (wt), P < 0.001 for both. Similarly, each MC1R RHC genotype (R/R, R/r, R/wt, r/r and r/wt) increased melanoma risk relative to wt/wt (P < 0.001 for all). In E318K+ cases, R alleles increased melanoma risk relative to the wt allele [odds ratio (OR) 2.04 (95% confidence interval 1.67-2.49); P = 0.01], while the r allele risk was comparable with the wt allele [OR 0.78 (0.54-1.14) vs. 1.00, respectively]. E318K+ cases with the r/r genotype had a lower but not significant melanoma risk relative to wt/wt [OR 0.52 (0.20-1.38)]. Within the E318K+ cohort, R genotypes (R/R, R/r and R/wt) conferred a significantly higher risk compared with non-R genotypes (r/r, r/wt and wt/wt) (P < 0.001). UK Biobank data supported our findings that r did not increase melanoma risk in E318K+ individuals. CONCLUSIONS: RHC alleles/genotypes modify melanoma risk differently in MITF E318K- and E318K+ individuals. Specifically, although all RHC alleles increase risk relative to wt in E318K- individuals, only MC1R R increases melanoma risk in E318K+ individuals. Importantly, in the E318K+ cohort the MC1R r allele risk is comparable with wt. These findings could inform counselling and management for MITF E318K+ individuals.

Head and neck cancer patient-derived tumouroid cultures: opportunities and challenges.

Head and neck cancers (HNC) are the seventh most prevalent cancer type globally. Despite their common categorisation, HNCs are a heterogeneous group of malignancies arising in various anatomical sites within the head and neck region. These cancers exhibit different clinical and biological manifestations, and this heterogeneity also contributes to the high rates of treatment failure and mortality. To evaluate patients who will respond to a particular treatment, there is a need to develop in vitro model systems that replicate in vivo tumour status. Among the methods developed, patient-derived cancer organoids, also known as tumouroids, recapitulate in vivo tumour characteristics including tumour architecture. Tumouroids have been used for general disease modelling and genetic instability studies in pan-cancer research. However, a limited number of studies have thus far been conducted using tumouroid-based drug screening. Studies have concluded that tumouroids can play an essential role in bringing precision medicine for highly heterogenous cancer types such as HNC.

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.

Rapid Optical Clearing for Semi-High-Throughput Analysis of Tumor Spheroids.

Tumor spheroids are fast becoming commonplace in basic cancer research and drug development. Obtaining data regarding protein expression within the spheroid at the cellular level is important for analysis, yet existing techniques are often expensive, laborious, use non-standard equipment, cause significant size distortion, or are limited to relatively small spheroids. This protocol presents a new method of mounting and clearing spheroids that address these issues while allowing for confocal analysis of the inner structure of spheroids. In contrast to existing approaches, this protocol provides for rapid mounting and clearing of a large number of spheroids using standard equipment and laboratory supplies. Mounting spheroids in a pH-neutral agarose-PBS gel solution before introducing a refractive-index-matched clearing solution minimizes size distortion common to other similar techniques. This allows for detailed quantitative and statistical analysis where the accuracy of size measurements is paramount. Furthermore, compared to liquid clearing solutions, the agarose gel technique keeps spheroids fixed in place, allowing for the collection of three-dimensional (3D) confocal images. The present article elaborates how the method yields high-quality two- and 3D images that provide information about inter-cell variability and inner spheroid structure.

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 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.

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.

PTEN: A novel target for vitamin D in melanoma.

Melanoma is the most dangerous form of skin cancer, with poor prognosis in advanced stages. Vitamin D, also produced by ultraviolet radiation, is known for its anti-proliferative properties in some cancers including melanoma. While vitamin D deficiency has been associated with advanced melanoma stage and higher levels of vitamin D have been associated with better outcomes, the role for vitamin D in melanoma remains unclear. Vitamin D synthesis is initiated upon UVB exposure of skin cells and results in formation of the active metabolite 1,25-dihydroxyvitamin D3 (1,25D). We have previously demonstrated that 1,25D plays a role in protection against ultraviolet radiation-induced DNA damage, immune suppression, and skin carcinogenesis. In this study 1,25D significantly reduced cell viability and increased caspase levels in human melanoma cell lines. This effect was not present in cells that lacked both phosphatase and tensin homolog deleted on chromosome 10 (PTEN), a well-known tumour suppressor, and the vitamin D receptor (VDR). PTEN is frequently lost or mutated in melanoma. Incubation of selected melanoma cell lines with 1,25D resulted in significant increases in PTEN levels and downregulation of the AKT pathway and its downstream effectors. This suggests that 1,25D may act to reduce melanoma cell viability by targeting PTEN.

Reciprocal Regulation of BRN2 and NOTCH1/2 Signaling Synergistically Drives Melanoma Cell Migration and Invasion.

Phenotypic plasticity drives cancer progression, impacts treatment response, and is a major driver of therapeutic resistance. In melanoma, a regulatory axis between the MITF and BRN2 transcription factors has been reported to promote tumor heterogeneity by mediating switching between proliferative and invasive phenotypes, respectively. Despite strong evidence that subpopulations of cells that exhibit a BRN2high/MITFlow expression profile switch to a predominantly invasive phenotype, the mechanisms by which this switch is propagated and promotes invasion remain poorly defined. We have found that a reciprocal relationship between BRN2 and NOTCH1/2 signaling exists in melanoma cells in vitro, within patient datasets, and in in vivo primary and metastatic human tumors that bolsters acquisition of invasiveness. Working through the epigenetic modulator EZH2, the BRN2‒NOTCH1/2 axis is potentially a key mechanism by which the invasive phenotype is maintained. Given the emergence of agents targeting both EZH2 and NOTCH, understanding the mechanism through which BRN2 promotes heterogeneity may provide crucial biomarkers to predict treatment response to prevent metastasis.

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.

Targeting Replication Stress Using CHK1 Inhibitor Promotes Innate and NKT Cell Immune Responses and Tumour Regression.

Drugs selectively targeting replication stress have demonstrated significant preclinical activity, but this has not yet translated into an effective clinical treatment. Here we report that targeting increased replication stress with a combination of Checkpoint kinase 1 inhibitor (CHK1i) with a subclinical dose of hydroxyurea targets also promotes pro-inflammatory cytokine/chemokine expression that is independent of cGAS-STING pathway activation and immunogenic cell death in human and murine melanoma cells. In vivo, this drug combination induces tumour regression which is dependent on an adaptive immune response. It increases cytotoxic CD8+ T cell activity, but the major adaptive immune response is a pronounced NKT cell tumour infiltration. Treatment also promotes an immunosuppressive tumour microenvironment through CD4+ Treg and FoxP3+ NKT cells. The number of these accumulated during treatment, the increase in FoxP3+ NKT cells numbers correlates with the decrease in activated NKT cells, suggesting they are a consequence of the conversion of effector to suppressive NKT cells. Whereas tumour infiltrating CD8+ T cell PD-1 and tumour PD-L1 expression was increased with treatment, peripheral CD4+ and CD8+ T cells retained strong anti-tumour activity. Despite increased CD8+ T cell PD-1, combination with anti-PD-1 did not improve response, indicating that immunosuppression from Tregs and FoxP3+ NKT cells are major contributors to the immunosuppressive tumour microenvironment. This demonstrates that therapies targeting replication stress can be well tolerated, not adversely affect immune responses, and trigger an effective anti-tumour immune response.

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.

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.