Single-cell meta-analyses reveal responses of tumor-reactive CXCL13+ T cells to immune-checkpoint blockade

Cancer-Associated Fibroblast-Induced Remodeling of Tumor Microenvironment in Recurrent Bladder Cancer

Bladder carcinoma (BC) recurrence is a major clinical challenge, and targeting the tumor microenvironment (TME) is a promising therapy. However, the relationship between individual TME components, particularly cancer-associated fibroblasts (CAFs), and tumor recurrence is unclear. Here, TME heterogeneity in primary and recurrent BC is investigated using single-cell RNA sequence profiling of 62 460 cells. Two cancer stem cell (CSC) subtypes are identified in recurrent BC. An inflammatory CAF subtype, ICAM1+ iCAFs, specifically associated with BC recurrence is also identified. iCAFs are found to secrete FGF2, which acts on the CD44 receptor of rCSC-M, thereby maintaining tumor stemness and epithelial-mesenchymal transition. Additionally, THBS1+ monocytes, a group of myeloid-derived suppressor cells (MDSCs), are enriched in recurrent BC and interacted with CAFs. ICAM1+ iCAFs are found to secrete CCL2, which binds to CCR2 in MDSCs. Moreover, elevated STAT3, NFKB2, VEGFA, and CTGF levels in iCAFs reshape the TME in recurrent tumors. CCL2 inhibition in an in situ BC mouse model suppressed tumor growth, decreased MDSCs and Tregs, and fostered tumor immune suppression. The study results highlight the role of iCAFs in TME cell-cell crosstalk during recurrent BC. The identification of pivotal signaling factors driving BC relapse is promising for the development of novel therapies.

Immune landscape of tertiary lymphoid structures in hepatocellular carcinoma (HCC) treated with neoadjuvant immune checkpoint blockade

Neoadjuvant immunotherapy is thought to produce long-term remissions through induction of antitumor immune responses before removal of the primary tumor. Tertiary lymphoid structures (TLS), germinal center-like structures that can arise within tumors, may contribute to the establishment of immunological memory in this setting, but understanding of their role remains limited. Here, we investigated the contribution of TLS to antitumor immunity in hepatocellular carcinoma (HCC) treated with neoadjuvant immunotherapy. We found that neoadjuvant immunotherapy induced the formation of TLS, which were associated with superior pathologic response, improved relapse free survival, and expansion of the intratumoral T and B cell repertoire. While TLS in viable tumor displayed a highly active mature morphology, in areas of tumor regression we identified an involuted TLS morphology, which was characterized by dispersion of the B cell follicle and persistence of a T cell zone enriched for ongoing antigen presentation and T cell-mature dendritic cell interactions. Involuted TLS showed increased expression of T cell memory markers and expansion of CD8+ cytotoxic and tissue resident memory clonotypes. Collectively, these data reveal the circumstances of TLS dissolution and suggest a functional role for late-stage TLS as sites of T cell memory formation after elimination of viable tumor.

The tumor microenvironment state associates with response to HPV therapeutic vaccination in patients with respiratory papillomatosis

Recurrent respiratory papillomatosis (RRP) is a rare, debilitating neoplastic disorder caused by chronic infection with human papillomavirus (HPV) type 6 or 11 and characterized by growth of papillomas in the upper aerodigestive tract. There is no approved medical therapy, and patients require repeated debulking procedures to maintain voice and airway function. PRGN-2012 is a gorilla adenovirus immune-therapeutic capable of enhancing HPV 6/11-specific T cell immunity. This first-in-human, phase 1 study (NCT04724980) of adjuvant PRGN-2012 treatment in adult patients with severe, aggressive RRP demonstrates the overall safety and clinically meaningful benefit observed with PRGN-2012, with a 50% complete response rate in patients treated at the highest dose. Responders demonstrate greater expansion of peripheral HPV-specific T cells compared with nonresponders. Additional correlative studies identify an association between reduced baseline papilloma HPV gene expression, greater interferon responses and expression of CXCL9 and CXCL10, and greater papilloma T cell infiltration in responders. Conversely, nonresponders were characterized by greater HPV and CXCL8 gene expression, increased neutrophilic cell infiltration, and reduced T cell papilloma infiltration. These results suggest that papilloma HPV gene expression may regulate interferon signaling and chemokine expression profiles within the tumor microenvironment that cooperate to govern clinical response to therapeutic HPV vaccination in patients with respiratory papillomatosis.

An Inflammatory Checkpoint Generated by IL1RN Splicing Offers Therapeutic Opportunity for KRAS-Mutant Intrahepatic Cholangiocarcinoma

KRAS mutations are causally linked to protumor inflammation and are identified as driving factors in tumorigenesis. Here, using multiomics data gathered from a large set of patients, we showed that KRAS mutation was associated with a specific landscape of alternative mRNA splicing that connected to myeloid inflammation in intrahepatic cholangiocarcinoma (iCCA). Then, we identified a negative feedback mechanism in which the upregulation of interleukin 1 receptor antagonist (IL1RN)-201/203 due to alternative splicing confers vital anti-inflammatory effects in KRAS-mutant iCCA. In KRAS-mutant iCCA mice, both IL1RN-201/203 upregulation and anakinra treatment ignited a significant antitumor immune response by altering neutrophil recruitment and phenotypes. Furthermore, anakinra treatment synergistically enhanced anti-PD-1 therapy to activate intratumoral GZMB+ CD8+ T cells in KRAS-mutant iCCA mice. Clinically, we found that high IL1RN-201/203 levels in patients with KRAS-mutant iCCA were significantly associated with superior response to anti-PD-1 immunotherapy.

SIGNIFICANCE

This work describes a novel inflammatory checkpoint mediated by IL1RN alternative splicing variants that may serve as a promising basis to develop therapeutic options for KRAS-mutant iCCA and other cancers. This article is featured in Selected Articles from This Issue, p. 2109.

