Intra-tumoral T cells in pediatric brain tumors display clonal expansion and effector properties.

Brain tumors in children are a devastating disease in a high proportion of patients. Owing to inconsistent results in clinical trials in unstratified patients, the role of immunotherapy remains unclear. We performed an in-depth survey of the single-cell transcriptomes and clonal relationship of intra-tumoral T cells from children with brain tumors. Our results demonstrate that a large fraction of T cells in the tumor tissue are clonally expanded with the potential to recognize tumor antigens. Such clonally expanded T cells display enrichment of transcripts linked to effector function, tissue residency, immune checkpoints and signatures of neoantigen-specific T cells and immunotherapy response. We identify neoantigens in pediatric brain tumors and show that neoantigen-specific T cell gene signatures are linked to better survival outcomes. Notably, among the patients in our cohort, we observe substantial heterogeneity in the degree of clonal expansion and magnitude of T cell response. Our findings suggest that characterization of intra-tumoral T cell responses may enable selection of patients for immunotherapy, an approach that requires prospective validation in clinical trials.

JAML immunotherapy targets recently activated tumor-infiltrating CD8+ T cells.

Junctional adhesion molecule-like protein (JAML) serves as a co-stimulatory molecule in γδ T cells. While it has recently been described as a cancer immunotherapy target in mice, its potential to cause toxicity, specific mode of action with regard to its cellular targets, and whether it can be targeted in humans remain unknown. Here, we show that JAML is induced by T cell receptor engagement, reveal that this induction is linked to cis-regulatory interactions between the CD3D and JAML gene loci. When compared with other immunotherapy targets plagued by low target specificity and end-organ toxicity, we find JAML to be mostly restricted to and highly expressed by tissue-resident memory CD8+ T cells in multiple cancer types. By delineating the key cellular targets and functional consequences of agonistic anti-JAML therapy in a murine melanoma model, we show its specific mode of action and the reason for its synergistic effects with anti-PD-1.

Single-cell analysis reveals prognostic fibroblast subpopulations linked to molecular and immunological subtypes of lung cancer.

Fibroblasts are poorly characterised cells that variably impact tumour progression. Here, we use single cell RNA-sequencing, multiplexed immunohistochemistry and digital cytometry (CIBERSORTx) to identify and characterise three major fibroblast subpopulations in human non-small cell lung cancer: adventitial, alveolar and myofibroblasts. Alveolar and adventitial fibroblasts (enriched in control tissue samples) localise to discrete spatial niches in histologically normal lung tissue and indicate improved overall survival rates when present in lung adenocarcinomas (LUAD). Trajectory inference identifies three phases of control tissue fibroblast activation, leading to myofibroblast enrichment in tumour samples: initial upregulation of inflammatory cytokines, followed by stress-response signalling and ultimately increased expression of fibrillar collagens. Myofibroblasts correlate with poor overall survival rates in LUAD, associated with loss of epithelial differentiation, TP53 mutations, proximal molecular subtypes and myeloid cell recruitment. In squamous carcinomas myofibroblasts were not prognostic despite being transcriptomically equivalent. These findings have important implications for developing fibroblast-targeting strategies for cancer therapy.

Single-cell immune profiling reveals long-term changes in myeloid cells and identifies a novel subset of CD9+ monocytes associated with COVID-19 hospitalization.

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection can result in severe immune dysfunction, hospitalization, and death. Many patients also develop long-COVID-19, experiencing symptoms months after infection. Although significant progress has been made in understanding the immune response to acute SARS-CoV-2 infection, gaps remain in our knowledge of how innate immunity influences disease kinetics and severity. We hypothesized that cytometry by time-of-flight analysis of PBMCs from healthy and infected subjects would identify novel cell surface markers and innate immune cell subsets associated with COVID-19 severity. In this pursuit, we identified monocyte and dendritic cell subsets that changed in frequency during acute SARS-CoV-2 infection and correlated with clinical parameters of disease severity. Subsets of nonclassical monocytes decreased in frequency in hospitalized subjects, yet increased in the most severe patients and positively correlated with clinical values associated with worse disease severity. CD9, CD163, PDL1, and PDL2 expression significantly increased in hospitalized subjects, and CD9 and 6-Sulfo LacNac emerged as the markers that best distinguished monocyte subsets amongst all subjects. CD9+ monocytes remained elevated, whereas nonclassical monocytes remained decreased, in the blood of hospitalized subjects at 3-4 months postinfection. Finally, we found that CD9+ monocytes functionally released more IL-8 and MCP-1 after LPS stimulation. This study identifies new monocyte subsets present in the blood of COVID-19 patients that correlate with disease severity, and links CD9+ monocytes to COVID-19 progression.

