Discovering dominant tumor immune archetypes in a pan-cancer census.

Cancers display significant heterogeneity with respect to tissue of origin, driver mutations, and other features of the surrounding tissue. It is likely that individual tumors engage common patterns of the immune system-here "archetypes"-creating prototypical non-destructive tumor immune microenvironments (TMEs) and modulating tumor-targeting. To discover the dominant immune system archetypes, the University of California, San Francisco (UCSF) Immunoprofiler Initiative (IPI) processed 364 individual tumors across 12 cancer types using standardized protocols. Computational clustering of flow cytometry and transcriptomic data obtained from cell sub-compartments uncovered dominant patterns of immune composition across cancers. These archetypes were profound insofar as they also differentiated tumors based upon unique immune and tumor gene-expression patterns. They also partitioned well-established classifications of tumor biology. The IPI resource provides a template for understanding cancer immunity as a collection of dominant patterns of immune organization and provides a rational path forward to learn how to modulate these to improve therapy.

HiTSelect: a comprehensive tool for high-complexity-pooled screen analysis.

Genetic screens of an unprecedented scale have recently been made possible by the availability of high-complexity libraries of synthetic oligonucleotides designed to mediate either gene knockdown or gene knockout, coupled with next-generation sequencing. However, several sources of random noise and statistical biases complicate the interpretation of the resulting high-throughput data. We developed HiTSelect, a comprehensive analysis pipeline for rigorously selecting screen hits and identifying functionally relevant genes and pathways by addressing off-target effects, controlling for variance in both gene silencing efficiency and sequencing depth of coverage and integrating relevant metadata. We document the superior performance of HiTSelect using data from both genome-wide RNAi and CRISPR/Cas9 screens. HiTSelect is implemented as an open-source package, with a user-friendly interface for data visualization and pathway exploration. Binary executables are available at http://sourceforge.net/projects/hitselect/, and the source code is available at https://github.com/diazlab/HiTSelect.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma.

Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter- and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of antigens. In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface antigens that could be targeted by antibodies and chimeric antigen receptor-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas) and 9166 normal tissue samples (from the Genotype-Tissue Expression project), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN, which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative antigens. Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

VC-resist glioblastoma cell state: vessel co-option as a key driver of chemoradiation resistance.

Glioblastoma (GBM) is a highly lethal type of cancer. GBM recurrence following chemoradiation is typically attributed to the regrowth of invasive and resistant cells. Therefore, there is a pressing need to gain a deeper understanding of the mechanisms underlying GBM resistance to chemoradiation and its ability to infiltrate. Using a combination of transcriptomic, proteomic, and phosphoproteomic analyses, longitudinal imaging, organotypic cultures, functional assays, animal studies, and clinical data analyses, we demonstrate that chemoradiation and brain vasculature induce cell transition to a functional state named VC-Resist (vessel co-opting and resistant cell state). This cell state is midway along the transcriptomic axis between proneural and mesenchymal GBM cells and is closer to the AC/MES1-like state. VC-Resist GBM cells are highly vessel co-opting, allowing significant infiltration into the surrounding brain tissue and homing to the perivascular niche, which in turn induces even more VC-Resist transition. The molecular and functional characteristics of this FGFR1-YAP1-dependent GBM cell state, including resistance to DNA damage, enrichment in the G2M phase, and induction of senescence/stemness pathways, contribute to its enhanced resistance to chemoradiation. These findings demonstrate how vessel co-option, perivascular niche, and GBM cell plasticity jointly drive resistance to therapy during GBM recurrence.

Challenges in the discovery of tumor-specific alternative splicing-derived cell-surface antigens in glioma.