Intratumoral CD103+ CD8+ T cells predict response to neoadjuvant chemoimmunotherapy in advanced head and neck squamous cell carcinoma

BACKGROUND

Immune cell heterogenicity is known to determine the therapeutic response to cancer progression. Neoadjuvant chemoimmunotherapy (NACI) has shown clinical benefits in some patients with advanced head and neck squamous cell carcinoma (HNSCC), but the underlying mechanism behind this clinical response is unknown. The efficacy of NACI needs to be potentiated by identifying accurate biomarkers to predict clinical responses. Here, we attempted to identify molecules predicting NACI response in advanced HNSCC.

METHODS

We performed combined single-cell RNA sequencing (scRNA-seq) and multiplex immunofluorescence (mIHC) staining with tumor samples derived from NACI-treated HNSCC patients to identify a new tumor-infiltrating cell (TIL) subtype, CD103+ CD8+ TILs, associated with clinical response, while both in vitro and in vivo assays were carried out to determine its antitumor efficiency. The regulatory mechanism of the CD103+ CD8+ TILs population was examined by performing cell-cell interaction analysis of the scRNA-seq data and spatial analysis of the mIHC images.

RESULTS

We established intratumoral CD103+ CD8+ TILs density as a determinant of NACI efficacy in cancers. Our scRNA-seq results indicated that the population of CD103+ CD8+ TILs was dramatically increased in the responders of NACI-treated HNSCC patients, while mIHC analysis confirmed the correlation between intratumoral CD103+ CD8+ TILs density and NACI efficacy in HNSCC patients. Further receiver operating characteristic curve analysis defined this TIL subset as a potent marker to predict patient response to NACI. Functional assays showed that CD103+ CD8+ TILs were tumor-reactive T cells, while programmed cell death protein-1 (PD-1) blockade enhanced CD103+ CD8+ TILs cytotoxicity against tumor growth in vivo. Mechanistically, targeting the triggering receptor expressed on myeloid cells 2-positive (TREM2+ ) macrophages might enhance the population of CD103+ CD8+ TILs and facilitate antitumor immunity during NACI treatment.

CONCLUSIONS

Our study highlights the impact of intratumoral CD103+ CD8+ TILs density on NACI efficacy in different cancers, while the efforts to elevate its population warrant further clinical investigation.

Decoding Human Biology and Disease Using Single-cell Omics Technologies

Over the past decade, advances in single-cell omics (SCO) technologies have enabled the investigation of cellular heterogeneity at an unprecedented resolution and scale, opening a new avenue for understanding human biology and disease. In this review, we summarize the developments of sequencing-based SCO technologies and computational methods, and focus on considerable insights acquired from SCO sequencing studies to understand normal and diseased properties, with a particular emphasis on cancer research. We also discuss the technological improvements of SCO and its possible contribution to fundamental research of the human, as well as its great potential in clinical diagnoses and personalized therapies of human disease.

Single-cell dissection of tumor microenvironmental response and resistance to cancer therapy

Cancer treatment strategies have evolved significantly over the years, with chemotherapy, targeted therapy, and immunotherapy as major pillars. Each modality leads to unique treatment outcomes by interacting with the tumor microenvironment (TME), which imposes a fundamental selective pressure on cancer progression. The advent of single-cell profiling technologies has revolutionized our understanding of the intricate and heterogeneous nature of the TME at an unprecedented resolution. This review delves into the commonalities and differential manifestations of how cancer therapies reshape the microenvironment in diverse cancer types. We highlight how groundbreaking immune checkpoint blockade (ICB) strategies alone or in combination with tumor-targeting treatments are endowed with comprehensive mechanistic insights when decoded at the single-cell level, aiming to drive forward future research directions on personalized treatments.

The technological landscape and applications of single-cell multi-omics

Single-cell multi-omics technologies and methods characterize cell states and activities by simultaneously integrating various single-modality omics methods that profile the transcriptome, genome, epigenome, epitranscriptome, proteome, metabolome and other (emerging) omics. Collectively, these methods are revolutionizing molecular cell biology research. In this comprehensive Review, we discuss established multi-omics technologies as well as cutting-edge and state-of-the-art methods in the field. We discuss how multi-omics technologies have been adapted and improved over the past decade using a framework characterized by optimization of throughput and resolution, modality integration, uniqueness and accuracy, and we also discuss multi-omics limitations. We highlight the impact that single-cell multi-omics technologies have had in cell lineage tracing, tissue-specific and cell-specific atlas production, tumour immunology and cancer genetics, and in mapping of cellular spatial information in fundamental and translational research. Finally, we discuss bioinformatics tools that have been developed to link different omics modalities and elucidate functionality through the use of better mathematical modelling and computational methods.

TCL1A+ B cells predict prognosis in triple-negative breast cancer through integrative analysis of single-cell and bulk transcriptomic data

Triple-negative breast cancer (TNBC) is an aggressive subtype with limited treatment options and high mortality rates. It remains a prevailing clinical need to distinguish whether the patient can benefit from therapy, such as chemotherapy. By integrating single-cell and global transcriptome data, we have for the first time identified TCL1A+ B cell functions that are prognostically relevant in TNBC. This finding broadens the perspective of traditional tumor-infiltrating lymphocytes in predicting survival, especially the potential value of B cells in TNBC. Single-cell RNA-seq data from five TNBC patients were collected to identify the association between immune cell populations and clinical outcomes. Functional analysis was according to gene set enrichment analysis using pathways from MsigDB. Subsequently, the gene signature of TCL1A+ B cells based on differential expression genes of TCL1A+ B cells versus other immune cells was used to explore the correlation with tumor microenvironment (TME) and construct a prognostic signature using a non-parametric and unsupervised method. We identified TCL1A+ B cells as a cluster of B cells associated with clinical outcomes in TNBC. Functional analysis demonstrated its function in B cell activation and regulation of immune response. The highly enriched TCL1A+ B cell population was found to be associated with a thermal TME with anti-tumor effects. A high abundance of TCL1A+ B cell population is positively correlated with a favorable therapeutic outcome, as indicated by longer overall survival. The present study suggests that TCL1A+ B cells play a key role in the treatment and prognostic prediction of TNBC, although further studies are needed to validate our findings. Moreover, the integration of transcriptome data at various resolutions provides a viable approach for the discovery of novel prognostic markers.

