RESUMEN
Cancer progression is an evolutionary process driven by the selection of cells adapted to gain growth advantage. We present the first formal study on the adaptation of gene expression in subclonal evolution. We model evolutionary changes in gene expression as stochastic Ornstein-Uhlenbeck processes, jointly leveraging the evolutionary history of subclones and single-cell expression data. Applying our model to sublines derived from single cells of a mouse melanoma revealed that sublines with distinct phenotypes are underlined by different patterns of gene expression adaptation, indicating non-genetic mechanisms of cancer evolution. Interestingly, sublines previously observed to be resistant to anti-CTLA-4 treatment showed adaptive expression of genes related to invasion and non-canonical Wnt signaling, whereas sublines that responded to treatment showed adaptive expression of genes related to proliferation and canonical Wnt signaling. Our results suggest that clonal phenotypes emerge as the result of specific adaptivity patterns of gene expression.
RESUMEN
Despite the promising results of immune checkpoint blockade (ICB) therapy, outcomes for patients with brain metastasis (BrM) remain poor. Identifying resistance mechanisms has been hindered by limited access to patient samples and relevant preclinical models. Here, we developed two mouse melanoma BrM models that recapitulate the disparate responses to ICB seen in patients. We demonstrate that these models capture the cellular and molecular complexity of human disease and reveal key factors shaping the tumor microenvironment and influencing ICB response. BR1-responsive tumor cells express inflammatory programs that polarize microglia into reactive states, eliciting robust T cell recruitment. In contrast, BR3-resistant melanoma cells are enriched in neurological programs and exploit tolerance mechanisms to maintain microglia homeostasis and limit T cell infiltration. In humans, BR1 and BR3 expression signatures correlate positively or negatively with T cell infiltration and BrM patient outcomes, respectively. Our study provides clinically relevant models and uncovers mechanistic insights into BrM ICB responses, offering potential biomarkers and therapeutic targets to improve therapy efficacy.
RESUMEN
Intratumoral heterogeneity (ITH) can promote cancer progression and treatment failure, but the complexity of the regulatory programs and contextual factors involved complicates its study. To understand the specific contribution of ITH to immune checkpoint blockade (ICB) response, we generated single cell-derived clonal sublines from an ICB-sensitive and genetically and phenotypically heterogeneous mouse melanoma model, M4. Genomic and single cell transcriptomic analyses uncovered the diversity of the sublines and evidenced their plasticity. Moreover, a wide range of tumor growth kinetics were observed in vivo , in part associated with mutational profiles and dependent on T cell-response. Further inquiry into melanoma differentiation states and tumor microenvironment (TME) subtypes of untreated tumors from the clonal sublines demonstrated correlations between highly inflamed and differentiated phenotypes with the response to anti-CTLA-4 treatment. Our results demonstrate that M4 sublines generate intratumoral heterogeneity at both levels of intrinsic differentiation status and extrinsic TME profiles, thereby impacting tumor evolution during therapeutic treatment. These clonal sublines proved to be a valuable resource to study the complex determinants of response to ICB, and specifically the role of melanoma plasticity in immune evasion mechanisms.