PhISCS-BnB: a fast branch and bound algorithm for the perfect tumor phylogeny reconstruction problem.

MOTIVATION: Recent advances in single-cell sequencing (SCS) offer an unprecedented insight into tumor emergence and evolution. Principled approaches to tumor phylogeny reconstruction via SCS data are typically based on general computational methods for solving an integer linear program, or a constraint satisfaction program, which, although guaranteeing convergence to the most likely solution, are very slow. Others based on Monte Carlo Markov Chain or alternative heuristics not only offer no such guarantee, but also are not faster in practice. As a result, novel methods that can scale up to handle the size and noise characteristics of emerging SCS data are highly desirable to fully utilize this technology. RESULTS: We introduce PhISCS-BnB (phylogeny inference using SCS via branch and bound), a branch and bound algorithm to compute the most likely perfect phylogeny on an input genotype matrix extracted from an SCS dataset. PhISCS-BnB not only offers an optimality guarantee, but is also 10-100 times faster than the best available methods on simulated tumor SCS data. We also applied PhISCS-BnB on a recently published large melanoma dataset derived from the sublineages of a cell line involving 20 clones with 2367 mutations, which returned the optimal tumor phylogeny in <4 h. The resulting phylogeny agrees with and extends the published results by providing a more detailed picture on the clonal evolution of the tumor. AVAILABILITY AND IMPLEMENTATION: https://github.com/algo-cancer/PhISCS-BnB. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Noncoding RNA circuitry in melanoma onset, plasticity, and therapeutic response.

Melanoma, the cancer of the melanocyte, is the deadliest form of skin cancer with an aggressive nature, propensity to metastasize and tendency to resist therapeutic intervention. Studies have identified that the re-emergence of developmental pathways in melanoma contributes to melanoma onset, plasticity, and therapeutic response. Notably, it is well known that noncoding RNAs play a critical role in the development and stress response of tissues. In this review, we focus on the noncoding RNAs, including microRNAs, long non-coding RNAs, circular RNAs, and other small RNAs, for their functions in developmental mechanisms and plasticity, which drive onset, progression, therapeutic response and resistance in melanoma. Going forward, elucidation of noncoding RNA-mediated mechanisms may provide insights that accelerate development of novel melanoma therapies.

Melanoma clonal subline analysis uncovers heterogeneity-driven immunotherapy resistance mechanisms.

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.

PARP14 inhibition restores PD-1 immune checkpoint inhibitor response following IFNγ-driven acquired resistance in preclinical cancer models.

Resistance mechanisms to immune checkpoint blockade therapy (ICBT) limit its response duration and magnitude. Paradoxically, Interferon γ (IFNγ), a key cytokine for cellular immunity, can promote ICBT resistance. Using syngeneic mouse tumour models, we confirm that chronic IFNγ exposure confers resistance to immunotherapy targeting PD-1 (α-PD-1) in immunocompetent female mice. We observe upregulation of poly-ADP ribosyl polymerase 14 (PARP14) in chronic IFNγ-treated cancer cell models, in patient melanoma with elevated IFNG expression, and in melanoma cell cultures from ICBT-progressing lesions characterised by elevated IFNγ signalling. Effector T cell infiltration is enhanced in tumours derived from cells pre-treated with IFNγ in immunocompetent female mice when PARP14 is pharmacologically inhibited or knocked down, while the presence of regulatory T cells is decreased, leading to restoration of α-PD-1 sensitivity. Finally, we determine that tumours which spontaneously relapse in immunocompetent female mice following α-PD-1 therapy upregulate IFNγ signalling and can also be re-sensitised upon receiving PARP14 inhibitor treatment, establishing PARP14 as an actionable target to reverse IFNγ-driven ICBT resistance.

The Hitchhiker's Guide across a Waddington's landscape of melanoma.

Melanoma evolution may recapitulate the embryonic development of its progenitor tissue, neural crest (NC), but the exact process is unclear. In a recent issue of Nature, Karras et al. (2022) demonstrate that melanoma expansion mirrors the hierarchic process of NC differentiation, generating cell subpopulations, each with distinct function, including growth and metastasis.

Fast intratumor heterogeneity inference from single-cell sequencing data.

We introduce HUNTRESS, a computational method for mutational intratumor heterogeneity inference from noisy genotype matrices derived from single-cell sequencing data, the running time of which is linear with the number of cells and quadratic with the number of mutations. We prove that, under reasonable conditions, HUNTRESS computes the true progression history of a tumor with high probability. On simulated and real tumor sequencing data, HUNTRESS is demonstrated to be faster than available alternatives with comparable or better accuracy. Additionally, the progression histories of tumors inferred by HUNTRESS on real single-cell sequencing datasets agree with the best known evolution scenarios for the associated tumors.

Melanoma models for the next generation of therapies.

