Leveraging transcriptional dynamics to improve BRAF inhibitor responses in melanoma

BACKGROUND: Melanoma is a heterogeneous tumour, but the impact of this heterogeneity upon therapeutic response is not well understood. METHODS: Single cell mRNA analysis was used to define the transcriptional heterogeneity of melanoma and its dynamic response to BRAF inhibitor therapy and treatment holidays. Discrete transcriptional states were defined in cell lines and melanoma patient specimens that predicted initial sensitivity to BRAF inhibition and the potential for effective re-challenge following resistance. A mathematical model was developed to maintain competition between the drug-sensitive and resistant states, which was validated in vivo. FINDINGS: Our analyses showed melanoma cell lines and patient specimens to be composed of >3 transcriptionally distinct states. The cell state composition was dynamically regulated in response to BRAF inhibitor therapy and drug holidays. Transcriptional state composition predicted for therapy response. The differences in fitness between the different transcriptional states were leveraged to develop a mathematical model that optimized therapy schedules to retain the drug sensitive population. In vivo validation demonstrated that the personalized adaptive dosing schedules outperformed continuous or fixed intermittent BRAF inhibitor schedules. INTERPRETATION: Our study provides the first evidence that transcriptional heterogeneity at the single cell level predicts for initial BRAF inhibitor sensitivity. We further demonstrate that manipulating transcriptional heterogeneity through personalized adaptive therapy schedules can delay the time to resistance. FUNDING: This work was funded by the National Institutes of Health. The funder played no role in assembly of the manuscript.