Purinergic Ca2+ signaling as a novel mechanism of drug tolerance in BRAF mutant melanoma.

Drug tolerance is a major cause of relapse after cancer treatment. In spite of intensive efforts1-9, its molecular basis remains poorly understood, hampering actionable intervention. We report a previously unrecognized signaling mechanism supporting drug tolerance in BRAF-mutant melanoma treated with BRAF inhibitors that could be of general relevance to other cancers. Its key features are cell-intrinsic intracellular Ca2+ signaling initiated by P2X7 receptors (purinergic ligand-gated cation channels), and an enhanced ability for these Ca2+ signals to reactivate ERK1/2 in the drug-tolerant state. Extracellular ATP, virtually ubiquitous in living systems, is the ligand that can initiate Ca2+ spikes via P2X7 channels. ATP is abundant in the tumor microenvironment and is released by dying cells, ironically implicating treatment-initiated cancer cell death as a source of trophic stimuli that leads to ERK reactivation and drug tolerance. Such a mechanism immediately offers an explanation of the inevitable relapse after BRAFi treatment in BRAF-mutant melanoma, and points to actionable strategies to overcome it.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics.

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

Unified tumor growth mechanisms from multimodel inference and dataset integration.

Mechanistic models of biological processes can explain observed phenomena and predict responses to a perturbation. A mathematical model is typically constructed using expert knowledge and informal reasoning to generate a mechanistic explanation for a given observation. Although this approach works well for simple systems with abundant data and well-established principles, quantitative biology is often faced with a dearth of both data and knowledge about a process, thus making it challenging to identify and validate all possible mechanistic hypothesis underlying a system behavior. To overcome these limitations, we introduce a Bayesian multimodel inference (Bayes-MMI) methodology, which quantifies how mechanistic hypotheses can explain a given experimental datasets, and concurrently, how each dataset informs a given model hypothesis, thus enabling hypothesis space exploration in the context of available data. We demonstrate this approach to probe standing questions about heterogeneity, lineage plasticity, and cell-cell interactions in tumor growth mechanisms of small cell lung cancer (SCLC). We integrate three datasets that each formulated different explanations for tumor growth mechanisms in SCLC, apply Bayes-MMI and find that the data supports model predictions for tumor evolution promoted by high lineage plasticity, rather than through expanding rare stem-like populations. In addition, the models predict that in the presence of cells associated with the SCLC-N or SCLC-A2 subtypes, the transition from the SCLC-A subtype to the SCLC-Y subtype through an intermediate is decelerated. Together, these predictions provide a testable hypothesis for observed juxtaposed results in SCLC growth and a mechanistic interpretation for tumor treatment resistance.

Involvement of Epithelial-Mesenchymal Transition Genes in Small Cell Lung Cancer Phenotypic Plasticity.

Small cell lung cancer (SCLC) is an aggressive cancer recalcitrant to treatment, arising predominantly from epithelial pulmonary neuroendocrine (NE) cells. Intratumor heterogeneity plays critical roles in SCLC disease progression, metastasis, and treatment resistance. At least five transcriptional SCLC NE and non-NE cell subtypes were recently defined by gene expression signatures. Transition from NE to non-NE cell states and cooperation between subtypes within a tumor likely contribute to SCLC progression by mechanisms of adaptation to perturbations. Therefore, gene regulatory programs distinguishing SCLC subtypes or promoting transitions are of great interest. Here, we systematically analyze the relationship between SCLC NE/non-NE transition and epithelial to mesenchymal transition (EMT)-a well-studied cellular process contributing to cancer invasiveness and resistance-using multiple transcriptome datasets from SCLC mouse tumor models, human cancer cell lines, and tumor samples. The NE SCLC-A2 subtype maps to the epithelial state. In contrast, SCLC-A and SCLC-N (NE) map to a partial mesenchymal state (M1) that is distinct from the non-NE, partial mesenchymal state (M2). The correspondence between SCLC subtypes and the EMT program paves the way for further work to understand gene regulatory mechanisms of SCLC tumor plasticity with applicability to other cancer types.