Background: Despite advancements in cancer immunotherapy, solid tumors remain formidable challenges. In glioma, profound inter-and intra-tumoral heterogeneity of antigen landscape hampers therapeutic development. Therefore, it is critical to consider alternative sources to expand the repertoire of targetable (neo-)antigens and improve therapeutic outcomes. Accumulating evidence suggests that tumor-specific alternative splicing (AS) could be an untapped reservoir of neoantigens. Results: In this study, we investigated tumor-specific AS events in glioma, focusing on those predicted to generate major histocompatibility complex (MHC)-presentation-independent, cell-surface neoantigens that could be targeted by antibodies and chimeric antigen receptor (CAR)-T cells. We systematically analyzed bulk RNA-sequencing datasets comparing 429 tumor samples (from The Cancer Genome Atlas [TCGA]) and 9,166 normal tissue samples (from the Genotype-Tissue Expression project [GTEx]), and identified 13 AS events in 7 genes predicted to be expressed in more than 10% of the patients, including PTPRZ1 and BCAN , which were corroborated by an external RNA-sequencing dataset. Subsequently, we validated our predictions and elucidated the complexity of the isoforms using full-length transcript amplicon sequencing on patient-derived glioblastoma cells. However, analyses of the RNA-sequencing datasets of spatially mapped and longitudinally collected clinical tumor samples unveiled remarkable spatiotemporal heterogeneity of the candidate AS events. Furthermore, proteomics analysis did not reveal any peptide spectra matching the putative neoantigens. Conclusions: Our investigation illustrated the diverse characteristics of the tumor-specific AS events and the challenges of antigen exploration due to their notable spatiotemporal heterogeneity and elusive nature at the protein levels. Redirecting future efforts toward intracellular, MHC-presented antigens could offer a more viable avenue.

A single-cell atlas of glioblastoma evolution under therapy reveals cell-intrinsic and cell-extrinsic therapeutic targets.

Recent longitudinal studies of glioblastoma (GBM) have demonstrated a lack of apparent selection pressure for specific DNA mutations in recurrent disease. Single-cell lineage tracing has shown that GBM cells possess a high degree of plasticity. Together this suggests that phenotype switching, as opposed to genetic evolution, may be the escape mechanism that explains the failure of precision therapies to date. We profiled 86 primary-recurrent patient-matched paired GBM specimens with single-nucleus RNA, single-cell open-chromatin, DNA and spatial transcriptomic/proteomic assays. We found that recurrent GBMs are characterized by a shift to a mesenchymal phenotype. We show that the mesenchymal state is mediated by activator protein 1. Increased T-cell abundance at recurrence was prognostic and correlated with hypermutation status. We identified tumor-supportive networks of paracrine and autocrine signals between GBM cells, nonmalignant neuroglia and immune cells. We present cell-intrinsic and cell-extrinsic targets and a single-cell multiomics atlas of GBM under therapy.

The evolution of alternative splicing in glioblastoma under therapy.

BACKGROUND: Alternative splicing is a rich source of tumor-specific neoantigen targets for immunotherapy. This holds promise for glioblastomas (GBMs), the most common primary tumors of the adult brain, which are resistant to standard-of-care therapy. Although most clinical trials enroll patients at recurrence, most preclinical studies have been done with specimens from primary disease. There are limited expression data from GBMs at recurrence and surprisingly little is known about the evolution of splicing patterns under therapy. RESULT: We profile 37 primary-recurrent paired human GBM specimens via RNA sequencing. We describe the landscape of alternative splicing in GBM at recurrence and contrast that to primary and non-malignant brain-tissue specimens. By screening single-cell atlases, we identify cell-type-specific splicing patterns and novel splicing events in cell-surface proteins that are suitable targets for engineered T cell therapies. We identify recurrent-specific isoforms of mitogen-activated kinase pathway genes that enhance invasiveness and are preferentially expressed by stem-like cells. CONCLUSION: These studies shed light on gene expression in recurrent GBM and identify novel targets for therapeutic development.

ATRX regulates glial identity and the tumor microenvironment in IDH-mutant glioma.

