Single-cell spatial multiomics reveals tumor microenvironment vulnerabilities in cancer resistance to immunotherapy.

Heterogeneous resistance to immunotherapy remains a major challenge in cancer treatment, often leading to disease progression and death. Using CITE-seq and matched 40-plex PhenoCycler tissue imaging, we performed longitudinal multimodal single-cell analysis of tumors from metastatic melanoma patients with innate resistance, acquired resistance, or response to immunotherapy. We established the multimodal integration toolkit to align transcriptomic features, cellular epitopes, and spatial information to provide deeper insights into the tumors. With longitudinal analysis, we identified an "immune-striving" tumor microenvironment marked by peri-tumor lymphoid aggregates and low infiltration of T cells in the tumor and the emergence of MITF+SPARCL1+ and CENPF+ melanoma subclones after therapy. The enrichment of B cell-associated signatures in the molecular composition of lymphoid aggregates was associated with better survival. These findings provide further insights into the establishment of microenvironmental cell interactions and molecular composition of spatial structures that could inform therapeutic intervention.

Tumor Biomechanics Alters Metastatic Dissemination of Triple Negative Breast Cancer via Rewiring Fatty Acid Metabolism.

In recent decades, the role of tumor biomechanics on cancer cell behavior at the primary site has been increasingly appreciated. However, the effect of primary tumor biomechanics on the latter stages of the metastatic cascade, such as metastatic seeding of secondary sites and outgrowth remains underappreciated. This work sought to address this in the context of triple negative breast cancer (TNBC), a cancer type known to aggressively disseminate at all stages of disease progression. Using mechanically tuneable model systems, mimicking the range of stiffness's typically found within breast tumors, it is found that, contrary to expectations, cancer cells exposed to softer microenvironments are more able to colonize secondary tissues. It is shown that heightened cell survival is driven by enhanced metabolism of fatty acids within TNBC cells exposed to softer microenvironments. It is demonstrated that uncoupling cellular mechanosensing through integrin ß1 blocking antibody effectively causes stiff primed TNBC cells to behave like their soft counterparts, both in vitro and in vivo. This work is the first to show that softer tumor microenvironments may be contributing to changes in disease outcome by imprinting on TNBC cells a greater metabolic flexibility and conferring discrete cell survival advantages.

Detailed DNA methylation characterisation of phyllodes tumours identifies a signature of malignancy and distinguishes phyllodes from metaplastic breast carcinoma.

Phyllodes tumours (PTs) are rare fibroepithelial lesions of the breast that are classified as benign, borderline, or malignant. As little is known about the molecular underpinnings of PTs, current diagnosis relies on histological examination. However, accurate classification is often difficult, particularly for distinguishing borderline from malignant PTs. Furthermore, PTs can be misdiagnosed as other tumour types with shared histological features, such as fibroadenoma and metaplastic breast cancers. As DNA methylation is a recognised hallmark of many cancers, we hypothesised that DNA methylation could provide novel biomarkers for diagnosis and tumour stratification in PTs, whilst also allowing insight into the molecular aetiology of this otherwise understudied tumour. We generated whole-genome methylation data using the Illumina EPIC microarray in a novel PT cohort (n = 33) and curated methylation microarray data from published datasets including PTs and other potentially histopathologically similar tumours (total n = 817 samples). Analyses revealed that PTs have a unique methylome compared to normal breast tissue and to potentially histopathologically similar tumours (metaplastic breast cancer, fibroadenoma and sarcomas), with PT-specific methylation changes enriched in gene sets involved in KRAS signalling and epithelial-mesenchymal transition. Next, we identified 53 differentially methylated regions (DMRs) (false discovery rate < 0.05) that specifically delineated malignant from non-malignant PTs. The top DMR in both discovery and validation cohorts was hypermethylation at the HSD17B8 CpG island promoter. Matched PT single-cell expression data showed that HSD17B8 had minimal expression in fibroblast (putative tumour) cells. Finally, we created a methylation classifier to distinguish PTs from metaplastic breast cancer samples, where we revealed a likely misdiagnosis for two TCGA metaplastic breast cancer samples. In conclusion, DNA methylation alterations are associated with PT histopathology and hold the potential to improve our understanding of PT molecular aetiology, diagnostics, and risk stratification. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Cancer Evolution: A Multifaceted Affair.

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery. SIGNIFICANCE: Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Gene-Expression Profiling to Decipher Breast Cancer Inter- and Intratumor Heterogeneity.

