An evolutionary framework for treating pediatric sarcomas

The effect of tumor composition on the success of adaptive therapy: The case of metastatic Castrate-Resistant Prostate Cancer

Prostate-specific antigen (PSA) is the most commonly used serum marker for prostate cancer. It plays a role in cancer detection, treatment monitoring, and more recently, in guiding adaptive therapy protocols, where treatment is alternated based on PSA levels. However, the relationship between PSA levels and tumor volume remains poorly understood. Empirical evidence suggests that different cancer cell types produce varying amounts of PSA. Despite this, current mathematical cancer models often assume either that all cell types contribute equally to PSA levels or that only certain subpopulations produce PSA at fixed rates. In this study, we compare Zhang et al.'s classical adaptive therapy protocol with the standard of care, which involves continuous maximum tolerable dose treatment, under different assumptions regarding PSA production. Specifically, we explore the possibility that testosterone-dependent, testosterone-producing, and testosterone-independent cells contribute to PSA production to varying degrees. We use the time to competitive release as a proxy for the time to disease progression. Our findings indicate that adaptive therapy consistently results in a longer time to competitive release compared to the standard of care, regardless of the assumptions about PSA production. However, when testosterone-independent cells are the sole PSA producers, Zhang et al.'s adaptive therapy protocol becomes inapplicable, as PSA levels never fall to half of their initial value, preventing therapy discontinuation. Additionally, we observe that the number and duration of treatment cycles in adaptive therapy are highly sensitive to assumptions about how much each cell type contributes to PSA production. Overall, our results emphasize the need for a deeper understanding of patient-specific PSA dynamics, which could enhance the effectiveness of adaptive therapy in prostate cancer treatment.

Pediatric Leukemia Roadmaps Are a Guide for Positive Metastatic Bone Sarcoma Trials

ALL cures require many MRD therapies. This strategy should drive experiments and trials in metastatic bone sarcomas.

Ecological and evolutionary dynamics to design and improve ovarian cancer treatment

Ovarian cancer ecosystems are exceedingly complex, consisting of a high heterogeneity of cancer cells. Development of drugs such as poly ADP-ribose polymerase (PARP) inhibitors, targeted therapies and immunotherapies offer more options for sequential or combined treatments. Nevertheless, mortality in metastatic ovarian cancer patients remains high because cancer cells consistently develop resistance to single and combination therapies, urging a need for treatment designs that target the evolvability of cancer cells. The evolutionary dynamics that lead to resistance emerge from the complex tumour microenvironment, the heterogeneous populations, and the individual cancer cell's plasticity. We propose that successful management of ovarian cancer requires consideration of the ecological and evolutionary dynamics of the disease. Here, we review current options and challenges in ovarian cancer treatment and discuss principles of tumour evolution. We conclude by proposing evolutionarily designed strategies for ovarian cancer, with the goal of integrating such principles with longitudinal, quantitative data to improve the treatment design and management of drug resistance. KEY POINTS/HIGHLIGHTS: Tumours are ecosystems in which cancer and non-cancer cells interact and evolve in complex and dynamic ways. Conventional therapies for ovarian cancer inevitably lead to the development of resistance because they fail to consider tumours' heterogeneity and cellular plasticity. Eco-evolutionarily designed therapies should consider cancer cell plasticity and patient-specific characteristics to improve clinical outcome and prevent relapse.

Preventing evolutionary rescue in cancer

First-line cancer treatment frequently fails due to initially rare therapeutic resistance. An important clinical question is then how to schedule subsequent treatments to maximize the probability of tumour eradication. Here, we provide a theoretical solution to this problem by using mathematical analysis and extensive stochastic simulations within the framework of evolutionary rescue theory to determine how best to exploit the vulnerability of small tumours to stochastic extinction. Whereas standard clinical practice is to wait for evidence of relapse, we confirm a recent hypothesis that the optimal time to switch to a second treatment is when the tumour is close to its minimum size before relapse, when it is likely undetectable. This optimum can lie slightly before or slightly after the nadir, depending on tumour parameters. Given that this exact time point may be difficult to determine in practice, we study windows of high extinction probability that lie around the optimal switching point, showing that switching after the relapse has begun is typically better than switching too early. We further reveal how treatment dose and tumour demographic and evolutionary parameters influence the predicted clinical outcome, and we determine how best to schedule drugs of unequal efficacy. Our work establishes a foundation for further experimental and clinical investigation of this evolutionarily-informed "extinction therapy" strategy.

An evolutionary differential game for regulating the role of monoclonal antibodies in treating signalling pathways in oesophageal cancer