BACKGROUND: Recent single-cell transcriptomic studies report that IDH-mutant gliomas share a common hierarchy of cellular phenotypes, independent of genetic subtype. However, the genetic differences between IDH-mutant glioma subtypes are prognostic, predictive of response to chemotherapy, and correlate with distinct tumor microenvironments. RESULTS: To reconcile these findings, we profile 22 human IDH-mutant gliomas using scATAC-seq and scRNA-seq. We determine the cell-type-specific differences in transcription factor expression and associated regulatory grammars between IDH-mutant glioma subtypes. We find that while IDH-mutant gliomas do share a common distribution of cell types, there are significant differences in the expression and targeting of transcription factors that regulate glial identity and cytokine elaboration. We knock out the chromatin remodeler ATRX, which suffers loss-of-function alterations in most IDH-mutant astrocytomas, in an IDH-mutant immunocompetent intracranial murine model. We find that both human ATRX-mutant gliomas and murine ATRX-knockout gliomas are more heavily infiltrated by immunosuppressive monocytic-lineage cells derived from circulation than ATRX-intact gliomas, in an IDH-mutant background. ATRX knockout in murine glioma recapitulates gene expression and open chromatin signatures that are specific to human ATRX-mutant astrocytomas, including drivers of astrocytic lineage and immune-cell chemotaxis. Through single-cell cleavage under targets and tagmentation assays and meta-analysis of public data, we show that ATRX loss leads to a global depletion in CCCTC-binding factor association with DNA, gene dysregulation along associated chromatin loops, and protection from therapy-induced senescence. CONCLUSIONS: These studies explain how IDH-mutant gliomas from different subtypes maintain distinct phenotypes and tumor microenvironments despite a common lineage hierarchy.

TIGIT and PD-1 Immune Checkpoint Pathways Are Associated With Patient Outcome and Anti-Tumor Immunity in Glioblastoma.

Glioblastoma (GBM) remains an aggressive brain tumor with a high rate of mortality. Immune checkpoint (IC) molecules are expressed on tumor infiltrating lymphocytes (TILs) and promote T cell exhaustion upon binding to IC ligands expressed by the tumor cells. Interfering with IC pathways with immunotherapy has promoted reactivation of anti-tumor immunity and led to success in several malignancies. However, IC inhibitors have achieved limited success in GBM patients, suggesting that other checkpoint molecules may be involved with suppressing TIL responses. Numerous IC pathways have been described, with current testing of inhibitors underway in multiple clinical trials. Identification of the most promising checkpoint pathways may be useful to guide the future trials for GBM. Here, we analyzed the The Cancer Genome Atlas (TCGA) transcriptomic database and identified PD1 and TIGIT as top putative targets for GBM immunotherapy. Additionally, dual blockade of PD1 and TIGIT improved survival and augmented CD8+ TIL accumulation and functions in a murine GBM model compared with either single agent alone. Furthermore, we demonstrated that this combination immunotherapy affected granulocytic/polymorphonuclear (PMN) myeloid derived suppressor cells (MDSCs) but not monocytic (Mo) MDSCs in in our murine gliomas. Importantly, we showed that suppressive myeloid cells express PD1, PD-L1, and TIGIT-ligands in human GBM tissue, and demonstrated that antigen specific T cell proliferation that is inhibited by immunosuppressive myeloid cells can be restored by TIGIT/PD1 blockade. Our data provide new insights into mechanisms of GBM αPD1/αTIGIT immunotherapy.

Clonal ZEB1-Driven Mesenchymal Transition Promotes Targetable Oncologic Antiangiogenic Therapy Resistance.

