Human lung cancer harbors spatially organized stem-immunity hubs associated with response to immunotherapy.

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens and found an association with beneficial response to PD-1 blockade. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcome. This hub is distinct from mature tertiary lymphoid structures and is enriched for stem-like TCF7+PD-1+CD8+ T cells, activated CCR7+LAMP3+ dendritic cells and CCL19+ fibroblasts as well as chemokines that organize these cells. Within the stem-immunity hub, we find preferential interactions between CXCL10+ macrophages and TCF7-CD8+ T cells as well as between mature regulatory dendritic cells and TCF7+CD4+ and regulatory T cells. These results provide a picture of the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Disrupting CD38-driven T cell dysfunction restores sensitivity to cancer immunotherapy.

A central problem in cancer immunotherapy with immune checkpoint blockade (ICB) is the development of resistance, which affects 50% of patients with metastatic melanoma1,2. T cell exhaustion, resulting from chronic antigen exposure in the tumour microenvironment, is a major driver of ICB resistance3. Here, we show that CD38, an ecto-enzyme involved in nicotinamide adenine dinucleotide (NAD+) catabolism, is highly expressed in exhausted CD8+ T cells in melanoma and is associated with ICB resistance. Tumour-derived CD38hiCD8+ T cells are dysfunctional, characterised by impaired proliferative capacity, effector function, and dysregulated mitochondrial bioenergetics. Genetic and pharmacological blockade of CD38 in murine and patient-derived organotypic tumour models (MDOTS/PDOTS) enhanced tumour immunity and overcame ICB resistance. Mechanistically, disrupting CD38 activity in T cells restored cellular NAD+ pools, improved mitochondrial function, increased proliferation, augmented effector function, and restored ICB sensitivity. Taken together, these data demonstrate a role for the CD38-NAD+ axis in promoting T cell exhaustion and ICB resistance, and establish the efficacy of CD38 directed therapeutic strategies to overcome ICB resistance using clinically relevant, patient-derived 3D tumour models.

Metabolic predictors of response to immune checkpoint blockade therapy.

Metabolism of immune cells in the tumor microenvironment (TME) plays a critical role in cancer patient response to immune checkpoint inhibitors (ICI). Yet, a metabolic characterization of immune cells in the TME of patients treated with ICI is lacking. To bridge this gap we performed a semi-supervised analysis of ∼1700 metabolic genes using single-cell RNA-seq data of > 1 million immune cells from ∼230 samples of cancer patients treated with ICI. When clustering cells based on their metabolic gene expression, we found that similar immunological cellular states are found in different metabolic states. Most importantly, we found metabolic states that are significantly associated with patient response. We then built a metabolic predictor based on a dozen gene signature, which significantly differentiates between responding and non-responding patients across different cancer types (AUC = 0.8-0.92). Taken together, our results demonstrate the power of metabolism in predicting patient response to ICI.

Spatial analysis of human lung cancer reveals organized immune hubs enriched for stem-like CD8 T cells and associated with immunotherapy response.

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially-localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens, and found that they were associated with beneficial responses to PD-1-blockade. Immunity hubs were enriched for many interferon-stimulated genes, T cells in multiple differentiation states, and CXCL9/10/11 + macrophages that preferentially interact with CD8 T cells. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcomes, distinct from mature tertiary lymphoid structures, and enriched for stem-like TCF7+PD-1+ CD8 T cells and activated CCR7 + LAMP3 + dendritic cells, as well as chemokines that organize these cells. These results elucidate the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Targeting TBK1 to overcome resistance to cancer immunotherapy.

Despite the success of PD-1 blockade in melanoma and other cancers, effective treatment strategies to overcome resistance to cancer immunotherapy are lacking1,2. Here we identify the innate immune kinase TANK-binding kinase 1 (TBK1)3 as a candidate immune-evasion gene in a pooled genetic screen4. Using a suite of genetic and pharmacological tools across multiple experimental model systems, we confirm a role for TBK1 as an immune-evasion gene. Targeting TBK1 enhances responses to PD-1 blockade by decreasing the cytotoxicity threshold to effector cytokines (TNF and IFNγ). TBK1 inhibition in combination with PD-1 blockade also demonstrated efficacy using patient-derived tumour models, with concordant findings in matched patient-derived organotypic tumour spheroids and matched patient-derived organoids. Tumour cells lacking TBK1 are primed to undergo RIPK- and caspase-dependent cell death in response to TNF and IFNγ in a JAK-STAT-dependent manner. Taken together, our results demonstrate that targeting TBK1 is an effective strategy to overcome resistance to cancer immunotherapy.