This work presents a new framework for a competitive evolutionary game between monoclonal antibodies and signalling pathways in oesophageal cancer. The framework is based on a novel dynamical model that takes into account the dynamic progression of signalling pathways, resistance mechanisms and monoclonal antibody therapies. This game involves a scenario in which signalling pathways and monoclonal antibodies are the players competing against each other, where monoclonal antibodies use Brentuximab and Pembrolizumab dosages as strategies to counter the evolutionary resistance strategy implemented by the signalling pathways. Their interactions are described by the dynamical model, which serves as the game's playground. The analysis and computation of two game-theoretic strategies, Stackelberg and Nash equilibria, are conducted within this framework to ascertain the most favourable outcome for the patient. By comparing Stackelberg equilibria with Nash equilibria, numerical experiments show that the Stackelberg equilibria are superior for treating signalling pathways and are critical for the success of monoclonal antibodies in improving oesophageal cancer patient outcomes.

Roadmap for the next generation of Children's Oncology Group rhabdomyosarcoma trials

Clinical trials conducted by the Intergroup Rhabdomyosarcoma (RMS) Study Group and the Children's Oncology Group have been pivotal to establishing current standards for diagnosis and therapy for RMS. Recent advancements in understanding the biology and clinical behavior of RMS have led to more nuanced approaches to diagnosis, risk stratification, and treatment. The complexities introduced by these advancements, coupled with the rarity of RMS, pose challenges to conducting large-scale phase 3 clinical trials to evaluate new treatment strategies for RMS. Given these challenges, systematic planning of future clinical trials in RMS is paramount to address pertinent questions regarding the therapeutic efficacy of drugs, biomarkers of response, treatment-related toxicity, and patient quality of life. Herein, the authors outline the proposed strategic approach of the Children's Oncology Group Soft Tissue Sarcoma Committee to the next generation of RMS clinical trials, focusing on five themes: improved novel agent identification and preclinical to clinical translation, more efficient trial development and implementation, expanded opportunities for knowledge generation during trials, therapeutic toxicity reduction and quality of life, and patient engagement.

A war on many fronts: cross disciplinary approaches for novel cancer treatment strategies

Cancer is a disease characterized by uncontrolled cellular growth where cancer cells take advantage of surrounding cellular populations to obtain resources and promote invasion. Carcinomas are the most common type of cancer accounting for almost 90% of cancer cases. One of the major subtypes of carcinomas are adenocarcinomas, which originate from glandular cells that line certain internal organs. Cancers such as breast, prostate, lung, pancreas, colon, esophageal, kidney are often adenocarcinomas. Current treatment strategies include surgery, chemotherapy, radiation, targeted therapy, and more recently immunotherapy. However, patients with adenocarcinomas often develop resistance or recur after the first line of treatment. Understanding how networks of tumor cells interact with each other and the tumor microenvironment is crucial to avoid recurrence, resistance, and high-dose therapy toxicities. In this review, we explore how mathematical modeling tools from different disciplines can aid in the development of effective and personalized cancer treatment strategies. Here, we describe how concepts from the disciplines of ecology and evolution, economics, and control engineering have been applied to mathematically model cancer dynamics and enhance treatment strategies.

Seclidemstat blocks the transcriptional function of multiple FET-fusion oncoproteins

Genes encoding the RNA-binding proteins FUS, EWSR1, and TAF15 (FET proteins) are involved in chromosomal translocations in rare sarcomas. FET-rearranged sarcomas are often aggressive malignancies affecting patients of all ages. New therapies are needed. These translocations fuse the 5' portion of the FET gene with a 3' partner gene encoding a transcription factor (TF). The resulting fusion proteins are oncogenic TFs with a FET protein low complexity domain (LCD) and a DNA binding domain. FET fusion proteins have proven stubbornly difficult to target directly and promising strategies target critical co-regulators. One candidate is lysine specific demethylase 1 (LSD1). LSD1 is recruited by multiple FET fusions, including EWSR1::FLI1. LSD1 promotes EWSR1::FLI1 activity and treatment with the noncompetitive inhibitor SP-2509 blocks EWSR1::FLI1 transcriptional function. A similar molecule, seclidemstat (SP-2577), is currently in clinical trials for FET-rearranged sarcomas (NCT03600649). However, whether seclidemstat has pharmacological activity against FET fusions has not been demonstrated. Here, we evaluate the in vitro potency of seclidemstat against multiple FET-rearranged sarcoma cell lines, including Ewing sarcoma, desmoplastic small round cell tumor, clear cell sarcoma, and myxoid liposarcoma. We also define the transcriptomic effects of seclidemstat treatment and evaluated the activity of seclidemstat against FET fusion transcriptional regulation. Seclidemstat showed potent activity in cell viability assays across FET-rearranged sarcomas and disrupted the transcriptional function of all tested fusions. Though epigenetic and targeted inhibitors are unlikely to be effective as a single agents in the clinic, these data suggest seclidemstat remains a promising new treatment strategy for patients with FET-rearranged sarcomas.

