Next Article in Journal
Genome-Wide Association Study Identifies Genetic Polymorphisms Associated with Estimated Minimum Effective Concentration of Fentanyl in Patients Undergoing Laparoscopic-Assisted Colectomy
Previous Article in Journal
Kallikrein-Related Peptidase 12 (KLK12) in Breast Cancer as a Favorable Prognostic Marker
Previous Article in Special Issue
Serum Amyloid A in Stable Patients with Chronic Obstructive Pulmonary Disease Does Not Reflect the Clinical Course of the Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Single-Cell Analysis in Immuno-Oncology

by
Maria-Ioanna Christodoulou
1,2 and
Apostolos Zaravinos
2,3,*
1
Tumor Immunology and Biomarkers Group, Basic and Translational Cancer Research Center (BTCRC), 1516 Nicosia, Cyprus
2
Department of Life Sciences, School of Sciences, European University Cyprus, 2404 Nicosia, Cyprus
3
Cancer Genetics, Genomics and Systems Biology Group, Basic and Translational Cancer Research Center (BTCRC), 1516 Nicosia, Cyprus
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(9), 8422; https://doi.org/10.3390/ijms24098422
Submission received: 31 March 2023 / Revised: 1 May 2023 / Accepted: 4 May 2023 / Published: 8 May 2023
(This article belongs to the Special Issue Immunogenic Cell Death, Immunogenic Surrender and Antitumor Immunity)

Abstract

:
The complexity of the cellular and non-cellular milieu surrounding human tumors plays a decisive role in the course and outcome of disease. The high variability in the distribution of the immune and non-immune compartments within the tumor microenvironments (TME) among different patients governs the mode of their response or resistance to current immunotherapeutic approaches. Through deciphering this diversity, one can tailor patients’ management to meet an individual’s needs. Single-cell (sc) omics technologies have given a great boost towards this direction. This review gathers recent data about how multi-omics profiling, including the utilization of single-cell RNA sequencing (scRNA-seq), assay for transposase-accessible chromatin with sequencing (scATAC-seq), T-cell receptor sequencing (scTCR-seq), mass, tissue-based, or microfluidics cytometry, and related bioinformatics tools, contributes to the high-throughput assessment of a large number of analytes at single-cell resolution. Unravelling the exact TCR clonotype of the infiltrating T cells or pinpointing the classical or novel immune checkpoints across various cell subsets of the TME provide a boost to our comprehension of adaptive immune responses, their antigen specificity and dynamics, and grant suggestions for possible therapeutic targets. Future steps are expected to merge high-dimensional data with tissue localization data, which can serve the investigation of novel multi-modal biomarkers for the selection and/or monitoring of the optimal treatment from the current anti-cancer immunotherapeutic armamentarium.

1. Introduction

The accumulating evidence over the last few decades supports the critical role of the immune system in the development and persistence of human cancers. The avoidance of immune destruction and tumor-promoting inflammation were added to Hanahan and Weinberg’s updated list of cancer hallmarks in 2011 [1]. The immune system appears to be a critical barrier to tumor formation and progression, according to research on genetically engineered murine models as well as clinical observations. Immune cell infiltrates surrounding tumors indicate an activated immune system that is targeting cancer cells. Indeed, the higher the rate of infiltration of cytotoxic T lymphocytes (CTLs) and natural killer (NK) cells, the better the patient’s prognosis [2,3]. Burnet introduced the concept of immune surveillance in the 1970s [4], and since then, cancer cells’ ability to avoid immune responses has become a focus for immune-oncologists seeking therapeutic interventions.
The tumor microenvironment (TME) is the diverse ecosystem that surrounds a tumor and plays an important role in its development and maintenance, as well as its response to treatment [5]. It consists of immune and non-immune cell subsets, an extracellular matrix (ECM), blood vessels, growth factors, chemokines, and cytokines, shaping the networks of interactions that impact the tumor’s phenotype [5,6]. These interplays are moderated by immune cells of diverse origin, some of them acting as “friends” and others as “foes” of cancer cells; T helper (Th) 1 and CTLs, NK and B cells, and M1 macrophages and mature dendritic cells (DC) act against tumor cells, while Th2, regulatory T cells (Tregs), M2 macrophages, neutrophils, cancer-associated fibroblasts (CAFs), immature DCs, and myeloid-derived suppressor cells (MDSCs) are in favor of their escape from immune surveillance [5]. Essentially, specific TME components serve as therapeutic targets; these may be immune checkpoints, such as PD-1/PD-L1, targeted by antibodies, or angiogenic factors, such as VEGF, targeted by inhibitory molecules, both of which are currently used in routine clinical practice and associated with improved clinical outcomes [7].
Despite the useful information that bulk next-generation sequencing (NGS) data provide, the signal produced by distinct, sometimes also rare, cell types, including subsets of circulating tumor cells (CTCs), immune cells, or cancer stem cells (CSCs), is often hindered due to the averaging of the signal stemming from the large and heterogenous populations of cells. During recent years, single-cell sequencing (sc-seq) techniques, i.e., genomic, transcriptomic, epigenomic, and/or proteomic analyses, have focused on the comparison of the certain -omes (i.e., genome, transcriptome, epigenome, and proteome, respectively) of individual cells. These techniques aim at unraveling the similarities and differences within a cell population, thus assessing its heterogeneity [8,9,10]. The evaluation of such differences among individual cells has provided the chance of detecting rare cell populations, including hyper-responsive immune cells within the TME [11], or tracing lineage and developmental relationships, including lymphocyte fate diversification [12]. Additionally, especially for sc-RNAseq, exploring gene co-expression patterns at the single-cell level can provide us with the opportunity to detect co-regulated gene modules and regulatory networks, which are possibly associated with the heterogeneity in cell function and specification [13,14]. Single-cell sequencing can have many other significant applications in the field of immuno-oncology, including immune infiltration, trajectory influence, functional enrichment and TCR repertoire analysis, allowing us to provide personalized medicine (Figure 1).
Applying sc-seq approaches to immune-oncology has allowed us to delve deeper into the mechanisms that drive an immune response or resistance to cancer immunotherapy, thus broadening its clinical use [15,16]. Indeed, a single immunotherapeutic scheme does not fit all patients, and even if overcoming primary resistance, patients who first respond to immunotherapy may finally develop secondary resistance and relapse [17]. Certain distribution patterns of tumor-infiltrating lymphocytes (TILs) and MDSCs, genetic alterations in IFNGR and JAK1/2, increased VEGF signaling, the absence of PD-L1 expression, or aberrant levels of cancer-derived microRNA molecules (miRNAs) are regarded as markers for resistance to immunotherapy [17,18]. The unravelling of the mechanisms underlying this unresponsiveness is of the utmost importance for the development of individualized approaches towards the improved survival rates of patients.
Other single-cell analysis techniques, including in situ hybridization (ISH) and flow cytometry, are also used as standard methods of detection, but these can detect only a few analytes; therefore, they have limited ability to distinguish between different cellular subtypes in the tumor microenvironment. Due to recent advances in sequencing techniques, cell-isolating microfluidics, and bioinformatic tools, we are now able to investigate big datasets and detect a large number of analytes at single-cell resolution across different cells of the TME [19,20,21,22].
In immuno-oncology, sc-seq can help us to better understand the interactions between tumor and immune cells, their heterogeneity, and how they are linked with the patient response to immunotherapy. Sc-seq can also give us a better understanding of the clonotypes and phenotypes of T cells, as well as of how immune checkpoints are distributed among them within the TME. The multi-omics profiling of the immune and non-immune cellular and non-cellular partners of the TME are essential for counteracting the limitations of current therapies and providing clinical benefits for a higher number of human cancer patients with malignancies [7].

2. Single-Cell Transcriptomics

2.1. Single-Cell RNA Sequencing (scRNA-Seq)

