Mutant IDH inhibitors induce lineage differentiation in IDH-mutant oligodendroglioma.

A subset of patients with IDH-mutant glioma respond to inhibitors of mutant IDH (IDHi), yet the molecular underpinnings of such responses are not understood. Here, we profiled by single-cell or single-nucleus RNA-sequencing three IDH-mutant oligodendrogliomas from patients who derived clinical benefit from IDHi. Importantly, the tissues were sampled on-drug, four weeks from treatment initiation. We further integrate our findings with analysis of single-cell and bulk transcriptomes from independent cohorts and experimental models. We find that IDHi treatment induces a robust differentiation toward the astrocytic lineage, accompanied by a depletion of stem-like cells and a reduction of cell proliferation. Furthermore, mutations in NOTCH1 are associated with decreased astrocytic differentiation and may limit the response to IDHi. Our study highlights the differentiating potential of IDHi on the cellular hierarchies that drive oligodendrogliomas and suggests a genetic modifier that may improve patient stratification.

Integrative spatial analysis reveals a multi-layered organization of glioblastoma.

Glioma contains malignant cells in diverse states. Here, we combine spatial transcriptomics, spatial proteomics, and computational approaches to define glioma cellular states and uncover their organization. We find three prominent modes of organization. First, gliomas are composed of small local environments, each typically enriched with one major cellular state. Second, specific pairs of states preferentially reside in proximity across multiple scales. This pairing of states is consistent across tumors. Third, these pairwise interactions collectively define a global architecture composed of five layers. Hypoxia appears to drive the layers, as it is associated with a long-range organization that includes all cancer cell states. Accordingly, tumor regions distant from any hypoxic/necrotic foci and tumors that lack hypoxia such as low-grade IDH-mutant glioma are less organized. In summary, we provide a conceptual framework for the organization of cellular states in glioma, highlighting hypoxia as a long-range tissue organizer.

Glioma synapses recruit mechanisms of adaptive plasticity.

The role of the nervous system in the regulation of cancer is increasingly appreciated. In gliomas, neuronal activity drives tumour progression through paracrine signalling factors such as neuroligin-3 and brain-derived neurotrophic factor1-3 (BDNF), and also through electrophysiologically functional neuron-to-glioma synapses mediated by AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptors4,5. The consequent glioma cell membrane depolarization drives tumour proliferation4,6. In the healthy brain, activity-regulated secretion of BDNF promotes adaptive plasticity of synaptic connectivity7,8 and strength9-15. Here we show that malignant synapses exhibit similar plasticity regulated by BDNF. Signalling through the receptor tropomyosin-related kinase B16 (TrkB) to CAMKII, BDNF promotes AMPA receptor trafficking to the glioma cell membrane, resulting in increased amplitude of glutamate-evoked currents in the malignant cells. Linking plasticity of glioma synaptic strength to tumour growth, graded optogenetic control of glioma membrane potential demonstrates that greater depolarizing current amplitude promotes increased glioma proliferation. This potentiation of malignant synaptic strength shares mechanistic features with synaptic plasticity17-22 that contributes to memory and learning in the healthy brain23-26. BDNF-TrkB signalling also regulates the number of neuron-to-glioma synapses. Abrogation of activity-regulated BDNF secretion from the brain microenvironment or loss of glioma TrkB expression robustly inhibits tumour progression. Blocking TrkB genetically or pharmacologically abrogates these effects of BDNF on glioma synapses and substantially prolongs survival in xenograft models of paediatric glioblastoma and diffuse intrinsic pontine glioma. Together, these findings indicate that BDNF-TrkB signalling promotes malignant synaptic plasticity and augments tumour progression.

Spectrum of Response to Platinum and PARP Inhibitors in Germline BRCA-Associated Pancreatic Cancer in the Clinical and Preclinical Setting.

