Unravelling the heterogeneity and dynamic relationships of tumor-infiltrating T cells by single-cell RNA sequencing analysis

Targeting the chromatin structural changes of antitumor immunity

Epigenomic imbalance drives abnormal transcriptional processes, promoting the onset and progression of cancer. Although defective gene regulation generally affects carcinogenesis and tumor suppression networks, tumor immunogenicity and immune cells involved in antitumor responses may also be affected by epigenomic changes, which may have significant implications for the development and application of epigenetic therapy, cancer immunotherapy, and their combinations. Herein, we focus on the impact of epigenetic regulation on tumor immune cell function and the role of key abnormal epigenetic processes, DNA methylation, histone post-translational modification, and chromatin structure in tumor immunogenicity, and introduce these epigenetic research methods. We emphasize the value of small-molecule inhibitors of epigenetic modulators in enhancing antitumor immune responses and discuss the challenges of developing treatment plans that combine epigenetic therapy and immunotherapy through the complex interaction between cancer epigenetics and cancer immunology.

Coordinated single-cell tumor microenvironment dynamics reinforce pancreatic cancer subtype

Bulk analyses of pancreatic ductal adenocarcinoma (PDAC) samples are complicated by the tumor microenvironment (TME), i.e. signals from fibroblasts, endocrine, exocrine, and immune cells. Despite this, we and others have established tumor and stroma subtypes with prognostic significance. However, understanding of underlying signals driving distinct immune and stromal landscapes is still incomplete. Here we integrate 92 single cell RNA-seq samples from seven independent studies to build a reproducible PDAC atlas with a focus on tumor-TME interdependence. Patients with activated stroma are synonymous with higher myofibroblastic and immunogenic fibroblasts, and furthermore show increased M2-like macrophages and regulatory T-cells. Contrastingly, patients with 'normal' stroma show M1-like recruitment, elevated effector and exhausted T-cells. To aid interoperability of future studies, we provide a pretrained cell type classifier and an atlas of subtype-based signaling factors that we also validate in mouse data. Ultimately, this work leverages the heterogeneity among single-cell studies to create a comprehensive view of the orchestra of signaling interactions governing PDAC.

The Contribution of Multiplexing Single Cell RNA Sequencing in Acute Myeloid Leukemia

Decades ago, the treatment for acute myeloid leukemia relied on cytarabine and anthracycline. However, advancements in medical research have introduced targeted therapies, initially employing monoclonal antibodies such as ant-CD52 and anti-CD123, and subsequently utilizing specific inhibitors that target molecular mutations like anti-IDH1, IDH2, or FLT3. The challenge lies in determining the role of these therapeutic options, considering the inherent tumor heterogeneity associated with leukemia diagnosis and the clonal drift that this type of tumor can undergo. Targeted drugs necessitate an examination of various therapeutic targets at the individual cell level rather than assessing the entire population. It is crucial to differentiate between the prognostic value and therapeutic potential of a specific molecular target, depending on whether it is found in a terminally differentiated cell with limited proliferative potential or a stem cell with robust capabilities for both proliferation and self-renewal. However, this cell-by-cell analysis is accompanied by several challenges. Firstly, the scientific aspect poses difficulties in comparing different single cell analysis experiments despite efforts to standardize the results through various techniques. Secondly, there are practical obstacles as each individual cell experiment incurs significant financial costs and consumes a substantial amount of time. A viable solution lies in the ability to process multiple samples simultaneously, which is a distinctive feature of the cell hashing technique. In this study, we demonstrate the applicability of the cell hashing technique for analyzing acute myeloid leukemia cells. By comparing it to standard single cell analysis, we establish a strong correlation in various parameters such as quality control, gene expression, and the analysis of leukemic blast markers in patients. Consequently, this technique holds the potential to become an integral part of the biological assessment of acute myeloid leukemia, contributing to the personalized and optimized management of the disease, particularly in the context of employing targeted therapies.

The Effect of Lung Resection for NSCLC on Circulating Immune Cells: A Pilot Study