HIF activation enhances FcγRIIb expression on mononuclear phagocytes impeding tumor targeting antibody immunotherapy.

BACKGROUND: Hypoxia is a hallmark of the tumor microenvironment (TME) and in addition to altering metabolism in cancer cells, it transforms tumor-associated stromal cells. Within the tumor stromal cell compartment, tumor-associated macrophages (TAMs) provide potent pro-tumoral support. However, TAMs can also be harnessed to destroy tumor cells by monoclonal antibody (mAb) immunotherapy, through antibody dependent cellular phagocytosis (ADCP). This is mediated via antibody-binding activating Fc gamma receptors (FcγR) and impaired by the single inhibitory FcγR, FcγRIIb. METHODS: We applied a multi-OMIC approach coupled with in vitro functional assays and murine tumor models to assess the effects of hypoxia inducible factor (HIF) activation on mAb mediated depletion of human and murine cancer cells. For mechanistic assessments, siRNA-mediated gene silencing, Western blotting and chromatin immune precipitation were utilized to assess the impact of identified regulators on FCGR2B gene transcription. RESULTS: We report that TAMs are FcγRIIbbright relative to healthy tissue counterparts and under hypoxic conditions, mononuclear phagocytes markedly upregulate FcγRIIb. This enhanced FcγRIIb expression is transcriptionally driven through HIFs and Activator protein 1 (AP-1). Importantly, this phenotype reduces the ability of macrophages to eliminate anti-CD20 monoclonal antibody (mAb) opsonized human chronic lymphocytic leukemia cells in vitro and EL4 lymphoma cells in vivo in human FcγRIIb+/+ transgenic mice. Furthermore, post-HIF activation, mAb mediated blockade of FcγRIIb can partially restore phagocytic function in human monocytes. CONCLUSION: Our findings provide a detailed molecular and cellular basis for hypoxia driven resistance to antitumor mAb immunotherapy, unveiling a hitherto unexplored aspect of the TME. These findings provide a mechanistic rationale for the modulation of FcγRIIb expression or its blockade as a promising strategy to enhance approved and novel mAb immunotherapies.

Analysis of Immune Landscape in Pancreatic and Ileal Neuroendocrine Tumours Demonstrates an Immune Cold Tumour Microenvironment.

INTRODUCTION: Neuroendocrine tumours (NETs) are rare tumours with an increasing incidence. While low- and intermediate-grade pancreatic NET (PanNET) and small intestinal NET (siNET) are slow growing, they have a relatively high rate of metastasizing to the liver, leading to substantially worse outcomes. In many solid tumours, the outcome is determined by the quality of the antitumour immune response. However, the quality and significance of antitumour responses in NETs are incompletely understood. This study provides clinico-pathological analyses of the tumour immune microenvironment in PanNET and siNETs. METHODS: Formalin-fixed paraffin-embedded tissue from consecutive resected PanNETs (61) and siNETs (131) was used to construct tissue microarrays (TMAs); 1-mm cores were taken from the tumour centre, stroma, tumour edge, and adjacent healthy tissue. TMAs were stained with antibodies against CD8, CD4, CD68, FoxP3, CD20, and NCR1. T-cell counts were compared with counts from lung cancers. RESULTS: For PanNET, median counts were CD8+ 35.4 cells/mm2, CD4+ 7.6 cells/mm2, and CD68+ macrophages 117.7 cells/mm2. For siNET, there were CD8+ 39.2 cells/mm2, CD4+ 24.1 cells/mm2, and CD68+ 139.2 cells/mm2. The CD8+ cell density in the tumour and liver metastases were significantly lower than in the adjacent normal tissues, without evidence of a cell-rich area at the tumour edge that might have suggested immune exclusion. T-cell counts in lung cancer were significantly higher than those in PanNET and siNETs: CD8+ 541 cells/mm2 and CD4+ 861 cells/mm2 (p ≤ 0.0001). CONCLUSION: PanNETs and siNETs are immune cold with no evidence of T cell exclusion; the low density of immune infiltrates indicates poor antitumour immune responses.

Intratumoral follicular regulatory T cells curtail anti-PD-1 treatment efficacy.