Breast cancer is heterogeneous and differs substantially across different tumors (intertumor heterogeneity) and even within an individual tumor (intratumor heterogeneity). Gene-expression profiling has considerably impacted our understanding of breast cancer biology. Four main "intrinsic subtypes" of breast cancer (i.e., luminal A, luminal B, HER2-enriched, and basal-like) have been consistently identified by gene expression, showing significant prognostic and predictive value in multiple clinical scenarios. Thanks to the molecular profiling of breast tumors, breast cancer is a paradigm of treatment personalization. Several standardized prognostic gene-expression assays are presently being used in the clinic to guide treatment decisions. Moreover, the development of single-cell-level resolution molecular profiling has allowed us to appreciate that breast cancer is also heterogeneous within a single tumor. There is an evident functional heterogeneity within the neoplastic and tumor microenvironment cells. Finally, emerging insights from these studies suggest a substantial cellular organization of neoplastic and tumor microenvironment cells, thus defining breast cancer ecosystems and highlighting the importance of spatial localizations.

Mitotic Dysregulation at Tumor Initiation Creates a Therapeutic Vulnerability to Combination Anti-Mitotic and Pro-Apoptotic Agents for MYCN-Driven Neuroblastoma.

MYCN amplification occurs in approximately 20-30% of neuroblastoma patients and correlates with poor prognosis. The TH-MYCN transgenic mouse model mimics the development of human high-risk neuroblastoma and provides strong evidence for the oncogenic function of MYCN. In this study, we identified mitotic dysregulation as a hallmark of tumor initiation in the pre-cancerous ganglia from TH-MYCN mice that persists through tumor progression. Single-cell quantitative-PCR of coeliac ganglia from 10-day-old TH-MYCN mice revealed overexpression of mitotic genes in a subpopulation of premalignant neuroblasts at a level similar to single cells derived from established tumors. Prophylactic treatment using antimitotic agents barasertib and vincristine significantly delayed the onset of tumor formation, reduced pre-malignant neuroblast hyperplasia, and prolonged survival in TH-MYCN mice. Analysis of human neuroblastoma tumor cohorts showed a strong correlation between dysregulated mitosis and features of MYCN amplification, such as MYC(N) transcriptional activity, poor overall survival, and other clinical predictors of aggressive disease. To explore the therapeutic potential of targeting mitotic dysregulation, we showed that genetic and chemical inhibition of mitosis led to selective cell death in neuroblastoma cell lines with MYCN over-expression. Moreover, combination therapy with antimitotic compounds and BCL2 inhibitors exploited mitotic stress induced by antimitotics and was synergistically toxic to neuroblastoma cell lines. These results collectively suggest that mitotic dysregulation is a key component of tumorigenesis in early neuroblasts, which can be inhibited by the combination of antimitotic compounds and pro-apoptotic compounds in MYCN-driven neuroblastoma.

Luminal breast cancer identity is determined by loss of glucocorticoid receptor activity.

Glucocorticoid receptor (GR) is a transcription factor that plays a crucial role in cancer biology. In this study, we utilized an in silico-designed GR activity signature to demonstrate that GR relates to the proliferative capacity of numerous primary cancer types. In breast cancer, the GR activity status determines luminal subtype identity and has implications for patient outcomes. We reveal that GR engages with estrogen receptor (ER), leading to redistribution of ER on the chromatin. Notably, GR activation leads to upregulation of the ZBTB16 gene, encoding for a transcriptional repressor, which controls growth in ER-positive breast cancer and associates with prognosis in luminal A patients. In relation to ZBTB16's repressive nature, GR activation leads to epigenetic remodeling and loss of histone acetylation at sites proximal to cancer-driving genes. Based on these findings, epigenetic inhibitors reduce viability of ER-positive breast cancer cells that display absence of GR activity. Our findings provide insights into how GR controls ER-positive breast cancer growth and may have implications for patients' prognostication and provide novel therapeutic candidates for breast cancer treatment.

Phenotypic heterogeneity drives differential disease outcome in a mouse model of triple negative breast cancer.

