The determination of individual cell trajectories through a high-dimensional cell-state space is an outstanding challenge for understanding biological changes ranging from cellular differentiation to epigenetic responses of diseased cells upon drugging. We integrate experiments and theory to determine the trajectories that single BRAFV600E mutant melanoma cancer cells take between drug-naive and drug-tolerant states. Although single-cell omics tools can yield snapshots of the cell-state landscape, the determination of individual cell trajectories through that space can be confounded by stochastic cell-state switching. We assayed for a panel of signaling, phenotypic, and metabolic regulators at points across 5 days of drug treatment to uncover a cell-state landscape with two paths connecting drug-naive and drug-tolerant states. The trajectory a given cell takes depends upon the drug-naive level of a lineage-restricted transcription factor. Each trajectory exhibits unique druggable susceptibilities, thus updating the paradigm of adaptive resistance development in an isogenic cell population.

Abstract Integrated proteomic and metabolic single-cell assays reveal multiple independent adaptive responses to drug tolerance in a BRAF-mutant melanoma cell line Cancers commonly develop resistance against chemotherapeutics or targeted therapies through various types of genetic or non-genetic mechanisms. Non-genetic mechanisms have been shown to occur early on and can provide a latent reservoir of cells for the emergence of various different type of mechanisms, yet very limited understanding of process were resolve main from bulk analysis. Considering the heterogeneous nature of the tumor cells, a single-cell level characterization of the process worth detailed further investigation. Using MAPK inhibition of BRAF-mutant melanomas as a model system, we resolved that cells take different paths to go from drug-sensitive to drug-resistant state. Using a microfludic-based single-cell integrated proteomic and metabolic assay, we assayed for a panel of signaling, phenotypic, and metabolic regulators at four time points during the first five days of drug treatment. Dimensional reduction of the resultant data set, coupled with information theoretic analysis, uncovered a complex cell state landscape and identified two distinct paths connecting drug-naïve and drug-tolerant states. Cells are shown to exclusively traverse one of the two pathways depending on the level of the lineage restricted transcription factor MITF in the drug-naïve cells. The two trajectories are associated with distinct signaling and metabolic susceptibilities, and are independently druggable. Our results update the paradigm of adaptive resistance development in an isogenic cell population and offer insight into the design of more effective combination therapies. Citation Format: Yapeng Su, Guideng Li, Melissa Ko, Hanjun Cheng, Ronghui Zhu, Min Xue, Lidia Robert, Raphael Levine, Antoni Ribas, Garry Nolan, Wei Wei, Sylvia Plevritis, David Baltimore, James R. Heath. Systems biology for investigating drug resistance mechanism of melanoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 6585.

We resolved a mechanism connecting tumor epigenetic plasticity with non-genetic adaptive resistance to therapy, with MAPK inhibition of BRAF-mutant melanomas providing the model. These cancer cells undergo multiple, reversible drug-induced cell-state transitions, ultimately yielding a drug-resistant mesenchymal-like phenotype. A kinetic series of transcriptome and epigenome data, collected over two months of drug treatment and release, revealed changing levels of thousands of genes and extensive chromatin remodeling. However, a 3-step computational algorithm greatly simplified the interpretation of these changes, and revealed that the whole adaptive process was controlled by a gene module activated within just three days of treatment, with RelA driving chromatin remodeling to establish an epigenetic program encoding long-term phenotype changes. These findings were confirmed across several patient-derived cell lines and in melanoma patients under MAPK inhibitor treatment. Co-targeting BRAF and histone-modifying enzymes arrests adaptive transitions towards drug tolerance in epigenetically plastic melanoma cells and may be exploited therapeutically.

The determination of individual cell trajectories through a high-dimensional cell-state space is an outstanding challenge, with relevance towards understanding biological changes ranging from cellular differentiation to epigenetic (adaptive) responses of diseased cells to drugging. We report on a combined experimental and theoretic method for determining the trajectories that specific highly plastic BRAF V600E mutant patient-derived melanoma cancer cells take between drug-naïve and drug-tolerant states. Recent studies have implicated non-genetic, fast-acting resistance mechanisms are activated in these cells following BRAF inhibition. While single-cell highly multiplex omics tools can yield snapshots of the cell state space landscape sampled at any given time point, individual cell trajectories must be inferred from a kinetic series of snapshots, and that inference can be confounded by stochastic cell state switching. Using a microfludic-based single-cell integrated proteomic and metabolic assay, we assayed for a panel of signaling, phenotypic, and metabolic regulators at four time points during the first five days of drug treatment. Dimensional reduction of the resultant data set, coupled with information theoretic analysis, uncovered a complex cell state landscape and identified two distinct paths connecting drug-naïve and drug-tolerant states. Cells are shown to exclusively traverse one of the two pathways depending on the level of the lineage restricted transcription factor MITF in the drug-naïve cells. The two trajectories are associated with distinct signaling and metabolic susceptibilities, and are independently druggable. Our results update the paradigm of adaptive resistance development in an isogenic cell population and offer insight into the design of more effective combination therapies.