Modeling phenotypic heterogeneity towards evolutionarily inspired osteosarcoma therapy

Osteosarcoma is the most common bone sarcoma in children and young adults. While universally delivered, chemotherapy only benefits roughly half of patients with localized disease. Increasingly, intratumoral heterogeneity is recognized as a source of therapeutic resistance. In this study, we develop and evaluate an in vitro model of osteosarcoma heterogeneity based on phenotype and genotype. Cancer cell populations vary in their environment-specific growth rates and in their sensitivity to chemotherapy. We present the genotypic and phenotypic characterization of an osteosarcoma cell line panel with a focus on co-cultures of the most phenotypically divergent cell lines, 143B and SAOS2. Modest environmental (pH, glutamine) or chemical perturbations dramatically shift the success and composition of cell lines. We demonstrate that in nutrient rich culture conditions 143B outcompetes SAOS2. But, under nutrient deprivation or conventional chemotherapy, SAOS2 growth can be favored in spheroids. Importantly, when the simplest heterogeneity state is evaluated, a two-cell line coculture, perturbations that affect the faster growing cell line have only a modest effect on final spheroid size. Thus the only evaluated therapies to eliminate the spheroids were by switching therapies from a first strike to a second strike. This extensively characterized, widely available system, can be modeled and scaled to allow for improved strategies to anticipate resistance in osteosarcoma due to heterogeneity.

Subclonal Somatic Copy-Number Alterations Emerge and Dominate in Recurrent Osteosarcoma

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease.

SIGNIFICANCE

The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Treatment of evolving cancers will require dynamic decision support

Cancer research has traditionally focused on developing new agents, but an underexplored question is that of the dose and frequency of existing drugs. Based on the modus operandi established in the early days of chemotherapies, most drugs are administered according to predetermined schedules that seek to deliver the maximum tolerated dose and are only adjusted for toxicity. However, we believe that the complex, evolving nature of cancer requires a more dynamic and personalized approach. Chronicling the milestones of the field, we show that the impact of schedule choice crucially depends on processes driving treatment response and failure. As such, cancer heterogeneity and evolution dictate that a one-size-fits-all solution is unlikely-instead, each patient should be mapped to the strategy that best matches their current disease characteristics and treatment objectives (i.e. their 'tumorscape'). To achieve this level of personalization, we need mathematical modeling. In this perspective, we propose a five-step 'Adaptive Dosing Adjusted for Personalized Tumorscapes (ADAPT)' paradigm to integrate data and understanding across scales and derive dynamic and personalized schedules. We conclude with promising examples of model-guided schedule personalization and a call to action to address key outstanding challenges surrounding data collection, model development, and integration.

Clinical trial designs for evaluating and exploiting cancer evolution

The evolution of drug-resistant cell subpopulations causes cancer treatment failure. Current preclinical evidence shows that it is possible to model herding of clonal evolution and collateral sensitivity where an initial treatment could favourably influence the response to a subsequent one. Novel therapy strategies exploiting this understanding are being considered, and clinical trial designs for steering cancer evolution are needed. Furthermore, preclinical evidence suggests that different subsets of drug-sensitive and resistant clones could compete between themselves for nutrients/blood supply, and clones that populate a tumour do so at the expense of other clones. Treatment paradigms based on this clinical application of exploiting cell-cell competition include intermittent dosing regimens or cycling different treatments before progression. This will require clinical trial designs different from the conventional practice of evaluating responses to individual therapy regimens. Next-generation sequencing to assess clonal dynamics longitudinally will improve current radiological assessment of clinical response/resistance and be incorporated into trials exploiting evolution. Furthermore, if understood, clonal evolution can be used to therapeutic advantage, improving patient outcomes based on a new generation of clinical trials.

Bad to the Bone: Emerging Approaches to Aggressive Bone Sarcomas

Bone sarcomas are rare heterogeneous tumors that affect patients of all ages including children, adolescent young adults, and older adults. They include many aggressive subtypes and patient groups with poor outcomes, poor access to clinical trials, and lack of defined standard therapeutic strategies. Conventional chondrosarcoma remains a surgical disease, with no defined role for cytotoxic therapy and no approved targeted systemic therapies. Here, we discuss promising novel targets and strategies undergoing evaluation in clinical trials. Multiagent chemotherapy has greatly improved outcomes for patients with Ewing sarcoma (ES) and osteosarcoma, but management of those with high-risk or recurrent disease remains challenging and controversial. We describe the impact of international collaborative trials, such as the rEECur study, that aim to define optimal treatment strategies for those with recurrent, refractory ES, and evidence for high-dose chemotherapy with stem-cell support. We also discuss current and emerging strategies for other small round cell sarcomas, such as CIC-rearranged, BCOR-rearranged tumors, and the evaluation of emerging novel therapeutics and trial designs that may offer a new paradigm to improve survival in these aggressive tumors with notoriously bad (to the bone) outcomes.