This pilot study sought to evaluate the circulating levels of immune cells, particularly regulatory T-cell (Treg) subsets, before and after lung resection for non-small cell lung cancer. Twenty-five patients consented and had specimens collected. Initially, peripheral blood of 21 patients was collected for circulating immune cell studies. Two of these patients were excluded due to technical issues, leaving 19 patients for the analyses of circulating immune cells. Standard gating and high-dimensional unsupervised clustering flow cytometry analyses were performed. The blood, tumors and lymph nodes were analyzed via single-cell RNA and TCR sequencing for Treg analyses in a total of five patients (including four additional patients from the initial 21 patients). Standard gating flow cytometry revealed a transient increase in neutrophils immediately following surgery, with a variable neutrophil-lymphocyte ratio and a stable CD4-CD8 ratio. Unexpectedly, the total Treg and Treg subsets did not change with surgery with standard gating in short- or long-term follow-up. Similarly, unsupervised clustering of Tregs revealed a dominant cluster that was stable perioperatively and long-term. Two small FoxP3hi clusters slightly increased following surgery. In the longer-term follow-up, these small FoxP3hi Treg clusters were not identified, indicating that they were likely a response to surgery. Single-cell sequencing demonstrated six CD4+FoxP3+ clusters among the blood, tumors and lymph nodes. These clusters had a variable expression of FoxP3, and several were mainly, or only, present in tumor and lymph node tissue. As such, serial monitoring of circulating Tregs may be informative, but not completely reflective of the Tregs present in the tumor microenvironment.

T-cell exhaustion and stemness in antitumor immunity: Characteristics, mechanisms, and implications

T cells play a crucial role in the regulation of immune response and are integral to the efficacy of cancer immunotherapy. Because immunotherapy has emerged as a promising treatment for cancer, increasing attention has been focused on the differentiation and function of T cells in immune response. In this review, we describe the research progress on T-cell exhaustion and stemness in the field of cancer immunotherapy and summarize advances in potential strategies to intervene and treat chronic infection and cancer by reversing T-cell exhaustion and maintaining and increasing T-cell stemness. Moreover, we discuss therapeutic strategies to overcome T-cell immunodeficiency in the tumor microenvironment and promote continuous breakthroughs in the anticancer activity of T cells.

Divergent outcomes of anti-PD-L1 treatment coupled with host-intrinsic differences in TCR repertoire and distinct T cell activation states in responding versus non-responding tumors

Differential responses to immune checkpoint inhibitors (ICI) may be attributed to tumor-intrinsic factors or environmental cues; however, these mechanisms cannot fully explain the variable ICI responses in different individuals. Here, we investigate the potential contribution of immunological heterogeneity with a focus on differences in T-cell receptor (TCR) repertoire to ICI responses, which has not been defined previously. To reveal additional factors underlying heterogeneous responses to ICI, we employed a squamous cell carcinoma (SCC) mouse model in which tumor-bearing recipients unambiguously diverged into responders (R) or non-responders (NR) upon anti-PD-L1 treatment. Treatment efficacy absolutely required CD8 T-cells and correlated positively with effector functions of CD8 tumor-infiltrating lymphocytes (TILs). We showed that TCR repertoires exhibited a similar magnitude of clonal expansion in R vs. NR CD8 TILs. However, the top expanded TCR clonotypes appeared to be mutually exclusive between R and NR CD8 TILs, which also occurred in a recipient-specific manner, demonstrating preferential expansion of distinct TCR clonotypes against the same SCC tumor. Unexpectedly, R vs. NR CD8 TILs reached all activation clusters and did not exhibit substantial global differences in transcriptomes. By linking single-cell transcriptomic data with unique TCR clonotypes, CD8 TILs harboring top TCR clonotypes were found to occupy distinct activation clusters and upregulate genes favoring anti-tumor immunity to different extents in R vs. NR. We conclude that stochastic differences in CD8 TIL TCR repertoire and distinct activation states of top TCR clonotypes may contribute to differential anti-PD-L1 responses. Our study suggests that host-intrinsic immunological heterogeneity may offer a new explanation for differential ICI responses in different individuals, which could impact on strategies for personalized cancer immunotherapy.