Germline BRCA-associated pancreatic ductal adenocarcinoma (glBRCA PDAC) tumors are susceptible to platinum and PARP inhibition. The clinical outcomes of 125 patients with glBRCA PDAC were stratified based on the spectrum of response to platinum/PARP inhibition: (i) refractory [overall survival (OS) <6 months], (ii) durable response followed by acquired resistance (OS <36 months), and (iii) long-term responders (OS >36 months). Patient-derived xenografts (PDX) were generated from 25 patients with glBRCA PDAC at different clinical time points. Response to platinum/PARP inhibition in vivo and ex vivo culture (EVOC) correlated with clinical response. We deciphered the mechanisms of resistance in glBRCA PDAC and identified homologous recombination (HR) proficiency and secondary mutations restoring partial functionality as the most dominant resistant mechanism. Yet, a subset of HR-deficient (HRD) patients demonstrated clinical resistance. Their tumors displayed basal-like molecular subtype and were more aneuploid. Tumor mutational burden was high in HRD PDAC and significantly higher in tumors with secondary mutations. Anti-PD-1 attenuated tumor growth in a novel humanized glBRCA PDAC PDX model. This work demonstrates the utility of preclinical models, including EVOC, to predict the response of glBRCA PDAC to treatment, which has the potential to inform time-sensitive medical decisions. SIGNIFICANCE: glBRCA PDAC has a favorable response to platinum/PARP inhibition. However, most patients develop resistance. Additional treatment options for this unique subpopulation are needed. We generated model systems in PDXs and an ex vivo system (EVOC) that faithfully recapitulate these specific clinical scenarios as a platform to investigate the mechanisms of resistance for further drug development. This article is highlighted in the In This Issue feature, p. 1749.

Hallmarks of transcriptional intratumour heterogeneity across a thousand tumours.

Each tumour contains diverse cellular states that underlie intratumour heterogeneity (ITH), a central challenge of cancer therapeutics1. Dozens of recent studies have begun to describe ITH by single-cell RNA sequencing, but each study typically profiled only a small number of tumours and provided a narrow view of transcriptional ITH2. Here we curate, annotate and integrate the data from 77 different studies to reveal the patterns of transcriptional ITH across 1,163 tumour samples covering 24 tumour types. Among the malignant cells, we identify 41 consensus meta-programs, each consisting of dozens of genes that are coordinately upregulated in subpopulations of cells within many tumours. The meta-programs cover diverse cellular processes including both generic (for example, cell cycle and stress) and lineage-specific patterns that we map into 11 hallmarks of transcriptional ITH. Most meta-programs of carcinoma cells are similar to those identified in non-malignant epithelial cells, suggesting that a large fraction of malignant ITH programs are variable even before oncogenesis, reflecting the biology of their cell of origin. We further extended the meta-program analysis to six common non-malignant cell types and utilize these to map cell-cell interactions within the tumour microenvironment. In summary, we have assembled a comprehensive pan-cancer single-cell RNA-sequencing dataset, which is available through the Curated Cancer Cell Atlas website, and leveraged this dataset to carry out a systematic characterization of transcriptional ITH.

Cellular states are coupled to genomic and viral heterogeneity in HPV-related oropharyngeal carcinoma.

Head and neck squamous cell carcinoma (HNSCC) includes a subset of cancers driven by human papillomavirus (HPV). Here we use single-cell RNA-seq to profile both HPV-positive and HPV-negative oropharyngeal tumors, uncovering a high level of cellular diversity within and between tumors. First, we detect diverse chromosomal aberrations within individual tumors, suggesting genomic instability and enabling the identification of malignant cells even at pathologically negative margins. Second, we uncover diversity with respect to HNSCC subtypes and other cellular states such as the cell cycle, senescence and epithelial-mesenchymal transitions. Third, we find heterogeneity in viral gene expression within HPV-positive tumors. HPV expression is lost or repressed in a subset of cells, which are associated with a decrease in HPV-associated cell cycle phenotypes, decreased response to treatment, increased invasion and poor prognosis. These findings suggest that HPV expression diversity must be considered during diagnosis and treatment of HPV-positive tumors, with important prognostic ramifications.