Estimating tumor mutational burden from RNA-sequencing without a matched-normal sample.

Detection of somatic mutations using patients sequencing data has many clinical applications, including the identification of cancer driver genes, detection of mutational signatures, and estimation of tumor mutational burden (TMB). We have previously developed a tool for detection of somatic mutations using tumor RNA and a matched-normal DNA. Here, we further extend it to detect somatic mutations from RNA sequencing data without a matched-normal sample. This is accomplished via a machine-learning approach that classifies mutations as either somatic or germline based on various features. When applied to RNA-sequencing of >450 melanoma samples high precision and recall are achieved, and both mutational signatures and driver genes are correctly identified. Finally, we show that RNA-based TMB is significantly associated with patient survival, showing similar or higher significance level as compared to DNA-based TMB. Our pipeline can be utilized in many future applications, analyzing novel and existing datasets where only RNA is available.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

Alternative Splicing of the Inhibitory Immune Checkpoint Receptor SLAMF6 Generates a Dominant Positive Form, Boosting T-cell Effector Functions.

SLAMF6 is a homotypic receptor of the Ig-superfamily associated with progenitor-exhausted T cells. Here we show that in humans, SLAMF6 has three splice isoforms involving its V-domain. Although the canonical receptor inhibited T-cell activation through SAP recruitment, the short isoform SLAMF6Δ17-65 had a strong agonistic effect. The costimulatory action depended on protein phosphatase SHP1 and led to a cytotoxic molecular profile mediated by the expression of TBX21 and RUNX3. Patients treated with immune checkpoint blockade showed a shift toward SLAMF6Δ17-65 in peripheral blood T cells. We developed splice-switching antisense oligonucleotides (ASO) designed to target the relevant SLAMF6 splice junction. Our ASOs enhanced SLAMF6Δ17-65 expression in human tumor-infiltrating lymphocytes and improved their capacity to inhibit human melanoma in mice. The yin-yang relationship of SLAMF6 splice isoforms may represent a balancing mechanism that could be exploited to improve cancer immunotherapy.

B cells and tertiary lymphoid structures promote immunotherapy response.

Treatment with immune checkpoint blockade (ICB) has revolutionized cancer therapy. Until now, predictive biomarkers1-10 and strategies to augment clinical response have largely focused on the T cell compartment. However, other immune subsets may also contribute to anti-tumour immunity11-15, although these have been less well-studied in ICB treatment16. A previously conducted neoadjuvant ICB trial in patients with melanoma showed via targeted expression profiling17 that B cell signatures were enriched in the tumours of patients who respond to treatment versus non-responding patients. To build on this, here we performed bulk RNA sequencing and found that B cell markers were the most differentially expressed genes in the tumours of responders versus non-responders. Our findings were corroborated using a computational method (MCP-counter18) to estimate the immune and stromal composition in this and two other ICB-treated cohorts (patients with melanoma and renal cell carcinoma). Histological evaluation highlighted the localization of B cells within tertiary lymphoid structures. We assessed the potential functional contributions of B cells via bulk and single-cell RNA sequencing, which demonstrate clonal expansion and unique functional states of B cells in responders. Mass cytometry showed that switched memory B cells were enriched in the tumours of responders. Together, these data provide insights into the potential role of B cells and tertiary lymphoid structures in the response to ICB treatment, with implications for the development of biomarkers and therapeutic targets.

PD-1 blockade in subprimed CD8 cells induces dysfunctional PD-1+CD38hi cells and anti-PD-1 resistance.

Understanding resistance to antibody to programmed cell death protein 1 (PD-1; anti-PD-1) is crucial for the development of reversal strategies. In anti-PD-1-resistant models, simultaneous anti-PD-1 and vaccine therapy reversed resistance, while PD-1 blockade before antigen priming abolished therapeutic outcomes. This was due to induction of dysfunctional PD-1+CD38hi CD8+ cells by PD-1 blockade in suboptimally primed CD8 cell conditions induced by tumors. This results in erroneous T cell receptor signaling and unresponsiveness to antigenic restimulation. On the other hand, PD-1 blockade of optimally primed CD8 cells prevented the induction of dysfunctional CD8 cells, reversing resistance. Depleting PD-1+CD38hi CD8+ cells enhanced therapeutic outcomes. Furthermore, non-responding patients showed more PD-1+CD38+CD8+ cells in tumor and blood than responders. In conclusion, the status of CD8+ T cell priming is a major contributor to anti-PD-1 therapeutic resistance. PD-1 blockade in unprimed or suboptimally primed CD8 cells induces resistance through the induction of PD-1+CD38hi CD8+ cells that is reversed by optimal priming. PD-1+CD38hi CD8+ cells serve as a predictive and therapeutic biomarker for anti-PD-1 treatment. Sequencing of anti-PD-1 and vaccine is crucial for successful therapy.