Understanding initiation and progression of hepatocellular carcinoma through single cell sequencing

Unsatisfied clinical outcome drives to better understand hepatic carcinogenesis, microenvironment and escape of immune surveillance in hepatocellular carcinoma (HCC). Single cell RNA sequencing (scRNA-Seq) has generated enormous data to pinpoint pathophysiologic alterations in tumor microenvironment (TME) or trace lineage development in cancer stem cells (CSCs), circulating tumor cells (CTCs), and subsets of immune cells, such as exhausting T cells, tumor-associated macrophages (TAMs), dendritic cells or other lineages. New insights have significantly advanced current understanding in progression, poor responses to molecular-targeted therapeutics or immune checkpoint inhibitors, metastasis in both basic research and clinical practice. The present review intends to cover a basic workflow of the scRNA-seq technology, existing limitations and improvement areas. Moreover, in-depth understanding in TME, exhausting T cells, CSCs, CTCs, tumor-associated macrophages, dendritic cells in HCC facilitates implementation of personalized and precise therapy in an era of availability with an array of systemic regimens.

Cyclophosphamide enhances the antitumor potency of GITR engagement by increasing oligoclonal cytotoxic T cell fitness

Only a subset of cancer patients responds to checkpoint blockade inhibition in the clinic. Strategies to overcome resistance are promising areas of investigation. Targeting glucocorticoid-induced tumor necrosis factor receptor-related protein (GITR) has shown efficacy in preclinical models, but GITR engagement is ineffective in controlling advanced, poorly immunogenic tumors, such as B16 melanoma, and has not yielded benefit in clinical trials. The alkylating agent cyclophosphamide (CTX) depletes regulatory T cells (Tregs), expands tumor-specific effector T cells (Teffs) via homeostatic proliferation, and induces immunogenic cell death. GITR agonism has an inhibitory effect on Tregs and activates Teffs. We therefore hypothesized that CTX and GITR agonism would promote effective antitumor immunity. Here we show that the combination of CTX and GITR agonism controlled tumor growth in clinically relevant mouse models. Mechanistically, we show that the combination therapy caused tumor cell death, clonal expansion of highly active CD8+ T cells, and depletion of Tregs by activation-induced cell death. Control of tumor growth was associated with the presence of an expanded population of highly activated, tumor-infiltrating, oligoclonal CD8+ T cells that led to a diminished TCR repertoire. Our studies show that the combination of CTX and GITR agonism is a rational chemoimmunotherapeutic approach that warrants further clinical investigation.

The Ratio of Exhausted to Resident Infiltrating Lymphocytes Is Prognostic for Colorectal Cancer Patient Outcome

Immunotherapy success in colorectal cancer is mainly limited to patients whose tumors exhibit high microsatellite instability (MSI). However, there is variability in treatment outcomes within this group, which is in part driven by the frequency and characteristics of tumor-infiltrating immune cells. Indeed, the presence of specific infiltrating immune-cell subsets has been shown to correlate with immunotherapy response and is in many cases prognostic of treatment outcome. Tumor-infiltrating lymphocytes (TIL) can undergo distinct differentiation programs, acquiring features of tissue-residency or exhaustion, a process during which T cells upregulate inhibitory receptors, such as PD-1, and lose functionality. Although residency and exhaustion programs of CD8⁺ T cells are relatively well studied, these programs have only recently been appreciated in CD4⁺ T cells and remain largely unknown in tumor-infiltrating natural killer (NK) cells. In this study, we used single-cell RNA sequencing (RNA-seq) data to identify signatures of residency and exhaustion in colorectal cancer-infiltrating lymphocytes, including CD8⁺, CD4⁺, and NK cells. We then tested these signatures in independent single-cell data from tumor and normal tissue-infiltrating immune cells. Furthermore, we used versions of these signatures designed for bulk RNA-seq data to explore tumor-intrinsic mutations associated with residency and exhaustion from TCGA data. Finally, using two independent transcriptomic datasets from patients with colon adenocarcinoma, we showed that combinations of these signatures, in particular combinations of NK-cell activity signatures, together with tumor-associated signatures, such as TGFβ signaling, were associated with distinct survival outcomes in patients with colon adenocarcinoma.