Elucidating the diversity of malignant mesenchymal states in glioblastoma by integrative analysis.

BACKGROUND: Multiple glioblastoma studies have described a mesenchymal (MES) state, with each study defining the MES program by distinct sets of genes and highlighting distinct functional associations, including both immune activation and suppression. These variable descriptions complicate our understanding of the MES state and its implications. Here, we hypothesize that there is a range of glioma MES states, possibly reflecting distinct prior states in which a MES program can be induced, and/or distinct mechanisms that induce the MES states in those cells. METHODS: We integrated multiple published single-cell and bulk RNA sequencing datasets and MES signatures to define a core MES program that recurs across studies, as well as multiple function-specific MES signatures that vary across MES cells. We then examined the co-occurrence of these signatures and their associations with genetic and microenvironmental features. RESULTS: Based on co-occurrence of MES signatures, we found three main variants of MES states: hypoxia-related (MES-Hyp), astrocyte-related (MES-Ast), and an intermediate state. Notably, the MES states are differentially associated with genetic and microenvironmental features. MES-Hyp is preferentially associated with NF1 deletion, overall macrophage abundance, a high macrophage/microglia ratio, and M2-related macrophages, consistent with previous studies that associated MES with immune suppression. In contrast, MES-Ast is associated with T cell abundance and cytotoxicity, consistent with immune activation through expression of MHC-I/II. CONCLUSIONS: Diverse MES states occur in glioblastoma. These states share a subset of core genes but differ primarily in their association with hypoxia vs. astrocytic expression programs, and with immune suppression vs. activation, respectively.

An activation to memory differentiation trajectory of tumor-infiltrating lymphocytes informs metastatic melanoma outcomes.

There is a need for better classification and understanding of tumor-infiltrating lymphocytes (TILs). Here, we applied advanced functional genomics to interrogate 9,000 human tumors and multiple single-cell sequencing sets using benchmarked T cell states, comprehensive T cell differentiation trajectories, human and mouse vaccine responses, and other human TILs. Compared with other T cell states, enrichment of T memory/resident memory programs was observed across solid tumors. Trajectory analysis of single-cell melanoma CD8+ TILs also identified a high fraction of memory/resident memory-scoring TILs in anti-PD-1 responders, which expanded post therapy. In contrast, TILs scoring highly for early T cell activation, but not exhaustion, associated with non-response. Late/persistent, but not early activation signatures, prognosticate melanoma survival, and co-express with dendritic cell and IFN-γ response programs. These data identify an activation-like state associated to poor response and suggest successful memory conversion, above resuscitation of exhaustion, is an under-appreciated aspect of successful anti-tumoral immunity.

Genome-wide CRISPR screen identifies PRC2 and KMT2D-COMPASS as regulators of distinct EMT trajectories that contribute differentially to metastasis.

Epithelial-mesenchymal transition (EMT) programs operate within carcinoma cells, where they generate phenotypes associated with malignant progression. In their various manifestations, EMT programs enable epithelial cells to enter into a series of intermediate states arrayed along the E-M phenotypic spectrum. At present, we lack a coherent understanding of how carcinoma cells control their entrance into and continued residence in these various states, and which of these states favour the process of metastasis. Here we characterize a layer of EMT-regulating machinery that governs E-M plasticity (EMP). This machinery consists of two chromatin-modifying complexes, PRC2 and KMT2D-COMPASS, which operate as critical regulators to maintain a stable epithelial state. Interestingly, loss of these two complexes unlocks two distinct EMT trajectories. Dysfunction of PRC2, but not KMT2D-COMPASS, yields a quasi-mesenchymal state that is associated with highly metastatic capabilities and poor survival of patients with breast cancer, suggesting that great caution should be applied when PRC2 inhibitors are evaluated clinically in certain patient cohorts. These observations identify epigenetic factors that regulate EMP, determine specific intermediate EMT states and, as a direct consequence, govern the metastatic ability of carcinoma cells.

