Parallel single-cell and bulk transcriptome analyses reveal key features of the gastric tumor microenvironment.

BACKGROUND: The tumor microenvironment (TME) has been shown to strongly influence treatment outcome for cancer patients in various indications and to influence the overall survival. However, the cells forming the TME in gastric cancer have not been extensively characterized. RESULTS: We combine bulk and single-cell RNA sequencing from tumors and matched normal tissue of 24 treatment-naïve GC patients to better understand which cell types and transcriptional programs are associated with malignant transformation of the stomach. Clustering 96,623 cells of non-epithelial origin reveals 81 well-defined TME cell types. We find that activated fibroblasts and endothelial cells are most prominently overrepresented in tumors. Intercellular network reconstruction and survival analysis of an independent cohort imply the importance of these cell types together with immunosuppressive myeloid cell subsets and regulatory T cells in establishing an immunosuppressive microenvironment that correlates with worsened prognosis and lack of response in anti-PD1-treated patients. In contrast, we find a subset of IFNγ activated T cells and HLA-II expressing macrophages that are linked to treatment response and increased overall survival. CONCLUSIONS: Our gastric cancer single-cell TME compendium together with the matched bulk transcriptome data provides a unique resource for the identification of new potential biomarkers for patient stratification. This study helps further to elucidate the mechanism of gastric cancer and provides insights for therapy.

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

Immune-checkpoint blockade (ICB) therapies represent a paradigm shift in the treatment of human cancers; however, it remains incompletely understood how tumor-reactive T cells respond to ICB across tumor types. Here, we demonstrate that measuring CXCL13 expression could effectively identify both precursor and terminally differentiated tumor-reactive CD8+ T cells within tumors. Applying this approach, we performed meta-analyses of published single-cell data for CXCL13+CD8+ T cells in 225 samples from 102 patients treated with ICB across five cancer types. We found that CXCL13+CD8+ T cells were correlated with favorable responses to ICB, and the treatment further increased such cells in responsive tumors. In addition, CXCL13+ tumor-reactive subsets exhibited variable responses to ICB in distinct contexts, likely due to different degrees of exhaustion-related immunosuppression. Our integrated analyses provide insights into mechanisms underlying ICB and suggest that bolstering precursor tumor-reactive CD8+ T cells might provide an effective therapeutic approach to improve cancer treatment.

Single-cell landscape reveals active cell subtypes and their interaction in the tumor microenvironment of gastric cancer.

Background: Gastric cancer remains the third most common cause of cancer-related death worldwide. The development of novel therapeutic strategies for gastric cancer requires a deep understanding of the tumor cells and microenvironment of gastric cancer. Methods: We performed the single-cell RNA sequencing (scRNA-seq) on nine untreated non-metastatic gastric cancer patients. The transcriptomic atlas and ligand-receptor-based intercellular communication networks of the single cells were characterized. Results: Here, we profiled the transcriptomes of 47,304 cells from nine patients with gastric cancer. Tregs cells were significantly enriched in the gastric tumor tissues with increased expression of immune suppression related genes, which suggest a more immunosuppressive microenvironment. We also observed the absence of separate exhausted CD8+ T cell cluster, and the low expression level of exhaustion markers PDCD1, CTLA4, HAVCR2, LAG-3, and TIGIT in those specific cells. These may serve as molecular-level evidence for the limited benefit of immunotherapy among gastric cancer patients. In addition, we found ACKR1 specifically expressed in tumor endothelial cells, associated with poor prognosis in the cohort data and potentially provided a novel target of gastric cancer treatment. Furthermore, the tight interaction between endothelial cells and fibroblast implied the important roles of fibroblast in tumor angiogenesis and the maintenance of tumor vasculature. Conclusions: In conclusion, this single-cell atlas provide understanding the cellular heterogeneity from molecular level in gastric cancer and will serve as a valuable resource for developing innovative early and companion diagnostics, as well as discovering novel targeted therapies for gastric cancer.

Immune phenotypic linkage between colorectal cancer and liver metastasis.

The tumor microenvironment (TME) is connected to immunotherapy responses, but it remains unclear how cancer cells and host tissues differentially influence the immune composition within TME. Here, we performed single-cell analyses for autologous samples from liver metastasized colorectal cancer to disentangle factors shaping TME. By aligning CD45+ cells across different tissues, we classified exhausted CD8+ T cells (Texs) and activated regulatory T cells as M-type, whose phenotypes were associated with the malignancy, while natural killer and mucosal-associated invariant T cells were defined as N-type, whose phenotypes were associated with the niche. T cell receptor sharing between Texs in primary and metastatic tumors implicated the presence of common peripheral non-exhausted precursors. For myeloid cells, a subset of dendritic cells (DC3s) and SPP1+ macrophages were M-type, and the latter were predominant in liver metastasis, indicating its pro-metastasis role. Our analyses bridge immune phenotypes of primary and metastatic tumors, thereby helping to understand the tumor-specific contexture and identify the pro-metastasis components.