Separate cellular entities of a tumor (cancer cells, stromal cells, or immune cells) can be isolated and then analyzed at the single-cell level using scRNA-seq, where their transcript levels can be quantified. Technological developments and protocol improvements led to the creation of the Human Cell Atlas project (https://www.humancellatlas.org/, accessed on 1 March 2023), in which different (rare and novel) cell types, cellular states, and phylogenies are classified based on data from scRNA-seq projects [23,24] (Figure 2). The Human Cell Atlas can serve as a reference for the molecular profiles of the diverse cell subsets of human tissues with a healthy status and their association with phenotypes and functions [24]. It can also facilitate the comprehension of the pathobiology of physiological cellular processes, developmental courses, and homeostatic interactions among cells, thus helping with the elucidation of cellular and molecular deregulations during the onset of a disorder [24].
The TME is an intricate environment containing different cellular types that can be separated at a single level and processed without further management. scRNA-seq allows us to evaluate specific cellular subsets, such as cancer cells, immune, and stromal cells, providing information on the tumor heterogeneity and an outline of the TME with regard to the patient’s clinical image. Patient stratification and sampling before and after immunotherapy allows for a more comprehensive understanding of their response or resistance. Monitoring populations of immune cells, particularly T cells via paired TCR α-β subunit sequencing, allows us to characterize the T cell clones involved in the patient response [25].

2.2. Immune Response to Immunotherapy and Tumor, Immune, and Stromal Cells

In general, a thorough understanding of the composition and state of immune cells is critical not only for predicting patient responses (or resistance) to immunotherapies, but also for developing novel therapeutic approaches. The first attempts at this were through the development of comprehensive maps of the immune cell composition in non-cancerous tissues, such as lungs, lymph nodes, bone marrow, and blood [26]. These scRNA-seq profiles were later compared to multiple tumor types (skin melanoma, lung, colorectal, and breast cancer) and helped in analyzing the state of the immune cells in them (for example, activated CD8+ but not CD4+ T cells in the tumor setting) [27,28,29,30,31,32].
In principal, an scRNA-seq analysis recapitulates the existence of intra- and inter-tumoral heterogeneity. The clustering of different types of cancer cells has revealed that the transcriptomic profiles of single cells, clustered mainly by patient of origin [22,23,24], attest to the different evolutionary path that each tumor has. Since heterogeneity among different individuals with tumors can complicate the analysis of the patients’ responses to therapy, an analysis of individual immune cells or stromal cells using scRNA-seq allows for the gaining of meaningful information regarding the patients’ responses to different types of immunotherapy. Despite this heterogeneity, cancer cells with common transcriptional profiles across various patients can be detected. The transcriptomic profiling of single cancer cells has allowed for tumor type classification [33,34,35,36,37,38,39,40] and identified the signatures linked with patient responses (or resistance) to immunotherapies. An scRNA-seq analysis of the transcriptomes of skin melanoma patients who received immune checkpoint inhibitors (ICIs) by Jerby-Arnon et al. revealed a panel of deregulated genes expressed by cancer cells, which are linked to the exclusion of CD8+ T cells and resistance to therapy [35]. The up-regulated genes included ones involved in antigen processing and presentation, the IFN-γ signaling pathway, the response to complement fragments, and the migration of immune cells at the inflammatory lesion, while the downregulated genes included CDK4, E2F targets, transcription factors, and the targets of CDK7 and Myc [35].
On the other hand, scRNA-seq on CD45+-sorted cells from the immune checkpoint therapy (anti-PD-1 and/or anti-CTLA-4) of melanoma responders and non-responders revealed two distinct CD8+ T cell states associated with tumor progression or regression [28]: the first state contained upregulated levels of genes associated with memory, activation, and cell survival (e.g., IL7R, TCF7, and STAT4) and downregulated levels of co-inhibitory molecules, while the second state was enriched for genes associated with cell exhaustion (e.g., CD38, PDCD1, LAG3, CTLA4, and PTPN6). The ratio between the two states was proposed by the authors as a marker of the response/resistance to immune checkpoint therapy in these patients. Furthermore, TCF7 expression was linked to patient response, whereas CD39 and TIM3 co-expression were linked with resistance to ICI therapy in melanoma patients from the same study [28]. CD39 inhibition, combined with TIM3 and PD-1 (alone or with CTLA-4 inhibition), improved tumor control, suggesting that scRNA-seq can predict new targets for combinational immunotherapies [28].
Both the above described studies, using scRNA-seq either on cancer or immune cells of the TME, provided new insights on the crucial involvement of CD8+ T cells in driving resistance to immunotherapeutic agents, as well as suggested that specific transcriptomic signatures can predict the responders to immunotherapy, according to the state of the CD8+ T cells [41]. In general, a clear view of a tumor’s immune landscape and microenvironment will provide a better understanding of how tumor and immune cells evolve together during immunotherapy and relapse at the single-cell resolution level.

2.3. Interaction Analysis

Although it is critical to know how tumor and immune cells interact in cancer, several groups have used scRNA-seq to infer information on either isolated tumor or immune cells and/or stromal cells [42,43]. One study explored the interactome and found a link between tumor-associated DCs and peripheral T cells in humans with hepatocellular carcinoma (HCC) [44]. In detail, using two scRNA-seq protocols (10× Genomics and SMART-seq2), the researchers developed the transcriptomes of the CD45+ cells from tumor, adjacent liver, hepatic lymph node (LN), blood, and ascites samples of patients with HCC. The cluster of the tumor-derived LAMP3+ DCs exhibited the potential to migrate to LN, expressed a series of ligands to immune receptors, and were found to be capable of regulating various subsets of lymphocytes. Additionally, TAMs presented a certain transcriptome profile associated with a poor patient prognosis. Within this signature, two specific genes, namely SLC40A1 and GPNMB, were also found to correlate with unfavorable prognoses and play an inflammatory role in these cells, as attested by in vitro assays. In addition, an analysis revealed the tumor and blood origin of the myeloid and lymphoid cells in the ascites of patients with HCC [44]. Such interaction approaches can facilitate our comprehensive profiling of immune cell distribution and understanding of the dynamic interplays among immune cells in different cancer-involved tissue sites. This information may prove valuable for a deeper investigation and detection of the novel therapeutic targets and possible biomarkers for response-to-immunotherapy monitoring.

2.4. Limitations of scRNA-Seq

Nevertheless, scRNA-seq has some drawbacks, such as the fact that it produces noisy data due to the non-constant but pulsing transcription in eukaryotes [45,46]. In addition, it is difficult to make correlations between genotypes and phenotypes using transcriptome data alone. For example, we cannot confidently detect the somatic mutations in individual tumor cells using scRNA-seq; therefore, we cannot assess their implication in therapeutic response or resistance. Getting and processing samples on different days can also lead to batch effects and cause significant transcriptome-wide changes [47]. Various tools can correct these batch effects [48], but they should be carefully applied in order to not mask the true biological differences.

3. Single-Cell Proteomics

3.1. Mass Cytometry

Traditionally, the analysis of single-cell proteomics in next-generation immunotherapeutics relies on antibodies. Additionally, the use of peptide aptamers [49,50], non-immunoglobulin scaffolds, or affibody molecules [51] is a new, promising alternative to antibodies that can be used as high-affinity agents in the biomedical field [52,53]. Flow cytometry has been the main technique for studying protein expression at the single-cell level, despite the limited number of different epitopes that can be simultaneously detected due to a spectral overlap between the different fluorophores. Nevertheless, this limitation could be overcome with high-throughput mass cytometry [54]. This is about a fusion of the two technologies (mass spectrometry and flow cytometry), where a large number (>50) of immunophenotypic parameters can be simultaneously analyzed at single-cell resolution using antibodies labelled with metal isotopes [55].
As a whole, single-cell proteomics methodologies outperform the current scRNA-seq techniques in terms of their analytical specificity. In contrast to untargeted scRNA-seq, which enables the synchronized detection of thousands of RNA transcripts, the targeted approach of mass cytometry analyses, akin to flow cytometry, entails the deliberate selection of specific analytes, potentially introducing a degree of bias into the identification of cellular subpopulations. Moreover, the quantification of proteins is less susceptible to the variations seen in scRNA-seq data, as proteins generally exhibit a longer half-life. One can further study the post-translational modifications or assess the morphological features via cell morphometry [56,57]. Cell morphometry could be characterized as an advanced version of mass cytometry with “pathologists” properties. “Morphometric” markers can distinguish the major types of human lymphomas and leukemias based on morphologic differences in the principal hematopoietic cell populations of the healthy bone marrow [57]. Importantly, sole morphometric parameters have a unique diagnostic capacity, e.g., the expression of lamin A/C distinguishes healthy from neoplastic T cells, whereas lamin B1 suggests the presence of acute myeloid leukemia (AML) [57]. The incorporation of traditional surface markers into morphometric frames paves the way for more precise, automated diagnoses of complex haematological disorders.
Mass cytometry was shown to be very effective in the characterization of the specific cellular phenotypes involved in the response to anti-PD-1 immunotherapy at the single-cell level. Krieg et al. used high-dimensional single-cell mass cytometry to show that the frequency of CD14+CD16-HLA-DRhi monocytes in the peripheral blood of patients with metastatic melanoma has a strong potential for predicting progression-free and overall survival in response to anti-PD-1 immunotherapy and may serve in clinical decisions [58]. A similar study by Gide et al. [59] used mass cytometry to study the immune cell clusters in skin melanoma patients who had received either single (anti-PD-1) or combinatorial immunotherapy (anti-PD-1 plus anti-CTLA-4). They identified three distinct T cell clusters, characterized by specific markers for T-cell activation, differentiation, exhaustion, recent activation, tissue residency, and other functions of T cells, and associated them with the response to mono or dual immunotherapy, highlighting the significance of assessing immune cell phenotypes at the single-cell level. Specifically, they named the EOMES+CD69+CD45RO+ effector memory T cell subpopulation as a predictor of the response to the aforementioned immunotherapeutic schemes. Additionally, they revealed that the mRNA expression profiles of these cells are associated with better progression-free survival in anti-PD-1-treated patients and a reduction in the tumor size in anti-PD-1/CTLA-4-treated patients [59].

3.2. Emerging Single-Cell Proteome Analysis Technologies

Next-generation spectral analyzers have the potential to bring fluorescence-based flow cytometry up to the level of mass cytometry. These systems have an advantage over mass cytometry, in that they can use many fluorescence-labeled antibodies to sort and isolate cells for ensuing experimentations, without destroying the cells. One more emergent technology is antibody tagging using barcode oligonucleotides as opposed to fluorophores or heavy metal ions.
Tissue-based cytometry and microfluidic cytometry are two supplementary antibody-targeted single-cell approaches [60,61]. The former approach provides an enhanced level of spatial resolution, while the latter involves the immobilization of cells on a microfluidic chip, as opposed to their analysis in a flow stream.
Multiple markers can be quantified on the same cells using chip cytometry due to the iterative staining that the technique allows [62,63]. Small volumes of samples with a low number of living cells that are self-immobilized within microfluidic chips are analyzed with an unlimited number of intracellular and surface molecules, using automated epifluorescence microscopy. When compared to flow cytometry, it was shown that chip-based cytometry offers a significantly increased sensitivity and a specificity of a similar level [60]. This approach may be also of a high diagnostic value since it can provide the detection and phenotyping of rare cell subsets in precious, small clinical samples [60,64,65]. IsoLight is also a new system that detects the secreted proteins from single cells based on single-cell barcode chip technology [66]. Using microfluidic chips and sandwich-ELISA-like assays for the quantification of the secreted ligands, the system captures data from several hundred to several thousand individual viable cells over a period of multiple hours. This system is also versatile enough to accommodate additional analytes such as metabolites. Although at an early stage of development and under validation in the clinic, these platforms are promising tools in immuno-oncology, since they may provide solutions for patient profiling in association with immunotherapy response. Additionally, the single-cell, multiplex proteomic profiling of CD19 CAR-T products on microfluidics devices unravels a signature that is specific for the immune effector responses of CD19 CAR-T cells to specific antigens and offers a novel platform for capturing CAR-T pre-infusion products, and thus an insight into the efficacy of CAR-T cell therapy [66,67].

4. Single-Cell Genomics

There are different techniques for single-cell genome and exome sequencing [68], but their high costs do not allow for their application to high-dimensional analyses. As a result, these approaches hinder the scalability of analyzing a vast number of cells per patient. Instead, the deep sequencing of bulk tumor DNA samples to identify single-nucleotide variants (sSNVs) would be a more appropriate path, followed by the targeted DNA sequencing of single cells to identify putative driver mutations. One such system used for the mutational analysis of thousands of single cells was developed by Pellegrino et al. [69]. The authors devised a novel microfluidic strategy, which involved the barcoding of amplified genomic DNA extracted from thousands of AML cells that were restricted within droplets. Subsequently, these barcodes were utilized to reconstruct the genetic profiles of the cells from the NGS data, thereby enabling the routine investigation of AML heterogeneity and facilitating improved patient stratification and selection of therapy [70,71,72,73]. Indeed, certain DNA alterations have been already associated with resistance against certain therapies in AML patients: receptor tyrosine kinase (RTK) pathway mutations contribute to primary resistance against ivosidenib, RTK pathway mutations and 2-HG (hydroxyglutarate)-restoring mutations are linked to acquired resistance [73], while FLT3 internal tandem duplication gain or TP53 loss are associated with cross-resistance to both venetoclax and cytotoxic-based therapies [72].

5. Single-Cell Epigenomics

Various NGS methodologies have been developed for the analysis of the epigenome at the single-cell level, including transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), which assesses genome-wide chromatin accessibility, bisulfite-based DNA methylation sequencing, chromatin immunoprecipitation sequencing (ChIP-seq), which identifies the binding sites of DNA-associated proteins and can be used to map the global binding sites for a given protein, and chromosome conformation capture (3C) technologies (or Hi-C), which analyze the spatial organization of chromatin in a cell [74,75,76]. Of these, single-cell (sc)ATAC-seq is mostly used due to its high throughput. For example, Satpathy et al. used scATAC-seq to analyze the chromatin profiles of thousands of single cells from the serial biopsies of basal cell carcinoma patients [77]. The study revealed the certain chromatin regulators of T cell populations that were responsive to PD-1 blockade therapy and detected increased numbers of exhausted CD8+ T cells and CD4+ T follicular helper cells post-treatment [77].
The simultaneous profiling of chromatin accessibility and gene expression is expected to further advance the field of single-cell analysis in oncology [78,79]. In this context, multiplexed single-cell reduced-representation bisulfite sequencing was applied on B cells from healthy donors and patients with CLL, showing that the common clonal origin of CLL had a consistently increased epimutation rate, with a low variability in the cell-to-cell epimutation rate [80]. In addition, Cleavage Under Targets and Tagmentation (CUT&Tag) is a recently developed method that involves enzyme tethering, facilitating the generation of high-quality sequencing libraries suitable for the detailed analysis of a wide range of chromatin components. Kaya-Okur et al. [81] showed that this strategy is efficient in profiling histone modifications, RNA Polymerase II, and transcription factors on low cell numbers and single cells, and can thus broaden our ability to evaluate epigenetic modifications.