Interactions between cancer cells and immune cells drive transitions to mesenchymal-like states in glioblastoma.

The mesenchymal subtype of glioblastoma is thought to be determined by both cancer cell-intrinsic alterations and extrinsic cellular interactions, but remains poorly understood. Here, we dissect glioblastoma-to-microenvironment interactions by single-cell RNA sequencing analysis of human tumors and model systems, combined with functional experiments. We demonstrate that macrophages induce a transition of glioblastoma cells into mesenchymal-like (MES-like) states. This effect is mediated, both in vitro and in vivo, by macrophage-derived oncostatin M (OSM) that interacts with its receptors (OSMR or LIFR) in complex with GP130 on glioblastoma cells and activates STAT3. We show that MES-like glioblastoma states are also associated with increased expression of a mesenchymal program in macrophages and with increased cytotoxicity of T cells, highlighting extensive alterations of the immune microenvironment with potential therapeutic implications.

Decoupling epithelial-mesenchymal transitions from stromal profiles by integrative expression analysis.

Epithelial-to-mesenchymal transition (EMT) is the most commonly cited mechanism for cancer metastasis, but it is difficult to distinguish from profiles of normal stromal cells in the tumour microenvironment. In this study we use published single cell RNA-seq data to directly compare mesenchymal signatures from cancer and stromal cells. Informed by these comparisons, we developed a computational framework to decouple these two sources of mesenchymal expression profiles using bulk RNA-seq datasets. This deconvolution offers the opportunity to characterise EMT across hundreds of tumours and examine its association with metastasis and other clinical features. With this approach, we find three distinct patterns of EMT, associated with squamous, gynaecological and gastrointestinal cancer types. Surprisingly, in most cancer types, EMT patterns are not associated with increased chance of metastasis, suggesting that other steps in the metastatic cascade may represent the main bottleneck. This work provides a comprehensive evaluation of EMT profiles and their functional significance across hundreds of tumours while circumventing the confounding effect of stromal cells.

Single-Cell Deconvolution of Head and Neck Squamous Cell Carcinoma.

Complexities in cell-type composition have rightfully led to skepticism and caution in the interpretation of bulk transcriptomic analyses. Recent studies have shown that deconvolution algorithms can be utilized to computationally estimate cell-type proportions from the gene expression data of bulk blood samples, but their performance when applied to tumor tissues, including those from head and neck, remains poorly characterized. Here, we use single-cell data (~6000 single cells) collected from 21 head and neck squamous cell carcinoma (HNSCC) samples to generate cell-type-specific gene expression signatures. We leverage bulk RNA-seq data from >500 HNSCC samples profiled by The Cancer Genome Atlas (TCGA), and using single-cell data as a reference, apply two newly developed deconvolution algorithms (CIBERSORTx and MuSiC) to the bulk transcriptome data to quantitatively estimate cell-type proportions for each tumor in TCGA. We show that these two algorithms produce similar estimates of constituent/major cell-type proportions and that a high T-cell fraction correlates with improved survival. By further characterizing T-cell subpopulations, we identify that regulatory T-cells (Tregs) were the major contributor to this improved survival. Lastly, we assessed gene expression, specifically in the Treg population, and found that TNFRSF4 (Tumor Necrosis Factor Receptor Superfamily Member 4) was differentially expressed in the core Treg subpopulation. Moreover, higher TNFRSF4 expression was associated with greater survival, suggesting that TNFRSF4 could play a key role in mechanisms underlying the contribution of Treg in HNSCC outcomes.

Decoding Cancer Biology One Cell at a Time.