Glioblastoma (GBM) responses to bevacizumab are invariably transient with acquired resistance. We profiled paired patient specimens and bevacizumab-resistant xenograft models pre- and post-resistance toward the primary goal of identifying regulators whose targeting could prolong the therapeutic window, and the secondary goal of identifying biomarkers of therapeutic window closure. Bevacizumab-resistant patient specimens and xenografts exhibited decreased vessel density and increased hypoxia versus pre-resistance, suggesting that resistance occurs despite effective therapeutic devascularization. Microarray analysis revealed upregulated mesenchymal genes in resistant tumors correlating with bevacizumab treatment duration and causing three changes enabling resistant tumor growth in hypoxia. First, perivascular invasiveness along remaining blood vessels, which co-opts vessels in a VEGF-independent and neoangiogenesis-independent manner, was upregulated in novel biomimetic 3D bioengineered platforms modeling the bevacizumab-resistant microenvironment. Second, tumor-initiating stem cells housed in the perivascular niche close to remaining blood vessels were enriched. Third, metabolic reprogramming assessed through real-time bioenergetic measurement and metabolomics upregulated glycolysis and suppressed oxidative phosphorylation. Single-cell sequencing of bevacizumab-resistant patient GBMs confirmed upregulated mesenchymal genes, particularly glycoprotein YKL-40 and transcription factor ZEB1, in later clones, implicating these changes as treatment-induced. Serum YKL-40 was elevated in bevacizumab-resistant versus bevacizumab-naïve patients. CRISPR and pharmacologic targeting of ZEB1 with honokiol reversed the mesenchymal gene expression and associated stem cell, invasion, and metabolic changes defining resistance. Honokiol caused greater cell death in bevacizumab-resistant than bevacizumab-responsive tumor cells, with surviving cells losing mesenchymal morphology. Employing YKL-40 as a resistance biomarker and ZEB1 as a target to prevent resistance could fulfill the promise of antiangiogenic therapy. SIGNIFICANCE: Bevacizumab resistance in GBM is associated with mesenchymal/glycolytic shifts involving YKL-40 and ZEB1. Targeting ZEB1 reduces bevacizumab-resistant GBM phenotypes. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/7/1498/F1.large.jpg.

The Phenotypes of Proliferating Glioblastoma Cells Reside on a Single Axis of Variation.

Although tumor-propagating cells can be derived from glioblastomas (GBM) of the proneural and mesenchymal subtypes, a glioma stem-like cell (GSC) of the classic subtype has not been identified. It is unclear whether mesenchymal GSCs (mGSC) and/or proneural GSCs (pGSC) alone are sufficient to generate the heterogeneity observed in GBM. We performed single-cell/single-nucleus RNA sequencing of 28 gliomas, and single-cell ATAC sequencing for 8 cases. We found that GBM GSCs reside on a single axis of variation, ranging from proneural to mesenchymal. In silico lineage tracing using both transcriptomics and genetics supports mGSCs as the progenitors of pGSCs. Dual inhibition of pGSC-enriched and mGSC-enriched growth and survival pathways provides a more complete treatment than combinations targeting one GSC phenotype alone. This study sheds light on a long-standing debate regarding lineage relationships among GSCs and presents a paradigm by which personalized combination therapies can be derived from single-cell RNA signatures, to overcome intratumor heterogeneity. SIGNIFICANCE: Tumor-propagating cells can be derived from mesenchymal and proneural glioblastomas. However, a stem cell of the classic subtype has yet to be demonstrated. We show that classic-subtype gliomas are comprised of proneural and mesenchymal cells. This study sheds light on a long-standing debate regarding lineage relationships between glioma cell types.See related commentary by Fine, p. 1650.This article is highlighted in the In This Issue feature, p. 1631.

Single-cell analysis of long non-coding RNAs in the developing human neocortex.

BACKGROUND: Long non-coding RNAs (lncRNAs) comprise a diverse class of transcripts that can regulate molecular and cellular processes in brain development and disease. LncRNAs exhibit cell type- and tissue-specific expression, but little is known about the expression and function of lncRNAs in the developing human brain. Furthermore, it has been unclear whether lncRNAs are highly expressed in subsets of cells within tissues, despite appearing lowly expressed in bulk populations. RESULTS: We use strand-specific RNA-seq to deeply profile lncRNAs from polyadenylated and total RNA obtained from human neocortex at different stages of development, and we apply this reference to analyze the transcriptomes of single cells. While lncRNAs are generally detected at low levels in bulk tissues, single-cell transcriptomics of hundreds of neocortex cells reveal that many lncRNAs are abundantly expressed in individual cells and are cell type-specific. Notably, LOC646329 is a lncRNA enriched in single radial glia cells but is detected at low abundance in tissues. CRISPRi knockdown of LOC646329 indicates that this lncRNA regulates cell proliferation. CONCLUSION: The discrete and abundant expression of lncRNAs among individual cells has important implications for both their biological function and utility for distinguishing neural cell types.