6. Single-Cell T Cell Receptor (TCR) Analysis

Most FDA-approved cancer immunotherapies focus on T cells, which are the main effectors of an anti-tumor immune response. As a consequence, there has been an increase in the interest in the analysis of individual T cell clonotypes. The TCR is made up of two ligand-binding chains (α and β or γ and δ) and six CD3 chains. The CD3 subunits form three dimers (δε, γε, and ζζ) and are essential for cell-surface expression and intracellular signaling (Malissen et al., 1999). Genetic recombination events lead to a highly diverse TCR repertoire. Most T cells contain individual TCRs and subunits that recognize a cognate peptide-MHC antigen [82]. Lymphocytes infiltrating immunogenic tumors (TILs), spearheaded by CD8+ T cells, can recognize, target, and kill cancer cells based on the interaction between their TCRs and the MHC-peptide complex [83]. Their “license to kill” can be determined using the single-cell sequencing of paired TCRs and subunits for TCR reconstruction and specificity testing. Moreover, for an immunotherapy to yield a tumor antigen-specific response, such an antigen should be presented as a complex with an MHC molecule, which will be recognized by the cognate TCRs on TILs.
Not all TILs can recognize and kill their target cells, though. Recent evidence has shown the existence of abundant bystander CD8+ T cells, which are phenotypically distinct and ineffective. Recent data suggest that we can quantify or isolate such bystander T cells through CD39 expression [84]. Therefore, many tumors are termed as “cold” due to the lack of proficient effector T cells, despite their large numbers in the tumor microenvironment [85].
Using single-cell sequencing, we can also identify paired TCR subunits to determine TIL specificity. This information can shed light on resistance to immunotherapy, but it can also help in discovering high-avidity anti-tumor TCRs for therapeutic purposes. Researchers initially sequenced individual TCR subunit combinations via TIL extraction and subcloning, but over the next years, they managed to sequence individual TCR complementarity-determining regions (CDR3) and antigen-specificity-allocating TCR subunits to individual T cells in a tumor [86,87,88,89]. With RNase H-dependent PCR-enabled T-cell receptor sequencing (rhTCRseq), for example, one can now define the alpha/beta TCR clonotypes in single cells or assess the alpha and beta TCR repertoire in bulk RNA samples on sorted tiny cell populations [86].
Droplet-based sequencing techniques now enable the analysis of paired TCR sequences at a higher throughput, which can subsequently be assessed and reconstructed within native T cells or reporter cell lines to determine their antigen specificity [90,91,92]. When scRNA-seq is combined with single-cell TCR sequencing (scTCR-seq), T cell phenotypes such as activation, memory, and exhaustion are linked to individual T cell TCR clonotypes. In order to exploit the capability to unambiguously assign single TCR sequences, microfluidic- and flow-cytometry-based approaches have been developed to directly identify the antigen-specific T cells in TIL, tissue, and peripheral blood samples [93,94,95]. Segaliny et al. developed a droplet microfluidics assay, where the functional screening and monitoring of single TCR T cell activation following the recognition of antigen-specific tumor cells in real-time can be performed [93]. It is important that each cell clone-of-interest can be sorted and used in downstream methods, laying the foundations for future novel cell-based immunotherapies. Tetramer-associated T-cell receptor sequencing (TetTCR-seq) is a high-throughput method that provides the ability to link TCR sequences to their cognate antigens in single cells [95]. Zhang et al. utilized this method to identify the patterns of TCR cross-reactivity with tumor neoantigens, followed by a rapid isolation of neoantigen-specific TCRs that do not cross-react with the wild-type antigen. This approach may potentially be applied in cancer immunotherapy for individual patients on the basis of an analysis of the antigen-binding T cells [95].
Recent clinical studies have suggested the usefulness of such approaches in the study of anti-tumor responses in certain patient cohorts. Using scTCR-seq, Keskin et al. analyzed glioblastoma patients treated with personalized neoantigen-targeting vaccines and found that the neoantigen-specific T cells from the peripheral blood of these patients could migrate to the intracranial glioblastoma tumor [96]. These findings show that neoantigen-targeting vaccines can favorably alter the immune milieu of glioblastomas.
A combination of scRNA-seq and scTCR-seq may be a very useful method for analyzing immune cell phenotypic and functional characteristics during immunotherapy treatment. An examination of biopsy samples obtained before and after anti-PD-1 treatment was administered to patients with advanced-stage basal cell or squamous cell carcinoma, revealed that the number of their activated and chronically activated (exhausted) CD8+ TILs, increased [97]. The post-treatment exhausted CD8+ TILs had the highest level of TCR clonality, indicating the lowest diversity. Additionally, most of these clonotypes were newly expanded clones that had not been previously identified. Prior to treatment, the exhausted TIL clones did not transition to non-exhausted phenotypes or expand. Although these clonotypes were not identified in the pre-treatment TILs, an analysis of the pre-treatment peripheral blood samples’ mass repertoire showed that 11.8% of the post-treatment TILs were new exhausted T cell clonotypes. These findings suggest that pre-existing tumor-specific TILs can only be partially reactivated and that ICIs can induce an influx and expansion of novel T cell clones (clonal replacement).
Moreover, a recent study involving the analysis of 8 primary uveal melanomas and 3 metastases through scTCR-seq and scRNA-seq, showed the existence of expanded populations of specific T cell clonotypes [98]. The study also found that among the exhaustion-associated immune checkpoints, LAG-3 exhibited the highest expression of genes, while PDCD1 and CTLA-4 had the lowest expressions. This finding may explain the lack of effectiveness of anti-PD-1 and anti-CTLA-4 therapies in this particular setting. Single-cell analyses, such as the one described, suggest that the effectiveness of ICIs may be dependent on non-exhausted T cells and T cells recruited from the periphery, rather than solely on the revitalization of exhausted T cells within the tumor.
The integration of multiple techniques, such as scRNA-seq, for finding respective TCR clonotypes, allows for a more detailed characterization of antigen-specific T cells and their roles in the response to immunotherapy. These comprehensive approaches significantly enhance immune cell analyses, enabling a more refined phenotypic profiling of T cells. Additionally, the non-invasive analysis of the expanded TCR clones in blood samples obtained during immunotherapy treatment has the potential to offer valuable understandings into the identity and functionality of the therapeutically relevant TCRs in malignancies. A recent study describing a practical method for detecting paired TCR subunit clonotypes in peripheral blood [99] may increase the capacity for a longitudinal analysis of blood samples collected from immunotherapy patients in clinical trials.

7. Analysis of High-Dimensional Spatial Data

To address the need to better understand how spatial diversity in the tumor immune milieu affects the patient response to immunotherapy, several spatial approaches have been developed for the diagnosis, prognosis, classification, and treatment of patients. A few examples include PD-L1 quantification using immunohistochemistry (IHC), immunofluorescence, and the in-situ hybridization of HER2 amplification.
Multiplexing IHC and immunofluorescence can also provide meaningful information for preclinical immunotherapy research. A meta-analysis of a large patient cohort with different types of cancer treated with anti-PD-1 or anti-PD-L1 antibodies revealed that, when these two were multiplexed, they provided a better prediction response of the patients to immunotherapy, in contrast to tumor mutational burden (TMB), gene expression profiling, and PD-L1 IHC [100,101]. The integration of data on protein co-expression and spatial relationships at the level of individual cells was crucial in facilitating an improved diagnostic accuracy. Supporting these results, Johnson et al. [102] showed that tumors with stronger spatial contacts between their PD-L1+ cells and PD-1+ non-tumor cells, as well as larger proportions of cancer cells co-expressing HLA-DR and IDO-1, were more likely to react to immunotherapy and be linked with better survival. The above spatial data on a single-cell resolution support that the sufficient and proximal expression of both PD-1 in T cells and PD-L1 in tumor cells is essential for effective ICI therapy.
The utilization of high-dimensional protein detection technologies has the potential to greatly augment the multiplexing capabilities of spatial analyses, thereby enabling the identification of 20–100 markers and significantly expanding the number of cell types and states that can be defined. This enhanced dimensionality permits a remarkable level of detail in the examination of the intricate spatial interactions between the networks of individual cell phenotypes. The spatial recognition of cellular subgroups in the tumor and/or its microenvironment offers the potential to identify prognostic biomarkers that are more accurate than those obtained through a single-cell phenotype quantification alone [103,104].
Helmink et al. used mass cytometry and scRNA-seq on melanoma patients treated with ICB and showed that switched memory B cells were enriched in the tumors of the responders to ICB [105]. Additionally, the authors showed a B cell clonal expansion and unique states of function. The data from this study supported the crucial functional involvement of B cells and the significance of their localization pattern in tertiary lymphoid structures (TLSs) for the response to ICB treatment, providing material for the further exploration of biomarkers and/or therapeutic target discovery.
Cabrita et al. investigated the role of B cells in the anti-tumor responses to metastatic melanomas [106]. The authors’ investigation revealed that B cells interact with the CD4+ T cells within TLSs, which, in turn, are associated with increased numbers of TCF7+-naive and memory T cells. In contrast, tumors lacking TLSs exhibited molecular signatures indicative of immunological dysfunction and/or exhaustion in TILs. To predict improved responses to ICIs, the authors developed a TLS gene signature that combined spatial data with scRNA-seq data. Taken together, these findings suggest that TLSs play a crucial role in the immune microenvironment of melanomas by conferring distinct T cell phenotypes [106].
Even though spatial recognition multiplexed methods are still in their early stages, they are revolutionizing our ability to investigate the contacts and physical interactions between tumor cells and cells in the microenvironment that impact tumor immunology and patient responses to immunotherapy.

8. Multimodal Analysis

8.1. Integration of Different Single Cell Sequencing Methods

Lately, there has been a growing trend of integrating different single-cell analysis techniques into combined technologies [107]. These approaches work towards deciphering the current state of the cell, determining the cell lineage, or are computational methods dictating the cells along a pseudotemporal trajectory. So far, such multimodal approaches to investigating immunotherapy responses have been limited, but a few of them that are promising in terms of their applicability are presented below [107,108].

8.2. Single-Cell DNA/RNA Combined with Protein Sequencing