RNA sequence analysis reveals macroscopic somatic clonal expansion across normal tissues.

How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.

Defining T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma.

Treatment of cancer has been revolutionized by immune checkpoint blockade therapies. Despite the high rate of response in advanced melanoma, the majority of patients succumb to disease. To identify factors associated with success or failure of checkpoint therapy, we profiled transcriptomes of 16,291 individual immune cells from 48 tumor samples of melanoma patients treated with checkpoint inhibitors. Two distinct states of CD8+ T cells were defined by clustering and associated with patient tumor regression or progression. A single transcription factor, TCF7, was visualized within CD8+ T cells in fixed tumor samples and predicted positive clinical outcome in an independent cohort of checkpoint-treated patients. We delineated the epigenetic landscape and clonality of these T cell states and demonstrated enhanced antitumor immunity by targeting novel combinations of factors in exhausted cells. Our study of immune cell transcriptomes from tumors demonstrates a strategy for identifying predictors, mechanisms, and targets for enhancing checkpoint immunotherapy.

Developmental and oncogenic programs in H3K27M gliomas dissected by single-cell RNA-seq.

Gliomas with histone H3 lysine27-to-methionine mutations (H3K27M-glioma) arise primarily in the midline of the central nervous system of young children, suggesting a cooperation between genetics and cellular context in tumorigenesis. Although the genetics of H3K27M-glioma are well characterized, their cellular architecture remains uncharted. We performed single-cell RNA sequencing in 3321 cells from six primary H3K27M-glioma and matched models. We found that H3K27M-glioma primarily contain cells that resemble oligodendrocyte precursor cells (OPC-like), whereas more differentiated malignant cells are a minority. OPC-like cells exhibit greater proliferation and tumor-propagating potential than their more differentiated counterparts and are at least in part sustained by PDGFRA signaling. Our study characterizes oncogenic and developmental programs in H3K27M-glioma at single-cell resolution and across genetic subclones, suggesting potential therapeutic targets in this disease.

Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer.

The diverse malignant, stromal, and immune cells in tumors affect growth, metastasis, and response to therapy. We profiled transcriptomes of ∼6,000 single cells from 18 head and neck squamous cell carcinoma (HNSCC) patients, including five matched pairs of primary tumors and lymph node metastases. Stromal and immune cells had consistent expression programs across patients. Conversely, malignant cells varied within and between tumors in their expression of signatures related to cell cycle, stress, hypoxia, epithelial differentiation, and partial epithelial-to-mesenchymal transition (p-EMT). Cells expressing the p-EMT program spatially localized to the leading edge of primary tumors. By integrating single-cell transcriptomes with bulk expression profiles for hundreds of tumors, we refined HNSCC subtypes by their malignant and stromal composition and established p-EMT as an independent predictor of nodal metastasis, grade, and adverse pathologic features. Our results provide insight into the HNSCC ecosystem and define stromal interactions and a p-EMT program associated with metastasis.

An integrated computational and experimental study uncovers FUT9 as a metabolic driver of colorectal cancer.

Metabolic alterations play an important role in cancer and yet, few metabolic cancer driver genes are known. Here we perform a combined genomic and metabolic modeling analysis searching for metabolic drivers of colorectal cancer. Our analysis predicts FUT9, which catalyzes the biosynthesis of Ley glycolipids, as a driver of advanced-stage colon cancer. Experimental testing reveals FUT9's complex dual role; while its knockdown enhances proliferation and migration in monolayers, it suppresses colon cancer cells expansion in tumorspheres and inhibits tumor development in a mouse xenograft models. These results suggest that FUT9's inhibition may attenuate tumor-initiating cells (TICs) that are known to dominate tumorspheres and early tumor growth, but promote bulk tumor cells. In agreement, we find that FUT9 silencing decreases the expression of the colorectal cancer TIC marker CD44 and the level of the OCT4 transcription factor, which is known to support cancer stemness. Beyond its current application, this work presents a novel genomic and metabolic modeling computational approach that can facilitate the systematic discovery of metabolic driver genes in other types of cancer.