Human tumors are composed of diverse malignant and nonmalignant cells, generating a complex ecosystem that governs tumor biology and response to treatments. Recent technological advances have enabled the characterization of tumors at single-cell resolution, providing a compelling strategy to dissect their intricate biology. Here we describe recent developments in single-cell expression profiling and the studies applying them in clinical settings. We highlight some of the powerful insights gleaned from these studies for tumor classification, stem cell programs, tumor microenvironment, metastasis, and response to targeted and immune therapies. SIGNIFICANCE: Intratumor heterogeneity (ITH) has been a major barrier to our understanding of cancer. Single-cell genomics is leading a revolution in our ability to systematically dissect ITH. In this review, we focus on single-cell expression profiling and lessons learned in key aspects of human tumor biology.

Glioblastoma cell populations with distinct oncogenic programs release podoplanin as procoagulant extracellular vesicles.

Vascular anomalies, including local and peripheral thrombosis, are a hallmark of glioblastoma (GBM) and an aftermath of deregulation of the cancer cell genome and epigenome. Although the molecular effectors of these changes are poorly understood, the upregulation of podoplanin (PDPN) by cancer cells has recently been linked to an increased risk for venous thromboembolism (VTE) in GBM patients. Therefore, regulation of this platelet-activating protein by transforming events in cancer cells is of considerable interest. We used single-cell and bulk transcriptome data mining, as well as cellular and xenograft models in mice, to analyze the nature of cells expressing PDPN, as well as their impact on the activation of the coagulation system and platelets. We report that PDPN is expressed by distinct (mesenchymal) GBM cell subpopulations and downregulated by oncogenic mutations of EGFR and IDH1 genes, along with changes in chromatin modifications (enhancer of zeste homolog 2) and DNA methylation. Glioma cells exteriorize their PDPN and/or tissue factor (TF) as cargo of exosome-like extracellular vesicles (EVs) shed from cells in vitro and in vivo. Injection of glioma-derived podoplanin carrying extracelluar vesicles (PDPN-EVs) activates platelets, whereas tissue factor carrying extracellular vesicles (TF-EVs) activate the clotting cascade. Similarly, an increase in platelet activation (platelet factor 4) or coagulation (D-dimer) markers occurs in mice harboring the corresponding glioma xenografts expressing PDPN or TF, respectively. Coexpression of PDPN and TF by GBM cells cooperatively affects tumor microthrombosis. Thus, in GBM, distinct cellular subsets drive multiple facets of cancer-associated thrombosis and may represent targets for phenotype- and cell type-based diagnosis and antithrombotic intervention.

A unique subset of glycolytic tumour-propagating cells drives squamous cell carcinoma.

Head and neck squamous cell carcinoma (SCC) remains among the most aggressive human cancers. Tumour progression and aggressiveness in SCC are largely driven by tumour-propagating cells (TPCs). Aerobic glycolysis, also known as the Warburg effect, is a characteristic of many cancers; however, whether this adaptation is functionally important in SCC, and at which stage, remains poorly understood. Here, we show that the NAD+-dependent histone deacetylase sirtuin 6 is a robust tumour suppressor in SCC, acting as a modulator of glycolysis in these tumours. Remarkably, rather than a late adaptation, we find enhanced glycolysis specifically in TPCs. More importantly, using single-cell RNA sequencing of TPCs, we identify a subset of TPCs with higher glycolysis and enhanced pentose phosphate pathway and glutathione metabolism, characteristics that are strongly associated with a better antioxidant response. Together, our studies uncover enhanced glycolysis as a main driver in SCC, and, more importantly, identify a subset of TPCs as the cell of origin for the Warburg effect, defining metabolism as a key feature of intra-tumour heterogeneity.

Inhibitory CD161 receptor identified in glioma-infiltrating T cells by single-cell analysis.