Cooperative RNA degradation stabilizes intermediate epithelial-mesenchymal states and supports a phenotypic continuum.

Multiple intermediate epithelial-mesenchymal transition (EMT) states reflecting hybrid epithelial and mesenchymal phenotypes were observed in physiological and pathological conditions. Previous theoretical models explaining multiple EMT states rely on regulatory loops involving transcriptional feedback, which produce three or four attractors. This is incompatible with the observed continuum-like EMT spectrum. Here, we used mass-action-based models to describe post-transcriptional regulations, finding that cooperative RNA degradation via multiple microRNA binding sites can generate four-attractor systems without transcriptional feedback. Furthermore, the newly identified intermediates-enabling circuits are common in the EMT regulatory network, and they can synergize with transcriptional feedback to support phenotypic continuum. Finally, our model predicted a role of miR-101 in multistate EMT, and we identified evidence from single-cell RNA-sequencing data that support the prediction. Our work reveals a previously unknown role of cooperative RNA degradation and microRNAs in EMT, providing a framework that can bridge the gap between mechanistic models and single-cell experiments.

MuSyC dosing of adjuvanted cancer vaccines optimizes antitumor responses.

With the clinical approval of T-cell-dependent immune checkpoint inhibitors for many cancers, therapeutic cancer vaccines have re-emerged as a promising immunotherapy. Cancer vaccines require the addition of immunostimulatory adjuvants to increase vaccine immunogenicity, and increasingly multiple adjuvants are used in combination to bolster further and shape cellular immunity to tumor antigens. However, rigorous quantification of adjuvants' synergistic interactions is challenging due to partial redundancy in costimulatory molecules and cytokine production, leading to the common assumption that combining both adjuvants at the maximum tolerated dose results in optimal efficacy. Herein, we examine this maximum dose assumption and find combinations of these doses are suboptimal. Instead, we optimized dendritic cell activation by extending the Multidimensional Synergy of Combinations (MuSyC) framework that measures the synergy of efficacy and potency between two vaccine adjuvants. Initially, we performed a preliminary in vitro screening of clinically translatable adjuvant receptor targets (TLR, STING, NLL, and RIG-I). We determined that STING agonist (CDN) plus TLR4 agonist (MPL-A) or TLR7/8 agonist (R848) as the best pairwise combinations for dendritic cell activation. In addition, we found that the combination of R848 and CDN is synergistically efficacious and potent in activating both murine and human antigen-presenting cells (APCs) in vitro. These two selected adjuvants were then used to estimate a MuSyC-dose optimized for in vivo T-cell priming using ovalbumin-based peptide vaccines. Finally, using B16 melanoma and MOC1 head and neck cancer models, MuSyC-dose-based adjuvating of cancer vaccines improved the antitumor response, increased tumor-infiltrating lymphocytes, and induced novel myeloid tumor infiltration changes. Further, the MuSyC-dose-based adjuvants approach did not cause additional weight changes or increased plasma cytokine levels compared to CDN alone. Collectively, our findings offer a proof of principle that our MuSyC-extended approach can be used to optimize cancer vaccine formulations for immunotherapy.

Archetype tasks link intratumoral heterogeneity to plasticity and cancer hallmarks in small cell lung cancer.

Small cell lung cancer (SCLC) tumors comprise heterogeneous mixtures of cell states, categorized into neuroendocrine (NE) and non-neuroendocrine (non-NE) transcriptional subtypes. NE to non-NE state transitions, fueled by plasticity, likely underlie adaptability to treatment and dismal survival rates. Here, we apply an archetypal analysis to model plasticity by recasting SCLC phenotypic heterogeneity through multi-task evolutionary theory. Cell line and tumor transcriptomics data fit well in a five-dimensional convex polytope whose vertices optimize tasks reminiscent of pulmonary NE cells, the SCLC normal counterparts. These tasks, supported by knowledge and experimental data, include proliferation, slithering, metabolism, secretion, and injury repair, reflecting cancer hallmarks. SCLC subtypes, either at the population or single-cell level, can be positioned in archetypal space by bulk or single-cell transcriptomics, respectively, and characterized as task specialists or multi-task generalists by the distance from archetype vertex signatures. In the archetype space, modeling single-cell plasticity as a Markovian process along an underlying state manifold indicates that task trade-offs, in response to microenvironmental perturbations or treatment, may drive cell plasticity. Stifling phenotypic transitions and plasticity may provide new targets for much-needed translational advances in SCLC. A record of this paper's Transparent Peer Review process is included in the supplemental information.