Resistance to checkpoint blockade therapy through inactivation of antigen presentation.

Treatment with immune checkpoint blockade (CPB) therapies often leads to prolonged responses in patients with metastatic melanoma, but the common mechanisms of primary and acquired resistance to these agents remain incompletely characterized and have yet to be validated in large cohorts. By analyzing longitudinal tumor biopsies from 17 metastatic melanoma patients treated with CPB therapies, we observed point mutations, deletions or loss of heterozygosity (LOH) in beta-2-microglobulin (B2M), an essential component of MHC class I antigen presentation, in 29.4% of patients with progressing disease. In two independent cohorts of melanoma patients treated with anti-CTLA4 and anti-PD1, respectively, we find that B2M LOH is enriched threefold in non-responders (~30%) compared to responders (~10%) and associated with poorer overall survival. Loss of both copies of B2M is found only in non-responders. B2M loss is likely a common mechanism of resistance to therapies targeting CTLA4 or PD1.

Decoupling genetics, lineages, and microenvironment in IDH-mutant gliomas by single-cell RNA-seq.

Tumor subclasses differ according to the genotypes and phenotypes of malignant cells as well as the composition of the tumor microenvironment (TME). We dissected these influences in isocitrate dehydrogenase (IDH)-mutant gliomas by combining 14,226 single-cell RNA sequencing (RNA-seq) profiles from 16 patient samples with bulk RNA-seq profiles from 165 patient samples. Differences in bulk profiles between IDH-mutant astrocytoma and oligodendroglioma can be primarily explained by distinct TME and signature genetic events, whereas both tumor types share similar developmental hierarchies and lineages of glial differentiation. As tumor grade increases, we find enhanced proliferation of malignant cells, larger pools of undifferentiated glioma cells, and an increase in macrophage over microglia expression programs in TME. Our work provides a unifying model for IDH-mutant gliomas and a general framework for dissecting the differences among human tumor subclasses.

Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma.

Although human tumours are shaped by the genetic evolution of cancer cells, evidence also suggests that they display hierarchies related to developmental pathways and epigenetic programs in which cancer stem cells (CSCs) can drive tumour growth and give rise to differentiated progeny. Yet, unbiased evidence for CSCs in solid human malignancies remains elusive. Here we profile 4,347 single cells from six IDH1 or IDH2 mutant human oligodendrogliomas by RNA sequencing (RNA-seq) and reconstruct their developmental programs from genome-wide expression signatures. We infer that most cancer cells are differentiated along two specialized glial programs, whereas a rare subpopulation of cells is undifferentiated and associated with a neural stem cell expression program. Cells with expression signatures for proliferation are highly enriched in this rare subpopulation, consistent with a model in which CSCs are primarily responsible for fuelling the growth of oligodendroglioma in humans. Analysis of copy number variation (CNV) shows that distinct CNV sub-clones within tumours display similar cellular hierarchies, suggesting that the architecture of oligodendroglioma is primarily dictated by developmental programs. Subclonal point mutation analysis supports a similar model, although a full phylogenetic tree would be required to definitively determine the effect of genetic evolution on the inferred hierarchies. Our single-cell analyses provide insight into the cellular architecture of oligodendrogliomas at single-cell resolution and support the cancer stem cell model, with substantial implications for disease management.

A joint analysis of transcriptomic and metabolomic data uncovers enhanced enzyme-metabolite coupling in breast cancer.

Disrupted regulation of cellular processes is considered one of the hallmarks of cancer. We analyze metabolomic and transcriptomic profiles jointly collected from breast cancer and hepatocellular carcinoma patients to explore the associations between the expression of metabolic enzymes and the levels of the metabolites participating in the reactions they catalyze. Surprisingly, both breast cancer and hepatocellular tumors exhibit an increase in their gene-metabolites associations compared to noncancerous adjacent tissues. Following, we build predictors of metabolite levels from the expression of the enzyme genes catalyzing them. Applying these predictors to a large cohort of breast cancer samples we find that depleted levels of key cancer-related metabolites including glucose, glycine, serine and acetate are significantly associated with improved patient survival. Thus, we show that the levels of a wide range of metabolites in breast cancer can be successfully predicted from the transcriptome, going beyond the limited set of those measured.