Therapies that boost the anti-tumor responses of cytotoxic T lymphocytes (CTLs) have shown promise; however, clinical responses to the immunotherapeutic agents currently available vary considerably, and the molecular basis of this is unclear. We performed transcriptomic profiling of tumor-infiltrating CTLs from treatment-naive patients with lung cancer to define the molecular features associated with the robustness of anti-tumor immune responses. We observed considerable heterogeneity in the expression of molecules associated with activation of the T cell antigen receptor (TCR) and of immunological-checkpoint molecules such as 4-1BB, PD-1 and TIM-3. Tumors with a high density of CTLs showed enrichment for transcripts linked to tissue-resident memory cells (TRM cells), such as CD103, and CTLs from CD103hi tumors displayed features of enhanced cytotoxicity. A greater density of TRM cells in tumors was predictive of a better survival outcome in lung cancer, and this effect was independent of that conferred by CTL density. Here we define the ‘molecular fingerprint’ of tumor-infiltrating CTLs and identify potentially new targets for immunotherapy.

Determining mechanisms of resistance to aPD-1/PD-L1 immune-checkpoint immunotherapy is key to developing new treatment strategies. Cancer-associated fibroblasts (CAF) have many tumor-promoting functions and promote immune evasion through multiple mechanisms, but as yet, no CAF-specific inhibitors are clinically available. Here we generated CAF-rich murine tumor models (TC1, MC38, and 4T1) to investigate how CAFs influence the immune microenvironment and affect response to different immunotherapy modalities [anticancer vaccination, TC1 (HPV E7 DNA vaccine), aPD-1, and MC38] and found that CAFs broadly suppressed response by specifically excluding CD8 þ T cells from tumors (not CD4 þ T cells or macrophages); CD8 þ T-cell exclusion was similarly present in CAF-rich human tumors. RNA sequencing of CD8 þ T cells from CAF-rich murine tumors and immunochemistry analysis of human tumors identified significant upregulation of CTLA-4 in the absence of other exhaustion markers; inhibiting CTLA-4 with a nondepleting antibody overcame the CD8 þ T-cell exclusion effect without affecting Tregs. We then examined the potential for CAF targeting, focusing on the ROS-producing enzyme NOX4, which is upregulated by CAF in many human cancers, and compared this with TGFb1 inhibition, a key regulator of the CAF phenotype. siRNA knockdown or pharmacologic inhibition [GKT137831 (Setanaxib)] of NOX4 "normalized" CAF to a quiescent phenotype and promoted intratumoral CD8 þ T-cell infiltration, overcoming the exclusion effect; TGFb1 inhibition could prevent, but not reverse, CAF differentiation. Finally, NOX4 inhibition restored immunotherapy response in CAFrich tumors. These findings demonstrate that CAF-mediated immunotherapy resistance can be effectively overcome through NOX4 inhibition and could improve outcome in a broad range of cancers. Significance: NOX4 is critical for maintaining the immunesuppressive CAF phenotype in tumors. Pharmacologic inhibition of NOX4 potentiates immunotherapy by overcoming CAFmediated CD8 þ T-cell exclusion.

High numbers of tissue-resident memory T (TRM) cells are associated with better clinical outcomes in cancer patients. However, the molecular characteristics that drive their efficient immune response to tumors are poorly understood. Here, single-cell and bulk transcriptomic analysis of TRM and non-TRM cells present in tumor and normal lung tissue from patients with lung cancer revealed that PD-1–expressing TRM cells in tumors were clonally expanded and enriched for transcripts linked to cell proliferation and cytotoxicity when compared with PD-1–expressing non-TRM cells. This feature was more prominent in the TRM cell subset coexpressing PD-1 and TIM-3, and it was validated by functional assays ex vivo and also reflected in their chromatin accessibility profile. This PD-1+TIM-3+ TRM cell subset was enriched in responders to PD-1 inhibitors and in tumors with a greater magnitude of CTL responses. These data highlight that not all CTLs expressing PD-1 are dysfunctional; on the contrary, TRM cells with PD-1 expression were enriched for features suggestive of superior functionality.

BackgroundThe role of tumor-associated macrophages (TAMs) in determining the outcome between the antitumor effects of the adaptive immune system and the tumor’s anti-immunity stratagems, is controversial. Macrophages modulate their activities and phenotypes by integration of signals in the tumor microenvironment. Depending on how macrophages are activated, they may adopt so-called M1-like, antitumor or M2-like, protumor profiles. In many solid tumors, a dominance of M2-like macrophages is associated with poor outcomes but in some tumor types, strong M1-like profiles are linked to better outcomes. We aimed to investigate the interrelationship of these TAM populations to establish how they modulate the efficacy of the adaptive immune system in early lung cancer.MethodsMacrophages from matched lung (non-tumor-associated macrophages (NTAMs)) and tumor samples (TAMs) from resected lung cancers were assessed by bulk and single-cell transcriptomic analysis. Protein expression of genes characteristic of M1-like (chemokine (C-X-C motif) ligand 9) or M2-like (matrix metallopeptidase 12) functions was confirmed by confocal microscopy. Immunohistochemistry related the distribution of TAM transcriptomic signatures to density of CD8+ tissue-resident memory T cells (TRM) in tumors and survival data from an independent cohort of 393 patients with lung cancer.ResultsTAMs have significantly different transcriptomic profiles from NTAMs with >1000 differentially expressed genes. TAMs displayed a strong M2-like signature with no significant variation between patients. However, single-cell RNA-sequencing supported by immuno-stained cells revealed that additionally, in 25% of patients the M2-like TAMs also co-expressed a strong/hot M1-like signature (M1hot). Importantly, there was a strong association between the density of M1hot TAMs and TRM cells in tumors that was in turn linked to better survival. Our data suggest a mechanism by which M1hot TAMs may recruit TRM cells via CXCL9 expression and sustain them by making available more of the essential fatty acids on which TRM depend.ConclusionsWe showed that in early lung cancer, expression of M1-like and M2-like gene signatures are not mutually exclusive since the same TAMs can simultaneously display both gene-expression profiles. The presence of M1hot TAMs was associated with a strong TRM tumor-infiltrate and better outcomes. Thus, therapeutic approaches to re-program TAMs to an M1hot phenotype are likely to augment the adaptive antitumor responses.