Activation of STAT3 through combined SRC and EGFR signaling drives resistance to a mitotic kinesin inhibitor in glioblastoma.

Inhibitors of the mitotic kinesin Kif11 are anti-mitotics that, unlike vinca alkaloids or taxanes, do not disrupt microtubules and are not neurotoxic. However, development of resistance has limited their clinical utility. While resistance to Kif11 inhibitors in other cell types is due to mechanisms that prevent these drugs from disrupting mitosis, we find that in glioblastoma (GBM), resistance to the Kif11 inhibitor ispinesib works instead through suppression of apoptosis driven by activation of STAT3. This form of resistance requires dual phosphorylation of STAT3 residues Y705 and S727, mediated by SRC and epidermal growth factor receptor (EGFR), respectively. Simultaneously inhibiting SRC and EGFR reverses this resistance, and combined targeting of these two kinases in vivo with clinically available inhibitors is synergistic and significantly prolongs survival in ispinesib-treated GBM-bearing mice. We thus identify a translationally actionable approach to overcoming Kif11 inhibitor resistance that may work to block STAT3-driven resistance against other anti-cancer therapies as well.

Organoids as a Systems Platform for SCLC Brain Metastasis.

Small Cell Lung Cancer (SCLC) is a highly aggressive, neuroendocrine tumor. Traditional reductionist approaches have proven ineffective to ameliorate the uniformly dismal outcomes for SCLC - survival at 5 years remains less than 5%. A major obstacle to improving treatment is that SCLC tumor cells disseminate early, with a strong propensity for metastasizing to the brain. Accumulating evidence indicates that, contrary to previous textbook knowledge, virtually every SCLC tumor is comprised of multiple subtypes. Important questions persist regarding the role that this intra-tumor subtype heterogeneity may play in supporting the invasive properties of SCLC. A recurrent hypothesis in the field is that subtype interactions and/or transition dynamics are major determinants of SCLC metastatic seeding and progression. Here, we review the advantages of cerebral organoids as an experimentally accessible platform for SCLC brain metastasis, amenable to genetic manipulations, drug perturbations, and assessment of subtype interactions when coupled, e.g., to temporal longitudinal monitoring by high-content imaging or high-throughput omics data generation. We then consider systems approaches that can produce mathematical and computational models useful to generalize lessons learned from ex vivo organoid cultures, and integrate them with in vivo observations. In summary, systems approaches combined with ex vivo SCLC cultures in brain organoids may effectively capture both tumor-tumor and host-tumor interactions that underlie general principles of brain metastasis.

Real-time luminescence enables continuous drug-response analysis in adherent and suspension cell lines.

The drug-induced proliferation (DIP) rate is a metric of in vitro drug response that avoids inherent biases in commonly used metrics such as 72 h viability. However, DIP rate measurements rely on direct cell counting over time, a laborious task that is subject to numerous challenges, including the need to fluorescently label cells and automatically segment nuclei. Moreover, it is incredibly difficult to directly count cells and accurately measure DIP rates for cell populations in suspension. As an alternative, we use real-time luminescence measurements derived from the cellular activity of NAD(P)H oxidoreductase to efficiently estimate drug response in both adherent and suspension cell populations to a panel of known anticancer agents. For the adherent cell lines, we collect both luminescence reads and direct cell counts over time simultaneously to assess their congruency. Our results demonstrate that the proposed approach significantly speeds up data collection, avoids the need for cellular labels and image segmentation, and opens the door to significant advances in high-throughput screening of anticancer drugs.