T cells are critical effectors of cancer immunotherapies, but little is known about their gene expression programs in diffuse gliomas. Here, we leverage single-cell RNA sequencing (RNA-seq) to chart the gene expression and clonal landscape of tumor-infiltrating T cells across 31 patients with isocitrate dehydrogenase (IDH) wild-type glioblastoma and IDH mutant glioma. We identify potential effectors of anti-tumor immunity in subsets of T cells that co-express cytotoxic programs and several natural killer (NK) cell genes. Analysis of clonally expanded tumor-infiltrating T cells further identifies the NK gene KLRB1 (encoding CD161) as a candidate inhibitory receptor. Accordingly, genetic inactivation of KLRB1 or antibody-mediated CD161 blockade enhances T cell-mediated killing of glioma cells in vitro and their anti-tumor function in vivo. KLRB1 and its associated transcriptional program are also expressed by substantial T cell populations in other human cancers. Our work provides an atlas of T cells in gliomas and highlights CD161 and other NK cell receptors as immunotherapy targets.

Pan-cancer single-cell RNA-seq identifies recurring programs of cellular heterogeneity.

Cultured cell lines are the workhorse of cancer research, but the extent to which they recapitulate the heterogeneity observed among malignant cells in tumors is unclear. Here we used multiplexed single-cell RNA-seq to profile 198 cancer cell lines from 22 cancer types. We identified 12 expression programs that are recurrently heterogeneous within multiple cancer cell lines. These programs are associated with diverse biological processes, including cell cycle, senescence, stress and interferon responses, epithelial-mesenchymal transition and protein metabolism. Most of these programs recapitulate those recently identified as heterogeneous within human tumors. We prioritized specific cell lines as models of cellular heterogeneity and used them to study subpopulations of senescence-related cells, demonstrating their dynamics, regulation and unique drug sensitivities, which were predictive of clinical response. Our work describes the landscape of heterogeneity within diverse cancer cell lines and identifies recurrent patterns of heterogeneity that are shared between tumors and specific cell lines.

Malignant cell-specific CXCL14 promotes tumor lymphocyte infiltration in oral cavity squamous cell carcinoma.

OBJECTIVES: To explore lymphocyte infiltration as a potential mechanism behind CXCL14-mediated tumor growth suppression in oral cavity squamous cell carcinoma (OSCC). METHODS: We analyzed single cell RNA-sequencing (scRNA-seq) data from OSCC to identify expression changes among malignant cells in lymph nodes (LN) versus primary tumors. CXCL14 expression in murine OSCC cell lines was quantified using qRT-PCR. Short hairpin RNA knockdown of CXCL14 was performed in mouse oral cavity (MOC)1 cells, and CXCL14 overexpression was performed in MOC2 cells. Cells in each condition were injected into C57BL/6 mice with and without T cell depletion, and tumor volume was measured. At 30 days, tumors were dissociated and analyzed by flow cytometry for CD45+CD3+ T cells. CXCL14 expression was correlated with gene expression signatures of tumor infiltrating lymphocytes (TIL) in scRNA-seq data, as well as TCGA tumors. RESULTS: scRNA-seq revealed CXCL14 as the most significantly downregulated gene among malignant cells in LNs relative to primary tumor, supporting a role in preventing invasion and/or metastasis. In a murine immunocompetent model, CXCL14 expression was higher in indolent MOC1 cells than in more aggressive MOC2 cells. Tumor growth in vivo was significantly increased by CXCL14 knockdown in MOC1 cells relative to control, with a corresponding decrease in TIL. In MOC2 cells, tumor growth was significantly reduced by CXCL14 overexpression relative to control and TIL were increased. Both effects were lost with T cell depletion. In a human tumor scRNA-seq cohort, we found that only malignant cell CXCL14, but not non-malignant cell or fibroblast CXCL14, was associated with TIL. Bulk CXCL14 from the TCGA cohort had no association with TIL. CONCLUSIONS: Higher CXCL14 expression by tumor cells is associated with reduced tumor growth and increased TIL, supporting immune-mediated suppression of tumor growth in OSCC. Given that CXCL14 is downregulated in LN metastases compared with primary tumors, our data raise the possibility that CXCL14-mediated immune infiltration may discourage invasion and metastasis. In human scRNA-seq data, only malignant cell-specific CXCL14 was associated with TIL, suggesting a critical context-dependent effect of CXCL14 expression.