Human papilloma virus (HPV)-associated head and neck squamous cell carcinoma (HNSCC) has a better prognosis than it's HPV negative (HPV(−)) counterpart. This may be due to the higher numbers of tumor-infiltrating lymphocytes (TILs) in HPV positive (HPV(+)) tumors. RNA-Sequencing (RNA-Seq) was used to evaluate whether the differences in clinical behaviour simply reflect a numerical difference in TILs or whether there is a fundamental behavioural difference between TILs in these two settings. Thirty-nine HNSCC tumors were scored for TIL density by immunohistochemistry. After the removal of 16 TILlow tumors, RNA-Seq analysis was performed on 23 TILhigh/med tumors (HPV(+) n=10 and HPV(−) n=13). Using EdgeR, differentially expressed genes (DEG) were identified. Immune subset analysis was performed using Functional Analysis of Individual RNA-Seq/ Microarray Expression (FAIME) and immune gene RNA transcript count analysis. In total, 1,634 DEGs were identified, with a dominant immune signature observed in HPV(+) tumors. After normalizing the expression profiles to account for differences in B- and T-cell number, 437 significantly DEGs remained. A B-cell associated signature distinguished HPV(+) from HPV(−) tumors, and included the DEGs CD200, GGA2, ADAM28, STAG3, SPIB, VCAM1, BCL2 and ICOSLG; the immune signal relative to T-cells was qualitatively similar between TILs of both tumor cohorts. Our findings were validated and confirmed in two independent cohorts using TCGA data and tumor-infiltrating B-cells from additional HPV(+) HNSCC patients. A B-cell associated signal segregated tumors relative to HPV status. Our data suggests that the role of B-cells in the adaptive immune response to HPV(+) HNSCC requires re-assessment.

Antibodies that block T-cell-regulatory checkpoints have recently emerged as a transformative approach to cancer treatment. However, the clinical efficacy of checkpoint blockade depends upon inherent tumor immunogenicity, with variation in infiltrating T cells contributing to differences in objective response rates. Here, we sought to understand the molecular correlates of tumor-infiltrating T lymphocytes (TIL) in squamous cell carcinoma (SCC), using a systems biologic approach to integrate publicly available omics datasets with histopathologic features. We provide evidence that links TIL abundance and therapeutic outcome to the regulation of tumor glycolysis by EGFR and HIF, both of which are attractive molecular targets for use in combination with immunotherapeutics. Cancer Res; 76(14); 4136-48. Ó2016 AACR.

Background Metabolic changes in tumour cells are used in clinical imaging and may provide potential therapeutic targets. Human papillomavirus (HPV) status is important in classifying head and neck cancers (HNSCC), identifying a distinct clinical phenotype; metabolic differences between these HNSCC subtypes remain poorly understood. Methods We used RNA sequencing to classify the metabolic expression profiles of HPV +ve and HPV −ve HNSCC, performed a meta-analysis on FDG-PET imaging characteristics and correlated results with in vitro extracellular flux analysis of HPV −ve and HPV +ve HNSCC cell lines. The monocarboxylic acid transporter-1 (MCT1) was identified as a potential metabolic target and tested in functional assays. Results Specific metabolic profiles were associated with HPV status, not limited to carbohydrate metabolism. There was dominance of all energy pathways in HPV-negative disease, with elevated expression of genes associated with glycolysis and oxidative phosphorylation. In vitro analysis confirmed comparative increased rates of oxidative phosphorylation and glycolysis in HPV-negative cell lines. PET SUV(max) scores however were unable to reliably differentiate between HPV-positive and HPV-negative tumours. MCT1 expression was significantly increased in HPV-negative tumours, and inhibition suppressed tumour cell invasion, colony formation and promoted radiosensitivity. Conclusion HPV-positive and negative HNSCC have different metabolic profiles which may have potential therapeutic applications.

Single-cell RNA sequencing (scRNA-Seq) provides a valuable platform for characterising multicellular ecosystems. Fibroblasts are a heterogeneous cell type involved in many physiological and pathological processes, but remain poorly-characterised. Analysis of fibroblasts is challenging: these cells are difficult to isolate from tissues, and are therefore commonly under-represented in scRNA-seq datasets. Here, we describe an optimised approach for fibroblast isolation from human lung tissues. We demonstrate the potential for this procedure in characterising stromal cell phenotypes using scRNA-Seq, analyse the effect of tissue disaggregation on gene expression, and optimise data processing to improve clustering quality. We also assess the impact of in vitro culture conditions on stromal cell gene expression and proliferation, showing that altering these conditions can skew phenotypes.