Comprehensive immunophenotyping of solid tumor-infiltrating immune cells reveals the expression characteristics of LAG-3 and its ligands

Background

Immune cell expression profiling from patient samples is critical for the successful development of immuno-oncology agents and is useful to understand mechanism-of-action, to identify exploratory biomarkers predictive of response, and to guide treatment selection and combination therapy strategies. LAG-3 is an inhibitory immune checkpoint that can suppress antitumor T-cell responses and targeting LAG-3, in combination with PD-1, is a rational approach to enhance antitumor immunity that has recently demonstrated clinical success. Here, we sought to identify human immune cell subsets that express LAG-3 and its ligands, to characterize the marker expression profile of these subsets, and to investigate the potential relationship between LAG-3 expressing subsets and clinical outcomes to immuno-oncology therapies.

Methods

Comprehensive high-parameter immunophenotyping was performed using mass and flow cytometry of tumor-infiltrating lymphocytes (TILs) and peripheral blood mononuclear cells (PBMCs) from two independent cohorts of samples from patients with various solid tumor types. Profiling of circulating immune cells by single cell RNA-seq was conducted on samples from a clinical trial cohort of melanoma patients treated with immunotherapy.

Results

LAG-3 was most highly expressed by subsets of tumor-infiltrating CD8 T central memory (TCM) and effector memory (TEM) cells and was frequently co-expressed with PD-1. We determined that these PD-1+ LAG-3+ CD8 memory T cells exhibited a unique marker profile, with greater expression of activation (CD69, HLA-DR), inhibitory (TIM-3, TIGIT, CTLA-4) and stimulatory (4-1BB, ICOS) markers compared to cells that expressed only PD-1 or LAG-3, or that were negative for both checkpoints. In contrast to tumors, LAG-3 expression was more limited in circulating immune cells from healthy donors and solid tumor patients. Additionally, we found abundant expression of the LAG-3 ligands MHC-II and galectin-3 in diverse immune cell types, whereas FGL1 and LSECtin were minimally expressed by immune cells in the tumor microenvironment (TME). Lastly, we found an inverse relationship between baseline and on-treatment levels of circulating LAG3 transcript-expressing CD8 memory T cells and response to combination PD-1 and CTLA-4 blockade in a clinical trial cohort of melanoma patients profiled by scRNAseq.

Conclusions

These results provide insights into the nature of LAG-3- and ligand-expressing immune cells within the TME, and suggest a biological basis for informing mechanistic hypotheses, treatment selection strategies, and combination immunotherapy approaches to support continued development of dual PD-1 and LAG-3 blockade.

Ligand-based adoptive T cell targeting CA125 in ovarian cancer

BACKGROUND

Ovarian cancer (OC) is a highly aggressive gynecological malignancy prevalent worldwide. Most OC cases are typically diagnosed at advanced stages, which has led to a 5-year overall survival rate of less than 35% following conventional treatment. Furthermore, immune checkpoint inhibitor therapy has shown limited efficacy in the treatment of patients with OC, and CAR-T therapy has also demonstrated modest results owing to inadequate T cell infiltration. Therefore, novel strategies must be developed to enhance T cell persistence and trafficking within the OC tumor microenvironment.

METHODS

In this study, we developed a novel adoptive T-cell therapy for ovarian cancer based on a chimeric antigen receptor structure. We used a ligand-receptor binding motif to enhance the therapeutic effect of targeting CA125. Since mesothelin can naturally bind to CA125 with high affinity, we concatenated the core-binding fragment of mesothelin with the 4-1BB and CD3ζ signal fragments to assemble a novel CA125-targeting chimeric receptor (CR). The CAR structure targeting CA125 derived from the 4H11 antibody was also constructed. CR- and CAR-encoding RNA were electroporated into T cells to evaluate their antitumor activity both in vitro and in vivo.

RESULTS

While CR-T or CAR-T cells exhibited moderate activity against two ovarian cancer cell lines, T cells co-expressing CR and CAR exhibited a superior killing effect compared to T cells expressing either CR or CAR alone. Furthermore, upon interaction with ovarian tumors, the ability of CR and CAR T cells to release activation markers and functional cytokines increased significantly. Similarly, CR and CAR co-expressing T cells persistently controlled the growth of transplanted ovarian cancer tumors in NSG mice and significantly prolonged the overall survival of tumor-challenged mice. Transcriptome sequencing revealed that the survival and cytotoxicity of T cells co-expressing CR and CAR were significantly altered compared with those of T cells expressing either CR or CAR.

CONCLUSION

Our findings demonstrate that CA125 targeting CR and CAR can synergistically kill ovarian cancer cells, indicating that CA125 targeting by the two binding motifs simultaneously in tumors may improve the therapeutic outcomes of ovarian cancer treatment.

Cellular and molecular waypoints along the path of T cell exhaustion

Thirty years of foundational research investigating molecular and cellular mechanisms promoting T cell exhaustion are now enabling rational design of T cell-based therapies for the treatment of chronic infections and cancer. Once described as a static cell fate, it is now well appreciated that the developmental path toward exhaustion is composed of a heterogeneous pool of cells with varying degrees of effector potential that ultimately converge on a terminally differentiated state. Recent description of the developmental stages along the differentiation trajectory of T cell exhaustion has provided insight into past immunotherapeutic success and future opportunities. Here, we discuss the hallmarks of distinct developmental stages occurring along the path to T cell dysfunction and the impact of these discrete CD8+ T cell fates on cancer immunotherapy.