Immune-checkpoint blockade (ICB) has shown remarkable clinical success in boosting antitumor immunity. However, the breadth of its cellular targets and specific mode of action remain elusive. We find that tumor-infiltrating follicular regulatory T (TFR) cells are prevalent in tumor tissues of several cancer types. They are primarily located within tertiary lymphoid structures and exhibit superior suppressive capacity and in vivo persistence as compared with regulatory T cells, with which they share a clonal and developmental relationship. In syngeneic tumor models, anti-PD-1 treatment increases the number of tumor-infiltrating TFR cells. Both TFR cell deficiency and the depletion of TFR cells with anti-CTLA-4 before anti-PD-1 treatment improve tumor control in mice. Notably, in a cohort of 271 patients with melanoma, treatment with anti-CTLA-4 followed by anti-PD-1 at progression was associated with better a survival outcome than monotherapy with anti-PD-1 or anti-CTLA-4, anti-PD-1 followed by anti-CTLA-4 at progression or concomitant combination therapy.

Severely ill COVID-19 patients display impaired exhaustion features in SARS-CoV-2-reactive CD8+ T cells.

The molecular properties of CD8+ T cells that respond to SARS-CoV-2 infection are not fully known. Here, we report on the single-cell transcriptomes of >80,000 virus-reactive CD8+ T cells, obtained using a modified Antigen-Reactive T cell Enrichment (ARTE) assay, from 39 COVID-19 patients and 10 healthy subjects. COVID-19 patients segregated into two groups based on whether the dominant CD8+ T cell response to SARS-CoV-2 was 'exhausted' or not. SARS-CoV-2-reactive cells in the exhausted subset were increased in frequency and displayed lesser cytotoxicity and inflammatory features in COVID-19 patients with mild compared to severe illness. In contrast, SARS-CoV-2-reactive cells in the dominant non-exhausted subset from patients with severe disease showed enrichment of transcripts linked to co-stimulation, pro-survival NF-κB signaling, and anti-apoptotic pathways, suggesting the generation of robust CD8+ T cell memory responses in patients with severe COVID-19 illness. CD8+ T cells reactive to influenza and respiratory syncytial virus from healthy subjects displayed polyfunctional features and enhanced glycolysis. Cells with such features were largely absent in SARS-CoV-2-reactive cells from both COVID-19 patients and healthy controls non-exposed to SARS-CoV-2. Overall, our single-cell analysis revealed substantial diversity in the nature of CD8+ T cells responding to SARS-CoV-2.

Imbalance of Regulatory and Cytotoxic SARS-CoV-2-Reactive CD4+ T Cells in COVID-19.

The contribution of CD4+ T cells to protective or pathogenic immune responses to SARS-CoV-2 infection remains unknown. Here, we present single-cell transcriptomic analysis of >100,000 viral antigen-reactive CD4+ T cells from 40 COVID-19 patients. In hospitalized patients compared to non-hospitalized patients, we found increased proportions of cytotoxic follicular helper cells and cytotoxic T helper (TH) cells (CD4-CTLs) responding to SARS-CoV-2 and reduced proportion of SARS-CoV-2-reactive regulatory T cells (TREG). Importantly, in hospitalized COVID-19 patients, a strong cytotoxic TFH response was observed early in the illness, which correlated negatively with antibody levels to SARS-CoV-2 spike protein. Polyfunctional TH1 and TH17 cell subsets were underrepresented in the repertoire of SARS-CoV-2-reactive CD4+ T cells compared to influenza-reactive CD4+ T cells. Together, our analyses provide insights into the gene expression patterns of SARS-CoV-2-reactive CD4+ T cells in distinct disease severities.

Single-Cell Transcriptomic Analysis of SARS-CoV-2 Reactive CD4 + T Cells.