A survey of open questions in adaptive therapy: Bridging mathematics and clinical translation

Adaptive therapy is a dynamic cancer treatment protocol that updates (or 'adapts') treatment decisions in anticipation of evolving tumor dynamics. This broad term encompasses many possible dynamic treatment protocols of patient-specific dose modulation or dose timing. Adaptive therapy maintains high levels of tumor burden to benefit from the competitive suppression of treatment-sensitive subpopulations on treatment-resistant subpopulations. This evolution-based approach to cancer treatment has been integrated into several ongoing or planned clinical trials, including treatment of metastatic castrate resistant prostate cancer, ovarian cancer, and BRAF-mutant melanoma. In the previous few decades, experimental and clinical investigation of adaptive therapy has progressed synergistically with mathematical and computational modeling. In this work, we discuss 11 open questions in cancer adaptive therapy mathematical modeling. The questions are split into three sections: (1) integrating the appropriate components into mathematical models (2) design and validation of dosing protocols, and (3) challenges and opportunities in clinical translation.

Modeling cancer's ecological and evolutionary dynamics

In this didactic paper, we present a theoretical modeling framework, called the G-function, that integrates both the ecology and evolution of cancer to understand oncogenesis. The G-function has been used in evolutionary ecology, but has not been widely applied to problems in cancer. Here, we build the G-function framework from fundamental Darwinian principles and discuss how cancer can be seen through the lens of ecology, evolution, and game theory. We begin with a simple model of cancer growth and add on components of cancer cell competition and drug resistance. To aid in exploration of eco-evolutionary modeling with this approach, we also present a user-friendly software tool. By the end of this paper, we hope that readers will be able to construct basic G function models and grasp the usefulness of the framework to understand the games cancer plays in a biologically mechanistic fashion.

Three unconventional maxims in the natural selection of cancer cells: Generation of induced tumor-suppressing cells (iTSCs)

Induced tumor-suppressing cells (iTSCs) can be generated from cancer and non-cancer cells. Here, three paradoxical maxims for the action of iTSCs are reviewed: the secretion of tumor-suppressing proteins, their role as a "double-edged" sword, and the elimination of lesser-fit cancer cells. "Super-fit" cancer cells secrete an array of proteins, most of which contribute to enhancing their growth and removing "lesser-fit" cancer cells. These maxims explain the potential dilemma with therapeutic agents since the inhibitory agents tend to promote the synthesis of tumor-promoting proteins. The maxims suggest the possibility of a novel treatment option using cancer-guided evolutionary-fit iTSCs.

The CMG helicase and cancer: a tumor "engine" and weakness with missing mutations

The replicative Cdc45-MCM-GINS (CMG) helicase is a large protein complex that functions in the DNA melting and unwinding steps as a component of replisomes during DNA replication in mammalian cells. Although the CMG performs this important role in cell growth, the CMG is not a simple bystander in cell cycle events. Components of the CMG, specifically the MCM precursors, are also involved in maintaining genomic stability by regulating DNA replication fork speeds, facilitating recovery from replicative stresses, and preventing consequential DNA damage. Given these important functions, MCM/CMG complexes are highly regulated by growth factors such as TGF-ß1 and by signaling factors such as Myc, Cyclin E, and the retinoblastoma protein. Mismanagement of MCM/CMG complexes when these signaling mediators are deregulated, and in the absence of the tumor suppressor protein p53, leads to increased genomic instability and is a contributor to tumorigenic transformation and tumor heterogeneity. The goal of this review is to provide insight into the mechanisms and dynamics by which the CMG is regulated during its assembly and activation in mammalian genomes, and how errors in CMG regulation due to oncogenic changes promote tumorigenesis. Finally, and most importantly, we highlight the emerging understanding of the CMG helicase as an exploitable vulnerability and novel target for therapeutic intervention in cancer.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

A Systems Biology Approach for Addressing Cisplatin Resistance in Non-Small Cell Lung Cancer