Immune checkpoint blockade induces distinct alterations in the microenvironments of primary and metastatic brain tumors

In comparison with responses in recurrent glioblastoma (rGBM), the intracranial response of brain metastases (BrM) to immune checkpoint blockade (ICB) is less well studied. Here, we present an integrated single-cell RNA-Seq (scRNA-Seq) study of 19 ICB-naive and 9 ICB-treated BrM samples from our own and published data sets. We compared them with our previously published scRNA-Seq data from rGBM and found that ICB led to more prominent T cell infiltration into BrM than rGBM. These BrM-infiltrating T cells exhibited a tumor-specific phenotype and displayed greater activated/exhausted features. We also used multiplex immunofluorescence and spatial transcriptomics to reveal that ICB reduced a distinct CD206+ macrophage population in the perivascular space, which may modulate T cell entry into BrM. Furthermore, we identified a subset of progenitor exhausted T cells that correlated with longer overall survival in BrM patients. Our study provides a comprehensive immune cellular landscape of ICB's effect on metastatic brain tumors and offers insights into potential strategies for improving ICB efficacy for brain tumor patients.

The screening, identification, design and clinical application of tumor-specific neoantigens for TCR-T cells

Recent advances in neoantigen research have accelerated the development of tumor immunotherapies, including adoptive cell therapies (ACTs), cancer vaccines and antibody-based therapies, particularly for solid tumors. With the development of next-generation sequencing and bioinformatics technology, the rapid identification and prediction of tumor-specific antigens (TSAs) has become possible. Compared with tumor-associated antigens (TAAs), highly immunogenic TSAs provide new targets for personalized tumor immunotherapy and can be used as prospective indicators for predicting tumor patient survival, prognosis, and immune checkpoint blockade response. Here, the identification and characterization of neoantigens and the clinical application of neoantigen-based TCR-T immunotherapy strategies are summarized, and the current status, inherent challenges, and clinical translational potential of these strategies are discussed.

scGEM: Unveiling the Nested Tree-Structured Gene Co-Expressing Modules in Single Cell Transcriptome Data

BACKGROUND

Single-cell transcriptome analysis has fundamentally changed biological research by allowing higher-resolution computational analysis of individual cells and subsets of cell types. However, few methods have met the need to recognize and quantify the underlying cellular programs that determine the specialization and differentiation of the cell types.

METHODS

In this study, we present scGEM, a nested tree-structured nonparametric Bayesian model, to reveal the gene co-expression modules (GEMs) reflecting transcriptome processes in single cells.

RESULTS

We show that scGEM can discover shared and specialized transcriptome signals across different cell types using peripheral blood mononuclear single cells and early brain development single cells. scGEM outperformed other methods in perplexity and topic coherence (p < 0.001) on our simulation data. Larger datasets, deeper trees and pre-trained models are shown to be positively associated with better scGEM performance. The GEMs obtained from triple-negative breast cancer single cells exhibited better correlations with lymphocyte infiltration (p = 0.009) and the cell cycle (p < 0.001) than other methods in additional validation on the bulk RNAseq dataset.

CONCLUSIONS

Altogether, we demonstrate that scGEM can be used to model the hidden cellular functions of single cells, thereby unveiling the specialization and generalization of transcriptomic programs across different types of cells.

Prognostic significance and mechanisms of CXCL genes in clear cell renal cell carcinoma

This study aimed to investigate the clinical significance, biological functions, and underlying mechanisms of CXCL genes in clear cell renal cell carcinoma (ccRcc) based on patient datasets and pan-cancer analysis. The interaction between CXCL genes in ccRcc and immune components, particularly in relation to neutrophil recruitment and polarization mechanisms, was also evaluated. Furthermore, a risk score was developed using a signature for neutrophil polarization. The role of CXCL2 was assessed through in vitro experiments. Results showed that five CXCL genes (CXCL 2, 5, 9, 10, and 11) were upregulated in renal cancer tissue, while seven genes (CXCL 1, 2, 3, 5, 8, 13, and 14) significantly impacted patient survival. Moreover, CXCL 1, 5, and 13 affected progression-free survival. Besides, differences in mRNA expression and immune components affected renal cancer outcomes. Furthermore, three pairs of CXCL gene-immune cell interactions (CXCL13-CD8+ T cells, CXCL9/10-M1 cells, CXCL1/2/3/8-neutrophils) were identified through single-cell and pan-cancer analysis. A TAN risk score with prognostic value for KIRC patients was constructed using 11 genes and a TAN signature. Neutrophil polarization significantly impacted survival. Notably, CXCL2 was involved in neutrophil recruitment and polarization, thus promoting ccRcc progression. In conclusion, seven prognostic CXCL genes (CXCL 1/2/3/5/8/13/14) for ccRcc patients and three pairs of CXCL gene-immune cell interactions were identified. Furthermore, results showed that CXCL 2 promotes ccRcc progression through neutrophil recruitment and polarization.

Molecular subtyping and immune score system by a novel pyroptosis-based gene signature precisely predict immune infiltrating, survival and response to immune-checkpoint blockade in breast cancer

Breast cancer is one of the most prevalent and lethal types of cancer affecting women globally. Pyroptosis is a recently elucidated form of inflammatory cell death mediated by the gasdermin family and is considered to be associated with the tumor immune microenvironment. However, the impact of pyroptosis on breast cancer patients remains unclear. In this study, we identified 31 Pyroptosis-Related-Genes (PRGs) and investigated their association with patient survival using data from the TCGA database. We then established a gene signature comprising 6 PRGs that were significantly correlated with prognosis, and used these genes to classify breast cancer into 2 molecular subtypes. We investigated the characteristics of these two subtypes and found that our molecular subtyping accurately separated the samples into two groups with distinct immune infiltration and prognosis. Patients with higher expression of these genes had significantly greater immune infiltrating, T cell signaling, and better prognosis. Moreover, we developed an immune score system based on these 6 genes that accurately predicted the immune infiltrating of patients and their response to immune-checkpoint blockade, which was difficult to achieve with previous models. Additionally, through single-cell analyses, we found that patients with higher immune scores had stronger cytotoxic immune cells. In summary, our study identified a novel gene set and developed an immune scoring system based on this gene set that can precisely predict the immune microenvironment and responses to immunotherapy of breast cancer (BRCA) patients, which could be useful in clinical trials.