There are various methodologies that integrate single-cell protein identification into scRNA-seq protocols (e.g., CITE-seq [109] or REAP-seq [110]). These can detect protein-based markers via scRNA-cDNA sequencing steps, allowing for the protein identification to progress without the restrictions set by spectral overlap or the number of resolvable metal ions. One such assay is the proximity extension assay [111], which allows for the creation of 96-plex panels for a protein detection using barcoded antibodies. The ability to detect a large number of proteins through such antibodies and incorporate this analysis into that of the transcriptome will significantly expand our capability to evaluate the TME. Furthermore, the assay’s performance is high in complex samples such as plasma and serum, using only 1 μL of specimen per test. It is expected that this multiplex technique will serve in the discovery of novel protein markers and their further exploration with regard to disease pathogenesis and responses to therapies. The Tapestri workflow also uses barcoded antibodies, combining protein expression with a mutational analysis at a single-cell resolution, enabling us to map somatic genotypes and clonal architectures with immunophenotypes [112]. This approach was applied for the analysis of clonal evolution in AML and unraveled the mechanisms of myeloid transformation, elucidating how the dynamics of clonal complexity accompany disease progression [112].
ScATAC-seq allows us to detect mitochondrial DNA (mtDNA) mutations, through which, we can trace cell lineage and cell fate (via the nuclear epigenome). One can reconstruct cell lineage relationships at a single-cell resolution with a high sensitivity and specificity using somatic mtDNA mutations [113,114]. scRNA-seq can also detect mtDNA mutations, but with a lower sensitivity due to its lower coverage [115]. The analysis of cell lineage is useful in cancer biology and immune-oncology, in understanding which cell clones expand or contract in response to immunotherapy, and are thus associated with resistance to the drug.

9. Future Directions

ScRNA-seq and its multiple iterations are groundbreaking new techniques that are helping us to understand the immunology of cancer cells and the processes driving patient responses or resistance to immunotherapy. While immunotherapies have revolutionized the treatment of various advanced-stage malignancies, the majority of patients do not respond or eventually develop a resistance to these therapies [16]. To tackle this issue, thousands of clinical studies have been conducted or are presently ongoing, which are exploring the efficacy of combination immunotherapy approaches that involve other immunotherapy medications, non-immunotherapy treatments, or a combination of both [116]. Nevertheless, the development of effective immunotherapy-containing combinations to overcome resistance requires a full characterization of the tumor immune microenvironment specific to each cancer type. Single-cell methods allow for complete TME profiling, making them ideal for researching tumor and immune cell heterogeneity, as well as how such changes in tumor biology contribute to treatment resistance [117,118].
In order to leverage single-cell technologies for enhancing our comprehension of immunotherapy response and resistance mechanisms and aiding in the development of innovative treatments, it is essential to meticulously design clinical trials that incorporate single-cell techniques in a practical and feasible manner. The current methods for single-cell analyses are expensive and require the separation of viable single cells from disaggregated tumor tissues, posing significant financial and logistical challenges that limit their widespread use in clinical settings. Nonetheless, there is now a steady decline in these sequencing costs, enabling the inclusion of single-cell analyses for a subset of patients achievable, especially if clinical trials are suitably planned in partnership with translational scientists. Valuable data from such studies can be furnished and extended further to more significant patient cohorts. The key factors to consider comprise the type and structure of the samples to be assembled (such as tumor tissue, non-malignant adjacent tissues, and/or peripheral blood), the optimal timing for collection (i.e., pretreatment, on-treatment, and/or post-treatment), and how the sample material will be quickly processed for an immediate single-cell analysis or cryopreserved for later examination [119]. Among immune TME partners, the isolation and assessment of the rare populations with key involvement in the regulation of local immune responses, such as regulatory lymphocytes, NK cells, or alternatively activated M2 macrophages, would unravel the key therapeutic targets associated with the progression of the disease. Importantly, as the management of cancer patients requires close monitoring for years after diagnosis and treatment, longitudinal approaches for the evaluation of the dynamic changes in sc-omics throughout the course of the disease and/or the response to therapy need to also be performed for the identification of promising biomarkers. In addition, the development of advanced computational tools, the standardization of data processing, and the analysis and involvement of specialized professionals are essential for universal and comparable approaches that could be integrated into routine clinical practices.
Although challenging, obtaining insights from a comprehensive evaluation of a limited number of patients (5–50) can have a significant impact on the analysis of a larger dataset, utilizing a less expensive and logistically simpler methodology (Figure 3). Using available scRNA-seq data rather than generating these data from scratch is an increasingly common method, as seen in a number of publications reported in this perspective. The Human Tumor Atlas Network’s formation has now institutionalized this sort of resource [120].
The exploration of high-dimensional single-cell data to identify the low-dimensional groups of genes and/or proteins that possess biomarker potential provides a clear roadmap towards therapeutic application. The application of bioinformatic methods to map RNA-seq data onto the cell types established by scRNA-seq and its associated approaches opens up a new route for the exploitation of data-rich single-cell information. These techniques for in silico dissection are less prone to the constraints that occur with employing only low-dimensional panels. CIBERSORTx is a machine-learning-based computational tool that can be used to deconvolute bulk gene expression data and estimate the cell type composition of complex tissues, including tumors, using a set of reference gene expression profiles [121,122]. This suggests that the highly detailed cell types identified through an RNA-seq data analysis can be measured in clinical samples that may be of a low quality, fixed, or challenging to isolate into individual cells. This could potentially remove the barriers to further research in the field of immuno-oncology. High-dimensional single-cell techniques that extend the information richness of established procedures such as IHC and flow cytometry might be more easily implemented in a clinic. Mass cytometry and flow cytometry panels with a higher number of fluorescent markers (such as 20–50) have previously been effectively utilized in clinical trials for samples from hematological malignancies and other diseases that may be readily disseminated into single cells (NCT03734198, NCT03915379, and NCT04181827). Spatially resolved multiplexed technologies have the potential to broaden our comprehension of the intercellular interactions that occur in the tumor immune environment by adding genotypic and phenotypic dimensions, and they represent the next frontier in explaining the processes driving immunotherapy resistance. Devices such as the in situ sequencing microfluidic tissue processor, which automates these technologies while staying compatible with conventional clinical laboratory procedures, indicate the potential future developments that will assist the implementation of these techniques [123].

10. Conclusions

The employment of scRNA-seq and comparable techniques is gaining prominence as an essential investigative tool for exploring immunotherapy response and resistance. Moreover, the analysis of scTCR is pivotal in assessing T cell antigen specificity, providing valuable insights into the adaptive immune system’s role in cancer progression, development, and treatment. The forthcoming horizon in single-cell analyses lies in merging high-dimensional content with precise tissue localization. High-dimensional analyses of single cells have various clinical implications, such as the detection of high-dimensional biomarkers for immunotherapy response, the precise design of low-dimensional biomarkers, an enhanced analysis of bulk sequencing data, and the advancement of our comprehension of the intricate spatial cell–cell interactions in the TME that underlie immunotherapy responses. Combination therapy is a promising treatment approach against cancer, given the rising diversity of immunotherapies accessible, as well as the customized nature of tumor heterogeneity and therapeutic resistance [124,125,126]. As a consequence, multimodal biomarker profiles will become more important in identifying the optimum therapy options for each cancer patient. The resolution and data richness necessary to develop these immuno-oncology signatures will be provided by high-dimensional single-cell technologies.

Funding