Translational research in medicine, defined as the transfer of knowledge and discovery from the basic sciences to the clinic, is typically achieved through interactions between members across scientific disciplines to overcome the traditional silos within the community. Thus, translational medicine underscores 'Team Medicine', the partnership between basic science researchers and clinicians focused on addressing a specific goal in medicine. Here, we highlight this concept from a City of Hope perspective. Using cisplatin resistance in non-small cell lung cancer (NSCLC) as a paradigm, we describe how basic research scientists, clinical research scientists, and medical oncologists, in true 'Team Science' spirit, addressed cisplatin resistance in NSCLC and identified a previously approved compound that is able to alleviate cisplatin resistance in NSCLC. Furthermore, we discuss how a 'Team Medicine' approach can help to elucidate the mechanisms of innate and acquired resistance in NSCLC and develop alternative strategies to overcome drug resistance.

Addressing Drug Resistance in Cancer: A Team Medicine Approach

Drug resistance remains one of the major impediments to treating cancer. Although many patients respond well initially, resistance to therapy typically ensues. Several confounding factors appear to contribute to this challenge. Here, we first discuss some of the challenges associated with drug resistance. We then discuss how a 'Team Medicine' approach, involving an interdisciplinary team of basic scientists working together with clinicians, has uncovered new therapeutic strategies. These strategies, referred to as intermittent or 'adaptive' therapy, which are based on eco-evolutionary principles, have met with remarkable success in potentially precluding or delaying the emergence of drug resistance in several cancers. Incorporating such treatment strategies into clinical protocols could potentially enhance the precision of delivering personalized medicine to patients. Furthermore, reaching out to patients in the network of hospitals affiliated with leading academic centers could help them benefit from such innovative treatment options. Finally, lowering the dose of the drug and its frequency (because of intermittent rather than continuous therapy) can also have a significant impact on lowering the toxicity and undesirable side effects of the drugs while lowering the financial burden carried by the patient and insurance providers.

The role of evolutionary game theory in spatial and non-spatial models of the survival of cooperation in cancer: a review

Evolutionary game theory (EGT) is a branch of mathematics which considers populations of individuals interacting with each other to receive pay-offs. An individual’s pay-off is dependent on the strategy of its opponent(s) as well as on its own, and the higher its pay-off, the higher its reproductive fitness. Its offspring generally inherit its interaction strategy, subject to random mutation. Over time, the composition of the population shifts as different strategies spread or are driven extinct. In the last 25 years there has been a flood of interest in applying EGT to cancer modelling, with the aim of explaining how cancerous mutations spread through healthy tissue and how intercellular cooperation persists in tumour-cell populations. This review traces this body of work from theoretical analyses of well-mixed infinite populations through to more realistic spatial models of the development of cooperation between epithelial cells. We also consider work in which EGT has been used to make experimental predictions about the evolution of cancer, and discuss work that remains to be done before EGT can make large-scale contributions to clinical treatment and patient outcomes.

Evolution-based mathematical models significantly prolong response to abiraterone in metastatic castrate-resistant prostate cancer and identify strategies to further improve outcomes

Background

Abiraterone acetate is an effective treatment for metastatic castrate-resistant prostate cancer (mCRPC), but evolution of resistance inevitably leads to progression. We present a pilot study in which abiraterone dosing is guided by evolution-informed mathematical models to delay onset of resistance.

Methods

In the study cohort, abiraterone was stopped when PSA was <50% of pretreatment value and resumed when PSA returned to baseline. Results are compared to a contemporaneous cohort who had >50% PSA decline after initial abiraterone administration and met trial eligibility requirements but chose standard of care (SOC) dosing.

Results

17 subjects were enrolled in the adaptive therapy group and 16 in the SOC group. All SOC subjects have progressed, but four patients in the study cohort remain stably cycling (range 53-70 months). The study cohort had significantly improved median time to progression (TTP; 33.5 months; p<0.001) and median overall survival (OS; 58.5 months; hazard ratio, 0.41, 95% confidence interval (CI), 0.20-0.83, p<0.001) compared to 14.3 and 31.3 months in the SOC cohort. On average, study subjects received no abiraterone during 46% of time on trial. Longitudinal trial data demonstrated the competition coefficient ratio (αRS/αSR) of sensitive and resistant populations, a critical factor in intratumoral evolution, was two- to threefold higher than pre-trial estimates. Computer simulations of intratumoral evolutionary dynamics in the four long-term survivors found that, due to the larger value for αRS/αSR, cycled therapy significantly decreased the resistant population. Simulations in subjects who progressed predicted further increases in OS could be achieved with prompt abiraterone withdrawal after achieving 50% PSA reduction.

Conclusions

Incorporation of evolution-based mathematical models into abiraterone monotherapy for mCRPC significantly increases TTP and OS. Computer simulations with updated parameters from longitudinal trial data can estimate intratumoral evolutionary dynamics in each subject and identify strategies to improve outcomes.