The triple negative breast cancer (TNBC) subtype is one of the most aggressive forms of breast cancer that has poor clinical outcome and is an unmet clinical challenge. Accumulating evidence suggests that intratumoral heterogeneity or the presence of phenotypically distinct cell populations within a tumor play a crucial role in chemoresistance, tumor progression and metastasis. An increased understanding of the molecular regulators of intratumoral heterogeneity is crucial to the development of effective therapeutic strategies in TNBC. To this end, we used an unbiased approach to identify a molecular mediator of intratumoral heterogeneity in breast cancer by isolating two tumor cell populations (T1 and T2) from the 4T1 TNBC model. Phenotypic characterization revealed that the cells are different in terms of their morphology, proliferation and self-renewal ability in vitro as well as primary tumor formation and metastatic potential in vivo. Bioinformatic analysis followed by Kaplan Meier survival analysis in TNBC patients identified Metastasis associated colon cancer 1 (Macc1) as one of the top candidate genes mediating the aggressive phenotype in the T1 tumor cells. The role of Macc1 in regulating the proliferative phenotype was validated and taken forward in a therapeutic context with Lovastatin, a small molecule transcriptional inhibitor of Macc1 to target the T1 cell population. This study increases our understanding of the molecular underpinnings of intratumoral heterogeneity in breast cancer that is critical to improve the treatment of women currently living with the highly aggressive TNBC subtype.

Performance of tumour microenvironment deconvolution methods in breast cancer using single-cell simulated bulk mixtures.

Cells within the tumour microenvironment (TME) can impact tumour development and influence treatment response. Computational approaches have been developed to deconvolve the TME from bulk RNA-seq. Using scRNA-seq profiling from breast tumours we simulate thousands of bulk mixtures, representing tumour purities and cell lineages, to compare the performance of nine TME deconvolution methods (BayesPrism, Scaden, CIBERSORTx, MuSiC, DWLS, hspe, CPM, Bisque, and EPIC). Some methods are more robust in deconvolving mixtures with high tumour purity levels. Most methods tend to mis-predict normal epithelial for cancer epithelial as tumour purity increases, a finding that is validated in two independent datasets. The breast cancer molecular subtype influences this mis-prediction. BayesPrism and DWLS have the lowest combined numbers of false positives and false negatives, and have the best performance when deconvolving granular immune lineages. Our findings highlight the need for more single-cell characterisation of rarer cell types, and suggest that tumour cell compositions should be considered when deconvolving the TME.

Inhibition of EphA3 Expression in Tumour Stromal Cells Suppresses Tumour Growth and Progression.

Tumour progression relies on interactions with untransformed cells in the tumour microenvironment (TME), including cancer-associated fibroblasts (CAFs), which promote blood supply, tumour progression, and immune evasion. Eph receptor tyrosine kinases are cell guidance receptors that are most active during development but re-emerge in cancer and are recognised drug targets. EphA3 is overexpressed in a wide range of tumour types, and we previously found expression particularly in stromal and vascular tissues of the TME. To investigate its role in the TME, we generated transgenic mice with inducible shRNA-mediated knockdown of EphA3 expression. EphA3 knockdown was confirmed in aortic mesenchymal stem cells (MSCs), which displayed reduced angiogenic capacity. In mice with syngeneic lung tumours, EphA3 knockdown reduced vasculature and CAF/MSC-like cells in tumours, and inhibited tumour growth, which was confirmed also in a melanoma model. Single cell RNA sequencing analysis of multiple human tumour types confirmed EphA3 expression in CAFs, including in breast cancer, where EphA3 was particularly prominent in perivascular- and myofibroblast-like CAFs. Our results thus indicate expression of the cell guidance receptor EphA3 in distinct CAF subpopulations is important in supporting tumour angiogenesis and tumour growth, highlighting its potential as a therapeutic target.

ß-cell function is regulated by metabolic and epigenetic programming of islet-associated macrophages, involving Axl, Mertk, and TGFß receptor signaling.

We have exploited islet-associated macrophages (IAMs) as a model of resident macrophage function, focusing on more physiological conditions than the commonly used extremes of M1 (inflammation) versus M2 (tissue remodeling) polarization. Under steady state, murine IAMs are metabolically poised between aerobic glycolysis and oxidative phosphorylation, and thereby exert a brake on glucose-stimulated insulin secretion (GSIS). This is underpinned by epigenetic remodeling via the metabolically regulated histone demethylase Kdm5a. Conversely, GSIS is enhanced by engaging Axl receptors on IAMs, or by augmenting their oxidation of glucose. Following high-fat feeding, efferocytosis is stimulated in IAMs in conjunction with Mertk and TGFß receptor signaling. This impairs GSIS and potentially contributes to ß-cell failure in pre-diabetes. Thus, IAMs serve as relays in many more settings than currently appreciated, fine-tuning insulin secretion in response to dynamic changes in the external environment. Intervening in this nexus might represent a means of preserving ß-cell function during metabolic disease.

DNA barcoding reveals ongoing immunoediting of clonal cancer populations during metastatic progression and immunotherapy response.