Synergism Between IL21 and Anti-PD-1 Combination Therapy is Underpinned by the Coordinated Reprogramming of the Immune Cellular Network in the Tumor Microenvironment

T cell-stimulating cytokines and immune checkpoint inhibitors (ICI) are an ideal combination for increasing response rates of cancer immunotherapy. However, the results of clinical trials have not been satisfying. It is important to understand the mechanism of synergy between these two therapeutic modalities. Here, through integrated analysis of multiple single-cell RNA sequencing (scRNA-seq) datasets of human tumor-infiltrating immune cells, we demonstrate that IL21 is produced by tumor-associated T follicular helper cells and hyperactivated/exhausted CXCL13+CD4+ T cells in the human tumor microenvironment (TME). In the mouse model, the hyperactivated/exhausted CD4+ T cell-derived IL21 enhances the helper function of CD4+ T cells that boost CD8+ T cell-mediated immune responses during PD-1 blockade immunotherapy. In addition, we demonstrated that IL21's antitumor activity did not require T-cell trafficking. Using scRNA-seq analysis of the whole tumor-infiltrating immune cells, we demonstrated that IL21 treatment in combination with anti-PD-1 blockade synergistically drives tumor antigen-specific CD8+ T cells to undergo clonal expansion and differentiate toward the hyperactive/exhausted functional state in the TME. In addition, IL21 treatment and anti-PD-1 blockade synergistically promote dendritic cell (DC) activation and maturation to mature DC as well as monocyte to type 1 macrophage (M1) differentiation in the TME. Furthermore, the combined treatment reprograms the immune cellular network by reshaping cell-cell communication in the TME. Our study establishes unique mechanisms of synergy between IL21 and PD-1-based ICI in the TME through the coordinated promotion of type 1 immune responses.

Significance

This study reveals how cytokine and checkpoint inhibitor therapy can be combined to increase the efficacy of cancer immunotherapy.

Integrated single cell and bulk sequencing analysis identifies tumor reactive CXCR6+ CD8 T cells as a predictor of immune infiltration and immunotherapy outcomes in hepatocellular carcinoma

Background

Various immune cell types in the tumor microenvironment (TME) of hepatocellular carcinoma (HCC) have been identified as important parameters associated with prognosis and responsiveness to immunotherapy. However, how various factors influence immune cell infiltration remains incompletely understood. Hence, we investigated the single cell multi-omics landscape of immune infiltration in HCC, particularly key gene and cell subsets that influence immune infiltration, thus potentially linking the immunotherapy response and immune cell infiltration.

Methods

We grouped patients with HCC according to immune cell infiltration scores calculated by single sample gene set enrichment analysis (ssGSEA). Differential expression analysis, functional enrichment, clinical trait association, gene mutation analysis, tumor immune dysfunction and exclusion (TIDE) and prognostic model construction were used to investigate the immune infiltration landscape through multi-omics. Stepwise regression was further used to identify key genes regulating immune infiltration. Single cell analysis was performed to explore expression patterns of candidate genes and investigate associated cellular populations. Correlation analysis, ROC analysis, Immunotherapy cohorts were used to explore and confirm the role of key gene and cellular population in predicting immune infiltration state and immunotherapy response. Immunohistochemistry and multiplexed fluorescence staining were used to further validated our results.

Results

Patients with HCC were clustered into high and low immune infiltration groups. Mutations of CTNNB1 and TTN were significantly associated with immune infiltration and altered enrichment of cell populations in the TME. TIDE analysis demonstrated that T cell dysfunction and the T cell exclusion score were elevated in the high and low infiltration groups, respectively. Six risk genes and five risk immune cell types were identified and used to construct risk scores and a nomogram model. CXCR6 and LTA, identified by stepwise regression, were highly associated with immune infiltration. Single cell analysis revealed that LTA was expressed primarily in tumor infiltrating T lymphocytes and partial B lymphocytes, whereas CXCR6 was enriched predominantly in T and NK cells. Notably, CXCR6+ CD8 T cells were characterized as tumor enriched cells that may be potential predictors of high immune infiltration and the immune-checkpoint blockade response, and may serve as therapeutic targets.

Conclusion

We constructed a comprehensive single cell and multi-omics landscape of immune infiltration in HCC, and delineated key genes and cellular populations regulating immune infiltration and immunotherapy response, thus providing insights into the mechanisms of immune infiltration and future therapeutic control.

Location of CD39+ T cell subpopulations within tumors predict differential outcomes in non-small cell lung cancer

PURPOSE

An improved mechanistic understanding of immunosuppressive pathways in non-small cell lung cancer (NSCLC) is important to develop novel diagnostic and therapeutic approaches. Here, we investigate the prognostic significance of the ectonucleotidases CD39 and CD73 in NSCLC.

EXPERIMENTAL DESIGN

The expression and localization of CD39, CD73 and CD103 was digitally quantified in a cohort of 162 early treatment naïve NSCLC patients using multiplex-immunofluorescence and related to patient outcome. Expression among different cell-populations was assessed via flow cytometry. Targeted RNA-Seq was performed on CD4+ and CD8+ T cells from digested NSCLC tumor tissue and single-cell RNA-Seq data was analyzed to investigate the functional significance of CD39+ T cell populations.

RESULTS