There is a lack of appropriate melanoma models that can be used to evaluate the efficacy of novel therapeutic modalities. Here, we discuss the current state of the art of melanoma models including genetically engineered mouse, patient-derived xenograft, zebrafish, and ex vivo and in vitro models. We also identify five major challenges that can be addressed using such models, including metastasis and tumor dormancy, drug resistance, the melanoma immune response, and the impact of aging and environmental exposures on melanoma progression and drug resistance. Additionally, we discuss the opportunity for building models for rare subtypes of melanomas, which represent an unmet critical need. Finally, we identify key recommendations for melanoma models that may improve accuracy of preclinical testing and predict efficacy in clinical trials, to help usher in the next generation of melanoma therapies.

PRL3-DDX21 Transcriptional Control of Endolysosomal Genes Restricts Melanocyte Stem Cell Differentiation.

Melanocytes, replenished throughout life by melanocyte stem cells (MSCs), play a critical role in pigmentation and melanoma. Here, we reveal a function for the metastasis-associated phosphatase of regenerating liver 3 (PRL3) in MSC regeneration. We show that PRL3 binds to the RNA helicase DDX21, thereby restricting productive transcription by RNAPII at master transcription factor (MITF)-regulated endolysosomal vesicle genes. In zebrafish, this mechanism controls premature melanoblast expansion and differentiation from MSCs. In melanoma patients, restricted transcription of this endolysosomal vesicle pathway is a hallmark of PRL3-high melanomas. Our work presents the conceptual advance that PRL3-mediated control of transcriptional elongation is a differentiation checkpoint mechanism for activated MSCs and has clinical relevance for the activity of PRL3 in regenerating tissue and cancer.

Melanoblast transcriptome analysis reveals pathways promoting melanoma metastasis.

Cutaneous malignant melanoma is an aggressive cancer of melanocytes with a strong propensity to metastasize. We posit that melanoma cells acquire metastatic capability by adopting an embryonic-like phenotype, and that a lineage approach would uncover metastatic melanoma biology. Using a genetically engineered mouse model to generate a rich melanoblast transcriptome dataset, we identify melanoblast-specific genes whose expression contribute to metastatic competence and derive a 43-gene signature that predicts patient survival. We identify a melanoblast gene, KDELR3, whose loss impairs experimental metastasis. In contrast, KDELR1 deficiency enhances metastasis, providing the first example of different disease etiologies within the KDELR-family of retrograde transporters. We show that KDELR3 regulates the metastasis suppressor, KAI1, and report an interaction with the E3 ubiquitin-protein ligase gp78, a regulator of KAI1 degradation. Our work demonstrates that the melanoblast transcriptome can be mined to uncover targetable pathways for melanoma therapy.

Multimodel preclinical platform predicts clinical response of melanoma to immunotherapy.

Although immunotherapy has revolutionized cancer treatment, only a subset of patients demonstrate durable clinical benefit. Definitive predictive biomarkers and targets to overcome resistance remain unidentified, underscoring the urgency to develop reliable immunocompetent models for mechanistic assessment. Here we characterize a panel of syngeneic mouse models, representing a variety of molecular and phenotypic subtypes of human melanomas and exhibiting their diverse range of responses to immune checkpoint blockade (ICB). Comparative analysis of genomic, transcriptomic and tumor-infiltrating immune cell profiles demonstrated alignment with clinical observations and validated the correlation of T cell dysfunction and exclusion programs with resistance. Notably, genome-wide expression analysis uncovered a melanocytic plasticity signature predictive of patient outcome in response to ICB, suggesting that the multipotency and differentiation status of melanoma can determine ICB benefit. Our comparative preclinical platform recapitulates melanoma clinical behavior and can be employed to identify mechanisms and treatment strategies to improve patient care.

Something Old, Something New: The Tumor Microenvironment Comes of Age.

In this issue, Weeraratna and colleagues demonstrate that observed differences in melanoma aggressiveness in younger versus older patients can be explained not just by cell-intrinsic alterations over time, but by age-dependent changes in fibroblasts and the extracellular matrix they help create. Their findings identify novel cellular targets for melanoma therapy, as well as candidate prognostic biomarkers to better inform clinical decisions for patients with melanoma.See related article by Kaur et al., p. 64.See related article by Ecker et al., p. 82.

Gut microbiota dependent anti-tumor immunity restricts melanoma growth in Rnf5-/- mice.

Accumulating evidence points to an important role for the gut microbiome in anti-tumor immunity. Here, we show that altered intestinal microbiota contributes to anti-tumor immunity, limiting tumor expansion. Mice lacking the ubiquitin ligase RNF5 exhibit attenuated activation of the unfolded protein response (UPR) components, which coincides with increased expression of inflammasome components, recruitment and activation of dendritic cells and reduced expression of antimicrobial peptides in intestinal epithelial cells. Reduced UPR expression is also seen in murine and human melanoma tumor specimens that responded to immune checkpoint therapy. Co-housing of Rnf5-/- and WT mice abolishes the anti-tumor immunity and tumor inhibition phenotype, whereas transfer of 11 bacterial strains, including B. rodentium, enriched in Rnf5-/- mice, establishes anti-tumor immunity and restricts melanoma growth in germ-free WT mice. Altered UPR signaling, exemplified in Rnf5-/- mice, coincides with altered gut microbiota composition and anti-tumor immunity to control melanoma growth.