MuSyC is a consensus framework that unifies multi-drug synergy metrics for combinatorial drug discovery.

Drug combination discovery depends on reliable synergy metrics but no consensus exists on the correct synergy criterion to characterize combined interactions. The fragmented state of the field confounds analysis, impedes reproducibility, and delays clinical translation of potential combination treatments. Here we present a mass-action based formalism to quantify synergy. With this formalism, we clarify the relationship between the dominant drug synergy principles, and present a mapping of commonly used frameworks onto a unified synergy landscape. From this, we show how biases emerge due to intrinsic assumptions which hinder their broad applicability and impact the interpretation of synergy in discovery efforts. Specifically, we describe how traditional metrics mask consequential synergistic interactions, and contain biases dependent on the Hill-slope and maximal effect of single-drugs. We show how these biases systematically impact synergy classification in large combination screens, potentially misleading discovery efforts. Thus the proposed formalism can provide a consistent, unbiased interpretation of drug synergy, and accelerate the translatability of synergy studies.

An in vitro model of tumor heterogeneity resolves genetic, epigenetic, and stochastic sources of cell state variability.

Tumor heterogeneity is a primary cause of treatment failure and acquired resistance in cancer patients. Even in cancers driven by a single mutated oncogene, variability in response to targeted therapies is well known. The existence of additional genomic alterations among tumor cells can only partially explain this variability. As such, nongenetic factors are increasingly seen as critical contributors to tumor relapse and acquired resistance in cancer. Here, we show that both genetic and nongenetic factors contribute to targeted drug response variability in an experimental model of tumor heterogeneity. We observe significant variability to epidermal growth factor receptor (EGFR) inhibition among and within multiple versions and clonal sublines of PC9, a commonly used EGFR mutant nonsmall cell lung cancer (NSCLC) cell line. We resolve genetic, epigenetic, and stochastic components of this variability using a theoretical framework in which distinct genetic states give rise to multiple epigenetic "basins of attraction," across which cells can transition driven by stochastic noise. Using mutational impact analysis, single-cell differential gene expression, and correlations among Gene Ontology (GO) terms to connect genomics to transcriptomics, we establish a baseline for genetic differences driving drug response variability among PC9 cell line versions. Applying the same approach to clonal sublines, we conclude that drug response variability in all but one of the sublines is due to epigenetic differences; in the other, it is due to genetic alterations. Finally, using a clonal drug response assay together with stochastic simulations, we attribute subclonal drug response variability within sublines to stochastic cell fate decisions and confirm that one subline likely contains genetic resistance mutations that emerged in the absence of drug treatment.

ASCL1 represses a SOX9+ neural crest stem-like state in small cell lung cancer.

ASCL1 is a neuroendocrine lineage-specific oncogenic driver of small cell lung cancer (SCLC), highly expressed in a significant fraction of tumors. However, ∼25% of human SCLC are ASCL1-low and associated with low neuroendocrine fate and high MYC expression. Using genetically engineered mouse models (GEMMs), we show that alterations in Rb1/Trp53/Myc in the mouse lung induce an ASCL1+ state of SCLC in multiple cells of origin. Genetic depletion of ASCL1 in MYC-driven SCLC dramatically inhibits tumor initiation and progression to the NEUROD1+ subtype of SCLC. Surprisingly, ASCL1 loss promotes a SOX9+ mesenchymal/neural crest stem-like state and the emergence of osteosarcoma and chondroid tumors, whose propensity is impacted by cell of origin. ASCL1 is critical for expression of key lineage-related transcription factors NKX2-1, FOXA2, and INSM1 and represses genes involved in the Hippo/Wnt/Notch developmental pathways in vivo. Importantly, ASCL1 represses a SOX9/RUNX1/RUNX2 program in vivo and SOX9 expression in human SCLC cells, suggesting a conserved function for ASCL1. Together, in a MYC-driven SCLC model, ASCL1 promotes neuroendocrine fate and represses the emergence of a SOX9+ nonendodermal stem-like fate that resembles neural crest.