Single-cell RNA sequencing in cancer: Applications, advances, and emerging challenges

Cancer has become one of the greatest threats to human health, and new technologies are urgently needed to further clarify the mechanisms of cancer so that better detection and treatment strategies can be developed. At present, extensive genomic analysis and testing of clinical specimens shape the insights into carcinoma. Nevertheless, carcinoma of humans is a complex ecosystem of cells, including carcinoma cells and immunity-related and stroma-related subsets, with accurate characteristics obscured by extensive genome-related approaches. A growing body of research shows that sequencing of single-cell RNA (scRNA-seq) is emerging to be an effective way for dissecting human tumor tissue at single-cell resolution, presenting one prominent way for explaining carcinoma biology. This review summarizes the research progress of scRNA-seq in the field of tumors, focusing on the application of scRNA-seq in tumor circulating cells, tumor stem cells, tumor drug resistance, the tumor microenvironment, and so on, which provides a new perspective for tumor research.

Applications of single-cell sequencing in cancer research: progress and perspectives

Single-cell sequencing, including genomics, transcriptomics, epigenomics, proteomics and metabolomics sequencing, is a powerful tool to decipher the cellular and molecular landscape at a single-cell resolution, unlike bulk sequencing, which provides averaged data. The use of single-cell sequencing in cancer research has revolutionized our understanding of the biological characteristics and dynamics within cancer lesions. In this review, we summarize emerging single-cell sequencing technologies and recent cancer research progress obtained by single-cell sequencing, including information related to the landscapes of malignant cells and immune cells, tumor heterogeneity, circulating tumor cells and the underlying mechanisms of tumor biological behaviors. Overall, the prospects of single-cell sequencing in facilitating diagnosis, targeted therapy and prognostic prediction among a spectrum of tumors are bright. In the near future, advances in single-cell sequencing will undoubtedly improve our understanding of the biological characteristics of tumors and highlight potential precise therapeutic targets for patients.

The Immune Atlas of Human Deciduas With Unexplained Recurrent Pregnancy Loss

Recurrent pregnancy loss (RPL) is a common fertility problem that affects 1%-2% of couples all over the world. Despite exciting discoveries regarding the important roles of the decidual natural killer cell (dNK) and regulatory T cell in pregnancy, the immune heterogeneity in patients with unexplained recurrent pregnancy loss (URPL) remains elusive. Here, we profiled the transcriptomes of 13,953 CD45⁺ cells from three normal and three URPL deciduas. Based on our data, the cellular composition revealed three major populations of immune cells including dNK cell, T cell, and macrophage, and four minor populations including monocytes, dendritic cell (DC), mast cell, and B cell. Especially, we identified a subpopulation of CSF1+ CD59+ KIRs-expressing dNK cells in normal deciduas, while the proportion of this subpopulation was decreased in URPL deciduas. We also identified a small subpopulation of activated dDCs that were accumulated mainly in URPL deciduas. Furthermore, our data revealed that in decidua at early pregnancy, CD8⁺ T cells exhibited cytotoxic properties. The decidual macrophages expressed high levels of both M1 and M2 feature genes, which made them unique to the conventional M1/M2 classification. Our single-cell data revealed the immune heterogeneity in decidua and the potentially pathogenic immune variations in URPL.

Evaluation of tumor immune contexture among intrinsic molecular subtypes helps to predict outcome in early breast cancer

Background

The prognosis of early breast cancer is linked to clinic-pathological stage and the molecular characteristics of intrinsic tumor cells. In some patients, the amount and quality of tumor-infiltrating immune cells appear to affect long term outcome. We aimed to propose a new tool to estimate immune infiltrate, and link these factors to patient prognosis according to breast cancer molecular subtypes.