The APC was partly funded by Institutional funds.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hanahan, D.; Weinberg, R.A. Hallmarks of Cancer: The next Generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed]
  2. Nelson, B.H. The Impact of T-Cell Immunity on Ovarian Cancer Outcomes. Immunol. Rev. 2008, 222, 101–116. [Google Scholar] [CrossRef] [PubMed]
  3. Pagès, F.; Galon, J.; Dieu-Nosjean, M.-C.; Tartour, E.; Sautès-Fridman, C.; Fridman, W.-H. Immune Infiltration in Human Tumors: A Prognostic Factor That Should Not Be Ignored. Oncogene 2010, 29, 1093–1102. [Google Scholar] [CrossRef] [PubMed]
  4. Burnet, F.M. The Concept of Immunological Surveillance. Prog. Exp. Tumor Res. 1970, 13, 1–27. [Google Scholar] [CrossRef]
  5. Nascimento, C.; Ferreira, F. Tumor Microenvironment of Human Breast Cancer, and Feline Mammary Carcinoma as a Potential Study Model. Biochim. Biophys. Acta Rev. Cancer 2021, 1876, 188587. [Google Scholar] [CrossRef]
  6. Villanueva, J.; Herlyn, M. Melanoma and the Tumor Microenvironment. Curr. Oncol. Rep. 2008, 10, 439–446. [Google Scholar] [CrossRef]
  7. Christodoulou, M.-I.; Zaravinos, A. New Clinical Approaches and Emerging Evidence on Immune-Checkpoint Inhibitors as Anti-Cancer Therapeutics: CTLA-4 and PD-1 Pathways and Beyond. Crit. Rev. Immunol. 2019, 39, 379–408. [Google Scholar] [CrossRef]
  8. Yan, L.; Yang, M.; Guo, H.; Yang, L.; Wu, J.; Li, R.; Liu, P.; Lian, Y.; Zheng, X.; Yan, J.; et al. Single-Cell RNA-Seq Profiling of Human Preimplantation Embryos and Embryonic Stem Cells. Nat. Struct. Mol. Biol. 2013, 20, 1131–1139. [Google Scholar] [CrossRef]
  9. Deng, Q.; Ramsköld, D.; Reinius, B.; Sandberg, R. Single-Cell RNA-Seq Reveals Dynamic, Random Monoallelic Gene Expression in Mammalian Cells. Science 2014, 343, 193–196. [Google Scholar] [CrossRef]
  10. Jaitin, D.A.; Kenigsberg, E.; Keren-Shaul, H.; Elefant, N.; Paul, F.; Zaretsky, I.; Mildner, A.; Cohen, N.; Jung, S.; Tanay, A.; et al. Massively Parallel Single-Cell RNA-Seq for Marker-Free Decomposition of Tissues into Cell Types. Science 2014, 343, 776–779. [Google Scholar] [CrossRef]
  11. Shalek, A.K.; Satija, R.; Shuga, J.; Trombetta, J.J.; Gennert, D.; Lu, D.; Chen, P.; Gertner, R.S.; Gaublomme, J.T.; Yosef, N.; et al. Single-Cell RNA-Seq Reveals Dynamic Paracrine Control of Cellular Variation. Nature 2014, 510, 363–369. [Google Scholar] [CrossRef]
  12. Lönnberg, T.; Svensson, V.; James, K.R.; Fernandez-Ruiz, D.; Sebina, I.; Montandon, R.; Soon, M.S.F.; Fogg, L.G.; Nair, A.S.; Liligeto, U.; et al. Single-Cell RNA-Seq and Computational Analysis Using Temporal Mixture Modelling Resolves Th1/Tfh Fate Bifurcation in Malaria. Sci. Immunol. 2017, 2, eaal2192. [Google Scholar] [CrossRef]
  13. Liu, S.; Trapnell, C. Single-Cell Transcriptome Sequencing: Recent Advances and Remaining Challenges. F1000Res 2016, 5, F1000 Faculty Rev-182. [Google Scholar] [CrossRef]
  14. Wagner, A.; Regev, A.; Yosef, N. Revealing the Vectors of Cellular Identity with Single-Cell Genomics. Nat. Biotechnol. 2016, 34, 1145–1160. [Google Scholar] [CrossRef]
  15. Iorgulescu, J.B.; Braun, D.; Oliveira, G.; Keskin, D.B.; Wu, C.J. Acquired Mechanisms of Immune Escape in Cancer Following Immunotherapy. Genome Med. 2018, 10, 87. [Google Scholar] [CrossRef]
  16. Sharma, P.; Hu-Lieskovan, S.; Wargo, J.A.; Ribas, A. Primary, Adaptive, and Acquired Resistance to Cancer Immunotherapy. Cell 2017, 168, 707–723. [Google Scholar] [CrossRef]
  17. Fares, C.M.; Van Allen, E.M.; Drake, C.G.; Allison, J.P.; Hu-Lieskovan, S. Mechanisms of Resistance to Immune Checkpoint Blockade: Why Does Checkpoint Inhibitor Immunotherapy Not Work for All Patients? Am. Soc. Clin. Oncol. Educ. Book 2019, 39, 147–164. [Google Scholar] [CrossRef]
  18. Yi, M.; Xu, L.; Jiao, Y.; Luo, S.; Li, A.; Wu, K. The role of cancer-derived microRNAs in cancer immune escape. J. Hematol. Oncol. 2020, 13, 25. [Google Scholar] [CrossRef]
  19. Wang, L.; Livak, K.J.; Wu, C.J. High-Dimension Single-Cell Analysis Applied to Cancer. Mol. Aspects Med. 2018, 59, 70–84. [Google Scholar] [CrossRef]
  20. Finotello, F.; Eduati, F. Multi-Omics Profiling of the Tumor Microenvironment: Paving the Way to Precision Immuno-Oncology. Front. Oncol. 2018, 8, 430. [Google Scholar] [CrossRef]
  21. Gomes, T.; Teichmann, S.A.; Talavera-López, C. Immunology Driven by Large-Scale Single-Cell Sequencing. Trends Immunol. 2019, 40, 1011–1021. [Google Scholar] [CrossRef] [PubMed]
  22. Tang, X.; Huang, Y.; Lei, J.; Luo, H.; Zhu, X. The Single-Cell Sequencing: New Developments and Medical Applications. Cell Biosci. 2019, 9, 53. [Google Scholar] [CrossRef] [PubMed]
  23. Svensson, V.; Vento-Tormo, R.; Teichmann, S.A. Exponential Scaling of Single-Cell RNA-Seq in the Past Decade. Nat. Protoc. 2018, 13, 599–604. [Google Scholar] [CrossRef] [PubMed]
  24. Regev, A.; Teichmann, S.A.; Lander, E.S.; Amit, I.; Benoist, C.; Birney, E.; Bodenmiller, B.; Campbell, P.; Carninci, P.; Clatworthy, M.; et al. The Human Cell Atlas. eLife 2017, 6, e27041. [Google Scholar] [CrossRef]
  25. Gohil, S.H.; Iorgulescu, J.B.; Braun, D.A.; Keskin, D.B.; Livak, K.J. Applying High-Dimensional Single-Cell Technologies to the Analysis of Cancer Immunotherapy. Nat. Rev. Clin. Oncol. 2021, 18, 244–256. [Google Scholar] [CrossRef]
  26. Szabo, P.A.; Levitin, H.M.; Miron, M.; Snyder, M.E.; Senda, T.; Yuan, J.; Cheng, Y.L.; Bush, E.C.; Dogra, P.; Thapa, P.; et al. Single-Cell Transcriptomics of Human T Cells Reveals Tissue and Activation Signatures in Health and Disease. Nat. Commun. 2019, 10, 4706. [Google Scholar] [CrossRef]
  27. Lambrechts, D.; Wauters, E.; Boeckx, B.; Aibar, S.; Nittner, D.; Burton, O.; Bassez, A.; Decaluwé, H.; Pircher, A.; Van den Eynde, K.; et al. Phenotype Molding of Stromal Cells in the Lung Tumor Microenvironment. Nat. Med. 2018, 24, 1277–1289. [Google Scholar] [CrossRef]
  28. Sade-Feldman, M.; Yizhak, K.; Bjorgaard, S.L.; Ray, J.P.; de Boer, C.G.; Jenkins, R.W.; Lieb, D.J.; Chen, J.H.; Frederick, D.T.; Barzily-Rokni, M.; et al. Defining T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma. Cell 2018, 175, 998–1013.e20. [Google Scholar] [CrossRef]
  29. Zhang, L.; Yu, X.; Zheng, L.; Zhang, Y.; Li, Y.; Fang, Q.; Gao, R.; Kang, B.; Zhang, Q.; Huang, J.Y.; et al. Lineage Tracking Reveals Dynamic Relationships of T Cells in Colorectal Cancer. Nature 2018, 564, 268–272. [Google Scholar] [CrossRef]
  30. Azizi, E.; Carr, A.J.; Plitas, G.; Cornish, A.E.; Konopacki, C.; Prabhakaran, S.; Nainys, J.; Wu, K.; Kiseliovas, V.; Setty, M.; et al. Single-Cell Map of Diverse Immune Phenotypes in the Breast Tumor Microenvironment. Cell 2018, 174, 1293–1308.e36. [Google Scholar] [CrossRef]
  31. Zhang, L.; Zhang, Z. Recharacterizing Tumor-Infiltrating Lymphocytes by Single-Cell RNA Sequencing. Cancer Immunol. Res. 2019, 7, 1040–1046. [Google Scholar] [CrossRef]
  32. Andreatta, M.; Corria-Osorio, J.; Müller, S.; Cubas, R.; Coukos, G.; Carmona, S.J. Interpretation of T Cell States from Single-Cell Transcriptomics Data Using Reference Atlases. Nat. Commun. 2021, 12, 2965. [Google Scholar] [CrossRef]
  33. Filbin, M.G.; Tirosh, I.; Hovestadt, V.; Shaw, M.L.; Escalante, L.E.; Mathewson, N.D.; Neftel, C.; Frank, N.; Pelton, K.; Hebert, C.M.; et al. Developmental and Oncogenic Programs in H3K27M Gliomas Dissected by Single-Cell RNA-Seq. Science 2018, 360, 331–335. [Google Scholar] [CrossRef]
  34. Tirosh, I.; Izar, B.; Prakadan, S.M.; Wadsworth, M.H.; Treacy, D.; Trombetta, J.J.; Rotem, A.; Rodman, C.; Lian, C.; Murphy, G.; et al. Dissecting the Multicellular Ecosystem of Metastatic Melanoma by Single-Cell RNA-Seq. Science 2016, 352, 189–196. [Google Scholar] [CrossRef]
  35. Jerby-Arnon, L.; Shah, P.; Cuoco, M.S.; Rodman, C.; Su, M.-J.; Melms, J.C.; Leeson, R.; Kanodia, A.; Mei, S.; Lin, J.-R.; et al. A Cancer Cell Program Promotes T Cell Exclusion and Resistance to Checkpoint Blockade. Cell 2018, 175, 984–997.e24. [Google Scholar] [CrossRef]
  36. Izar, B.; Tirosh, I.; Stover, E.H.; Wakiro, I.; Cuoco, M.S.; Alter, I.; Rodman, C.; Leeson, R.; Su, M.-J.; Shah, P.; et al. A Single-Cell Landscape of High-Grade Serous Ovarian Cancer. Nat. Med. 2020, 26, 1271–1279. [Google Scholar] [CrossRef]
  37. Tirosh, I.; Venteicher, A.S.; Hebert, C.; Escalante, L.E.; Patel, A.P.; Yizhak, K.; Fisher, J.M.; Rodman, C.; Mount, C.; Filbin, M.G.; et al. Single-Cell RNA-Seq Supports a Developmental Hierarchy in Human Oligodendroglioma. Nature 2016, 539, 309–313. [Google Scholar] [CrossRef]
  38. Venteicher, A.S.; Tirosh, I.; Hebert, C.; Yizhak, K.; Neftel, C.; Filbin, M.G.; Hovestadt, V.; Escalante, L.E.; Shaw, M.L.; Rodman, C.; et al. Decoupling Genetics, Lineages, and Microenvironment in IDH-Mutant Gliomas by Single-Cell RNA-Seq. Science 2017, 355, eaai8478. [Google Scholar] [CrossRef]
  39. Neftel, C.; Laffy, J.; Filbin, M.G.; Hara, T.; Shore, M.E.; Rahme, G.J.; Richman, A.R.; Silverbush, D.; Shaw, M.L.; Hebert, C.M.; et al. An Integrative Model of Cellular States, Plasticity, and Genetics for Glioblastoma. Cell 2019, 178, 835–849.e21. [Google Scholar] [CrossRef]
  40. Young, M.D.; Mitchell, T.J.; Vieira Braga, F.A.; Tran, M.G.B.; Stewart, B.J.; Ferdinand, J.R.; Collord, G.; Botting, R.A.; Popescu, D.-M.; Loudon, K.W.; et al. Single-Cell Transcriptomes from Human Kidneys Reveal the Cellular Identity of Renal Tumors. Science 2018, 361, 594–599. [Google Scholar] [CrossRef]
  41. A Deep Dive into Immunotherapy Resistance. Cancer Discov. 2019, 9, 10. [CrossRef] [PubMed]
  42. Efremova, M.; Vento-Tormo, M.; Teichmann, S.A.; Vento-Tormo, R. CellPhoneDB: Inferring Cell-Cell Communication from Combined Expression of Multi-Subunit Ligand-Receptor Complexes. Nat. Protoc. 2020, 15, 1484–1506. [Google Scholar] [CrossRef] [PubMed]
  43. Cabello-Aguilar, S.; Alame, M.; Kon-Sun-Tack, F.; Fau, C.; Lacroix, M.; Colinge, J. SingleCellSignalR: Inference of Intercellular Networks from Single-Cell Transcriptomics. Nucleic Acids Res. 2020, 48, e55. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, Q.; He, Y.; Luo, N.; Patel, S.J.; Han, Y.; Gao, R.; Modak, M.; Carotta, S.; Haslinger, C.; Kind, D.; et al. Landscape and Dynamics of Single Immune Cells in Hepatocellular Carcinoma. Cell 2019, 179, 829–845.e20. [Google Scholar] [CrossRef]
  45. Chubb, J.R.; Trcek, T.; Shenoy, S.M.; Singer, R.H. Transcriptional Pulsing of a Developmental Gene. Curr. Biol. 2006, 16, 1018–1025. [Google Scholar] [CrossRef]
  46. Dar, R.D.; Razooky, B.S.; Singh, A.; Trimeloni, T.V.; McCollum, J.M.; Cox, C.D.; Simpson, M.L.; Weinberger, L.S. Transcriptional Burst Frequency and Burst Size Are Equally Modulated across the Human Genome. Proc. Natl. Acad. Sci. USA 2012, 109, 17454–17459. [Google Scholar] [CrossRef]
  47. van den Brink, S.C.; Sage, F.; Vértesy, Á.; Spanjaard, B.; Peterson-Maduro, J.; Baron, C.S.; Robin, C.; van Oudenaarden, A. Single-Cell Sequencing Reveals Dissociation-Induced Gene Expression in Tissue Subpopulations. Nat. Methods 2017, 14, 935–936. [Google Scholar] [CrossRef]
  48. Tran, H.T.N.; Ang, K.S.; Chevrier, M.; Zhang, X.; Lee, N.Y.S.; Goh, M.; Chen, J. A Benchmark of Batch-Effect Correction Methods for Single-Cell RNA Sequencing Data. Genome Biol. 2020, 21, 12. [Google Scholar] [CrossRef]
  49. Reverdatto, S.; Burz, D.S.; Shekhtman, A. Peptide Aptamers: Development and Applications. Curr. Top. Med. Chem. 2015, 15, 1082–1101. [Google Scholar] [CrossRef]
  50. Zhou, Z.; Liu, M.; Jiang, J. The Potential of Aptamers for Cancer Research. Anal. Biochem. 2018, 549, 91–95. [Google Scholar] [CrossRef]
  51. Ståhl, S.; Gräslund, T.; Eriksson Karlström, A.; Frejd, F.Y.; Nygren, P.-Å.; Löfblom, J. Affibody Molecules in Biotechnological and Medical Applications. Trends Biotechnol. 2017, 35, 691–712. [Google Scholar] [CrossRef]
  52. Löfblom, J.; Frejd, F.Y.; Ståhl, S. Non-Immunoglobulin Based Protein Scaffolds. Curr. Opin. Biotechnol. 2011, 22, 843–848. [Google Scholar] [CrossRef]
  53. Banta, S.; Dooley, K.; Shur, O. Replacing Antibodies: Engineering New Binding Proteins. Annu. Rev. Biomed. Eng. 2013, 15, 93–113. [Google Scholar] [CrossRef]
  54. Spitzer, M.H.; Nolan, G.P. Mass Cytometry: Single Cells, Many Features. Cell 2016, 165, 780–791. [Google Scholar] [CrossRef]
  55. Behbehani, G.K. Immunophenotyping by Mass Cytometry. Methods Mol. Biol. 2019, 2032, 31–51. [Google Scholar] [CrossRef]
  56. Fisher, J.; Sharma, R.; Don, D.W.; Barisa, M.; Hurtado, M.O.; Abramowski, P.; Porter, L.; Day, W.; Borea, R.; Inglott, S.; et al. Engineering ΓδT Cells Limits Tonic Signaling Associated with Chimeric Antigen Receptors. Sci. Signal. 2019, 12, eaax1872. [Google Scholar] [CrossRef]
  57. Tsai, A.G.; Glass, D.R.; Juntilla, M.; Hartmann, F.J.; Oak, J.S.; Fernandez-Pol, S.; Ohgami, R.S.; Bendall, S.C. Multiplexed Single-Cell Morphometry for Hematopathology Diagnostics. Nat. Med. 2020, 26, 408–417. [Google Scholar] [CrossRef]
  58. Krieg, C.; Nowicka, M.; Guglietta, S.; Schindler, S.; Hartmann, F.J.; Weber, L.M.; Dummer, R.; Robinson, M.D.; Levesque, M.P.; Becher, B. High-Dimensional Single-Cell Analysis Predicts Response to Anti-PD-1 Immunotherapy. Nat. Med. 2018, 24, 144–153. [Google Scholar] [CrossRef]
  59. Gide, T.N.; Quek, C.; Menzies, A.M.; Tasker, A.T.; Shang, P.; Holst, J.; Madore, J.; Lim, S.Y.; Velickovic, R.; Wongchenko, M.; et al. Distinct Immune Cell Populations Define Response to Anti-PD-1 Monotherapy and Anti-PD-1/Anti-CTLA-4 Combined Therapy. Cancer Cell 2019, 35, 238–255.e6. [Google Scholar] [CrossRef]
  60. Hennig, C.; Adams, N.; Hansen, G. A Versatile Platform for Comprehensive Chip-Based Explorative Cytometry. Cytometry A 2009, 75, 362–370. [Google Scholar] [CrossRef]
  61. Teo, J.; Mirenska, A.; Tan, M.; Lee, Y.; Oh, J.; Hong, L.Z.; Wnek, R.; Yap, Y.-S.; Shih, S.-J.; Bhagat, S.A.A.; et al. A Preliminary Study for the Assessment of PD-L1 and PD-L2 on Circulating Tumor Cells by Microfluidic-Based Chipcytometry. Future Sci. OA 2017, 3, FSO244. [Google Scholar] [CrossRef] [PubMed]
  62. Shi, Q.; Qin, L.; Wei, W.; Geng, F.; Fan, R.; Shin, Y.S.; Guo, D.; Hood, L.; Mischel, P.S.; Heath, J.R. Single-Cell Proteomic Chip for Profiling Intracellular Signaling Pathways in Single Tumor Cells. Proc. Natl. Acad. Sci. USA 2012, 109, 419–424. [Google Scholar] [CrossRef] [PubMed]
  63. Su, Y.; Ko, M.E.; Cheng, H.; Zhu, R.; Xue, M.; Wang, J.; Lee, J.W.; Frankiw, L.; Xu, A.; Wong, S.; et al. Multi-Omic Single-Cell Snapshots Reveal Multiple Independent Trajectories to Drug Tolerance in a Melanoma Cell Line. Nat. Commun. 2020, 11, 2345. [Google Scholar] [CrossRef] [PubMed]
  64. Ni, W.; Han, E.X.; Cyr, M.; Mackay, S.; Zhou, J. Applying Single-Cell Highly Multiplexed Secretome Proteomics to Characterize Immunotherapeutic Products and Predict Clinical Responses. Proteomics 2023, 2023, e2200242. [Google Scholar] [CrossRef]
  65. Li, L.; Yan, S.; Lin, B.; Shi, Q.; Lu, Y. Single-Cell Proteomics for Cancer Immunotherapy. Adv. Cancer Res. 2018, 139, 185–207. [Google Scholar] [CrossRef]
  66. Liu, D.; Paczkowski, P.; Mackay, S.; Ng, C.; Zhou, J. Single-Cell Multiplexed Proteomics on the IsoLight Resolves Cellular Functional Heterogeneity to Reveal Clinical Responses of Cancer Patients to Immunotherapies. Methods Mol. Biol. 2020, 2055, 413–431. [Google Scholar] [CrossRef]
  67. Xue, Q.; Bettini, E.; Paczkowski, P.; Ng, C.; Kaiser, A.; McConnell, T.; Kodrasi, O.; Quigley, M.F.; Heath, J.; Fan, R.; et al. Single-Cell Multiplexed Cytokine Profiling of CD19 CAR-T Cells Reveals a Diverse Landscape of Polyfunctional Antigen-Specific Response. J. Immunother. Cancer 2017, 5, 85. [Google Scholar] [CrossRef]
  68. Gawad, C.; Koh, W.; Quake, S.R. Single-Cell Genome Sequencing: Current State of the Science. Nat. Rev. Genet. 2016, 17, 175–188. [Google Scholar] [CrossRef]
  69. Pellegrino, M.; Sciambi, A.; Treusch, S.; Durruthy-Durruthy, R.; Gokhale, K.; Jacob, J.; Chen, T.X.; Geis, J.A.; Oldham, W.; Matthews, J.; et al. High-Throughput Single-Cell DNA Sequencing of Acute Myeloid Leukemia Tumors with Droplet Microfluidics. Genome Res. 2018, 28, 1345–1352. [Google Scholar] [CrossRef]
  70. Xu, L.; Durruthy-Durruthy, R.; Eastburn, D.J.; Pellegrino, M.; Shah, O.; Meyer, E.; Zehnder, J. Clonal Evolution and Changes in Two AML Patients Detected with A Novel Single-Cell DNA Sequencing Platform. Sci. Rep. 2019, 9, 11119. [Google Scholar] [CrossRef]
  71. Ediriwickrema, A.; Aleshin, A.; Reiter, J.G.; Corces, M.R.; Köhnke, T.; Stafford, M.; Liedtke, M.; Medeiros, B.C.; Majeti, R. Single-Cell Mutational Profiling Enhances the Clinical Evaluation of AML MRD. Blood Adv. 2020, 4, 943–952. [Google Scholar] [CrossRef]
  72. DiNardo, C.D.; Tiong, I.S.; Quaglieri, A.; MacRaild, S.; Loghavi, S.; Brown, F.C.; Thijssen, R.; Pomilio, G.; Ivey, A.; Salmon, J.M.; et al. Molecular Patterns of Response and Treatment Failure after Frontline Venetoclax Combinations in Older Patients with AML. Blood 2020, 135, 791–803. [Google Scholar] [CrossRef]
  73. Choe, S.; Wang, H.; DiNardo, C.D.; Stein, E.M.; de Botton, S.; Roboz, G.J.; Altman, J.K.; Mims, A.S.; Watts, J.M.; Pollyea, D.A.; et al. Molecular Mechanisms Mediating Relapse Following Ivosidenib Monotherapy in IDH1-Mutant Relapsed or Refractory AML. Blood Adv. 2020, 4, 1894–1905. [Google Scholar] [CrossRef]
  74. Kelsey, G.; Stegle, O.; Reik, W. Single-Cell Epigenomics: Recording the Past and Predicting the Future. Science 2017, 358, 69–75. [Google Scholar] [CrossRef]
  75. Shema, E.; Bernstein, B.E.; Buenrostro, J.D. Single-Cell and Single-Molecule Epigenomics to Uncover Genome Regulation at Unprecedented Resolution. Nat. Genet. 2019, 51, 19–25. [Google Scholar] [CrossRef]
  76. Ludwig, C.H.; Bintu, L. Mapping Chromatin Modifications at the Single Cell Level. Development 2019, 146, dev170217. [Google Scholar] [CrossRef]
  77. Satpathy, A.T.; Granja, J.M.; Yost, K.E.; Qi, Y.; Meschi, F.; McDermott, G.P.; Olsen, B.N.; Mumbach, M.R.; Pierce, S.E.; Corces, M.R.; et al. Massively Parallel Single-Cell Chromatin Landscapes of Human Immune Cell Development and Intratumoral T Cell Exhaustion. Nat. Biotechnol. 2019, 37, 925–936. [Google Scholar] [CrossRef]
  78. Cao, J.; Cusanovich, D.A.; Ramani, V.; Aghamirzaie, D.; Pliner, H.A.; Hill, A.J.; Daza, R.M.; McFaline-Figueroa, J.L.; Packer, J.S.; Christiansen, L.; et al. Joint Profiling of Chromatin Accessibility and Gene Expression in Thousands of Single Cells. Science 2018, 361, 1380–1385. [Google Scholar] [CrossRef]
  79. Reyes, M.; Billman, K.; Hacohen, N.; Blainey, P.C. Simultaneous Profiling of Gene Expression and Chromatin Accessibility in Single Cells. Adv. Biosyst. 2019, 3, 1900065. [Google Scholar] [CrossRef]
  80. Gaiti, F.; Chaligne, R.; Gu, H.; Brand, R.M.; Kothen-Hill, S.; Schulman, R.C.; Grigorev, K.; Risso, D.; Kim, K.-T.; Pastore, A.; et al. Epigenetic Evolution and Lineage Histories of Chronic Lymphocytic Leukaemia. Nature 2019, 569, 576–580. [Google Scholar] [CrossRef]
  81. Kaya-Okur, H.S.; Wu, S.J.; Codomo, C.A.; Pledger, E.S.; Bryson, T.D.; Henikoff, J.G.; Ahmad, K.; Henikoff, S. CUT&Tag for Efficient Epigenomic Profiling of Small Samples and Single Cells. Nat. Commun. 2019, 10, 1930. [Google Scholar] [CrossRef] [PubMed]
  82. Davis, M.M. T Cell Receptor Gene Diversity and Selection. Annu. Rev. Biochem. 1990, 59, 475–496. [Google Scholar] [CrossRef] [PubMed]
  83. Yang, J.C.; Rosenberg, S.A. Adoptive T-Cell Therapy for Cancer. Adv. Immunol. 2016, 130, 279–294. [Google Scholar] [CrossRef] [PubMed]
  84. Simoni, Y.; Becht, E.; Fehlings, M.; Loh, C.Y.; Koo, S.-L.; Teng, K.W.W.; Yeong, J.P.S.; Nahar, R.; Zhang, T.; Kared, H.; et al. Bystander CD8+ T Cells Are Abundant and Phenotypically Distinct in Human Tumour Infiltrates. Nature 2018, 557, 575–579. [Google Scholar] [CrossRef] [PubMed]
  85. Scheper, W.; Kelderman, S.; Fanchi, L.F.; Linnemann, C.; Bendle, G.; de Rooij, M.A.J.; Hirt, C.; Mezzadra, R.; Slagter, M.; Dijkstra, K.; et al. Low and Variable Tumor Reactivity of the Intratumoral TCR Repertoire in Human Cancers. Nat. Med. 2019, 25, 89–94. [Google Scholar] [CrossRef]
  86. Li, S.; Sun, J.; Allesøe, R.; Datta, K.; Bao, Y.; Oliveira, G.; Forman, J.; Jin, R.; Olsen, L.R.; Keskin, D.B.; et al. RNase H-Dependent PCR-Enabled T-Cell Receptor Sequencing for Highly Specific and Efficient Targeted Sequencing of T-Cell Receptor MRNA for Single-Cell and Repertoire Analysis. Nat. Protoc. 2019, 14, 2571–2594. [Google Scholar] [CrossRef]
  87. Dash, P.; Wang, G.C.; Thomas, P.G. Single-Cell Analysis of T-Cell Receptor Aβ Repertoire. Methods Mol. Biol. 2015, 1343, 181–197. [Google Scholar] [CrossRef]
  88. Han, A.; Glanville, J.; Hansmann, L.; Davis, M.M. Linking T-Cell Receptor Sequence to Functional Phenotype at the Single-Cell Level. Nat. Biotechnol. 2014, 32, 684–692. [Google Scholar] [CrossRef]
  89. Wang, G.C.; Dash, P.; McCullers, J.A.; Doherty, P.C.; Thomas, P.G. T Cell Receptor Aβ Diversity Inversely Correlates with Pathogen-Specific Antibody Levels in Human Cytomegalovirus Infection. Sci. Transl. Med. 2012, 4, 128ra42. [Google Scholar] [CrossRef]
  90. Hu, Z.; Anandappa, A.J.; Sun, J.; Kim, J.; Leet, D.E.; Bozym, D.J.; Chen, C.; Williams, L.; Shukla, S.A.; Zhang, W.; et al. A Cloning and Expression System to Probe T-Cell Receptor Specificity and Assess Functional Avidity to Neoantigens. Blood 2018, 132, 1911–1921. [Google Scholar] [CrossRef]
  91. Guo, X.-Z.J.; Dash, P.; Calverley, M.; Tomchuck, S.; Dallas, M.H.; Thomas, P.G. Rapid Cloning, Expression, and Functional Characterization of Paired Aβ and Γδ T-Cell Receptor Chains from Single-Cell Analysis. Mol. Ther. Methods Clin. Dev. 2016, 3, 15054. [Google Scholar] [CrossRef]
  92. Kobayashi, E.; Mizukoshi, E.; Kishi, H.; Ozawa, T.; Hamana, H.; Nagai, T.; Nakagawa, H.; Jin, A.; Kaneko, S.; Muraguchi, A. A New Cloning and Expression System Yields and Validates TCRs from Blood Lymphocytes of Patients with Cancer within 10 Days. Nat. Med. 2013, 19, 1542–1546. [Google Scholar] [CrossRef]
  93. Segaliny, A.I.; Li, G.; Kong, L.; Ren, C.; Chen, X.; Wang, J.K.; Baltimore, D.; Wu, G.; Zhao, W. Functional TCR T Cell Screening Using Single-Cell Droplet Microfluidics. Lab. Chip 2018, 18, 3733–3749. [Google Scholar] [CrossRef]
  94. Kula, T.; Dezfulian, M.H.; Wang, C.I.; Abdelfattah, N.S.; Hartman, Z.C.; Wucherpfennig, K.W.; Lyerly, H.K.; Elledge, S.J. T-Scan: A Genome-Wide Method for the Systematic Discovery of T Cell Epitopes. Cell 2019, 178, 1016.e13–1028.e13. [Google Scholar] [CrossRef]
  95. Zhang, S.-Q.; Ma, K.-Y.; Schonnesen, A.A.; Zhang, M.; He, C.; Sun, E.; Williams, C.M.; Jia, W.; Jiang, N. High-Throughput Determination of the Antigen Specificities of T Cell Receptors in Single Cells. Nat. Biotechnol. 2018, 1, 4282. [Google Scholar] [CrossRef]
  96. Keskin, D.B.; Anandappa, A.J.; Sun, J.; Tirosh, I.; Mathewson, N.D.; Li, S.; Oliveira, G.; Giobbie-Hurder, A.; Felt, K.; Gjini, E.; et al. Neoantigen Vaccine Generates Intratumoral T Cell Responses in Phase Ib Glioblastoma Trial. Nature 2019, 565, 234–239. [Google Scholar] [CrossRef]
  97. Yost, K.E.; Satpathy, A.T.; Wells, D.K.; Qi, Y.; Wang, C.; Kageyama, R.; McNamara, K.L.; Granja, J.M.; Sarin, K.Y.; Brown, R.A.; et al. Clonal Replacement of Tumor-Specific T Cells Following PD-1 Blockade. Nat. Med. 2019, 25, 1251–1259. [Google Scholar] [CrossRef]
  98. Durante, M.A.; Rodriguez, D.A.; Kurtenbach, S.; Kuznetsov, J.N.; Sanchez, M.I.; Decatur, C.L.; Snyder, H.; Feun, L.G.; Livingstone, A.S.; Harbour, J.W. Single-Cell Analysis Reveals New Evolutionary Complexity in Uveal Melanoma. Nat. Commun. 2020, 11, 496. [Google Scholar] [CrossRef]
  99. Tanno, H.; McDaniel, J.R.; Stevens, C.A.; Voss, W.N.; Li, J.; Durrett, R.; Lee, J.; Gollihar, J.; Tanno, Y.; Delidakis, G.; et al. A Facile Technology for the High-Throughput Sequencing of the Paired VH:VL and TCRβ:TCRα Repertoires. Sci. Adv. 2020, 6, eaay9093. [Google Scholar] [CrossRef]
  100. Lu, S.; Stein, J.E.; Rimm, D.L.; Wang, D.W.; Bell, J.M.; Johnson, D.B.; Sosman, J.A.; Schalper, K.A.; Anders, R.A.; Wang, H.; et al. Comparison of Biomarker Modalities for Predicting Response to PD-1/PD-L1 Checkpoint Blockade: A Systematic Review and Meta-Analysis. JAMA Oncol. 2019, 5, 1195–1204. [Google Scholar] [CrossRef]
  101. Lu, S.; Wang, H.; Taube, J.M. Different Biomarker Modalities and Response to Anti-PD-1/PD-L1 Therapies-Reply. JAMA Oncol. 2020, 6, 299. [Google Scholar] [CrossRef] [PubMed]
  102. Johnson, D.B.; Bordeaux, J.; Kim, J.Y.; Vaupel, C.; Rimm, D.L.; Ho, T.H.; Joseph, R.W.; Daud, A.I.; Conry, R.M.; Gaughan, E.M.; et al. Quantitative Spatial Profiling of PD-1/PD-L1 Interaction and HLA-DR/IDO-1 Predicts Improved Outcomes of Anti-PD-1 Therapies in Metastatic Melanoma. Clin. Cancer Res. 2018, 24, 5250–5260. [Google Scholar] [CrossRef] [PubMed]
  103. Jackson, H.W.; Fischer, J.R.; Zanotelli, V.R.T.; Ali, H.R.; Mechera, R.; Soysal, S.D.; Moch, H.; Muenst, S.; Varga, Z.; Weber, W.P.; et al. The Single-Cell Pathology Landscape of Breast Cancer. Nature 2020, 578, 615–620. [Google Scholar] [CrossRef] [PubMed]
  104. Bosisio, F.M.; Antoranz, A.; van Herck, Y.; Bolognesi, M.M.; Marcelis, L.; Chinello, C.; Wouters, J.; Magni, F.; Alexopoulos, L.; Stas, M.; et al. Functional Heterogeneity of Lymphocytic Patterns in Primary Melanoma Dissected through Single-Cell Multiplexing. Elife 2020, 9, e53008. [Google Scholar] [CrossRef]
  105. Helmink, B.A.; Reddy, S.M.; Gao, J.; Zhang, S.; Basar, R.; Thakur, R.; Yizhak, K.; Sade-Feldman, M.; Blando, J.; Han, G.; et al. B Cells and Tertiary Lymphoid Structures Promote Immunotherapy Response. Nature 2020, 577, 549–555. [Google Scholar] [CrossRef]
  106. Cabrita, R.; Lauss, M.; Sanna, A.; Donia, M.; Skaarup Larsen, M.; Mitra, S.; Johansson, I.; Phung, B.; Harbst, K.; Vallon-Christersson, J.; et al. Tertiary Lymphoid Structures Improve Immunotherapy and Survival in Melanoma. Nature 2020, 577, 561–565. [Google Scholar] [CrossRef]
  107. Stuart, T.; Satija, R. Integrative Single-Cell Analysis. Nat. Rev. Genet. 2019, 20, 257–272. [Google Scholar] [CrossRef]
  108. Optimizing Biological Inferences from Single-Cell Data. Nat. Rev. Genet. 2019, 20, 249. [CrossRef]
  109. Stoeckius, M.; Hafemeister, C.; Stephenson, W.; Houck-Loomis, B.; Chattopadhyay, P.K.; Swerdlow, H.; Satija, R.; Smibert, P. Simultaneous Epitope and Transcriptome Measurement in Single Cells. Nat. Methods 2017, 14, 865–868. [Google Scholar] [CrossRef]
  110. Peterson, V.M.; Zhang, K.X.; Kumar, N.; Wong, J.; Li, L.; Wilson, D.C.; Moore, R.; McClanahan, T.K.; Sadekova, S.; Klappenbach, J.A. Multiplexed Quantification of Proteins and Transcripts in Single Cells. Nat. Biotechnol. 2017, 35, 936–939. [Google Scholar] [CrossRef]
  111. Assarsson, E.; Lundberg, M.; Holmquist, G.; Björkesten, J.; Thorsen, S.B.; Ekman, D.; Eriksson, A.; Rennel Dickens, E.; Ohlsson, S.; Edfeldt, G.; et al. Homogenous 96-Plex PEA Immunoassay Exhibiting High Sensitivity, Specificity, and Excellent Scalability. PLoS ONE 2014, 9, e95192. [Google Scholar] [CrossRef]
  112. Miles, L.A.; Bowman, R.L.; Merlinsky, T.R.; Csete, I.S.; Ooi, A.T.; Durruthy-Durruthy, R.; Bowman, M.; Famulare, C.; Patel, M.A.; Mendez, P.; et al. Single-Cell Mutation Analysis of Clonal Evolution in Myeloid Malignancies. Nature 2020, 587, 477–482. [Google Scholar] [CrossRef]
  113. Ludwig, L.S.; Lareau, C.A.; Ulirsch, J.C.; Christian, E.; Muus, C.; Li, L.H.; Pelka, K.; Ge, W.; Oren, Y.; Brack, A.; et al. Lineage Tracing in Humans Enabled by Mitochondrial Mutations and Single-Cell Genomics. Cell 2019, 176, 1325–1339.e22. [Google Scholar] [CrossRef]
  114. Xu, J.; Nuno, K.; Litzenburger, U.M.; Qi, Y.; Corces, M.R.; Majeti, R.; Chang, H.Y. Single-Cell Lineage Tracing by Endogenous Mutations Enriched in Transposase Accessible Mitochondrial DNA. Elife 2019, 8, e45105. [Google Scholar] [CrossRef]
  115. Peng, A.; Mao, X.; Zhong, J.; Fan, S.; Hu, Y. Single-Cell Multi-Omics and Its Prospective Application in Cancer Biology. Proteomics 2020, 20, e1900271. [Google Scholar] [CrossRef]
  116. Xin Yu, J.; Hodge, J.P.; Oliva, C.; Neftelinov, S.T.; Hubbard-Lucey, V.M.; Tang, J. Trends in Clinical Development for PD-1/PD-L1 Inhibitors. Nat. Rev. Drug. Discov. 2020, 19, 163–164. [Google Scholar] [CrossRef]
  117. McGranahan, N.; Furness, A.J.S.; Rosenthal, R.; Ramskov, S.; Lyngaa, R.; Saini, S.K.; Jamal-Hanjani, M.; Wilson, G.A.; Birkbak, N.J.; Hiley, C.T.; et al. Clonal Neoantigens Elicit T Cell Immunoreactivity and Sensitivity to Immune Checkpoint Blockade. Science 2016, 351, 1463–1469. [Google Scholar] [CrossRef]
  118. Reuben, A.; Spencer, C.N.; Prieto, P.A.; Gopalakrishnan, V.; Reddy, S.M.; Miller, J.P.; Mao, X.; De Macedo, M.P.; Chen, J.; Song, X.; et al. Genomic and Immune Heterogeneity Are Associated with Differential Responses to Therapy in Melanoma. NPJ Genom. Med. 2017, 2, 10. [Google Scholar] [CrossRef]
  119. Valdes-Mora, F.; Handler, K.; Law, A.M.K.; Salomon, R.; Oakes, S.R.; Ormandy, C.J.; Gallego-Ortega, D. Single-Cell Transcriptomics in Cancer Immunobiology: The Future of Precision Oncology. Front. Immunol. 2018, 9, 2582. [Google Scholar] [CrossRef]
  120. Rozenblatt-Rosen, O.; Regev, A.; Oberdoerffer, P.; Nawy, T.; Hupalowska, A.; Rood, J.E.; Ashenberg, O.; Cerami, E.; Coffey, R.J.; Demir, E.; et al. The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution. Cell 2020, 181, 236–249. [Google Scholar] [CrossRef]
  121. Newman, A.M.; Steen, C.B.; Liu, C.L.; Gentles, A.J.; Chaudhuri, A.A.; Scherer, F.; Khodadoust, M.S.; Esfahani, M.S.; Luca, B.A.; Steiner, D.; et al. Determining Cell Type Abundance and Expression from Bulk Tissues with Digital Cytometry. Nat. Biotechnol. 2019, 37, 773–782. [Google Scholar] [CrossRef] [PubMed]
  122. Rusk, N. Expanded CIBERSORTx. Nat. Methods 2019, 16, 577. [Google Scholar] [CrossRef] [PubMed]
  123. Maïno, N.; Hauling, T.; Cappi, G.; Madaboosi, N.; Dupouy, D.G.; Nilsson, M. A Microfluidic Platform towards Automated Multiplexed in Situ Sequencing. Sci. Rep. 2019, 9, 3542. [Google Scholar] [CrossRef] [PubMed]
  124. Maus, M.V.; Fraietta, J.A.; Levine, B.L.; Kalos, M.; Zhao, Y.; June, C.H. Adoptive Immunotherapy for Cancer or Viruses. Annu. Rev. Immunol. 2014, 32, 189–225. [Google Scholar] [CrossRef]
  125. Martin-Liberal, J.; Ochoa de Olza, M.; Hierro, C.; Gros, A.; Rodon, J.; Tabernero, J. The Expanding Role of Immunotherapy. Cancer Treat. Rev. 2017, 54, 74–86. [Google Scholar] [CrossRef]
  126. Sanghera, C.; Sanghera, R. Immunotherapy—Strategies for Expanding Its Role in the Treatment of All Major Tumor Sites. Cureus 2019, 11, e5938. [Google Scholar] [CrossRef]