Temporal single-cell tracing reveals clonal revival and expansion of precursor exhausted T cells during anti-PD-1 therapy in lung cancer.

Anti-PD-1 treatment has shown unprecedented clinical success in the treatment of non-small-cell lung cancer (NSCLC), but the underlying mechanisms remain incompletely understood. Here, we performed temporal single-cell RNA and paired T-cell receptor sequencing on 47 tumor biopsies from 36 patients with NSCLC following PD-1-based therapies. We observed increased levels of precursor exhausted T (Texp) cells in responsive tumors after treatment, characterized by low expression of coinhibitory molecules and high expression of GZMK. By contrast, nonresponsive tumors failed to accumulate Texp cells. Our data suggested that Texp cells were unlikely to be derived from the reinvigoration of terminally exhausted cells; instead, they were accumulated by (1) local expansion and (2) replenishment by peripheral T cells with both new and pre-existing clonotypes, a phenomenon we named clonal revival. Our study provides insights into mechanisms underlying PD-1-based therapies, implicating clonal revival and expansion of Texp cells as steps to improve NSCLC treatment.

Pan-cancer single-cell landscape of tumor-infiltrating T cells.

T cells play a central role in cancer immunotherapy, but we lack systematic comparison of the heterogeneity and dynamics of tumor-infiltrating T cells across cancer types. We built a single-cell RNA-sequencing pan-cancer atlas of T cells for 316 donors across 21 cancer types and revealed distinct T cell composition patterns. We found multiple state-transition paths in the exhaustion of CD8+ T cells and the preference of those paths among different tumor types. Certain T cell populations showed specific correlation with patient properties such as mutation burden, shedding light on the possible determinants of the tumor microenvironment. T cell compositions within tumors alone could classify cancer patients into groups with clinical trait specificity, providing new insights into T cell immunity and precision immunotherapy targeting T cells.

Single-cell analyses reveal key immune cell subsets associated with response to PD-L1 blockade in triple-negative breast cancer.

In triple-negative breast cancer (TNBC), the benefit of combining chemotherapy with checkpoint inhibitors is still not very clear. We utilize single-cell RNA- and ATAC-sequencing to examine the immune cell dynamics in 22 patients with advanced TNBC treated with paclitaxel or its combination with the anti-PD-L1 atezolizumab. We demonstrate that high levels of baseline CXCL13+ T cells are linked to the proinflammatory features of macrophages and can predict effective responses to the combination therapy. In responsive patients, lymphoid tissue inducer (LTi) cells, follicular B (Bfoc) cells, CXCL13+ T cells, and conventional type 1 dendritic cells (cDC1) concertedly increase following the combination therapy, but instead decrease after paclitaxel monotherapy. Our data highlight the importance of CXCL13+ T cells in effective responses to anti-PD-L1 therapies and suggest that their reduction by paclitaxel regimen may compromise the clinical outcomes of accompanying atezolizumab for TNBC treatment.

Tumour heterogeneity and intercellular networks of nasopharyngeal carcinoma at single cell resolution.

The heterogeneous nature of tumour microenvironment (TME) underlying diverse treatment responses remains unclear in nasopharyngeal carcinoma (NPC). Here, we profile 176,447 cells from 10 NPC tumour-blood pairs, using single-cell transcriptome coupled with T cell receptor sequencing. Our analyses reveal 53 cell subtypes, including tumour-infiltrating CD8+ T, regulatory T (Treg), and dendritic cells (DCs), as well as malignant cells with different Epstein-Barr virus infection status. Trajectory analyses reveal exhausted CD8+ T and immune-suppressive TNFRSF4+ Treg cells in tumours might derive from peripheral CX3CR1+CD8+ T and naïve Treg cells, respectively. Moreover, we identify immune-regulatory and tolerogenic LAMP3+ DCs. Noteworthily, we observe intensive inter-cell interactions among LAMP3+ DCs, Treg, exhausted CD8+ T, and malignant cells, suggesting potential cross-talks to foster an immune-suppressive niche for the TME. Collectively, our study uncovers the heterogeneity and interacting molecules of the TME in NPC at single-cell resolution, which provide insights into the mechanisms underlying NPC progression and the development of precise therapies for NPC.

A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cells.

Tumor-infiltrating myeloid cells (TIMs) are key regulators in tumor progression, but the similarity and distinction of their fundamental properties across different tumors remain elusive. Here, by performing a pan-cancer analysis of single myeloid cells from 210 patients across 15 human cancer types, we identified distinct features of TIMs across cancer types. Mast cells in nasopharyngeal cancer were found to be associated with better prognosis and exhibited an anti-tumor phenotype with a high ratio of TNF+/VEGFA+ cells. Systematic comparison between cDC1- and cDC2-derived LAMP3+ cDCs revealed their differences in transcription factors and external stimulus. Additionally, pro-angiogenic tumor-associated macrophages (TAMs) were characterized with diverse markers across different cancer types, and the composition of TIMs appeared to be associated with certain features of somatic mutations and gene expressions. Our results provide a systematic view of the highly heterogeneous TIMs and suggest future avenues for rational, targeted immunotherapies.

Dynamics of peripheral T cell clones during PD-1 blockade in non-small cell lung cancer.

Understanding of the functional states and clonal dynamics of T cells after immune checkpoint blockade (ICB) is valuable for improving these therapeutic strategies. Here we performed Smart-seq2 single-cell RNA sequencing (scRNA-seq) analysis on 3,110 peripheral T cells of non-small cell lung cancer (NSCLC) patients before and after the initiation of programmed cell death protein 1 (PD-1) blockade. We identified individual peripheral T cell clones based on the full-length T cell receptor (TCR) sequences and monitored their dynamics during immunotherapy. We found a higher cytotoxic activity in the tumor-related CD4+ T cell clones than in the CD8+ T cell clones. Based on a large tumor-related CD4+ T cell clone, we observed a dramatically decreased abundance after progression, as well as a reduction in the percentage of PD-1+ T cells. We also detected 25 genes, such as CXCR4, DUSP2 and ZFP36, that were noticeably upregulated or downregulated following progression. In addition, the pseudotime trajectory of CD8+ T cell clones corresponded to the treatment time points, showing a decreased activity in the "cytokine and cytokine receptor interaction" pathway. These analyses provided an insight into the dynamics of peripheral T cell clones during PD-1 blockade in NSCLC.

Single-Cell Analyses Inform Mechanisms of Myeloid-Targeted Therapies in Colon Cancer.

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for defining cellular diversity in tumors, but its application toward dissecting mechanisms underlying immune-modulating therapies is scarce. We performed scRNA-seq analyses on immune and stromal populations from colorectal cancer patients, identifying specific macrophage and conventional dendritic cell (cDC) subsets as key mediators of cellular cross-talk in the tumor microenvironment. Defining comparable myeloid populations in mouse tumors enabled characterization of their response to myeloid-targeted immunotherapy. Treatment with anti-CSF1R preferentially depleted macrophages with an inflammatory signature but spared macrophage populations that in mouse and human expresses pro-angiogenic/tumorigenic genes. Treatment with a CD40 agonist antibody preferentially activated a cDC population and increased Bhlhe40+ Th1-like cells and CD8+ memory T cells. Our comprehensive analysis of key myeloid subsets in human and mouse identifies critical cellular interactions regulating tumor immunity and defines mechanisms underlying myeloid-targeted immunotherapies currently undergoing clinical testing.

Landscape and Dynamics of Single Immune Cells in Hepatocellular Carcinoma.

The immune microenvironment of hepatocellular carcinoma (HCC) is poorly characterized. Combining two single-cell RNA sequencing technologies, we produced transcriptomes of CD45+ immune cells for HCC patients from five immune-relevant sites: tumor, adjacent liver, hepatic lymph node (LN), blood, and ascites. A cluster of LAMP3+ dendritic cells (DCs) appeared to be the mature form of conventional DCs and possessed the potential to migrate from tumors to LNs. LAMP3+ DCs also expressed diverse immune-relevant ligands and exhibited potential to regulate multiple subtypes of lymphocytes. Of the macrophages in tumors that exhibited distinct transcriptional states, tumor-associated macrophages (TAMs) were associated with poor prognosis, and we established the inflammatory role of SLC40A1 and GPNMB in these cells. Further, myeloid and lymphoid cells in ascites were predominantly linked to tumor and blood origins, respectively. The dynamic properties of diverse CD45+ cell types revealed by this study add new dimensions to the immune landscape of HCC.

Deep single-cell RNA sequencing data of individual T cells from treatment-naïve colorectal cancer patients.

T cells, as a crucial compartment of the tumour microenvironment, play vital roles in cancer immunotherapy. However, the basic properties of tumour-infiltrating T cells (TILs) such as the functional state, migratory capability and clonal expansion remain elusive. Here, using Smart-seq2 protocol, we have generated a RNA sequencing dataset of 11,138 T cells isolated from peripheral blood, adjacent normal and tumour tissues of 12 colorectal cancer (CRC) patients, including 4 with microsatellite instability (MSI). The dataset contained an expression profile of 10,805 T cells, as well as the full-length T cell receptor (TCR) sequences of 9,878 cells after quality control. To facilitate data mining of our T cell dataset, we developed a web-based application to deliver systematic interrogations and customizable functionalities ( http://crctcell.cancer-pku.cn/ ). Functioning with our dataset, the web tool enables the characterization of TILs based on both transcriptome and assembled TCR sequences at the single cell level, which will help unleash the potential value of our CRC T cell data resource.

Distinct epigenetic features of tumor-reactive CD8+ T cells in colorectal cancer patients revealed by genome-wide DNA methylation analysis.

BACKGROUND: Tumor-reactive CD8+ tumor-infiltrating lymphocytes (TILs) represent a subtype of T cells that can recognize and destroy tumor specifically. Understanding the regulatory mechanism of tumor-reactive CD8+ T cells has important therapeutic implications. Yet the DNA methylation status of this T cell subtype has not been elucidated. RESULTS: In this study, we segregate tumor-reactive and bystander CD8+ TILs, as well as naïve and effector memory CD8+ T cell subtypes as controls from colorectal cancer patients, to compare their transcriptome and methylome characteristics. Transcriptome profiling confirms previous conclusions that tumor-reactive TILs have an exhausted tissue-resident memory signature. Whole-genome methylation profiling identifies a distinct methylome pattern of tumor-reactive CD8+ T cells, with tumor-reactive markers CD39 and CD103 being specifically demethylated. In addition, dynamic changes are observed during the transition of naïve T cells into tumor-reactive CD8+ T cells. Transcription factor binding motif enrichment analysis identifies several immune-related transcription factors, including three exhaustion-related genes (NR4A1, BATF, and EGR2) and VDR, which potentially play an important regulatory role in tumor-reactive CD8+ T cells. CONCLUSION: Our study supports the involvement of DNA methylation in shaping tumor-reactive and bystander CD8+ TILs, and provides a valuable resource for the development of novel DNA methylation markers and future therapeutics.

Lineage tracking reveals dynamic relationships of T cells in colorectal cancer.

T cells are key elements of cancer immunotherapy1 but certain fundamental properties, such as the development and migration of T cells within tumours, remain unknown. The enormous T cell receptor (TCR) repertoire, which is required for the recognition of foreign and self-antigens2, could serve as lineage tags to track these T cells in tumours3. Here we obtained transcriptomes of 11,138 single T cells from 12 patients with colorectal cancer, and developed single T cell analysis by RNA sequencing and TCR tracking (STARTRAC) indices to quantitatively analyse the dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. Although both CD8+ effector and 'exhausted' T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, most tumour-infiltrating T regulatory (Treg) cells showed clonal exclusivity, whereas certain Treg cell clones were developmentally linked to several T helper (TH) cell clones. Notably, we identified two IFNG+ TH1-like cell clusters in tumours that were associated with distinct IFNγ-regulating transcription factors -the GZMK+ effector memory T cells, which were associated with EOMES and RUNX3, and CXCL13+BHLHE40+ TH1-like cell clusters, which were associated with BHLHE40. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in patients with microsatellite-instable tumours, and this might explain their favourable responses to immune-checkpoint blockade. Furthermore, IGFLR1 was highly expressed in both CXCL13+BHLHE40+ TH1-like cells and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provide a powerful approach to dissect the T cell properties in colorectal cancer comprehensively, and could provide insights into the dynamic relationships of T cells in other cancers.

Publisher Correction: Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing.

In the version of this article originally published, the P statistic described in Fig. 3d was incorrect. It was described as "P < 22 × 10-16". It should have been "P < 2.2 × 10-16". Also, the "CD8+ Treg" label in Fig. 4f was incorrect. It should have been "CD4+ Treg". The errors have been corrected in the HTML and PDF versions of this article.

Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing.

Cancer immunotherapies have shown sustained clinical responses in treating non-small-cell lung cancer1-3, but efficacy varies and depends in part on the amount and properties of tumor infiltrating lymphocytes4-6. To depict the baseline landscape of the composition, lineage and functional states of tumor infiltrating lymphocytes, here we performed deep single-cell RNA sequencing for 12,346 T cells from 14 treatment-naïve non-small-cell lung cancer patients. Combined expression and T cell antigen receptor based lineage tracking revealed a significant proportion of inter-tissue effector T cells with a highly migratory nature. As well as tumor-infiltrating CD8+ T cells undergoing exhaustion, we observed two clusters of cells exhibiting states preceding exhaustion, and a high ratio of "pre-exhausted" to exhausted T cells was associated with better prognosis of lung adenocarcinoma. Additionally, we observed further heterogeneity within the tumor regulatory T cells (Tregs), characterized by the bimodal distribution of TNFRSF9, an activation marker for antigen-specific Tregs. The gene signature of those activated tumor Tregs, which included IL1R2, correlated with poor prognosis in lung adenocarcinoma. Our study provides a new approach for patient stratification and will help further understand the functional states and dynamics of T cells in lung cancer.

RNA editing of AZIN1 induces the malignant progression of non-small-cell lung cancers.

RNA editing is a widespread post-transcriptional mechanism that confers specific and reproducible nucleotide changes in selected RNA transcripts and plays a critical role in many human cancers. However, little is known about how RNA editing operates in non-small-cell lung cancers. Here, we measured the sequence and expression level of genes of antizyme inhibitor 1 and adenosine deaminase acting on RNA family in 30 non-small-cell lung cancer patient samples and 13 cell lines and revealed RNA editing S367G in antizyme inhibitor 1 is a high-frequent molecular events. We determined overexpression of antizyme inhibitor 1 with RNA editing, implying the oncogenic function of this alteration. We also detected the association of adenosine deaminase acting on RNA overexpression with RNA editing occurred in antizyme inhibitor 1. Furthermore, the RNA editing could cause a cytoplasmic-to-nuclear translocation of antizyme inhibitor 1 protein and conferred the malignant phenotype of non-small-cell lung cancer cells. The in vivo experiment confirmed that this RNA editing confers higher capacity of tumor migration as well. In conclusion, antizyme inhibitor 1 RNA editing and its involvement in tumorigenesis of non-small-cell lung cancer pave a new way for potential clinical management of non-small-cell lung cancer.

RNA editing is induced by type I interferon in esophageal squamous cell carcinoma.

In recent years, abnormal RNA editing has been shown to play an important role in the development of esophageal squamous cell carcinoma, as such abnormal editing is catalyzed by ADAR (adenosine deaminases acting on RNA). However, the regulatory mechanism of ADAR1 in esophageal squamous cell carcinomas remains largely unknown. In this study, we investigated ADAR1 expression and its association with RNA editing in esophageal squamous cell carcinomas. RNA sequencing applied to esophageal squamous cell carcinoma clinical samples showed that ADAR1 expression was correlated with the expression of STAT1, STAT2, and IRF9. In vitro experiments showed that the abundance of ADAR1 protein was associated with the induced activation of the JAK/STAT pathway by type I interferon. RNA sequencing results showed that treatment with type I interferon caused an increase in the number and degree of RNA editing in esophageal squamous cell carcinoma cell lines. In conclusion, the activation of the JAK/STAT pathway is a regulatory mechanism of ADAR1 expression and causes abnormal RNA editing profile in esophageal squamous cell carcinoma. This mechanism may serve as a new target for esophageal squamous cell carcinoma therapy.

Landscape of Infiltrating T Cells in Liver Cancer Revealed by Single-Cell Sequencing.

Systematic interrogation of tumor-infiltrating lymphocytes is key to the development of immunotherapies and the prediction of their clinical responses in cancers. Here, we perform deep single-cell RNA sequencing on 5,063 single T cells isolated from peripheral blood, tumor, and adjacent normal tissues from six hepatocellular carcinoma patients. The transcriptional profiles of these individual cells, coupled with assembled T cell receptor (TCR) sequences, enable us to identify 11 T cell subsets based on their molecular and functional properties and delineate their developmental trajectory. Specific subsets such as exhausted CD8+ T cells and Tregs are preferentially enriched and potentially clonally expanded in hepatocellular carcinoma (HCC), and we identified signature genes for each subset. One of the genes, layilin, is upregulated on activated CD8+ T cells and Tregs and represses the CD8+ T cell functions in vitro. This compendium of transcriptome data provides valuable insights and a rich resource for understanding the immune landscape in cancers.