We demonstrate that flow cytometry of early untreated NSCLC patients shows an upregulation of CD39 expression in the tumor tissue among natural killer (NK) cells, fibroblasts and T cells. CD73 expression is mainly found among fibroblasts and Epcam+cells in the tumor tissue. Multiplex Immunofluorescence in a cohort of 162 early untreated NSCLC patients demonstrates that CD39 expression is mainly localized in the tumor stroma while CD73 expression is equally distributed between tumor nest and stroma, and high expression of CD39 and CD73 in the tumor stroma is associated with poor recurrence-free survival (RFS) at 5 years. Additionally, we find that CD8+T cells located in the tumor nest express CD103 and the density of CD39+CD103+CD8+ T cells in the tumor nest predicts improved RFS at 5 years. Targeted RNA-Seq shows that the tumor microenvironment of NSCLC upregulates regulatory pathways in CD4+ T cells and exhaustion in CD8+ T cells, and analysis of a single cell RNA sequencing dataset shows that CD39+CD4+ cells are enriched in Treg signature gene-sets, and CD39+CD103+ cytotoxic T lymphocyte show gene signatures indicative of an exhausted cytotoxic phenotype with upregulated expression of CXCL13.

CONCLUSIONS

Knowledge of patterns of distribution and location are required to understand the prognostic impact of CD39+ T cell populations in NSCLC. This study provides an improved understanding of spatial and functional characteristics of CD39+ T cells and their significance to patient outcome.

Cancer-associated fibroblasts drive CXCL13 production in activated T cells via TGF-beta

Introduction

Tumour-reactive T cells producing the B-cell attractant chemokine CXCL13, in solid tumours, promote development of tertiary lymphoid structures (TLS) and are associated with improved prognosis and responsiveness to checkpoint immunotherapy. Cancer associated fibroblasts are the dominant stromal cell type in non-small cell lung cancer (NSCLC) where they co-localise with T cells and can influence T cell activation and exhaustion. We questioned whether CAF directly promote CXCL13-production during T cell activation.

Methods

We characterised surface markers, cytokine production and transcription factor expression in CXCL13-producing T cells in NSCLC tumours and paired non-cancerous lung samples using flow cytometry. We then assessed the influence of human NSCLC-derived primary CAF lines on T cells from healthy donors and NSCLC patients during activation in vitro measuring CXCL13 production and expression of cell-surface markers and transcription factors by flow cytometry.

Results

CAFs significantly increased the production of CXCL13 by both CD4+ and CD8+ T cells. CAF-induced CXCL13-producing cells lacked expression of CXCR5 and BCL6 and displayed a T peripheral helper cell phenotype. Furthermore, we demonstrate CXCL13 production by T cells is induced by TGF-β and limited by IL-2. CAF provide TGF-β during T cell activation and reduce availability of IL-2 both directly (by reducing the capacity for IL-2 production) and indirectly, by expanding a population of activated Treg. Inhibition of TGF-β signalling prevented both CAF-driven upregulation of CXCL13 and Treg expansion.

Discussion

Promoting CXCL13 production represents a newly described immune-regulatory function of CAF with the potential to shape the immune infiltrate of the tumour microenvironment both by altering the effector-function of tumour infiltrating T-cells and their capacity to attract B cells and promote TLS formation.

Defining and using immune archetypes to classify and treat cancer

Tumours are surrounded by a host immune system that can suppress or promote tumour growth. The tumour microenvironment (TME) has often been framed as a singular entity, suggesting a single type of immune state that is defective and in need of therapeutic intervention. By contrast, the past few years have highlighted a plurality of immune states that can surround tumours. In this Perspective, we suggest that different TMEs have 'archetypal' qualities across all cancers - characteristic and repeating collections of cells and gene-expression profiles at the level of the bulk tumour. We discuss many studies that together support a view that tumours typically draw from a finite number (around 12) of 'dominant' immune archetypes. In considering the likely evolutionary origin and roles of these archetypes, their associated TMEs can be predicted to have specific vulnerabilities that can be leveraged as targets for cancer treatment with expected and addressable adverse effects for patients.

Pan-cancer T cell atlas links a cellular stress response state to immunotherapy resistance

Tumor-infiltrating T cells offer a promising avenue for cancer treatment, yet their states remain to be fully characterized. Here we present a single-cell atlas of T cells from 308,048 transcriptomes across 16 cancer types, uncovering previously undescribed T cell states and heterogeneous subpopulations of follicular helper, regulatory and proliferative T cells. We identified a unique stress response state, TSTR, characterized by heat shock gene expression. TSTR cells are detectable in situ in the tumor microenvironment across various cancer types, mostly within lymphocyte aggregates or potential tertiary lymphoid structures in tumor beds or surrounding tumor edges. T cell states/compositions correlated with genomic, pathological and clinical features in 375 patients from 23 cohorts, including 171 patients who received immune checkpoint blockade therapy. We also found significantly upregulated heat shock gene expression in intratumoral CD4/CD8+ cells following immune checkpoint blockade treatment, particularly in nonresponsive tumors, suggesting a potential role of TSTR cells in immunotherapy resistance. Our well-annotated T cell reference maps, web portal and automatic alignment/annotation tool could provide valuable resources for T cell therapy optimization and biomarker discovery.

Intratumoral dendritic cell-CD4+ T helper cell niches enable CD8+ T cell differentiation following PD-1 blockade in hepatocellular carcinoma

Despite no apparent defects in T cell priming and recruitment to tumors, a large subset of T cell rich tumors fail to respond to immune checkpoint blockade (ICB). We leveraged a neoadjuvant anti-PD-1 trial in patients with hepatocellular carcinoma (HCC), as well as additional samples collected from patients treated off-label, to explore correlates of response to ICB within T cell-rich tumors. We show that ICB response correlated with the clonal expansion of intratumoral CXCL13+CH25H+IL-21+PD-1+CD4+ T helper cells ("CXCL13+ TH") and Granzyme K+ PD-1+ effector-like CD8+ T cells, whereas terminally exhausted CD39hiTOXhiPD-1hiCD8+ T cells dominated in nonresponders. CD4+ and CD8+ T cell clones that expanded post-treatment were found in pretreatment biopsies. Notably, PD-1+TCF-1+ (Progenitor-exhausted) CD8+ T cells shared clones mainly with effector-like cells in responders or terminally exhausted cells in nonresponders, suggesting that local CD8+ T cell differentiation occurs upon ICB. We found that these Progenitor CD8+ T cells interact with CXCL13+ TH within cellular triads around dendritic cells enriched in maturation and regulatory molecules, or "mregDC". These results suggest that discrete intratumoral niches that include mregDC and CXCL13+ TH control the differentiation of tumor-specific Progenitor exhasuted CD8+ T cells following ICB.

CD39 (ENTPD1) in tumors: a potential therapeutic target and prognostic biomarker

As a regulator of the dynamic balance between immune-activated extracellular ATP and immunosuppressive adenosine, CD39 ectonucleotidase impairs the ability of immune cells to exert anticancer immunity and plays an important role in the immune escape of tumor cells within the tumor microenvironment. In addition, CD39 has been studied in cancer patients to evaluate the prognosis, the efficacy of immunotherapy (e.g., PD-1 blockade) and the prediction of recurrence. This article reviews the importance of CD39 in tumor immunology, summarizes the preclinical evidence on targeting CD39 to treat tumors and focuses on the potential of CD39 as a biomarker to evaluate the prognosis and the response to immune checkpoint inhibitors in tumors.

Phenotypic diversity of T cells in human primary and metastatic brain tumors revealed by multiomic interrogation

The immune-specialized environment of the healthy brain is tightly regulated to prevent excessive neuroinflammation. However, after cancer development, a tissue-specific conflict between brain-preserving immune suppression and tumor-directed immune activation may ensue. To interrogate potential roles of T cells in this process, we profiled these cells from individuals with primary or metastatic brain cancers via integrated analyses on the single-cell and bulk population levels. Our analysis revealed similarities and differences in T cell biology between individuals, with the most pronounced differences observed in a subgroup of individuals with brain metastasis, characterized by accumulation of CXCL13-expressing CD39+ potentially tumor-reactive T (pTRT) cells. In this subgroup, high pTRT cell abundance was comparable to that in primary lung cancer, whereas all other brain tumors had low levels, similar to primary breast cancer. These findings indicate that T cell-mediated tumor reactivity can occur in certain brain metastases and may inform stratification for treatment with immunotherapy.

Transitioning single-cell genomics into the clinic

The use of genomics is firmly established in clinical practice, resulting in innovations across a wide range of disciplines such as genetic screening, rare disease diagnosis and molecularly guided therapy choice. This new field of genomic medicine has led to improvements in patient outcomes. However, most clinical applications of genomics rely on information generated from bulk approaches, which do not directly capture the genomic variation that underlies cellular heterogeneity. With the advent of single-cell technologies, research is rapidly uncovering how genomic data at cellular resolution can be used to understand disease pathology and mechanisms. Both DNA-based and RNA-based single-cell technologies have the potential to improve existing clinical applications and open new application spaces for genomics in clinical practice, with oncology, immunology and haematology poised for initial adoption. However, challenges in translating cellular genomics from research to a clinical setting must first be overcome.

PD-1 blockade potentiates neoadjuvant chemotherapy in NSCLC via increasing CD127+ and KLRG1+ CD8 T cells

The combination of PD-1 blockade with neoadjuvant chemotherapy (NAC) has achieved unprecedented clinical success in non-small cell lung cancer (NSCLC) compared to NAC alone, but the underlying mechanisms by which PD-1 blockade augments the effects of chemotherapy remain incompletely elucidated. Single-cell RNA sequencing was performed on CD45⁺ immune cells isolated from surgically resected fresh tumors of seven NSCLC patients receiving NAC or neoadjuvant pembrolizumab and chemotherapy (NAPC). Multiplex fluorescent immunohistochemistry was performed on FFPE tissues before and after NAC or NAPC from 65 resectable NSCLC patients, and results were validated with GEO dataset. NAC resulted in an increase only of CD20⁺ B cells, whereas NAPC increased the infiltration of CD20⁺ B cells, CD4⁺ T cells, CD4⁺CD127⁺ T cells, CD8⁺ T cells, CD8⁺CD127⁺ and CD8⁺KLRG1⁺ T cells. Synergistic increase in B and T cells promotes favorable therapeutic response after NAPC. Spatial distribution analysis discovered that CD8⁺ T cells and their CD127⁺ and KLRG1⁺ subsets were in closer proximity to CD4⁺ T/CD20⁺ B cells in NAPC versus NAC. GEO dataset validated that B-cell, CD4, memory, and effector CD8 signatures correlated with therapeutic responses and clinical outcomes. The addition of PD-1 blockade to NAC promoted anti-tumor immunity through T and B cells recruitment in the tumor microenvironment and induced tumor-infiltrating CD8⁺ T cells skewed toward CD127⁺ and KLRG1⁺ phenotypes, which may be assisted by CD4⁺ T cells and B cells. Our comprehensive study identified key immune cell subsets exerting anti-tumor responses during PD-1 blockade therapy and that may be therapeutically targeted to improve upon existing immunotherapies for NSCLC.

Single-cell sequencing technology in colorectal cancer: a new technology to disclose the tumor heterogeneity and target precise treatment

Colorectal Cancer (CRC) is one of the most common gastrointestinal tumors, and its high tumor heterogeneity makes traditional sequencing methods incapable of obtaining information about the heterogeneity of individual cancer cells in CRC. Therefore, single-cell sequencing technology can be applied to better analyze the differences in genetic and protein information between cells, to obtain genomic sequence information of single cells, and to more thoroughly analyze the cellular characteristics and interactions in the CRC microenvironment. This will provide a more comprehensive understanding of colorectal cancer development and metastasis and indicate the treatment plan and prognosis. In this study, we review the application of single-cell sequencing to analyze the tumor microenvironment of CRC, explore the mechanisms involved in CRC metastasis and progression, and provide a reference for potential treatment options.

Phenotypic plasticity and reduced tissue retention of exhausted tumor-infiltrating T cells following neoadjuvant immunotherapy in head and neck cancer

Neoadjuvant immunotherapies (NITs) have led to clinical benefits in several cancers. Characterization of the molecular mechanisms underlying responses to NIT may lead to improved treatment strategies. Here we show that exhausted, tumor-infiltrating CD8+ T (Tex) cells display local and systemic responses to concurrent neoadjuvant TGF-β and PD-L1 blockade. NIT induces a significant and selective increase in circulating Tex cells associated with reduced intratumoral expression of the tissue-retention marker CD103. TGF-β-driven CD103 expression on CD8+ T cells is reversed following TGF-β neutralization in vitro, implicating TGF-β in T cell tissue retention and impaired systemic immunity. Transcriptional changes implicate T cell receptor signaling and glutamine metabolism as important determinants of enhanced or reduced Tex treatment response, respectively. Our analysis illustrates physiological and metabolic changes underlying T cell responses to NIT, highlighting the interplay between immunosuppression, tissue retention, and systemic anti-tumor immunity and suggest antagonism of T cell tissue retention as a promising neoadjuvant treatment strategy.

Dynamics and specificities of T cells in cancer immunotherapy

Recent advances in cancer immunotherapy - ranging from immune-checkpoint blockade therapy to adoptive cellular therapy and vaccines - have revolutionized cancer treatment paradigms, yet the variability in clinical responses to these agents has motivated intense interest in understanding how the T cell landscape evolves with respect to response to immune intervention. Over the past decade, the advent of multidimensional single-cell technologies has provided the unprecedented ability to dissect the constellation of cell states of lymphocytes within a tumour microenvironment. In particular, the rapidly expanding capacity to definitively link intratumoural phenotypes with the antigen specificity of T cells provided by T cell receptors (TCRs) has now made it possible to focus on investigating the properties of T cells with tumour-specific reactivity. Moreover, the assessment of TCR clonality has enabled a molecular approach to track the trajectories, clonal dynamics and phenotypic changes of antitumour T cells over the course of immunotherapeutic intervention. Here, we review the current knowledge on the cellular states and antigen specificities of antitumour T cells and examine how fine characterization of T cell dynamics in patients has provided meaningful insights into the mechanisms underlying effective cancer immunotherapy. We highlight those T cell subsets associated with productive T cell responses and discuss how diverse immunotherapies might leverage the pre-existing tumour-reactive T cell pool or instruct de novo generation of antitumour specificities. Future studies aimed at elucidating the factors associated with the elicitation of productive antitumour T cell immunity are anticipated to instruct the design of more efficacious treatment strategies.

Accelerating the understanding of cancer biology through the lens of genomics

A core mission of cancer genomics is to comprehensively chart molecular underpinnings of cancer-driving events and to provide personalized therapeutic strategies. Primarily focused on cancer cells, cancer genomics studies have successfully uncovered many drivers for major cancer types. Since the emergence of cancer immune evasion as a critical cancer hallmark, the paradigm has been elevated to the holistic tumor ecosystem, with distinct cellular components and their functional states elucidated. We highlight the milestones of cancer genomics, depict the evolving path of the field, and discuss future directions in completing the understanding of the tumor ecosystem and in advancing therapeutic strategies.

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.

GITR and TIGIT immunotherapy provokes divergent multi-cellular responses in the tumor microenvironment of gastrointestinal cancers

Understanding the cellular mechanisms of novel immunotherapy agents in the human tumor microenvironment (TME) is critical to their clinical success. We examined GITR and TIGIT immunotherapy in gastric and colon cancer patients using ex vivo slice tumor slice cultures derived from cancer surgical resections. This primary culture system maintains the original TME in a near-native state. We applied paired single-cell RNA and TCR sequencing to identify cell type specific transcriptional reprogramming. The GITR agonist was limited to increasing effector gene expression only in cytotoxic CD8 T cells. The TIGIT antagonist increased TCR signaling and activated both cytotoxic and dysfunctional CD8 T cells, including clonotypes indicative of potential tumor antigen reactivity. The TIGIT antagonist also activated T follicular helper-like cells and dendritic cells, and reduced markers of immunosuppression in regulatory T cells. Overall, we identified cellular mechanisms of action of these two immunotherapy targets in the patients' TME.

CD83 expression characterizes precursor exhausted T cell population

T cell exhaustion is a main obstacle against effective cancer immunotherapy. Exhausted T cells include a subpopulation that maintains proliferative capacity, referred to as precursor exhausted T cells (T_PEX). While functionally distinct and important for antitumor immunity, T_PEX possess some overlapping phenotypic features with the other T-cell subsets within the heterogeneous tumor-infiltrating T-lymphocytes (TIL). Here we explore surface marker profiles unique to T_PEX using the tumor models treated by chimeric antigen receptor (CAR)-engineered T cells. We find that CD83 is predominantly expressed in the CCR7⁺PD1⁺ intratumoral CAR-T cells compared with the CCR7^-PD1⁺ (terminally differentiated) and CAR-negative (bystander) T cells. The CD83⁺CCR7⁺ CAR-T cells exhibit superior antigen-induced proliferation and IL-2 production compared with the CD83^- T cells. Moreover, we confirm selective expression of CD83 in the CCR7⁺PD1⁺ T-cell population in primary TIL samples. Our findings identify CD83 as a marker to discriminate T_PEX from terminally exhausted and bystander TIL.

Advancing beyond the twists and turns of T cell exhaustion in cancer

Chronic antigen stimulation leads to T cell exhaustion. Nutrient restrictions and other suppressive factors in the tumor microenvironment further exacerbate T cell dysfunction. Better understanding of heterogeneity and dynamics of exhausted CD8 T cells will guide novel therapies that modulate T cell differentiation to achieve more effective antitumor responses.