Methods

We performed in silico analyses in more than 2800 early breast cancer transcriptomes with corresponding clinical annotations. We first developed a new gene expression deconvolution algorithm that accurately estimates the quantity of immune cell populations (tumor immune contexture, TIC) in tumors. Then, we studied associations between these immune profiles and relapse-free and overall survival among the different intrinsic molecular subtypes of breast cancer defined by PAM50 classification.

Results

TIC estimates the abundance of 15 immune cell subsets. Both myeloid and lymphoid subpopulations show different spread among intrinsic molecular breast cancer subtypes. A high abundance of myeloid cells was associated with poor outcome, while lymphoid cells were associated with favorable prognosis. Unsupervised clustering describing the 15 immune cell subsets revealed four subgroups of breast tumors associated with distinct patient survival, but independent from PAM50. Adding this information to clinical stage and PAM50 strongly improves the prediction of relapse or death.

Conclusions

Our findings make it possible to refine the survival stratification of early patients with breast cancer by incorporating TIC in addition to PAM50 and clinical tumor burden in a prognostic model validated in training and validation cohorts.

STARTRAC analyses of scRNAseq data from tumor models reveal T cell dynamics and therapeutic targets

Single-cell RNA sequencing is a powerful tool to examine cellular heterogeneity, novel markers and target genes, and therapeutic mechanisms in human cancers and animal models. Here, we analyzed single-cell RNA sequencing data of T cells obtained from multiple mouse tumor models by PCA-based subclustering coupled with TCR tracking using the STARTRAC algorithm. This approach revealed various differentiated T cell subsets and activation states, and a correspondence of T cell subsets between human and mouse tumors. STARTRAC analyses demonstrated peripheral T cell subsets that were developmentally connected with tumor-infiltrating CD8⁺ cells, CD4⁺ Th1 cells, and T reg cells. In addition, large amounts of paired TCRα/β sequences enabled us to identify a specific enrichment of paired public TCR clones in tumor. Finally, we identified CCR8 as a tumor-associated T reg cell marker that could preferentially deplete tumor-associated T reg cells. We showed that CCR8-depleting antibody treatment provided therapeutic benefit in CT26 tumors and synergized with anti–PD-1 treatment in MC38 and B16F10 tumor models.

Differences in TCR repertoire and T cell activation underlie the divergent outcomes of antitumor immune responses in tumor-eradicating versus tumor-progressing hosts

Background

Antitumor immunity is highly heterogeneous between individuals; however, underlying mechanisms remain elusive, despite their potential to improve personalized cancer immunotherapy. Head and neck squamous cell carcinomas (HNSCCs) vary significantly in immune infiltration and therapeutic responses between patients, demanding a mouse model with appropriate heterogeneity to investigate mechanistic differences.

Methods

We developed a unique HNSCC mouse model to investigate underlying mechanisms of heterogeneous antitumor immunity. This model system may provide a better control for tumor-intrinsic and host-genetic variables, thereby uncovering the contribution of the adaptive immunity to tumor eradication. We employed single-cell T-cell receptor (TCR) sequencing coupled with single-cell RNA sequencing to identify the difference in TCR repertoire of CD8 tumor-infiltrating lymphocytes (TILs) and the unique activation states linked with different TCR clonotypes.

Results

We discovered that genetically identical wild-type recipient mice responded heterogeneously to the same squamous cell carcinoma tumors orthotopically transplanted into the buccal mucosa. While tumors initially grew in 100% of recipients and most developed aggressive tumors, ~25% of recipients reproducibly eradicated tumors without intervention. Heterogeneous antitumor responses were dependent on CD8 T cells. Consistently, CD8 TILs in regressing tumors were significantly increased and more activated. Single-cell TCR-sequencing revealed that CD8 TILs from both growing and regressing tumors displayed evidence of clonal expansion compared with splenic controls. However, top TCR clonotypes and TCR specificity groups appear to be mutually exclusive between regressing and growing TILs. Furthermore, many TCRα/TCRβ sequences only occur in one recipient. By coupling single-cell transcriptomic analysis with unique TCR clonotypes, we found that top TCR clonotypes clustered in distinct activation states in regressing versus growing TILs. Intriguingly, the few TCR clonotypes shared between regressors and progressors differed greatly in their activation states, suggesting a more dominant influence from tumor microenvironment than TCR itself on T cell activation status.