Cancers evade the immune system through the process of cancer immunoediting. While immune checkpoint inhibitors are effective for reactivating tumour immunity in some cancer types, many other solid cancers, including breast cancer, remain largely non-responsive. Understanding how non-responsive cancers evade immunity and whether this occurs at the clonal level will improve immunotherapeutic design. Here we use DNA barcoding to track murine mammary cancer cell clones during immunoediting and determine clonal transcriptional profiles that allow immune evasion following anti-PD1 plus anti-CTLA4 immunotherapy. Clonal diversity is significantly restricted by immunotherapy treatment in both primary tumours and metastases, demonstrating selection for pre-existing breast cancer cell populations and ongoing immunoediting during metastasis and treatment. Immunotherapy resistant clones express a common gene signature associated with poor survival of basal-like breast cancer patient cohorts. At least one of these genes has an existing small molecule that can potentially be used to improve immunotherapy response.

STAT3 gain-of-function mutations connect leukemia with autoimmune disease by pathological NKG2Dhi CD8+ T cell dysregulation and accumulation.

The association between cancer and autoimmune disease is unexplained, exemplified by T cell large granular lymphocytic leukemia (T-LGL) where gain-of-function (GOF) somatic STAT3 mutations correlate with co-existing autoimmunity. To investigate whether these mutations are the cause or consequence of CD8+ T cell clonal expansions and autoimmunity, we analyzed patients and mice with germline STAT3 GOF mutations. STAT3 GOF mutations drove the accumulation of effector CD8+ T cell clones highly expressing NKG2D, the receptor for stress-induced MHC-class-I-related molecules. This subset also expressed genes for granzymes, perforin, interferon-γ, and Ccl5/Rantes and required NKG2D and the IL-15/IL-2 receptor IL2RB for maximal accumulation. Leukocyte-restricted STAT3 GOF was sufficient and CD8+ T cells were essential for lethal pathology in mice. These results demonstrate that STAT3 GOF mutations cause effector CD8+ T cell oligoclonal accumulation and that these rogue cells contribute to autoimmune pathology, supporting the hypothesis that somatic mutations in leukemia/lymphoma driver genes contribute to autoimmune disease.

Spatial transcriptomics.

Spatial transcriptomics, with other spatial technologies, has enabled scientists to dissect the organization and interaction of different cell types within the tumor microenvironment. We asked experts to discuss some aspects of this technology from revealing the tumor microenvironment and heterogeneity, to tracking tumor evolution, to guiding tumor therapy, to current technical challenges.

Temporal profiling of the breast tumour microenvironment reveals collagen XII as a driver of metastasis.

The tumour stroma, and in particular the extracellular matrix (ECM), is a salient feature of solid tumours that plays a crucial role in shaping their progression. Many desmoplastic tumours including breast cancer involve the significant accumulation of type I collagen. However, recently it has become clear that the precise distribution and organisation of matrix molecules such as collagen I is equally as important in the tumour as their abundance. Cancer-associated fibroblasts (CAFs) coexist within breast cancer tissues and play both pro- and anti-tumourigenic roles through remodelling the ECM. Here, using temporal proteomic profiling of decellularized tumours, we interrogate the evolving matrisome during breast cancer progression. We identify 4 key matrisomal clusters, and pinpoint collagen type XII as a critical component that regulates collagen type I organisation. Through combining our proteomics with single-cell transcriptomics, and genetic manipulation models, we show how CAF-secreted collagen XII alters collagen I organisation to create a pro-invasive microenvironment supporting metastatic dissemination. Finally, we show in patient cohorts that collagen XII may represent an indicator of breast cancer patients at high risk of metastatic relapse.

Unchecked oxidative stress in skeletal muscle prevents outgrowth of disseminated tumour cells.

Skeletal muscle has long been recognized as an inhospitable site for disseminated tumour cells (DTCs). Yet its antimetastatic nature has eluded a thorough mechanistic examination. Here, we show that DTCs traffic to and persist within skeletal muscle in mice and in humans, which raises the question of how this tissue suppresses colonization. Results from mouse and organotypic culture models along with metabolomic profiling suggested that skeletal muscle imposes a sustained oxidative stress on DTCs that impairs their proliferation. Functional studies demonstrated that disrupting reduction-oxidation homeostasis via chemogenetic induction of reactive oxygen species slowed proliferation in a more fertile organ: the lung. Conversely, enhancement of the antioxidant potential of tumour cells through ectopic expression of catalase in the tumour or host mitochondria allowed robust colonization of skeletal muscle. These findings reveal a profound metabolic bottleneck imposed on DTCs and sustained by skeletal muscle. A thorough understanding of this biology could reveal previously undocumented DTC vulnerabilities that can be exploited to prevent metastasis in other more susceptible tissues.