The contribution of CD4+ T cells to protective or pathogenic immune responses to SARS-CoV-2 infection remains unknown. Here, we present large-scale single-cell transcriptomic analysis of viral antigen-reactive CD4+ T cells from 32 COVID-19 patients. In patients with severe disease compared to mild disease, we found increased proportions of cytotoxic follicular helper (TFH) cells and cytotoxic T helper cells (CD4-CTLs) responding to SARS-CoV-2, and reduced proportion of SARS-CoV-2 reactive regulatory T cells. Importantly, the CD4-CTLs were highly enriched for the expression of transcripts encoding chemokines that are involved in the recruitment of myeloid cells and dendritic cells to the sites of viral infection. Polyfunctional T helper (TH)1 cells and TH17 cell subsets were underrepresented in the repertoire of SARS-CoV-2-reactive CD4+ T cells compared to influenza-reactive CD4+ T cells. Together, our analyses provide so far unprecedented insights into the gene expression patterns of SARS-CoV-2 reactive CD4+ T cells in distinct disease severities. Funding: This work was funded by NIH grants U19AI142742 (P.V., A.S., C.H.O), U19AI118626 (P.V., A.S., G.S.), R01HL114093 (P.V., F.A., G.S.,), R35-GM128938 (F.A), S10RR027366 (BD FACSAria-II), S10OD025052 (Illumina Novaseq6000), the William K. Bowes Jr Foundation (P.V.), and Whittaker foundation (P.V., C.H.O.). Supported by the Wessex Clinical Research Network and National Institute of Health Research UK. Conflict of Interest: The authors declare no competing financial interests. Ethical Approval: Ethical approval for this study from the Berkshire Research Ethics Committee 20/SC/0155 and the Ethics Committee of La Jolla Institute for Immunology (LJI) was in place. Written consent was obtained from all subjects.

Severely ill COVID-19 patients display augmented functional properties in SARS-CoV-2-reactive CD8 + T cells.

The molecular properties of CD8 + T cells that respond to SARS-CoV-2 infection are not fully known. Here, we report on the single-cell transcriptomes of >80,000 virus-reactive CD8 + T cells from 39 COVID-19 patients and 10 healthy subjects. COVID-19 patients segregated into two groups based on whether the dominant CD8 + T cell response to SARS-CoV-2 was 'exhausted' or not. SARS-CoV-2-reactive cells in the exhausted subset were increased in frequency and displayed lesser cytotoxicity and inflammatory features in COVID-19 patients with mild compared to severe illness. In contrast, SARS-CoV-2-reactive cells in the non-exhausted subsets from patients with severe disease showed enrichment of transcripts linked to co-stimulation, pro-survival NF-κB signaling, and anti-apoptotic pathways, suggesting the generation of robust CD8 + T cell memory responses in patients with severe COVID-19 illness. CD8 + T cells reactive to influenza and respiratory syncytial virus from healthy subjects displayed polyfunctional features. Cells with such features were mostly absent in SARS-CoV-2 responsive cells from both COVID-19 patients and healthy controls non-exposed to SARS-CoV-2. Overall, our single-cell analysis revealed substantial diversity in the nature of CD8 + T cells responding to SARS-CoV-2.

M1hot tumor-associated macrophages boost tissue-resident memory T cells infiltration and survival in human lung cancer.

BACKGROUND: The 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. METHODS: Macrophages 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. RESULTS: TAMs 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. CONCLUSIONS: We 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.

Single-cell transcriptomic analysis of SARS-CoV-2 reactive CD4 + T cells.

The contribution of CD4 + T cells to protective or pathogenic immune responses to SARS-CoV-2 infection remains unknown. Here, we present large-scale single-cell transcriptomic analysis of viral antigen-reactive CD4 + T cells from 32 COVID-19 patients. In patients with severe disease compared to mild disease, we found increased proportions of cytotoxic follicular helper (T FH ) cells and cytotoxic T helper cells (CD4-CTLs) responding to SARS-CoV-2, and reduced proportion of SARS-CoV-2 reactive regulatory T cells. Importantly, the CD4-CTLs were highly enriched for the expression of transcripts encoding chemokines that are involved in the recruitment of myeloid cells and dendritic cells to the sites of viral infection. Polyfunctional T helper (T H )1 cells and T H 17 cell subsets were underrepresented in the repertoire of SARS-CoV-2-reactive CD4 + T cells compared to influenza-reactive CD4 + T cells. Together, our analyses provide so far unprecedented insights into the gene expression patterns of SARS-CoV-2 reactive CD4 + T cells in distinct disease severities.

Autophagy inhibition-mediated epithelial-mesenchymal transition augments local myofibroblast differentiation in pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF), the prototypic progressive fibrotic interstitial lung disease, is thought to be a consequence of repetitive micro-injuries to an ageing, susceptible alveolar epithelium. Ageing is a risk factor for IPF and incidence has been demonstrated to increase with age. Decreased (macro)autophagy with age has been reported extensively in a variety of systems and diseases, including IPF. However, it is undetermined whether the role of faulty autophagy is causal or coincidental in the context of IPF. Here, we report that in alveolar epithelial cells inhibition of autophagy promotes epithelial-mesenchymal transition (EMT), a process implicated in embryonic development, wound healing, cancer metastasis and fibrosis. We further demonstrate that this is attained, at least in part, by increased p62/SQSTM1 expression that promotes p65/RELA mediated-transactivation of an EMT transcription factor, Snail2 (SNAI2), which not only controls EMT but also regulates the production of locally acting profibrogenic mediators. Our data suggest that reduced autophagy induces EMT of alveolar epithelial cells and can contribute to fibrosis via aberrant epithelial-fibroblast crosstalk.