Funding

Moffitt internal grants and NIH/NCI U54CA143970-05 (Physical Science Oncology Network).

Prognostic value of p16, p53, and pcna in sarcoma and an evaluation of immune infiltration

Background

p16, p53, and proliferating cell nuclear antigen (pcna) genes play significant roles in many chromatin modifications and have been found to be highly expressed in a variety of tumor tissues. Therefore, they have been used as target genes for some tumor therapies. However, the differential expressions of the p16, p53, and pcna genes in human sarcomas and their effects on prognosis have not been widely reported.

Methods

The Oncomine dataset was used to analyze the transcription levels of p16, p53, and pcna genes, and the gene expression profile interactive analysis (GEPIA) dataset was used to analyze the differential expressions of p16, p53, and pcna. The expression levels of p16, p53, and pcna were further analyzed by Western Blotting. GEPIA and Kaplan–Meier analyses were used to analyze the prognostic value of p16, p53, and pcna. Furthermore, p16, p53, and pcna gene mutations and their association with overall survival (OS) and disease-free survival (DFS) were analyzed using cBioPortal datasets. In addition, genes co-expressed with p16, p53, and pcna were analyzed using Oncomine. The DAVID dataset was used to analyze the functional enrichment of p16, p53, pcna, and their co-expressed genes by Gene Ontology (GO) and Metascape were used to construct a network map. Finally, the immune cell infiltration of p16, p53, and pcna in patients with sarcoma was reported by Tumor Immune Estimation Resource (TIMER).

Results

p16, p53, and pcna were up-regulated in human sarcoma tissues and almost all sarcoma cell lines. Western Blotting showed that the expression of p16, p53, and pcna was elevated in osteosarcoma cell lines. The expression of pcna was correlated with OS, the expression of p16, p53, and pcna was correlated with relapse-free survival, and the genetic mutation of p16 was negatively correlated with OS and DFS. We also found that p16, p53, and pcna genes were positively/negatively correlated with immune cell infiltration in sarcoma.

Conclusions

The results of this study showed that p16, p53, and pcna can significantly affect the survival and immune status of sarcoma patients. Therefore, p16, p53, and pcna could be used as potential biomarkers of prognosis and immune infiltration in human sarcoma and provide a possible therapeutic target for sarcoma.

A Proposed Trial Design for the Treatment of Widely Metastatic Ewing Sarcoma Inspired by Evolutionary Dynamics

Metastatic Ewing sarcoma has dismal long-term survival despite multiple attempts to intensify standard therapy through the addition of new agents to the existing chemotherapy backbone. Here, based on the application of evolutionary dynamics to pediatric sarcoma, we propose an alternative treatment strategy that varies exposure to agents and dosing intensities, termed sequential second-strike therapy (SSST). We announce an upcoming clinical trial to apply these principles to patients with widely metastatic Ewing sarcoma, those with metastatic disease beyond the lungs.

Optimizing Adaptive Therapy Based on the Reachability to Tumor Resistant Subpopulation

Simple Summary

The intra-competition among tumor subpopulations is a promising target to modify and control the outgrowth of the resistant subpopulation. Adaptive therapy lives up to this principle well, but the gain of tumors with an aggressive resistant subpopulation is not superior to maximum tolerated dose therapy (MTD). How to integrate these two therapies to maximize the outcome? According to the model and system reachability, the ‘restore index’ is proposed to evaluate the timing of the transition from the treatment cycle of adaptive therapy to high-frequency administration, and to juggle the benefits of intra-competition and killing of the sensitive subpopulation. Based on the simulation and animal experiment, the effectiveness of this method in treating tumors with an aggressive resistant subpopulation has been confirmed.

Abstract

Adaptive therapy exploits the self-organization of tumor cells to delay the outgrowth of resistant subpopulations successfully. When the tumor has aggressive resistant subpopulations, the outcome of adaptive therapy was not superior to maximum tolerated dose therapy (MTD). To explore methods to improve the adaptive therapy’s performance of this case, the tumor system was constructed by osimertinib-sensitive and resistant cell lines and illustrated by the Lotka-Volterra model in this study. Restore index proposed to assess the system reachability can predict the duration of each treatment cycle. Then the threshold of the restore index was estimated to evaluate the timing of interrupting the treatment cycle and switching to high-frequency administration. The introduced reachability-based adaptive therapy and classic adaptive therapy were compared through simulation and animal experiments. The results suggested that reachability-based adaptive therapy showed advantages when the tumor has an aggressive resistant subpopulation. This study provides a feasible method for evaluating whether to continue the adaptive therapy treatment cycle or switch to high-frequency administration. This method improves the gain of adaptive therapy by taking into account the benefits of tumor intra-competition and the tumor control of killing sensitive subpopulation.

The Contribution of Evolutionary Game Theory to Understanding and Treating Cancer

Evolutionary game theory mathematically conceptualizes and analyzes biological interactions where one's fitness not only depends on one's own traits, but also on the traits of others. Typically, the individuals are not overtly rational and do not select, but rather inherit their traits. Cancer can be framed as such an evolutionary game, as it is composed of cells of heterogeneous types undergoing frequency-dependent selection. In this article, we first summarize existing works where evolutionary game theory has been employed in modeling cancer and improving its treatment. Some of these game-theoretic models suggest how one could anticipate and steer cancer's eco-evolutionary dynamics into states more desirable for the patient via evolutionary therapies. Such therapies offer great promise for increasing patient survival and decreasing drug toxicity, as demonstrated by some recent studies and clinical trials. We discuss clinical relevance of the existing game-theoretic models of cancer and its treatment, and opportunities for future applications. Moreover, we discuss the developments in cancer biology that are needed to better utilize the full potential of game-theoretic models. Ultimately, we demonstrate that viewing tumors with evolutionary game theory has medically useful implications that can inform and create a lockstep between empirical findings and mathematical modeling. We suggest that cancer progression is an evolutionary competition between different cell types and therefore needs to be viewed as an evolutionary game.

Evolutionary Dynamics of Treatment-Induced Resistance in Cancer Informs Understanding of Rapid Evolution in Natural Systems

Rapid evolution is ubiquitous in nature. We briefly review some of this quite broadly, particularly in the context of response to anthropogenic disturbances. Nowhere is this more evident, replicated and accessible to study than in cancer. Curiously cancer has been late - relative to fisheries, antibiotic resistance, pest management and evolution in human dominated landscapes - in recognizing the need for evolutionarily informed management strategies. The speed of evolution matters. Here, we employ game-theoretic modeling to compare time to progression with continuous maximum tolerable dose to that of adaptive therapy where treatment is discontinued when the population of cancer cells gets below half of its initial size and re-administered when the cancer cells recover, forming cycles with and without treatment. We show that the success of adaptive therapy relative to continuous maximum tolerable dose therapy is much higher if the population of cancer cells is defined by two cell types (sensitive vs. resistant in a polymorphic population). Additionally, the relative increase in time to progression increases with the speed of evolution. These results hold with and without a cost of resistance in cancer cells. On the other hand, treatment-induced resistance can be modeled as a quantitative trait in a monomorphic population of cancer cells. In that case, when evolution is rapid, there is no advantage to adaptive therapy. Initial responses to therapy are blunted by the cancer cells evolving too quickly. Our study emphasizes how cancer provides a unique system for studying rapid evolutionary changes within tumor ecosystems in response to human interventions; and allows us to contrast and compare this system to other human managed or dominated systems in nature.

Starvation and Climate Change—How to Constrain Cancer Cell Epigenetic Diversity and Adaptability to Enhance Treatment Efficacy

Advanced metastatic cancer is currently not curable and the major barrier to eliminating the disease in patients is the resistance of subpopulations of tumor cells to drug treatments. These resistant subpopulations can arise stochastically among the billions of tumor cells in a patient or emerge over time during therapy due to adaptive mechanisms and the selective pressures of drug therapies. Epigenetic mechanisms play important roles in tumor cell diversity and adaptability, and are regulated by metabolic pathways. Here, I discuss knowledge from ecology, evolution, infectious disease, species extinction, metabolism and epigenetics to synthesize a roadmap to a clinically feasible approach to help homogenize tumor cells and, in combination with drug treatments, drive their extinction. Specifically, cycles of starvation and hyperthermia could help synchronize tumor cells and constrain epigenetic diversity and adaptability by limiting substrates and impairing the activity of chromatin modifying enzymes. Hyperthermia could also help prevent cancer cells from entering dangerous hibernation-like states. I propose steps to a treatment paradigm to help drive cancer extinction that builds on the successes of fasting, hyperthermia and immunotherapy and is achievable in patients. Finally, I highlight the many unknowns, opportunities for discovery and that stochastic gene and allele level epigenetic mechanisms pose a major barrier to cancer extinction that warrants deeper investigation.

Treatment-induced evolutionary dynamics in nonmetastatic locally advanced rectal adenocarcinoma