Thunor: visualization and analysis of high-throughput dose-response datasets.

High-throughput cell proliferation assays to quantify drug-response are becoming increasingly common and powerful with the emergence of improved automation and multi-time point analysis methods. However, pipelines for analysis of these datasets that provide reproducible, efficient, and interactive visualization and interpretation are sorely lacking. To address this need, we introduce Thunor, an open-source software platform to manage, analyze, and visualize large, dose-dependent cell proliferation datasets. Thunor supports both end-point and time-based proliferation assays as input. It provides a simple, user-friendly interface with interactive plots and publication-quality images of cell proliferation time courses, dose-response curves, and derived dose-response metrics, e.g. IC50, including across datasets or grouped by tags. Tags are categorical labels for cell lines and drugs, used for aggregation, visualization and statistical analysis, e.g. cell line mutation or drug class/target pathway. A graphical plate map tool is included to facilitate plate annotation with cell lines, drugs and concentrations upon data upload. Datasets can be shared with other users via point-and-click access control. We demonstrate the utility of Thunor to examine and gain insight from two large drug response datasets: a large, publicly available cell viability database and an in-house, high-throughput proliferation rate dataset. Thunor is available from www.thunor.net.

Beyond Programmed Death-Ligand 1: B7-H6 Emerges as a Potential Immunotherapy Target in SCLC.

INTRODUCTION: The programmed death-ligand 1 (PD-L1) immune checkpoint inhibitors, atezolizumab and durvalumab, have received regulatory approval for the first-line treatment of patients with extensive-stage SCLC. Nevertheless, when used in combination with platinum-based chemotherapy, these PD-L1 inhibitors only improve overall survival by 2 to 3 months. This may be due to the observation that less than 20% of SCLC tumors express PD-L1 at greater than 1%. Evaluating the composition and abundance of checkpoint molecules in SCLC may identify molecules beyond PD-L1 that are amenable to therapeutic targeting. METHODS: We analyzed RNA-sequencing data from SCLC cell lines (n = 108) and primary tumor specimens (n = 81) for expression of 39 functionally validated inhibitory checkpoint ligands. Furthermore, we generated tissue microarrays containing SCLC cell lines and patient with SCLC specimens to confirm expression of these molecules by immunohistochemistry. We annotated patient outcomes data, including treatment response and overall survival. RESULTS: The checkpoint protein B7-H6 (NCR3LG1) exhibited increased protein expression relative to PD-L1 in cell lines and tumors (p < 0.05). Higher B7-H6 protein expression correlated with longer progression-free survival (p = 0.0368) and increased total immune infiltrates (CD45+) in patients. Furthermore, increased B7-H6 gene expression in SCLC tumors correlated with a decreased activated natural killer cell gene signature, suggesting a complex interplay between B7-H6 expression and immune signature in SCLC. CONCLUSIONS: We investigated 39 inhibitory checkpoint molecules in SCLC and found that B7-H6 is highly expressed and associated with progression-free survival. In addition, 26 of 39 immune checkpoint proteins in SCLC tumors were more abundantly expressed than PD-L1, indicating an urgent need to investigate additional checkpoint targets for therapy in addition to PD-L1.

Patterns of transcription factor programs and immune pathway activation define four major subtypes of SCLC with distinct therapeutic vulnerabilities.

Despite molecular and clinical heterogeneity, small cell lung cancer (SCLC) is treated as a single entity with predictably poor results. Using tumor expression data and non-negative matrix factorization, we identify four SCLC subtypes defined largely by differential expression of transcription factors ASCL1, NEUROD1, and POU2F3 or low expression of all three transcription factor signatures accompanied by an Inflamed gene signature (SCLC-A, N, P, and I, respectively). SCLC-I experiences the greatest benefit from the addition of immunotherapy to chemotherapy, while the other subtypes each have distinct vulnerabilities, including to inhibitors of PARP, Aurora kinases, or BCL-2. Cisplatin treatment of SCLC-A patient-derived xenografts induces intratumoral shifts toward SCLC-I, supporting subtype switching as a mechanism of acquired platinum resistance. We propose that matching baseline tumor subtype to therapy, as well as manipulating subtype switching on therapy, may enhance depth and duration of response for SCLC patients.