Single-cell transcriptomic analysis of tissue-resident memory T cells in human lung cancer.

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.

Paracrine signalling during ZEB1-mediated epithelial-mesenchymal transition augments local myofibroblast differentiation in lung fibrosis.

The contribution of epithelial-mesenchymal transition (EMT) to human lung fibrogenesis is controversial. Here we provide evidence that ZEB1-mediated EMT in human alveolar epithelial type II (ATII) cells contributes to the development of lung fibrosis by paracrine signalling to underlying fibroblasts. Activation of EGFR-RAS-ERK signalling in ATII cells induced EMT via ZEB1. ATII cells had extremely low extracellular matrix gene expression even after induction of EMT, however conditioned media from ATII cells undergoing RAS-induced EMT augmented TGFß-induced profibrogenic responses in lung fibroblasts. This epithelial-mesenchymal crosstalk was controlled by ZEB1 via the expression of tissue plasminogen activator (tPA). In human fibrotic lung tissue, nuclear ZEB1 expression was detected in alveolar epithelium adjacent to sites of extracellular matrix (ECM) deposition, suggesting that ZEB1-mediated paracrine signalling has the potential to contribute to early fibrotic changes in the lung interstitium. Targeting this novel ZEB1 regulatory axis may be a viable strategy for the treatment of pulmonary fibrosis.

Nanoscale dysregulation of collagen structure-function disrupts mechano-homeostasis and mediates pulmonary fibrosis.

Matrix stiffening with downstream activation of mechanosensitive pathways is strongly implicated in progressive fibrosis; however, pathologic changes in extracellular matrix (ECM) that initiate mechano-homeostasis dysregulation are not defined in human disease. By integrated multiscale biomechanical and biological analyses of idiopathic pulmonary fibrosis lung tissue, we identify that increased tissue stiffness is a function of dysregulated post-translational collagen cross-linking rather than any collagen concentration increase whilst at the nanometre-scale collagen fibrils are structurally and functionally abnormal with increased stiffness, reduced swelling ratio, and reduced diameter. In ex vivo and animal models of lung fibrosis, dual inhibition of lysyl oxidase-like (LOXL) 2 and LOXL3 was sufficient to normalise collagen fibrillogenesis, reduce tissue stiffness, and improve lung function in vivo. Thus, in human fibrosis, altered collagen architecture is a key determinant of abnormal ECM structure-function, and inhibition of pyridinoline cross-linking can maintain mechano-homeostasis to limit the self-sustaining effects of ECM on progressive fibrosis.

Evaluating the effect of immune cells on the outcome of patients with mesothelioma.

BACKGROUND: We systematically assessed the prognostic and predictive value of infiltrating adaptive and innate immune cells in a large cohort of patients with advanced mesothelioma. METHODS: A tissue microarray from 302 samples was constructed. Markers of adaptive immune response in T-cells (CD8+, FOXP3+, CD4+, CD45RO+, CD3+) and B-cells (CD20+), and of innate immune response; neutrophils (NP57+), natural killer cells (CD56+) and macrophages (CD68+) were evaluated. RESULTS: We found that in the epithelioid tumours, high CD4+ and CD20+ counts, and low FOXP3+, CD68+ and NP57+ counts linked to better outcome. In the non-epithelioid group low CD8+ and low FOXP3+ counts were beneficial.On multivariate analysis low FOXP3+ remained independently associated with survival in both groups. In the epithelioid group additionally high CD4+, high CD20+, and low NP57+ counts were prognostic. CONCLUSIONS: Our data demonstrate for the first time, in predominately advanced disease, the association of key markers of adaptive and innate immunity with survival and the differential effect of histology. A better understanding of the immunological drivers of the different subtypes of mesothelioma will assist prognostication and disease-specific clinical decision-making.

Tissue-resident memory features are linked to the magnitude of cytotoxic T cell responses in human lung cancer.

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.