Multi-modal treatment of non-metastatic locally advanced rectal adenocarcinoma (LARC) includes chemotherapy, radiation, and life-altering surgery. Although highly effective for local cancer control, metastatic failure remains significant and drives rectal cancer-related mortality. A consistent observation of this tri-modality treatment paradigm is that histologic response of the primary tumor to neoadjuvant treatment(s), which varies across patients, predicts overall oncologic outcome. In this chapter, we will examine this treatment response heterogeneity in the context of evolutionary dynamics. We hypothesize that improved understanding of eco-evolutionary pressures rendering small cancer cell populations vulnerable to extinction may influence treatment strategies and improve patient outcomes. Applying effective treatment(s) to cancer populations causes a "race to extinction." We explore principles of eco-evolutionary extinction in the context of these small cancer cell populations, evaluating how treatment(s) aim to eradicate the cancer populations to ultimately result in cure. In this chapter, we provide an evolutionary rationale for limiting continuous treatment(s) with the same agent or combination of agents to avoid selection of resistant cancer subpopulation phenotypes, allowing "evolutionary rescue." We draw upon evidence from nature demonstrating species extinction rarely occurring as a single event phenomenon, but rather a series of events in the slide to extinction. We posit that eradicating small cancer populations, similar to small populations in natural extinctions, will usually require a sequence of different external perturbations that produce negative, synergistic dynamics termed the "extinction vortex." By exploiting these unique extinction vulnerabilities of small cancer populations, the optimal therapeutic sequences may be informed by evolution-informed strategies for patients with LARC.

Deep-learning and MR images to target hypoxic habitats with evofosfamide in preclinical models of sarcoma

Rationale: Hypoxic regions (habitats) within tumors are heterogeneously distributed and can be widely variant. Hypoxic habitats are generally pan-therapy resistant. For this reason, hypoxia-activated prodrugs (HAPs) have been developed to target these resistant volumes. The HAP evofosfamide (TH-302) has shown promise in preclinical and early clinical trials of sarcoma. However, in a phase III clinical trial of non-resectable soft tissue sarcomas, TH-302 did not improve survival in combination with doxorubicin (Dox), possibly due to a lack of patient stratification based on hypoxic status. Therefore, we used magnetic resonance imaging (MRI) to identify hypoxic habitats and non-invasively follow therapies response in sarcoma mouse models. Methods: We developed deep-learning (DL) models to identify hypoxia, using multiparametric MRI and co-registered histology, and monitored response to TH-302 in a patient-derived xenograft (PDX) of rhabdomyosarcoma and a syngeneic model of fibrosarcoma (radiation-induced fibrosarcoma, RIF-1). Results: A DL convolutional neural network showed strong correlations (>0.76) between the true hypoxia fraction in histology and the predicted hypoxia fraction in multiparametric MRI. TH-302 monotherapy or in combination with Dox delayed tumor growth and increased survival in the hypoxic PDX model (p<0.05), but not in the RIF-1 model, which had a lower volume of hypoxic habitats. Control studies showed that RIF-1 resistance was due to hypoxia and not other causes. Notably, PDX tumors developed resistance to TH-302 under prolonged treatment that was not due to a reduction in hypoxic volumes. Conclusion: Artificial intelligence analysis of pre-therapy MR images can predict hypoxia and subsequent response to HAPs. This approach can be used to monitor therapy response and adapt schedules to forestall the emergence of resistance.

The Evolution and Ecology of Resistance in Cancer Therapy

Despite the continuous deployment of new treatment strategies and agents over many decades, most disseminated cancers remain fatal. Cancer cells, through their access to the vast information of the human genome, have a remarkable capacity to deploy adaptive strategies for even the most effective treatments. We note there are two critical steps in the clinical manifestation of treatment resistance. The first, which is widely investigated, requires molecular machinery necessary to eliminate the cytotoxic effect of the treatment. However, the emergence of a resistant phenotype is not in itself clinically significant. That is, resistant cells affect patient outcomes only when they succeed in the second step of resistance by proliferating into a sufficiently large population to allow tumor progression and treatment failure. Importantly, proliferation of the resistant phenotype is by no means certain and, in fact, depends on complex Darwinian dynamics governed by the costs and benefits of the resistance mechanisms in the context of the local environment and competing populations. Attempts to target the molecular machinery of resistance have had little clinical success largely because of the diversity within the human genome-therapeutic interruption of one mechanism simply results in its replacement by an alternative. Here we explore evolutionarily informed strategies (adaptive, double-bind, and extinction therapies) for overcoming treatment resistance that seek to understand and exploit the critical evolutionary dynamics that govern proliferation of the resistant phenotypes. In general, this approach has demonstrated that, while emergence of resistance mechanisms in cancer cells to every current therapy is inevitable, proliferation of the resistant phenotypes is not and can be delayed and even prevented with sufficient understanding of the underlying eco-evolutionary dynamics.