An Integrative Gene Expression and Mathematical Flux Balance Analysis Identifies Targetable Redox Vulnerabilities in Melanoma Cells.

Melanomas harboring BRAF mutations can be treated with BRAF inhibitors (BRAFi), but responses are varied and tumor recurrence is inevitable. Here we used an integrative approach of experimentation and mathematical flux balance analyses in BRAF-mutated melanoma cells to discover that elevated antioxidant capacity is linked to BRAFi sensitivity in melanoma cells. High levels of antioxidant metabolites in cells with reduced BRAFi sensitivity confirmed this conclusion. By extending our analyses to other melanoma subtypes in The Cancer Genome Atlas, we predict that elevated redox capacity is a general feature of melanomas, not previously observed. We propose that redox vulnerabilities could be exploited for therapeutic benefits and identify unsuspected combination targets to enhance the effects of BRAFi in any melanoma, regardless of mutational status. SIGNIFICANCE: An integrative bioinformatics, flux balance analysis, and experimental approach identify targetable redox vulnerabilities and show the potential for modulation of cancer antioxidant defense to augment the benefits of existing therapies in melanoma.

Drug-Tolerant Idling Melanoma Cells Exhibit Theory-Predicted Metabolic Low-Low Phenotype.

Cancer cells adjust their metabolic profiles to evade treatment. Metabolic adaptation is complex and hence better understood by an integrated theoretical-experimental approach. Using a minimal kinetic model, we predicted a previously undescribed Low/Low (L/L) phenotype, characterized by low oxidative phosphorylation (OXPHOS) and low glycolysis. Here, we report that L/L metabolism is observed in BRAF-mutated melanoma cells that enter a drug-tolerant "idling state" upon long-term MAPK inhibition (MAPKi). Consistently, using publicly available RNA-sequencing data of both cell lines and patient samples, we show that melanoma cells decrease their glycolysis and/or OXPHOS activity upon MAPKi and converge toward the L/L phenotype. L/L metabolism is unfavorable for tumor growth, yet supports successful cell division at ~50% rate. Thus, L/L drug-tolerant idling cells are a reservoir for accumulating mutations responsible for relapse, and it should be considered as a target subpopulation for improving MAPKi outcomes in melanoma treatment.

Charting the Fragmented Landscape of Drug Synergy.

Even as the clinical impact of drug combinations continues to accelerate, no consensus on how to quantify drug synergy has emerged. Rather, surveying the landscape of drug synergy reveals the persistence of historical fissures regarding the appropriate domains of conflicting synergy models - fissures impacting all aspects of combination therapy discovery and deployment. Herein we chronicle the impact of these divisions on: (i) the design, interpretation, and reproducibility of high-throughput combination screens; (ii) the performance of algorithms to predict synergistic mixtures; and (iii) the search for higher-order synergistic interactions. Further progress in each of these subfields hinges on reaching a consensus regarding the long-standing rifts in the field.

New Approaches to SCLC Therapy: From the Laboratory to the Clinic.

The outcomes of patients with SCLC have not yet been substantially impacted by the revolution in precision oncology, primarily owing to a paucity of genetic alterations in actionable driver oncogenes. Nevertheless, systemic therapies that include immunotherapy are beginning to show promise in the clinic. Although, these results are encouraging, many patients do not respond to, or rapidly recur after, current regimens, necessitating alternative or complementary therapeutic strategies. In this review, we discuss ongoing investigations into the pathobiology of this recalcitrant cancer and the therapeutic vulnerabilities that are exposed by the disease state. Included within this discussion, is a snapshot of the current biomarker and clinical trial landscapes for SCLC. Finally, we identify key knowledge gaps that should be addressed to advance the field in pursuit of reduced SCLC mortality. This review largely summarizes work presented at the Third Biennial International Association for the Study of Lung Cancer SCLC Meeting.