Conclusions

We reveal that intrinsic differences in the TCR repertoire of TILs and their different transcriptional trajectories may underlie the heterogeneous antitumor immune responses in different hosts. We suggest that antitumor immune responses are highly individualized and different hosts employ different TCR specificities against the same tumors, which may have important implications for developing personalized cancer immunotherapy.

The Evolving Knowledge on T and NK Cells in Classic Hodgkin Lymphoma: Insights into Novel Subsets Populating the Immune Microenvironment

Simple Summary

In classic Hodgkin lymphoma, T and NK cells constitute a significant fraction of the reactive microenvironment established by malignant Hodgkin Reed–Sternberg cells. Despite their abundance, T and NK cells remain largely ineffective because of two coordinated levels of immune evasion. The first is based on the acquisition of regulatory properties or exhausted phenotypes that cripple their antitumor activity. The second is represented by their peculiar spatial distribution, with the most immunosuppressive subpopulations lying in close proximity of neoplastic cells. Recent discoveries about the functional role and the spatial orientation of T and NK cells in classic Hodgkin lymphoma are the focus of this review.

Abstract

Classic Hodgkin lymphoma (cHL) is a unique lymphoid neoplasm characterized by extensive immune infiltrates surrounding rare malignant Hodgkin Reed–Sternberg (HRS) cells. Different subsets of T and NK cells have long been recognized in the cHL microenvironment, yet their distinct contribution to disease pathogenesis has remained enigmatic. Very recently, novel platforms for high dimensional analysis of immune cells, such as single-cell RNA sequencing and mass cytometry, have revealed unanticipated insights into the composition of T- and NK-cell compartments in cHL. Advances in imaging techniques have better defined specific T-helper subpopulations physically interacting with neoplastic cells. In addition, the identification of novel cytotoxic subsets with an exhausted phenotype, typically enriched in cHL milieu, is shedding light on previously unrecognized immune evasion mechanisms. This review examines the immunological features and the functional properties of T and NK subsets recently identified in the cHL microenvironment, highlighting their pathological interplay with HRS cells. We also discuss how this knowledge can be exploited to predict response to immunotherapy and to design novel strategies to improve PD-1 blockade efficacy.

Unravelling the heterogeneity and dynamic relationships of tumor-infiltrating T cells by single-cell RNA sequencing analysis

T cells are crucial for the success of immune-based cancer therapy. Reinvigorating antitumor T cell activity by blocking checkpoint inhibitory receptors has provided clinical benefits for many cancer patients. However, the efficacy of these treatments varies in cancer patients and the mechanisms underlying these diverse responses remain elusive. The density and status of tumor-infiltrating T cells have been shown to positively correlate with patient response to checkpoint blockades. Therefore, further understanding of the heterogeneity, clonal expansion, migration, and effector functions of tumor-infiltrating T cells will provide fundamental insights into antitumor immune responses. To this end, recent advances in single-cell RNA sequencing technology have enabled profound and extensive characterization of intratumoral immune cells and have improved our understanding of their dynamic relationships. Here, we summarize recent progress in single-cell RNA sequencing technology and current strategies to uncover heterogeneous tumor-infiltrating T cell subsets. In particular, we discuss how the coupling of deep transcriptome information with T cell receptor (TCR)-based lineage tracing has furthered our understanding of intratumoral T cell populations. We also discuss the functional implications of various T cell subsets in tumors and highlight the identification of novel T cell markers with therapeutic or prognostic potential.