Figure 1. (a) The principle of single-cell sequencing includes the isolation of a single cell suspension, extracting their genetic material, and amplifying it for sequencing in order to study cell heterogeneity. (b) Single-cell sequencing has different applications in immune-oncology, from studying cell heterogeneity, to molecular mapping, as well as studying cell classification, immune infiltration, epigenetic changes, differential gene expression and others. Three common single-cell analysis techniques that are applied in immuno-oncology, include scRNA-seq, mass cytometry and flow cytometry.
Figure 1. (a) The principle of single-cell sequencing includes the isolation of a single cell suspension, extracting their genetic material, and amplifying it for sequencing in order to study cell heterogeneity. (b) Single-cell sequencing has different applications in immune-oncology, from studying cell heterogeneity, to molecular mapping, as well as studying cell classification, immune infiltration, epigenetic changes, differential gene expression and others. Three common single-cell analysis techniques that are applied in immuno-oncology, include scRNA-seq, mass cytometry and flow cytometry.
Ijms 24 08422 g001
Figure 2. Application of single-cell (sc) sequencing in immuno-oncology. (a) The TME is an intricate environment containing different cellular types that can be separated at a single level and processed without further management. scRNA-seq allows us to evaluate specific cellular subsets, such as cancer cells, immune, and stromal cells, providing information on the tumor heterogeneity and an outline of the TME with regard to the patient’s clinical image. (b) Patient response and resistance to cancer treatment can vary widely depending on the type of therapy used. One approach that has gained significant attention in recent years is immunotherapy. Patient stratification and sampling before and after immunotherapy allows for a more comprehensive understanding of their response or resistance. Monitoring populations of immune cells, particularly T cells, via paired TCR α-β subunit sequencing, allows us to characterize T cell clones involved in patient response. (c) Different aspects of the tumor’s biology, such as cell communication and signaling pathways, the identification of biomarkers, using techniques such as flow cytometry and mass cytometry, as well as therapeutic targets (e.g., cancer stem cells) that can be applied using single-cell sequencing technologies. (Figure adapted from Gohil et al. (2021) [25]).
Figure 2. Application of single-cell (sc) sequencing in immuno-oncology. (a) The TME is an intricate environment containing different cellular types that can be separated at a single level and processed without further management. scRNA-seq allows us to evaluate specific cellular subsets, such as cancer cells, immune, and stromal cells, providing information on the tumor heterogeneity and an outline of the TME with regard to the patient’s clinical image. (b) Patient response and resistance to cancer treatment can vary widely depending on the type of therapy used. One approach that has gained significant attention in recent years is immunotherapy. Patient stratification and sampling before and after immunotherapy allows for a more comprehensive understanding of their response or resistance. Monitoring populations of immune cells, particularly T cells, via paired TCR α-β subunit sequencing, allows us to characterize T cell clones involved in patient response. (c) Different aspects of the tumor’s biology, such as cell communication and signaling pathways, the identification of biomarkers, using techniques such as flow cytometry and mass cytometry, as well as therapeutic targets (e.g., cancer stem cells) that can be applied using single-cell sequencing technologies. (Figure adapted from Gohil et al. (2021) [25]).
Ijms 24 08422 g002
Figure 3. A generic strategy for adding single-cell data into immunotherapy clinical trials. (a) The pursuit of novel therapies for a limited patient cohort is a crucial goal in oncology. Achieving this goal relies heavily on a comprehensive assessment of pre-treatment, on-treatment, and post-treatment (resistance) tumor biopsies and peripheral immune cells, which can be achieved through single-cell and bulk techniques. Such methods allow us to identify critical biological factors underlying the efficacy of therapeutic interventions and the development of resistance mechanisms. These analyses would only be conducted on a small number of individuals since they are technically challenging and expensive. (b) Validation using traditional methods on a broader patient population. For a higher number of participants in clinical trials, it is possible to collect more standard samples such as frozen or formalin-fixed paraffin-embedded (FFPE) tumor tissue at baseline and peripheral blood cells before and during therapy. These samples can be utilized to confirm potential biological determinants of response and resistance to immunotherapy using standardized assays. PBMCs, peripheral blood mononuclear cells; IHC, immunohistochemistry; Hi-Dim, high-dimensional TCR, T-cell receptor; OCT, optimal cutting temperature compound; and scRNA-seq, single-cell RNA sequencing. (Figure adapted from Gohil et al. (2021) [25]).
Figure 3. A generic strategy for adding single-cell data into immunotherapy clinical trials. (a) The pursuit of novel therapies for a limited patient cohort is a crucial goal in oncology. Achieving this goal relies heavily on a comprehensive assessment of pre-treatment, on-treatment, and post-treatment (resistance) tumor biopsies and peripheral immune cells, which can be achieved through single-cell and bulk techniques. Such methods allow us to identify critical biological factors underlying the efficacy of therapeutic interventions and the development of resistance mechanisms. These analyses would only be conducted on a small number of individuals since they are technically challenging and expensive. (b) Validation using traditional methods on a broader patient population. For a higher number of participants in clinical trials, it is possible to collect more standard samples such as frozen or formalin-fixed paraffin-embedded (FFPE) tumor tissue at baseline and peripheral blood cells before and during therapy. These samples can be utilized to confirm potential biological determinants of response and resistance to immunotherapy using standardized assays. PBMCs, peripheral blood mononuclear cells; IHC, immunohistochemistry; Hi-Dim, high-dimensional TCR, T-cell receptor; OCT, optimal cutting temperature compound; and scRNA-seq, single-cell RNA sequencing. (Figure adapted from Gohil et al. (2021) [25]).
Ijms 24 08422 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Christodoulou, M.-I.; Zaravinos, A. Single-Cell Analysis in Immuno-Oncology. Int. J. Mol. Sci. 2023, 24, 8422. https://doi.org/10.3390/ijms24098422

AMA Style

Christodoulou M-I, Zaravinos A. Single-Cell Analysis in Immuno-Oncology. International Journal of Molecular Sciences. 2023; 24(9):8422. https://doi.org/10.3390/ijms24098422

Chicago/Turabian Style

Christodoulou, Maria-Ioanna, and Apostolos Zaravinos. 2023. "Single-Cell Analysis in Immuno-Oncology" International Journal of Molecular Sciences 24, no. 9: 8422. https://doi.org/10.3390/ijms24098422

APA Style

Christodoulou, M. -I., & Zaravinos, A. (2023). Single-Cell Analysis in Immuno-Oncology. International Journal of Molecular Sciences, 24(9), 8422. https://doi.org/10.3390/ijms24098422

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop