Applications of Single-Cell Omics to Dissect Tumor Microenvironment
- 1Institute of Respiratory Health, Frontiers Science Center for Disease-related Molecular Network, West China Hospital, Sichuan University, Chengdu, China
- 2Precision Medicine Research Center, West China Hospital, Sichuan University, Chengdu, China
- 3Department of Respiratory and Critical Care Medicine, West China Hospital, Sichuan University, Chengdu, China
- 4The Research Units of West China, Chinese Academy of Medical Sciences, West China Hospital, Chengdu, China
- 5Precision Medicine Key Laboratory of Sichuan Province, Chengdu, China
The recent technical and computational advances in single-cell sequencing technologies have significantly broaden our toolkit to study tumor microenvironment (TME) directly from human specimens. The TME is the complex and dynamic ecosystem composed of multiple cell types, including tumor cells, immune cells, stromal cells, endothelial cells, and other non-cellular components such as the extracellular matrix and secreted signaling molecules. The great success on immune checkpoint blockade therapy has highlighted the importance of TME on anti-tumor immunity and has made it a prime target for further immunotherapy strategies. Applications of single-cell transcriptomics on studying TME has yielded unprecedented resolution of the cellular and molecular complexity of the TME, accelerating our understanding of the heterogeneity, plasticity, and complex cross-interaction between different cell types within the TME. In this review, we discuss the recent advances by single-cell sequencing on understanding the diversity of TME and its functional impact on tumor progression and immunotherapy response driven by single-cell sequencing. We primarily focus on the major immune cell types infiltrated in the human TME, including T cells, dendritic cells, and macrophages. We further discuss the limitations of the existing methodologies and the prospects on future studies utilizing single-cell multi-omics technologies. Since immune cells undergo continuous activation and differentiation within the TME in response to various environmental cues, we highlight the importance of integrating multimodal datasets to enable retrospective lineage tracing and epigenetic profiling of the tumor infiltrating immune cells. These novel technologies enable better characterization of the developmental lineages and differentiation states that are critical for the understanding of the underlying mechanisms driving the functional diversity of immune cells within the TME. We envision that with the continued accumulation of single-cell omics datasets, single-cell sequencing will become an indispensable aspect of the immune-oncology experimental toolkit. It will continue to drive the scientific innovations in precision immunotherapy and will be ultimately adopted by routine clinical practice in the foreseeable future.
A tumor grows within a highly complex and dynamic local environment composed of immune cells, stromal cells, endothelial cells, as well as other non-cellular components such as secreted cytokines, chemokines and extracellular matrix (ECM). This complex ecosystem is collectively termed as the tumor microenvironment (TME), characterized by its profound heterogeneity, plasticity, and complex cross-interaction (Fridman et al., 2017). Our growing understanding of the TME has laid scientific foundation for cancer immunotherapy, which is arguably one of the greatest breakthroughs in cancer therapeutics over the past decade (Hodi et al., 2010; Robert et al., 2014; Garon et al., 2015; Ribas and Wolchok, 2018). The clinical success on immune checkpoint blockade (ICB) therapies has shifted the cancer therapeutic paradigm by demonstrating the great importance of T-cell mediated anti-tumor immunity; therefore, the TME has become a prime target for future anti-tumor therapies (Sharma and Allison, 2015; Wei et al., 2018). The recent technical and computational advances on single-cell sequencing provide a powerful tool to profile the TME cellular landscape in an unbiased and comprehensive way, which is enabling researchers to further dissect the cellular and molecular mechanisms underlying tumor-specific immune responses, tumor-immune cell interactions and immune evasion. Recent studies applying single-cell RNA sequencing (scRNA-seq) on the TME have yielded unprecedented resolution of its cellular and molecular complexity, which is likely to form a key determining factor in tumor progression and therapeutic response (Giladi and Amit, 2018; Papalexi and Satija, 2018; Ren et al., 2018; Gomes et al., 2019). We envision that single-cell sequencing will become an indispensable aspect of the immune-oncology experimental toolkit. With the continued accumulation of single-cell transcriptomics datasets and further extension to single-cell genomics, epigenomics and proteomics, we are geared to resolve not only the cellular composition and functional states of immune cells, but also the developmental history, regulatory network and cellular interactions of tumor and immune cells. Such increasing knowledge of the TME will pave the way for future precision immunotherapy through the establishment of connection between the TME within each patient’s tumor and the corresponding response to immunotherapy, and thus facilitate the design of optimal therapeutic strategies tailored to each patient. Meanwhile, unknown immune cell subtypes and states in the TME could be identified with the help of single-cell sequencing to yield novel targets for the development of future immunotherapy. In the review, we will discuss the recent advances on understanding the diversity of TME and its functional impact on tumor progression and immunotherapy response driven by single-cell sequencing, with the primary focus on the major immune cell types identified within the TME.
Recent TME Studies Driven by Single-Cell RNA Sequencing
Tumor Infiltrating Lymphocytes
Single-cell RNA sequencing has been widely used to profile the intratumoral immune landscape of various types of tumors, including skin, breast, lung, liver, colorectal, head and neck, brain, bladder, renal, and endometrial cancers (Table 1; Tirosh et al., 2016; Lavin et al., 2017; Müller et al., 2017; Puram et al., 2017; Zheng et al., 2017; Azizi et al., 2018; Guo et al., 2018; Jerby-Arnon et al., 2018; Lambrechts et al., 2018; Sade-Feldman et al., 2018; Savas et al., 2018; Zhang et al., 2018, 2019, 2020; Clarke et al., 2019; Li H. et al., 2019; Yost et al., 2019; Zilionis et al., 2019; Cillo et al., 2020; Oh D.Y. et al., 2020; Wu et al., 2020). Tumor infiltrating lymphocytes (TILs) have been the primary focus in many early studies since they are considered as the major immune cell type responsible for tumor-specific immunity. Accordingly, the level of T cell infiltration has been most frequently correlated with good prognosis in multiple cancer types (Fridman et al., 2017; Binnewies et al., 2018; Galon and Bruni, 2020). Immune checkpoint blockade therapy functions to restore the T cell-mediated anti-tumor immune responses by blocking the checkpoint receptors, such as PD-1 and CTLA-4, which are expressed in T cells, or their ligands, such as PD-L1, which is upregulated in tumor and antigen presenting cells (APCs) upon chronic antigen stimulation within the TME (Chen and Mellman, 2017; Wei et al., 2018; Callahan and Wolchok, 2019). Although the ICB therapy has achieved remarkable and durable clinical responses on selected patients and cancer types, its overall response rate has been limited, and many patients with initial response suffer from refractory disease or acquired resistance. The observed variation in ICB efficacy has been linked to various aspects of the TME, particularly correlated with signatures of the intratumoral T cell states, including overall T cell infiltration, activation and exhaustion (Sharma et al., 2017; Binnewies et al., 2018).
Cytotoxic CD8+ T Cells
Cytotoxic CD8+ T cells are considered as the primary T cell subtype responsible for direct tumor killing. Multiple scRNA-seq studies have demonstrated that intratumoral cytotoxic CD8+ T cells are predominantly in the dysfunctional and exhausted state, characterized by their lack of classical effector cytokine secretion and cytolytic functions, as well as their sustained expression of T cell exhaustion markers (Zheng et al., 2017; Azizi et al., 2018; Guo et al., 2018; Zhang et al., 2018; Li H. et al., 2019; Yost et al., 2019; Cillo et al., 2020). Moreover, comprehensive profiling on tens of thousands cytotoxic CD8+ T cells has yielded an unprecedented resolution to uncover a continuous spectrum of cell states from the early effector “transitioning” into the fully dysfunctional and exhausted state (Figure 1; Azizi et al., 2018; Li H. et al., 2019; Yost et al., 2019). Therefore, the intratumoral cytotoxic CD8+ T cells are not composed of discrete cell subtypes, but a continuum of cell states driven by ongoing activation and differentiation in response to TME stimuli. This level of complexity was under appreciated until the recent technical advances on large-scale single-cell sequencing.
Figure 1. Intratumoral expansion and differentiation of CD8+ T cells within the human TME. Upon tumor infiltration, tumor-reactive T cells that recognize tumor-specific antigens undergo antigen-driven T cell activation and differentiation, coupled with antigen-driven clonal expansion. The freshly infiltrated “pre-dysfunctional” tumor-reactive T cells continuously differentiate into various “dysfunctional” states, transitioning from the “early dysfunctional”, “dysfunctional” into the “late dysfunctional” state, characterized by the lack of cytolytic effector capacity and high expression of co-inhibitory markers (e.g., PDCD1, LAG3 and TIM3) as well as T cell exhaustion markers (e.g., CD39 and LAYN). In addition, the tumor-reactive T cells are usually double-positive for tissue-resident T cell marker CD103 and T cell exhaustion marker CD39, whereas the bystander T cells are usually negative for both. Resident-memory-like T cells (Trm) positive for CD103 and PDCD1 also present in the TME and are likely to be tumor-reactive. It remains unclear whether these Trm cells can differentiation into cytolytic tumor-reactive T cells and whether they correspond to the anti-PD-1 responsive TCF7+ progenitor-like T cell subset identified in mouse models.
To elucidate the developmental lineages of T cells along the differentiation axis, lineage tracing at the single-cell level is required. Joint analysis of TCR repertoire and transcriptome within a single T cell serves as a powerful tool to connect the phenotype of a given cell with its developmental lineage. Such analyses have revealed that the “dysfunctional” CD8+ T cells have undergone the highest level of clonal expansion among the various tumor infiltrating T cell subpopulations (Figure 1; Zheng et al., 2017; Azizi et al., 2018; Guo et al., 2018; Sade-Feldman et al., 2018; Zhang et al., 2018; Li H. et al., 2019; Yost et al., 2019). CD8+ T cells can directly recognize tumor cells in an MHC-I-dependent manner to activate the TCR signaling pathway; and then trigger T cell clonal expansion and effector functions to kill the targeted tumor cells via cytolytic activity and secretion of inflammatory cytokines. Therefore, the clonal expansion of intratumoral CD8+ T cells reflects their tumor reactivity. On the other hand, several recent studies have unveiled that not all intratumoral CD8+ T cells are equally tumor-reactive: the tumor-reactive CD8+ T cells are the most clonally expanded subpopulations with a tissue-resident memory phenotype and exhausted phenotype marked by CD103 and CD39, respectively; whereas the bystander CD8+ T cells are less clonally expanded and negative for CD103 and CD39 (Figure 1; Duhen et al., 2018; Simoni et al., 2018; Li H. et al., 2019; Scheper et al., 2019). Collectively, the functional state and antigen-specificity are the two key determinants of the anti-tumor capacity of CD8+ T cells; and the combination of scRNA-seq and scTCR repertoire analyses enables efficient identification and lineage tracing of tumor-reactive T cells within the human cancer TME.
By applying joint scRNA-seq and scTCR-seq, researchers have generated novel insight into the mechanism underlying PD-1 blockade directly from human specimens. Yost et al. (2019) showed that the majority of CD8+ T cells identified from basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) after PD-1 blockade did not share TCRs with the pre-existing intratumoral CD8+ T cells, suggesting that the pre-existing intratumoral CD8+ T cells had limited capacity of reinvigoration; instead, a fresh population of tumor-reactive CD8+ T cells, presumably originated outside of the TME, were recruited to the tumor site in response to PD-1 blockade. Similar observation was made on melanoma, suggesting that the CD8+ T cells presented predominately in post-ICB treated tumor samples hardly shared TCRs with the pre-treatment samples (Sade-Feldman et al., 2018). Moreover, despite their lack of reinvigoration capacity, dysfunctional tumor-reactive CD8+ T cells can secrete cytokine CXCL13, the ligand for CXCR5 that is commonly expressed on B cells, T follicular helper (TFH) cells, and follicular CXCR5-expressing CD8+ T cells (Moser et al., 2002; He et al., 2016; Thommen et al., 2018). These findings suggest that in spite of their altered effector functions, these dysfunctional CD8+ T cells may play a novel role in recruitment of other immune cells to the TME. Taken together, the recent single-cell studies on human TME has yielded tremendous insight into the functional diversity and developmental lineages of CD8+ T cells, advanced our understanding on the mechanism underlying PD-1 blockade, such that the pre-existing tumor-reactive dysfunctional CD8+ TILs may have limited reinvigoration capacity following PD-1 blockade, but they may function together with other immune cells (i.e., CD4+ T cells, B cells, dendritic cells (DCs), natural killer cells (NKs), etc.) to promote the priming and recruitment of de novo tumor-reactive T cells to the tumor. Therefore, although CD8+ T cells are the primary tumor-reactive cell type, their infiltration and functionality are largely modulated by other immune cell types, which are also essential for anti-tumor immunity.
CD4+ Helper T Cells
CD4+ T cells are further classified into pro-inflammatory CD4+ helper T (TH) cells and immunosuppressive CD4+ regulatory T (TREG) cells. In addition to their well-established role on promoting B cell activation and differentiation, TH cells also promote priming and clonal expansion of antigen-specific CD8+ T cells and modulate CD8+ T cell-intrinsic effector functionality in the TME and secondary lymphoid organs (Borst et al., 2018). Expansion of an ICOS+ TH1-like CD4+ effector subset was observed in response to CTLA4 blockade; therefore, TH cells, particularly the TH1 cells, are also pivotal players in T-cell mediated anti-tumor immunity and can respond to immunotherapy (Wei et al., 2017). Nevertheless, the cellular composition and functional plasticity of intratumoral TH cells are much less well-characterized compared to that of CD8+ T cells. Individual scRNA-seq studies have correlated the presence of specific subtypes of TH cells to clinical outcomes in certain cancer types. For example, scRNA-seq on the TME of colorectal tumors identified that intratumoral CXCL13+ TH1-like cells were enriched in microsatellite-instable (MSI) colorectal tumors, which is associated with significantly higher responsiveness to ICB treatment compared to the microsatellite-stable (MSS) colorectal cancers (Zhang et al., 2018). Moreover, intratumoral CD4+ T follicular helper (TFH) cells were reported to enrich in HPV+ head and neck squamous cell carcinoma (HNSCC) and associate with longer progression-free survival (Cillo et al., 2020). Notably, the formation of tertiary lymphoid structure (TLS) that involves TFH cells and B cells has been recently shown to promote immunotherapy response and survival in multiple cancer types (Thommen et al., 2018; Sautès-Fridman et al., 2019; Cabrita et al., 2020; Helmink et al., 2020; Petitprez et al., 2020). Therefore, CXCL13 expressing CD8+ and TH1 T cells may work together with TFH cells and B cells on facilitating the priming and recruitment of fresh tumor-reactive T cells to the TME. Additionally, a recent scRNA-seq study on bladder cancer identified two cytotoxic CD4+ TH subsets that were clonally expanded and capable of killing analogous tumor cells in an MHC-II-dependent manner. Meanwhile, gene signature of the cytotoxic CD4+ T cell subsets predicted clinical response in metastatic bladder cancer patients treated with anti-PD-L1 (Oh D.Y. et al., 2020). Conversely, scRNA-seq analyses on CRC and NSCLC indicated that the most clonally expanded cytotoxic CD4+ T cell subpopulation was rather enriched in blood and normal adjacent tissues, suggestive of their migratory feature and potential function outside of the TME (Guo et al., 2018; Zhang et al., 2018). Therefore, the cytotoxicity of CD4+ TH cells toward tumor cells is likely cancer type specific, presumably depends on the MHC-II expression level of different tumor cell types.
CD4+ Regulatory T Cells
In contrast to CD4+ helper T cells, which mainly participate in anti-tumor immunity and prevent tumor progression, CD4+ regulatory T (TREG) cells, characterized by the expression of lineage-specific transcription factor FOXP3, are the immunosuppressive T cells correlated with tumor progression and poor clinical outcomes. TREG cells are the key component of immune homeostasis that maintains self-tolerance; and dysfunction in TREG cells often leads to autoimmune diseases (Wing and Sakaguchi, 2010; Togashi et al., 2019). TREG cells exert their immunosuppressive functions through several independent mechanisms including the CTLA-4-mediated suppression of antigen-presenting cells, competition for and consumption of IL-2, secretion of immunosuppressive cytokines such as IL-10 and TGF-β, CD39-dependent conversion of ATP into adenosine, and direct cytotoxic effect on effector cells via granzyme and/or perforin secretion (Togashi et al., 2019). Intratumoral TREG cells are noted as one of the major immune cell types that contribute to the immunosuppressive TME in human cancers. The depletion of intratumoral TREG cells was observed from patients responsive to CTLA4 blockade, whereas the persisted or elevated presence of TREG cells may contribute to ICB resistance (Sharma et al., 2017; Tanaka and Sakaguchi, 2017; Wei et al., 2018; Togashi et al., 2019). scRNA-seq analyses on human cancer TME unveiled phenotypical distinctions of intratumoral TREG cells from circulating and tissue-resident TREG cells, and demonstrated that the intratumoral TREG cells underwent significantly greater clonal expansion compared to the other two TREG populations (Zheng et al., 2017; Guo et al., 2018; Zhang et al., 2018). The clonally expanded TCRs were mostly exclusive to intratumoral TREG themselves, suggesting that they were the tumor-reactive cells locally expanded upon activation by tumor-specific antigens within the TME (Zheng et al., 2017; Guo et al., 2018; Zhang et al., 2018). Furthermore, a considerable level of heterozygosity in terms of functional activation states was observed within the intratumoral TREG population; and only the “activated” TREG expressing high level of both co-inhibitory and co-stimulatory genes correlated with poor prognosis in NSCLC (Guo et al., 2018; Li H. et al., 2019). These findings are in large agreement with another independent functional study demonstrating that the activation and clonal expansion of intratumoral TREG cells within the TME were driven by their tumor antigen reactivity (Ahmadzadeh et al., 2019). Furthermore, through side-by-side scRNA-seq analysis on mouse and human TREG cells, Zemmour et al. (2018) showed that TREG cells could be generally divided into activated/memory and resting/naïve states, and TCR signaling seemed to shape the different facets of activated TREG cells such as their differentiation and effector functions. Therefore, it is conceivable that the functionality of intratumoral TREG cells are shaped by tumor antigens in the TME, contributing to their proliferation and functional plasticity. Whether intratumoral TREG cells were recruited from peripheral blood and/or adjacent non-malignant tissue, or converted from intratumoral TH cells, has not been unequivocally determined. Shared TCR clones were identified from paired analysis of peripheral blood and tumor sample of human breast cancer (Wang et al., 2019), indicating the blood origin of intratumoral TREG cells, whereas bulk RNA transcriptomic analyses on multiple cancer types suggested a closer phenotypic similarity between tumor and tissue-resident TREG cells (Plitas et al., 2016; Magnuson et al., 2018). As compared to blood and tissue-resident TREG cells, intratumoral TREG cells were upregulated for chemokine receptor CCR8 and exhaustion markers such as LAG3, CD39, and LAYN, possibly reflective of their activated and highly immunosuppressive state (De Simone et al., 2016; Plitas et al., 2016; Miragaia et al., 2019; Wang et al., 2019). ScRNA-seq analysis also identified low level of TCR repertoire sharing of intratumoral TREG cells with TREG from both blood and normal mucosa, as well as from intratumoral TH cells, arguing for a potential multi-origin model of tumor TREG cells (Zheng et al., 2017; Guo et al., 2018; Zhang et al., 2018). More extensive single-cell based analyses are necessary to further address the immunosuppressive functions of intratumoral TREG cells derived from different origins to fully understand the mechanism underlying intratumoral TREG cell recruitment, differentiation, and activation. Furthermore, close interaction between intratumoral TREG cells and tumor-associated macrophages (TAMs) was observed from human breast cancer, and this interaction may be important for the development of an immunosuppressive TME (Wang et al., 2019).
Taken together, intratumoral T cells display high level of complexity both in terms of their intrinsic functional states and interactions between other cell subtypes. Single-cell sequencing allowed the identification of “transitional” T cell states, delineation of T cell differentiation trajectory, predication of intercellular interactions, and association of novel T cell subsets with clinical outcomes. These findings significantly enhanced our understanding on T-cell mediated anti-tumor immunity and the mechanism underlying ICB therapy response.
Tumor Infiltrating Myeloid Cells
In addition to T lymphocytes, tumor infiltrating myeloid cells (TIMs) are also critical mediators of tumor progression and immune evasion. Myeloid cells can modulate T cell functions through their ability to present tumor antigens to T cells and to secrete key cytokines and chemokines for T cell differentiation and recruitment (Fridman et al., 2017; Binnewies et al., 2018). Thus, TIM cell composition appears to modulate lymphocytes infiltration, activation, and their anti-tumor functions. Understanding of the TIM compartment will shed important insight into the regulation of T-cell mediated anti-tumor immunity. However, due to the lack of sufficiently distinctive cell markers and questioned conservation between human and mouse model, human TIM cell composition is much less well-characterized compared to TILs (Zilionis et al., 2019). Single-cell sequencing provides an opportunity to comprehensively profile the TIM compartment independent of pre-defined cell surface markers and assumptions on cross-species conservation, thereby allowing for identification of novel cell lineages and/or cell states, as well as comparison of cell compositions between different species. ScRNA-seq has first been used to identify novel subtypes of DCs, monocytes and progenitor cells from human peripheral blood, scratching the surface of the previously under-appreciated complexity of the myeloid lineage (Villani et al., 2017). Meanwhile, single-cell sequencing is becoming an emerging strategy to validate results from animal models for translation to human immunobiology. For example, a recent scRNA-seq study conducted side-by-side comparison between human and mouse myeloid cell populations from lung cancer TME, and concluded that the population structures of intratumoral DCs, monocytes and neutrophils were largely conserved between mouse and human, whereas that of macrophages exhibited significant difference across species (Zilionis et al., 2019). Knowledge of the human-mouse correspondences and discrepancies of the TIM is tremendously valuable in guiding functional studies on distinct TIM subpopulations using mouse models, and in connecting the abundance of certain TIM subsets with clinical outcomes and therapeutic responses in both human and mouse models.
Dendritic cells comprise heterogenous lineages of myeloid cells, including conventional DCs (cDCs), plasmacytoid DCs (pDCs), and monocytes-derived DCs (mo-DCs). The increased abundance of intratumoral cDCs is generally associated with good prognosis, suggesting an anti-tumor role of cDCs (Hashimoto et al., 2011; Veglia and Gabrilovich, 2017; Michea et al., 2018; Moussion and Mellman, 2018). The primary function of cDCs in anti-tumor immunity is to acquire tumor-specific antigens at the tumor site, then migrate to tumor-draining lymph nodes (tdLNs) to prime T cells for de novo tumor-specific response. Intratumoral DCs identified from human liver and lung cancers by scRNA-seq consisted of two major subpopulations corresponding to cDC1 and cDC2, a third “activated” cDC subpopulation marked by CCR7 and LAMP3, and a fourth pDC subpopulation derived from a distinct lineage (Figure 2; Lavin et al., 2017; Zhang et al., 2019; Zilionis et al., 2019). cDC1, marked by CD141+ in human and CD103+ in mouse, excels in cross-presentation and priming of tumor antigens to CD8+ T cells; cDC2, marked by CD1c+ in human and CD11b+ in mouse, mainly primes CD4+ T cells. Specific depletion of CD141+ cDC1s was observed from early stage lung cancer TME by scRNA-seq and CyToF analyses, suggesting that DC depletion is one of the early hallmark events leading to immune evasion and tumor progression (Lavin et al., 2017). This observation is in line with mouse model studies on CD103+ cDC1s, such that the expansion and activation of CD103+ DCs dramatically increased CD8+ T cell recruitment to the tumor site and transformed response to checkpoint blockade (Roberts et al., 2016; Salmon et al., 2016). More specifically, the CD103+CCR7+ migratory subset of cDC1 is required for tumor antigen trafficking and priming of CD8+ T cells in the tumor draining lymph nodes (tdLN; Roberts et al., 2016; Salmon et al., 2016). The identification of an “activated” cDC subpopulation with predicted migratory capacity in human lung, liver and colorectal cancers further confirmed the conservation between human and mouse in DC biology (Zhang et al., 2019, 2020; Zilionis et al., 2019). The activated CCR7+LAMP3+ cDCs was comprised of populations derived from both cDC1 and cDC2 cells and interact intensively with T cells and NK cells via a rich array of chemokines and co-stimulatory, co-inhibitory molecules (Zhang et al., 2019, 2020). Two recent studies revealed that NK cells were involved in the recruitment of mouse CD103+ cDC1s to the TME through NK cell derived chemoattractants CCL5 and XCL1, and cytokine FLT3LG (Barry et al., 2018; Bottcher et al., 2018). Accordingly, NK cell abundance positively correlated with the presence of cDC1 in human melanoma, increased patient responsiveness to anti-PD-1 therapy, and better clinical outcome (Barry et al., 2018; Bottcher et al., 2018). Notably, the depletion of cDC1 from early stage lung cancer TME was coupled with reduced and impaired NK cells, further supporting the critical role of NK cells in cDC1 recruitment (Lavin et al., 2017). In addition to its role in tumor antigen presentation and CD8+ T cell priming, cDC1 was also shown to directly regulate CD8+ T cell functions and promote anti-PD-1 therapeutic efficacy through alternative mechanisms. Based on scRNA-seq analysis on mouse models, Garris et al. (2018) showed that intratumoral cDC1s specifically secreted IL-12, which was required for restoring the effector functions of CD8+ TILs following PD-1 blockade. Likewise, scRNA-seq analysis on mouse colorectal cancer model MC38 suggested that anti-CD40 agonist specifically activated a subset of cDC1, leading to upregulation of co-stimulatory marker CD80 and CD86, as well as cytokine IL-12 (Zhang et al., 2020). Additionally, cDC1s was reported to secrete chemokines CXCL9 and CXCL10, which interact with intratumoral T-cell expressed receptor CXCR3, to regulate CD8+ T cell response to PD-1 blockade (Chow et al., 2019). In addition to studies focusing on cDC1s, a scRNA-seq study on mouse and human tdLN samples identified two subpopulations of migratory cDC2s that were directly responsible for the priming of tumor antigen-specific CD4+ TH cells. More importantly, these cDC2 subsets were subjected to suppression by TREG cells in the TME, and the depletion of tumor TREG cells led to enhanced cDC2 migration and generation of functionally matured TH1-like cells characterized by ICOShiPD-1lo (Binnewies et al., 2019). Altogether, cDCs play multifaced anti-tumor roles in priming, recruiting and regulating intratumoral T cell differentiation, survival and effector functions.
Figure 2. Dendritic cell subtypes identified in the human TME and their functions. Conventional DCs, including the CD141+ cDC1 and CD1c+ cDC2, and plasmacytoid DCs (pDCs) are identified from the human TME. Both cDC1 and cDC2 undergo intratumoral activation and maturation into “activated” cDC1s and cDC2s, marked by CCR7 and LAMP3. The activated cDCs can present tumor antigens and migrate from the TME to tumor-draining lymph nodes (tdLN) and/or the tertiary lymphoid structure (TLS) to prime T cells for tumor antigen-specific reactivity. cDC1s mainly prime for CD8+ T cells and cDC2 mainly prime for CD4+ T cells. Tumor-antigen primed T cells are then recruited back to the TME via chemokine axis secreted by the activated cDCs and other immune cells from the TME. DCs also interact extensively with various subtypes of lymphocytes, such as natural killer cells (NKs), regulatory T cells (TREG) and CD8+ T cells, via a rich array of chemokine and cytokine secretion. NKs can recruit and activate cDC1s by secreting XCL1, CCL5 and FLT3L; TREGs can directly bind to cDC2s and inhibit their migration and priming of CD4+ T cells; and cDC1 can secrete IL-12 to modulate CD8+ T cell response to PD-1 blockade.
In contrast to DCs, the population structure of intratumoral macrophages exhibited considerable human-specific pattern, highlighting the importance on direct analysis of human samples in understanding the pro- and anti-tumor roles of the monocyte/macrophage lineage in human TME (Zilionis et al., 2019). The major monocyte/macrophage lineage-derived cell populations identified in the human TME by scRNA-seq and CyToF analyses consist of classical monocytes (CD14+CD16–), inflammatory monocytes (CD14intCD16+), immigrant macrophages (HLA-DRintCD192+), tissue-resident macrophages (HLA-DRintCD206+), TAMs (HLA-DRhiCD68+CD64+), myeloid-derived suppressor cells (MDSCs, HLA-DR–/low), and monocyte-derived DCs (mo-DCs; Chevrier et al., 2017; Lavin et al., 2017; Azizi et al., 2018; Li H. et al., 2019; Wagner et al., 2019; Zhang et al., 2019; Zilionis et al., 2019; Cillo et al., 2020). Each of these populations can further diversify into a spectrum of activation and differentiation states in response to TME stimuli. Three convergent differentiation trajectories were identified across different cancer types by scRNA-seq analyses: (1) the first trajectory was an ongoing differentiation of monocytes from classical circulating monocytes derived from the tumor vasculature into activated intratumoral monocytes (Azizi et al., 2018; Li H. et al., 2019); (2) the second trajectory was the differentiation of intratumoral monocytes into TAMs (Azizi et al., 2018); and (3) and the third trajectory was the differentiation of classical circulating monocytes into early immigrant macrophages and then into mature tissue-resident macrophages (Figure 3; Chevrier et al., 2017; Wagner et al., 2019). Compared to juxta-tumoral stroma and normal adjacent tissue, tumors were generally enriched for TAMs, intratumoral monocytes and MDSCs, whereas depleted of tissue-resident macrophages, classical circulating monocytes, and inflammatory monocytes (Lavin et al., 2017; Azizi et al., 2018; Li H. et al., 2019; Wagner et al., 2019; Zhang et al., 2019). TAMs comprise a highly heterogenous population derived either from monocytes or tissue-resident macrophages; and they simultaneously express activated anti-tumor M1 markers and pro-tumor M2 markers, suggesting that the classical M1/M2 polarization paradigm may not be applicable to human TAMs (Lavin et al., 2017; Müller et al., 2017; Azizi et al., 2018; Binnewies et al., 2018; Wagner et al., 2019; Zhang et al., 2019). Although the infiltration of TAMs in the TME has been generally associated with tumor progression and poor prognosis, the functional role of TAMs has been controversial and may vary among different cancer types and stages (Keeley et al., 2019). A scRNA-seq study on paired early-stage lung adenocarcinoma and non-cancerous lung tissue identified a tumor-enriched TAM population preferentially expressing high level of PPARγ, TREM, CD81, MARCO, and APOE, and was negatively associated with clinical outcome, suggesting a pro-tumor role is this TAM population (Lavin et al., 2017). A separate scRNA-seq study on lung cancer TME identified three macrophage subpopulations strongly associated with poor clinical outcomes, and each of the three subsets expressed a distinct chemokine, presumably interacting with distinct cell types expressing the corresponding receptors within the microenvironment (Zilionis et al., 2019). A third study revealed that TAMs infiltrated in early stage lung cancer TME were highly plastic, co-expressing co-inhibitory and co-stimulatory receptors, suggesting their potential to both stimulate and inhibit T cell functions (Singhal et al., 2019). Although the widespread expression of co-inhibitory ligand PD-L1 in TAMs is thought to involve in the immunosuppression of T cells, two recent studies on early-stage lung cancer suggested otherwise. Singhal et al. (2019) showed that unlike PD-L1 expressing tumor cells, PD-L1+ TAMs did not directly suppress T cell functions in vitro, but rather protected themselves from T cell mediated cytotoxicity. On the other hand, PD-L1 was also highly expressed in tissue-resident macrophages, and these PD-L1+ macrophages accumulated in tight clusters at the tumor invasive margin and negatively correlated with T cell infiltration in tumor (Lavin et al., 2017). Two studies carried out on human breast cancer by scRNA-seq and CyToF analyses identified TAMs with heterogenous expression profiles, phenotypically differed from tissue-resident macrophages in juxta-tumoral regions (Azizi et al., 2018; Wagner et al., 2019). PD-L1 expression was common in these breast cancer infiltrated TAMs, and the presence of PD-L1+ TAMs was associated with high grade tumor and poor prognosis, suggesting their role in promoting breast cancer progression (Azizi et al., 2018; Wagner et al., 2019). A scRNA-seq study carried out on human liver cancer identified two distinct states of tumor-enriched macrophages, one of which expressing TAM-like signature genes such as APOE, TREM2, GPNMB, and SLC40A1, and the other expressing MDSC-like signatures. Only the TAM-like population was associated with poor prognosis in liver cancer, probably through production of pro-inflammatory cytokines such as TNFα and IL-6 (Zhang et al., 2019). Another recent scRNA-seq study on colorectal cancer also identified two distinct TAM subsets: the C1QC+ TAMs primarily involved in complement activation and antigen presentation, and the SPP1+ TAMs primarily involved in tumor angiogenesis and ECM interaction (Zhang et al., 2020). Intriguingly, the two TAM subtypes exhibited differential sensitivity toward CSF1R blockade; and the resistance of the SPP1+ TAMs against anti-CSF1R treatment may explain the minimal therapeutic benefit observed following anti-CSF1R monotherapy (Zhang et al., 2020). Lastly, a scRNA-seq carried out on human glioma identified two TAM populations with distinct lineages, one corresponding to blood-derived TAMs and the other corresponding to tissue-resident microglial TAMs; and only the blood-derived TAMs was associated with poor clinical outcomes (Müller et al., 2017).
Figure 3. Intratumoral differentiation of the monocyte/macrophage lineage. Three conserved monocyte/macrophage differentiation trajectories are identified from human TME. The first trajectory is the tumor infiltration and activation of conventional monocytes (CD14+CD16–) from the vasculature into activated intratumoral monocytes; the second trajectory is the ongoing intratumoral macrophage differentiation from the activated intratumoral monocytes into tumor-associated macrophages (TAMs), characterized by high expression of HLA-DR, CD68 and CD64; the third trajectory is the differentiation of circulating conventional monocytes (CD14+CD16–) into early immigrant macrophages (HLA-DRintCD192+) and then continue to tissue-resident macrophages (HLA-DRintCD206+) through stepwise tissue adaptation within the adjacent normal tissue. Tissue-resident macrophages (HLA-DRintCD206+) can then across the invasive margin to infiltrate tumor and potentially be converted to TAMs (HLA-DRhiCD68+CD64+). The intratumoral monocyte may also differentiate into myeloid-derived suppressor cells (MDSCs), whose presence often negatively associates with clinical outcomes.
Single-cell sequencing is emerging as an important research tool to directly study the human cancer infiltrated myeloid cell compartment, which heavily relied on mouse models and pre-defined surface markers in the past. These single-cell studies have revealed the expression heterogeneity and functional plasticity of myeloid lineage cells and highlighted their complex crosstalk with lymphocytes. Dendritic cells are crucial for T-cell mediated anti-tumor immunity: they function as the primary tumor antigen-presenting cells for de novo priming of tumor-reactive T cells; meanwhile, they secrete a rich array of chemokines, cytokines, and co-stimulatory molecules to modulate T cell functions and their response to ICB therapy (Roberts et al., 2016; Salmon et al., 2016; Garris et al., 2018; Chow et al., 2019). On the other hand, DCs are subjected to intensive regulation by lymphocytes, such as NK cells and TREG cells (Barry et al., 2018; Bottcher et al., 2018; Binnewies et al., 2019). The identification of the “activated” state of cDCs within the human TME by scRNA-seq further validated the capacity of single-cell sequencing in recovering rare transiting cell state (Zhang et al., 2019; Zilionis et al., 2019). Tumor-associated macrophages exhibited the highest cross-species discrepancy among the myeloid lineage, reflected by their divergence from the classical polarization model, and their unique population structure (Müller et al., 2017; Azizi et al., 2018; Zhang et al., 2019; Zilionis et al., 2019). Therefore, understanding of the TAM functions and their clinical relevance to immunotherapy heavily rely on analysis of human samples, attributed by the advances in single-cell sequencing technologies. Nevertheless, current studies on TIMs relying solely on scRNA-seq present several limitations: the lack of spatial resolution limits our understanding on how the spatial organization of immune and stromal cells may impact on immune cell functions and tumor cell invasion; the lack of tracible genetic variations limits the confidence on inferring the developmental lineages of infiltrating myeloid cells; and the lack of direct cell-to-cell interaction information limits our in depth understanding on intercellular crosstalk that modulates the cell localization and functions.
Tumor Cell Intrinsic Programs and Immunoediting
In addition to immune and stromal cell populations, the tumor cells also play direct and critical roles on regulating immune cell infiltration and their anti-tumor functions (Pitt et al., 2016; Nguyen and Spranger, 2020). Meanwhile, the tumor cells are constantly under selection pressure imposed by the host immune system, and thereby are driven to acquire somatic mutations enabling their escape from immune surveillance, a phenomenon termed as “immunoediting” (Mittal et al., 2014; Galon and Bruni, 2020). Therefore, the interplay between an evolving cancer and a dynamic immune microenvironment is one of the major contributors of tumor evolution. Through scRNA-seq analysis on melanoma patients pre- and post-ICB treatment, Jerby-Arnon et al. (2018) identified an ICB resistance program-upregulation of the CDK4/6 pathway by melanoma tumor cells-that led to T cell exclusion and immune evasion. The authors further demonstrated that ICB therapy shaped the tumor cell expression toward the resistance program, suggestive of immunoediting on tumor cells driven by the restored anti-tumor immunity following ICB treatment. Using similar approach, Puram et al. (2017) identified a partial epithelial-to-mesenchymal transition (p-EMT) program expressed by HNSCC malignant cells localized to the leading edge of primary tumor, driving efficient metastasis. These studies generated novel predictive biomarkers for clinical outcomes and therapeutic responses, highlighting the clinical utility of single-cell sequencing on TME profiling. Tumor-derived chemokines were suggested as another tumor-intrinsic mechanism to modulate the TME: epigenetic silencing of CXCL9/10 was shown to mediate T cell exclusion in mouse ovarian cancer model; and secretion of CXCL1 was reported to recruit myeloid cells and deplete T cells in mouse pancreatic ductal adenocarcinoma (PDAC) model (Peng et al., 2015; Li J. et al., 2018). On the other hand, a series of studies by TRACERx consortium reported a strong correlation between clonal neoantigen burden with T cell infiltration and corresponding intratumoral TCR expansion, suggestive of intensive interplay between the tumor-responsive T cells and their tumor cell targets (McGranahan et al., 2016; Joshi et al., 2019; Rosenthal et al., 2019). Infiltrated T cells exerted strong selection pressure on tumor cells, driving their ongoing immunoediting to achieve immune escape. The major mechanisms include: (1) disruption of the antigen presentation machinery by HLA loss of heterozygosity (LOH) or B2M loss-of-function mutations; (2) selection for neoantigen depletion by DNA copy number loss, transcriptional repression or epigenetic repression; and (3) mutations in IFN and IL-2 signaling components (Rooney et al., 2015; McGranahan et al., 2016; Joshi et al., 2019; Rosenthal et al., 2019). The aforementioned studies were mainly based on transcriptional profiles of bulk tumor samples, with limited resolution to dissect the cell type specific cancer-immune interplay. Conceivably, the future application of single-cell technologies on this topic will greatly facilitate us to further understand the tumor-intrinsic mechanisms that lead to immune escape and their respective impacts on clinical outcomes.
Implications of Different Immune Cell Types on Immunotherapy
It has been generally recognized that the abundance of CD8+ cytotoxic T cells in TME is associated with favorable clinical outcomes in various cancer types (Fridman et al., 2017); and corresponding biomarker utilizing the CD8+/CD3+ T cell density ratio has been validated to predict the risk of recurrence for colorectal cancer patients, independent from TNM staging information (Pagès et al., 2018). However, a robust immune-based biomarker, independent from tumor-based PD-L1 level and TMB, has yet been validated in predicting treatment response of immunotherapy. It underscores the importance of identifying immunotherapy-responsive immune cell subsets from the local TME, as well as from the secondary lymphoid organs and peripheral blood. A “progenitor-like” intratumoral CD8+ T cell subset that may contribute to the durable response to ICB treatment has been identified in a series of mouse model studies: this Tcf1+ CD8+ T cell subset is featured by its polyfunctionality to proliferate and differentiate into more cytolytic “terminally exhausted” CD8+ T cells in response to PD-1 blockade or TCR stimulation (Im et al., 2016; Utzschneider et al., 2016; Wu et al., 2016; Snell et al., 2018; Chen et al., 2019; Miller et al., 2019; Siddiqui et al., 2019; Yao et al., 2019). TCF1+ TILs and CXCR5+ TILs have also been observed in human cancer TME; however, whether the human TCF1+ or CXCR5+ TILs resemble their mouse counterparts remains highly controversial (Brummelman et al., 2018; Sade-Feldman et al., 2018). More interestingly, two mouse model studies utilizing scRNA-seq together with bulk ATAC-seq identified distinct epigenetic signatures between the Tcf1+ “progenitor exhausted” and the Tcf1– “terminally exhausted” T cells (Jadhav et al., 2019; Miller et al., 2019). These studies also demonstrated that the epigenetic signatures of different T cell subsets remained stable after PD-1 blockade, underscoring the importance of epigenetic state in regulating T cell fate and function (Jadhav et al., 2019; Miller et al., 2019).
Epigenetic mechanisms, acting in conjunction with transcriptional factors, play a pivotal role in regulation of immune cell differentiation and function. Most importantly, unlike transcriptional profiles that can be reversed upon PD-1 blockade, the epigenetic state of exhausted T cells remained stable, responsible for the failed memory development of the transiently rejuvenated exhausted T cells and their fast re-exhaustion (Youngblood et al., 2011; Ladle et al., 2016; Pauken et al., 2016; Sen et al., 2016; Ghoneim et al., 2017; Mognol et al., 2017). Therefore, the epigenetic imprinting on exhausted or dysfunctional T cells appears to be a major roadblock preventing them from sustainable rejuvenation by ICB treatment. Consistently, distinct epigenetic states corresponded to distinct dysfunctional states of tumor-specific CD8+ T cells in tumor mouse model: the “early dysfunctional” TILs obtained a plastic chromatin state that is reprogrammable, whereas the “late dysfunctional” TILs obtained a fixed chromatin state resistant to reprogramming (Schietinger et al., 2016; Philip et al., 2017). This mirrored observations on human TME that tumor-reactive CD8+ TILs exhibited various degrees of dysfunction at the transcriptome level, and thereby predicts that the degree of dysfunction might be further sculpted by the epigenetic state which has direct impact on the rejuvenation capacity of exhausted T cells in response to ICB treatment.
In addition to T cells, several recent studies have demonstrated that PD-L1 blockade reinvigorated DC function to enhance T cell priming, and thereby generated potent anticancer T cell immunity (Mayoux et al., 2020; Oh S.A. et al., 2020). A DC gene signature is associated with improved overall survival in patients with renal cancer and NSCLC treated with PD-L1 blockade (Mayoux et al., 2020). Consistent to this view, two recent scRNA-seq studies have observed the expansion of anti-tumor TCR reperitore in the peripheral blood and clonal replacement of intratumoral tumor-specific T cells following PD-1/PD-L1 blockade (Yost et al., 2019; Wu et al., 2020), suggesting the enhanced T cell priming of tumor neoantigens as a major mechanism underlying effective immunotherapy. Altogether, these results have demonstrated the dual-effect of PD-L1/PD-1 blockade on (1) enhancing T cell priming via activation or rejuvantation of DCs and (2) triggering proliferation and differentiation of the intratumoral “progenitor-like” T cells. It predicts that the presence of activated DC subsets and/or proliferative intratumoral T cell subsets would positively associate with the clinical outcome of immunotherapy; and developing biomarkers representing these specific immune cell subsets may yield promising predictive value on immunotherapy response.
Prospects on Applications of Single-Cell Multi-Omics on TME Studies
Recent applications of scRNA-seq has significantly advanced our understanding on the functional diversity of immune cells in human cancer TME. Moreover, an ongoing activation and differentiation trajectory within the TME was identified from multiple immune cell lineages, such as T cells and macrophages, highlighting the complex crosstalk between the microenvironment and immune cells to regulate their differentiation and function. Single-cell transcriptomics generate a static snapshot of the transcriptional phenotype of individual cells at a single time point. Computational algorithms can be used to further project all sampled cells into a differentiation trajectory based on the assumptions that cells spanning a continuum of transitional states are sufficiently sampled in the dataset, and cells with similar transcriptional profiles are developmentally related (Kester and van Oudenaarden, 2018; Baron and van Oudenaarden, 2019; Tritschler et al., 2019). The most commonly used computational algorithms in TME studies include: (1) Monocle/monocle 2, which projects a minimum spanning tree that connects cells with similar transcriptional profiles and then constructs a “pseudotime” that serves as the backbone of the predicted lineage trajectory (Trapnell et al., 2014; Qiu et al., 2017); (2) and RNA velocity, which predicts the future state of a cell based on the fraction of spliced and unspliced transcripts (La Manno et al., 2018). However, the differentiation trajectories predicted by single-cell transcriptomics data are purely phenotypic, not necessarily reflecting the true genetic relationship between lineages of cells (Kester and van Oudenaarden, 2018; Baron and van Oudenaarden, 2019; Tritschler et al., 2019). Moreover, several key aspects of cellular identity, including but not limited to the epigenetic state, protein profile and spatial location, are missing from the scRNA-seq derived maps. Therefore, moving from scRNA-seq snapshots to multimodal measurements and integration of the genome, epigenome, transcriptome, proteome, and spatial organization datasets will further extend the power of single-cell genomics in immunology studies (Ren et al., 2018; Nathan et al., 2019; Efremova and Teichmann, 2020; Schier, 2020; Zhu et al., 2020).
Growing efforts on developing more sophisticated computational platforms to integrate multimodal single-cell datasets led to more accurate definition of cell identity and state with multiple layers of information. These efforts will ultimately build access into the gene regulatory network that shapes the phenotype and behavior of a pure cell population within complex biological systems (Stuart et al., 2019; Welch et al., 2019). Emerging techniques that allow for simultaneous assessment of multi-omics information from the same cell have also been developed and continue to evolve for greater scalability and higher genomics coverage in the recent years (Table 2; Dixit et al., 2016; Hou et al., 2016; Guo et al., 2017; Peterson et al., 2017; Pott, 2017; Stoeckius et al., 2017; Alemany et al., 2018; Bian et al., 2018; Biddy et al., 2018; Cao et al., 2018; Chen et al., 2018; Clark et al., 2018; Li L. et al., 2018; Satpathy et al., 2018; Spanjaard et al., 2018; Argelaguet et al., 2019; Li G. et al., 2019; Liu et al., 2019; Ludwig et al., 2019; Mimitou et al., 2019; Rodriques et al., 2019; Rooijers et al., 2019; Rubin et al., 2019; Xu et al., 2019; Zhu et al., 2019; Weinreb et al., 2020). For instance, joint analysis of TCR repertoire and chromatin accessibly at the single-cell level was achieved by a novel method named as transcript-indexed ATAC-seq (T-ATAC-seq; Satpathy et al., 2018). By profiling of human peripheral T cells, T-ATAC-seq identified cis- and trans- regulators of naïve and memory CD4+ T cell states and substantial epigenomic heterogeneity within the surface-marker-defined T cell populations (Satpathy et al., 2018). This method enables analysis of epigenetic state of clonal T cells, which will be particularly useful for lineage tracing of clonal T cells undergoing differentiation to understand the epigenetic regulation of T cell differentiation and memory formation in response to different environmental stimulations.
The joint single-cell analysis of TCR repertoire and transcriptome has yielded tremendous insight into T cell differentiation process by connecting the genetic lineage of a given T cell reflected by its highly polymorphic TCR sequence to its functional phenotype defined by the transcriptome. A landmark study applied this approach to delineate the complete T cell differentiation roadmap of human, from early progenitors residing in the hematopoietic fetal liver and thymus into fully matured T cell types, via sampling of developing human thymus ranged from the early embryonic stage to the adulthood (Park et al., 2020). In the context of human TME, similar analyses have unveiled that a shared TCR clonotype can span a wide range of T cell differentiation states, supporting that these cell states are in fact developmentally linked (Guo et al., 2018; Zhang et al., 2018; Li H. et al., 2019). Conversely, certain TCR clonotypes are almost restricted to a specific cell subset, suggestive of its distinct developmental lineage (Guo et al., 2018; Zhang et al., 2018; Li H. et al., 2019; Yost et al., 2019). Despite of these important observations, our direct knowledge on the developmental lineages of immune cells in human TME remains very limited. For instance, the origin of intratumoral TREG cells remains elusive since vast majority of their TCR clonotypes are exclusive to themselves, with only a small fraction of TCR sequences shared with TREG cells from blood or adjacent normal tissue, or with intratumoral TH cells (Guo et al., 2018; Zhang et al., 2018). Therefore, it remains largely unknown whether intratumoral TREG cells are primarily derived from blood TREG, tissue-resident TREG or intratumoral TH (iTREG), and whether intratumoral TREG cells derived from different origins exert distinct functions. Another important but unanswered question is the cell origin of the expanded T cell population in response to PD-1 blockade. Two groups have independently observed the “clonotype replacement” phenomenon following responsive anti-PD-1 treatment in skin cancer patients: clonotypes expanded after PD-1 blockade were different from the pre-existing clonotypes identified from the pre-treatment TME, suggesting limited reinvigoration capacity of the majority pre-existing T cells in the TME (Sade-Feldman et al., 2018; Yost et al., 2019). However, it remains undetermined whether these expanded T cell clones are originated from newly primed T cells outside of the TME, or from clonal expansion and differentiation of pre-existing “progenitor-like” T cells that are previously undetectable due to their rarity.
Prospective lineage tracing enabled by simultaneous profiling of single-cell transcriptome and the unique DNA barcodes introduced by genetic manipulation allows for tracking the clonal dynamics of different cell lineages by coupling the transcriptionally defined cell differentiation states with the perspective clonal identifiers that record the past history of cells (Biddy et al., 2018; Spanjaard et al., 2018; Weinreb et al., 2020). Weinreb et al. (2020) applied this approach to study the fate determination in hematopoiesis and precisely identified a continuous spectrum of cell states with primed fate potential, demonstrating the superior performance of this multimodal system compared to the previously used single-modal scRNA-seq only or clonal barcoding only systems. Moreover, this study has drawn two important conclusions that are instrumental to future lineage tracing studies: (1) matured cells derived from distinct lineages are differentially imprinted and thereby phenotypically distinct; (2) transcriptome alone is not sufficient to define the fate potential of progenitor cells, suggesting other heritable properties are missing from the current system, presumably including the epigenetic states, protein abundances, cell organizations, and the microenvironment (Weinreb et al., 2020). Despite of its robustness, prospective lineage tracing can be only applied to in vitro and animal model systems due to the requirement of genetic manipulation. Instead, retrospective lineage tracing, taking advantage of the naturally occurring somatic mutations during development, has been widely used for analysis of human samples (Woodworth et al., 2017). A variety of inheritable genetic alternations have been used as the natural “DNA barcodes”, such as TCR/BCR sequences (Han et al., 2014; Stubbington et al., 2016), copy number variations (CNVs; McConnell et al., 2013; Cai et al., 2014), single nucleotide variations (SNVs; Lodato et al., 2015), retrotransposon elements (i.e., LINE-1; Evrony et al., 2012), microsatellite repeats (Evrony et al., 2012, 2015) and mitochondrial DNA (mtDNA; Ludwig et al., 2019; Xu et al., 2019). Lineage tracing based on mtDNA mutations has been recently applied to study the lineage of macrophages in human liver cancer TME, and generated consistent results with independent RNA velocity-based prediction (Zhang et al., 2019). Conceivably, with the emergence and further adoption of single-cell multimodal sequencing technologies, lineage tracing analysis on immune cell populations from mouse and human TME will generate important insight into the regulatory mechanisms underlying immune cell differentiation, exhaustion and memory formation in response to local TME stimuli and various types of immunotherapy.
Current Limitations and Challenges of Single-Cell Multiomic Technologies
Despite of its great potential, the current single-cell technologies suffer from their limited scalability, sparse coverage, allelic dropouts, PCR errors, and most importantly, lack of spatial information (Nam et al., 2020). For instance, the current scRNA-seq technologies typically generate hundreds to thousands of cells per biological sample, and thousands to hundred thousands of cells per study. Compared to FACS and CyToF analyses, which typically generate data from millions of cells, accurate quantification of cell subsets, especially for rare subsets, remains a challenge for current single-cell technologies. Fortunately, the recent development of combinatorial indexing-based single-cell technology holds promise to further extent the scalability of current single-cell technologies by 1 to 2 orders of magnitude without increase in cost (Cao et al., 2017, 2018; Rosenberg et al., 2018; Zhu et al., 2019), and thereby approaching the throughput of conventional FACS and CyToF methods. Furthermore, since current single-cell sequencing technologies mainly rely on front-end molecular amplification of the picogram scale genomic content from individual cells, they inevitably suffer from PCR errors and allelic dropouts. Meanwhile, due to the limited genomic amplification efficiency and high sequencing cost, the genomic coverage of each individual cell is often very sparse, resulting technological challenges in data processing and interpretation. Therefore, specialized bioinformatic algorithms are necessary to compensate the errors and missing information associated with the typical single-cell sequencing datasets, limiting the general accessibility of these technologies to the greater research community. Last but most importantly, the current single-cell technologies often require front-end tissue dissociation, which inevitably destroys the spatial architecture of the biological specimens, eliminating a critical layer of information that contributes to the biological identity of a cell. To address this limitation, spatial sequencing techniques, such as Slide-Seq (Rodriques et al., 2019), have been developed; and the data integration of spatial transcriptomics is hopeful to help dissect important cell interactions at single-cell resolution in the near future (Adey, 2019; Welch et al., 2019; Nam et al., 2020).
Recent technological advancements in single-cell genomics have greatly impacted the way immunologist to conduct their research, particularly in the immune-oncology field. Due to the highly heterogeneous nature of the immune system, single-cell based analysis such as FACS has always been the gold standard for immunologists to study the phenotypes of immune cells. However, conventional FACS analysis is limited by the number of markers that can be simultaneously analyzed and is largely restricted to analysis of surface markers. Although the development of CyToF largely increased the number of markers to be assessed in parallel, it remains constrained by the availability of known surface markers, which limits the discovery of previously unknown cell states. The emergence of scRNA-seq has revolutionized the way immunologists to “phenotype” their cellular population of interest: instead of assigning a cell into a particular cell type positive or negative for a set of pre-defined markers, scRNA-seq places a cell into a phenotypic continuum reflecting various aspects of its cellular identity including activation, differentiation, metabolic states and more. This leads to the identification of immune cells with much more functional diversity than previously appreciated. Moreover, scRNA-seq unveiled that immune cells differentiate along a continuous trajectory attributed by intrinsic regulatory network and extrinsic environmental stimuli, instead of differentiating step-wisely into discrete cell types or intermediates defined by selected markers. Additionally, the development of single-cell multimodal sequencing enables simultaneous profiling of genome, epigenome, transcriptome, proteome and spatial localization from the same cell (Guo et al., 2017; Peterson et al., 2017; Stoeckius et al., 2017; Chen et al., 2018; Clark et al., 2018; Li L. et al., 2018; Argelaguet et al., 2019; Mimitou et al., 2019). Combining genome, epigenome and transcriptome from the same cell allows for retrospective lineage tracing directly from human samples, which is critical for the understanding of immune cell differentiation, activation and exhaustion in response to TME stimuli. Combining epigenome, transcriptome and proteome allows for more accurately defining the cell identity from multiple layers to fully elucidate the functional plasticity of a cell, as well as its past history and future potential. Addition of the spatial localization to the molecular profiles further complements the definition of cell identity, which is also regulated by spatial position for its normal function. Single-cell multi-omics, in a spatially resolved context, is important for understanding the interactions between cells of the immune system and crosstalk between tumor and other cell types in the TME. For instance, the TCR spatial heterogeneity, directed by the physical interaction between TCR and its targeting neoantigen, reflects genomic intratumoral heterogeneity (Joshi et al., 2019); and the spatial distribution of intratumoral PD-L1+ macrophages impacts T cell infiltration (Lavin et al., 2017). Lastly, the integration of large-scale datasets across platforms, omics and species, and combination with relevant functional and clinical information, will ultimately transform our understanding of human anti-tumor immunity. This will hopefully allow for improved patient stratifications, biomarker discovery and druggable targets identification, leading to the achievement of precision immunotherapy with better efficacy and less toxicity (Giladi and Amit, 2018).
XC conceived the project. TG and XC wrote the manuscript with input from WL. All authors read and approved the final manuscript.
This work was supported by National Natural Science Foundation of China (81972911 to XC and 81871890 and 91859203 to WL), CAMS Innovation Fund for Medical Science (2019TX310002), and National Guided Science and Technology Development Project of Sichuan Province (2020ZYD009).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Ahmadzadeh, M., Pasetto, A., Jia, L., Deniger, D. C., Stevanović, S., Robbins, P. F., et al. (2019). Tumor-infiltrating human CD4+ regulatory T cells display a distinct TCR repertoire and exhibit tumor and neoantigen reactivity. Sci. Immunol. 4:eaao4310. doi: 10.1126/sciimmunol.aao4310
Argelaguet, R., Clark, S. J., Mohammed, H., Stapel, L. C., Krueger, C., Kapourani, C.-A., et al. (2019). Multi-omics profiling of mouse gastrulation at single-cell resolution. Nature 576, 487–491. doi: 10.1038/s41586-019-1825-8
Azizi, E., Carr, A. J., Plitas, G., Cornish, A. E., Konopacki, C., Prabhakaran, S., et al. (2018). Single-cell map of diverse immune phenotypes in the breast tumor microenvironment. Cell 174, 1293–1308.e1236. doi: 10.1016/j.cell.2018.05.060
Baron, C. S., and van Oudenaarden, A. (2019). Unravelling cellular relationships during development and regeneration using genetic lineage tracing. Nat. Rev. Mol. Cell Biol. 20, 753–765. doi: 10.1038/s41580-019-0186-3
Barry, K. C., Hsu, J., Broz, M. L., Cueto, F. J., Binnewies, M., Combes, A. J., et al. (2018). A natural killer-dendritic cell axis defines checkpoint therapy-responsive tumor microenvironments. Nat. Med. 24, 1178–1191. doi: 10.1038/s41591-018-0085-8
Biddy, B. A., Kong, W., Kamimoto, K., Guo, C., Waye, S. E., Sun, T., et al. (2018). Single-cell mapping of lineage and identity in direct reprogramming. Nature 564, 219–224. doi: 10.1038/s41586-018-0744-4
Binnewies, M., Mujal, A. M., Pollack, J. L., Combes, A. J., Hardison, E. A., Barry, K. C., et al. (2019). Unleashing type-2 dendritic cells to drive protective antitumor CD4(+) T cell immunity. Cell 177, 556–571.e516. doi: 10.1016/j.cell.2019.02.005
Binnewies, M., Roberts, E. W., Kersten, K., Chan, V., Fearon, D. F., Merad, M., et al. (2018). Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat. Med. 24, 541–550. doi: 10.1038/s41591-018-0014-x
Borst, J., Ahrends, T., Ba̧bała, N., Melief, C. J. M., and Kastenmüller, W. (2018). CD4+ T cell help in cancer immunology and immunotherapy. Nat. Rev. Immunol. 18, 635–647. doi: 10.1038/s41577-018-0044-0
Bottcher, J. P., Bonavita, E., Chakravarty, P., Blees, H., Cabeza-Cabrerizo, M., Sammicheli, S., et al. (2018). NK cells stimulate recruitment of cDC1 into the tumor microenvironment promoting cancer immune control. Cell 172, 1022–1037.e1014. doi: 10.1016/j.cell.2018.01.004
Brummelman, J., Mazza, E. M. C., Alvisi, G., Colombo, F. S., Grilli, A., Mikulak, J., et al. (2018). High-dimensional single cell analysis identifies stem-like cytotoxic CD8(+) T cells infiltrating human tumors. J. Exp. Med. 215, 2520–2535. doi: 10.1084/jem.20180684
Cabrita, R., Lauss, M., Sanna, A., Donia, M., Skaarup Larsen, M., Mitra, S., et al. (2020). Tertiary lymphoid structures improve immunotherapy and survival in melanoma. Nature 577, 561–565. doi: 10.1038/s41586-019-1914-8
Cai, X., Evrony, G. D., Lehmann, H. S., Elhosary, P. C., Mehta, B. K., Poduri, A., et al. (2014). Single-cell, genome-wide sequencing identifies clonal somatic copy-number variation in the human brain. Cell Reports 8, 1280–1289. doi: 10.1016/j.celrep.2014.07.043
Cao, J., Cusanovich, D. A., Ramani, V., Aghamirzaie, D., Pliner, H. A., Hill, A. J., et al. (2018). Joint profiling of chromatin accessibility and gene expression in thousands of single cells. Science 361, 1380–1385. doi: 10.1126/science.aau0730
Cao, J., Packer, J. S., Ramani, V., Cusanovich, D. A., Huynh, C., Daza, R., et al. (2017). Comprehensive single-cell transcriptional profiling of a multicellular organism. Science 357, 661–667. doi: 10.1126/science.aam8940
Chen, X., Litzenburger, U. M., Wei, Y., Schep, A. N., LaGory, E. L., Choudhry, H., et al. (2018). Joint single-cell DNA accessibility and protein epitope profiling reveals environmental regulation of epigenomic heterogeneity. Nat. Commun. 9:4590. doi: 10.1038/s41467-018-07115-y
Chen, Z., Ji, Z., Ngiow, S. F., Manne, S., Cai, Z., Huang, A. C., et al. (2019). TCF-1-centered transcriptional network drives an effector versus exhausted CD8 T cell-fate decision. Immunity 51, 840–855.e845. doi: 10.1016/j.immuni.2019.09.013
Chevrier, S., Levine, J. H., Zanotelli, V. R. T., Silina, K., Schulz, D., Bacac, M., et al. (2017). An immune atlas of clear cell renal cell carcinoma. Cell 169, 736–749.e718. doi: 10.1016/j.cell.2017.04.016
Chow, M. T., Ozga, A. J., Servis, R. L., Frederick, D. T., Lo, J. A., Fisher, D. E., et al. (2019). Intratumoral activity of the CXCR3 chemokine system is required for the efficacy of anti-PD-1 therapy. Immunity 50, 1498–1512.e1495. doi: 10.1016/j.immuni.2019.04.010
Cillo, A. R., Kürten, C. H. L., Tabib, T., Qi, Z., Onkar, S., Wang, T., et al. (2020). Immune landscape of viral- and carcinogen-driven head and neck cancer. Immunity 52, 183–199.e189. doi: 10.1016/j.immuni.2019.11.014
Clark, S. J., Argelaguet, R., Kapourani, C. A., Stubbs, T. M., Lee, H. J., Alda-Catalinas, C., et al. (2018). scNMT-seq enables joint profiling of chromatin accessibility DNA methylation and transcription in single cells. Nat. Commun. 9:781. doi: 10.1038/s41467-018-03149-4
Clarke, J., Panwar, B., Madrigal, A., Singh, D., Gujar, R., Wood, O., et al. (2019). Single-cell transcriptomic analysis of tissue-resident memory T cells in human lung cancer. J. Exp. Med. 216, 2128–2149. doi: 10.1084/jem.20190249
De Simone, M., Arrigoni, A., Rossetti, G., Gruarin, P., Ranzani, V., Politano, C., et al. (2016). Transcriptional landscape of human tissue lymphocytes unveils uniqueness of tumor-infiltrating T regulatory cells. Immunity 45, 1135–1147. doi: 10.1016/j.immuni.2016.10.021
Dixit, A., Parnas, O., Li, B., Chen, J., Fulco, C. P., Jerby-Arnon, L., et al. (2016). Perturb-Seq: dissecting molecular circuits with scalable single-cell RNA profiling of pooled genetic screens. Cell 167, 1853–1866.e1817. doi: 10.1016/j.cell.2016.11.038
Duhen, T., Duhen, R., Montler, R., Moses, J., Moudgil, T., de Miranda, N. F., et al. (2018). Co-expression of CD39 and CD103 identifies tumor-reactive CD8 T cells in human solid tumors. Nat. Commun. 9:2724. doi: 10.1038/s41467-018-05072-0
Evrony, G. D., Cai, X., Lee, E., Hills, L. B., Elhosary, P. C., Lehmann, H. S., et al. (2012). Single-neuron sequencing analysis of L1 retrotransposition and somatic mutation in the human brain. Cell 151, 483–496. doi: 10.1016/j.cell.2012.09.035
Evrony, G. D., Lee, E., Mehta, B. K., Benjamini, Y., Johnson, R. M., Cai, X., et al. (2015). Cell lineage analysis in human brain using endogenous retroelements. Neuron 85, 49–59. doi: 10.1016/j.neuron.2014.12.028
Garon, E. B., Rizvi, N. A., Hui, R., Leighl, N., Balmanoukian, A. S., Eder, J. P., et al. (2015). Pembrolizumab for the treatment of non–small-cell lung cancer. N. Engl. J. Med. 372, 2018–2028. doi: 10.1056/NEJMoa1501824
Garris, C. S., Arlauckas, S. P., Kohler, R. H., Trefny, M. P., Garren, S., Piot, C., et al. (2018). Successful anti-PD-1 cancer immunotherapy requires T cell-dendritic cell crosstalk involving the cytokines IFN-gamma and IL-12. Immunity 49, 1148–1161.e1147. doi: 10.1016/j.immuni.2018.09.024
Ghoneim, H. E., Fan, Y., Moustaki, A., Abdelsamed, H. A., Dash, P., Dogra, P., et al. (2017). De novo epigenetic programs inhibit PD-1 blockade-mediated T cell rejuvenation. Cell 170, 142–157.e119. doi: 10.1016/j.cell.2017.06.007
Guo, X., Zhang, Y., Zheng, L., Zheng, C., Song, J., Zhang, Q., et al. (2018). Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing. Nat. Med. 24, 978–985. doi: 10.1038/s41591-018-0045-3
Helmink, B. A., Reddy, S. M., Gao, J., Zhang, S., Basar, R., Thakur, R., et al. (2020). B cells and tertiary lymphoid structures promote immunotherapy response. Nature 577, 549–555. doi: 10.1038/s41586-019-1922-8
Hodi, F. S., O’Day, S. J., McDermott, D. F., Weber, R. W., Sosman, J. A., Haanen, J. B., et al. (2010). Improved survival with ipilimumab in patients with metastatic melanoma. N. Engl. J. Med. 363, 711–723. doi: 10.1056/NEJMoa1003466
Hou, Y., Guo, H., Cao, C., Li, X., Hu, B., Zhu, P., et al. (2016). Single-cell triple omics sequencing reveals genetic, epigenetic, and transcriptomic heterogeneity in hepatocellular carcinomas. Cell Res. 26, 304–319. doi: 10.1038/cr.2016.23
Im, S. J., Hashimoto, M., Gerner, M. Y., Lee, J., Kissick, H. T., Burger, M. C., et al. (2016). Defining CD8+ T cells that provide the proliferative burst after PD-1 therapy. Nature 537, 417–421. doi: 10.1038/nature19330
Jadhav, R. R., Im, S. J., Hu, B., Hashimoto, M., Li, P., Lin, J. X., et al. (2019). Epigenetic signature of PD-1+ TCF1+ CD8 T cells that act as resource cells during chronic viral infection and respond to PD-1 blockade. Proc. Natl. Acad. Sci. U.S.A. 116, 14113–14118. doi: 10.1073/pnas.1903520116
Jerby-Arnon, L., Shah, P., Cuoco, M. S., Rodman, C., Su, M. J., Melms, J. C., et al. (2018). A cancer cell program promotes T cell exclusion and resistance to checkpoint blockade. Cell 175, 984–997.e924. doi: 10.1016/j.cell.2018.09.006
Joshi, K., Robert de Massy, M., Ismail, M., Reading, J. L., Uddin, I., et al. (2019). Spatial heterogeneity of the T cell receptor repertoire reflects the mutational landscape in lung cancer. Nat. Med. 25, 1549–1559. doi: 10.1038/s41591-019-0592-2
Keeley, T., Costanzo-Garvey, D. L., and Cook, L. M. (2019). Unmasking the many faces of tumor-associated neutrophils and macrophages: considerations for targeting innate immune cells in cancer. Trends Cancer 5, 789–798. doi: 10.1016/j.trecan.2019.10.013
Ladle, B. H., Li, K. P., Phillips, M. J., Pucsek, A. B., Haile, A., Powell, J. D., et al. (2016). De novo DNA methylation by DNA methyltransferase 3a controls early effector CD8+ T-cell fate decisions following activation. Proc. Natl. Acad. Sci. U.S.A. 113, 10631–10636. doi: 10.1073/pnas.1524490113
Lambrechts, D., Wauters, E., Boeckx, B., Aibar, S., Nittner, D., Burton, O., et al. (2018). Phenotype molding of stromal cells in the lung tumor microenvironment. Nat. Med. 24, 1277–1289. doi: 10.1038/s41591-018-0096-5
Lavin, Y., Kobayashi, S., Leader, A., Amir, E.-A. D., Elefant, N., Bigenwald, C., et al. (2017). Innate immune landscape in early lung adenocarcinoma by paired single-cell analyses. Cell 169, 750–765.e717. doi: 10.1016/j.cell.2017.04.014
Li, G., Liu, Y., Zhang, Y., Kubo, N., Yu, M., Fang, R., et al. (2019). Joint profiling of DNA methylation and chromatin architecture in single cells. Nat. Methods 16, 991–993. doi: 10.1038/s41592-019-0502-z
Li, H., van der Leun, A. M., Yofe, I., Lubling, Y., Gelbard-Solodkin, D., van Akkooi, A. C. J., et al. (2019). Dysfunctional CD8 T cells form a proliferative, dynamically regulated compartment within human melanoma. Cell 176, 775–789.e718. doi: 10.1016/j.cell.2018.11.043
Li, J., Byrne, K. T., Yan, F., Yamazoe, T., Chen, Z., Baslan, T., et al. (2018). Tumor cell-intrinsic factors underlie heterogeneity of immune cell infiltration and response to immunotherapy. Immunity 49, 178–193.e177. doi: 10.1016/j.immuni.2018.06.006
Liu, L., Liu, C., Quintero, A., Wu, L., Yuan, Y., Wang, M., et al. (2019). Deconvolution of single-cell multi-omics layers reveals regulatory heterogeneity. Nat. Commun. 10:470. doi: 10.1038/s41467-018-08205-7
Lodato, M. A., Woodworth, M. B., Lee, S., Evrony, G. D., Mehta, B. K., Karger, A., et al. (2015). Somatic mutation in single human neurons tracks developmental and transcriptional history. Science 350, 94–98. doi: 10.1126/science.aab1785
Ludwig, L. S., Lareau, C. A., Ulirsch, J. C., Christian, E., Muus, C., Li, L. H., et al. (2019). Lineage tracing in humans enabled by mitochondrial mutations and single-cell genomics. Cell 176:e1322. doi: 10.1016/j.cell.2019.01.022
Magnuson, A. M., Kiner, E., Ergun, A., Park, J. S., Asinovski, N., Ortiz-Lopez, A., et al. (2018). Identification and validation of a tumor-infiltrating Treg transcriptional signature conserved across species and tumor types. Proc. Natl. Acad. Sci. U.S.A. 115, E10672–E10681. doi: 10.1073/pnas.1810580115
Mayoux, M., Roller, A., Pulko, V., Sammicheli, S., Chen, S., Sum, E., et al. (2020). Dendritic cells dictate responses to PD-L1 blockade cancer immunotherapy. Sci. Transl. Med. 12:eaav7431. doi: 10.1126/scitranslmed.aav7431
McConnell, M. J., Lindberg, M. R., Brennand, K. J., Piper, J. C., Voet, T., Cowing-Zitron, C., et al. (2013). Mosaic copy number variation in human neurons. Science 342, 632–637. doi: 10.1126/science.1243472
McGranahan, N., Furness, A. J., Rosenthal, R., Ramskov, S., Lyngaa, R., Saini, S. K., et al. (2016). Clonal neoantigens elicit T cell immunoreactivity and sensitivity to immune checkpoint blockade. Science 351, 1463–1469. doi: 10.1126/science.aaf1490
Michea, P., Noel, F., Zakine, E., Czerwinska, U., Sirven, P., Abouzid, O., et al. (2018). Adjustment of dendritic cells to the breast-cancer microenvironment is subset specific. Nat. Immunol. 19, 885–897. doi: 10.1038/s41590-018-0145-8
Miller, B. C., Sen, D. R., Al Abosy, R., Bi, K., Virkud, Y. V., LaFleur, M. W., et al. (2019). Subsets of exhausted CD8+ T cells differentially mediate tumor control and respond to checkpoint blockade. Nat. Immunol. 20, 326–336. doi: 10.1038/s41590-019-0312-6
Mimitou, E. P., Cheng, A., Montalbano, A., Hao, S., Stoeckius, M., Legut, M., et al. (2019). Multiplexed detection of proteins, transcriptomes, clonotypes and CRISPR perturbations in single cells. Nat. Methods 16, 409–412. doi: 10.1038/s41592-019-0392-0
Miragaia, R. J., Gomes, T., Chomka, A., Jardine, L., Riedel, A., Hegazy, A. N., et al. (2019). Single-cell transcriptomics of regulatory T cells reveals trajectories of tissue adaptation. Immunity 50, 493–504.e497. doi: 10.1016/j.immuni.2019.01.001
Mittal, D., Gubin, M. M., Schreiber, R. D., and Smyth, M. J. (2014). New insights into cancer immunoediting and its three component phases–elimination, equilibrium and escape. Curr. Opin. Immunol. 27, 16–25. doi: 10.1016/j.coi.2014.01.004
Mognol, G. P., Spreafico, R., Wong, V., Scott-Browne, J. P., Togher, S., Hoffmann, A., et al. (2017). Exhaustion-associated regulatory regions in CD8+ tumor-infiltrating T cells. Proc. Natl. Acad. Sci. U.S.A. 114, E2776–E2785. doi: 10.1073/pnas.1620498114
Müller, S., Kohanbash, G., Liu, S. J., Alvarado, B., Carrera, D., Bhaduri, A., et al. (2017). Single-cell profiling of human gliomas reveals macrophage ontogeny as a basis for regional differences in macrophage activation in the tumor microenvironment. Genome Biol. 18:234. doi: 10.1186/s13059-017-1362-4
Nathan, A., Baglaenko, Y., Fonseka, C. Y., Beynor, J. I., and Raychaudhuri, S. (2019). Multimodal single-cell approaches shed light on T cell heterogeneity. Curr. Opin. Immunol. 61, 17–25. doi: 10.1016/j.coi.2019.07.002
Oh, D. Y., Kwek, S. S., Raju, S. S., Li, T., McCarthy, E., Chow, E., et al. (2020). Intratumoral CD4(+) T cells mediate anti-tumor cytotoxicity in human bladder cancer. Cell 181, 1612–1625.e1613. doi: 10.1016/j.cell.2020.05.017
Oh, S. A., Wu, D.-C., Cheung, J., Navarro, A., Xiong, H., Cubas, R., et al. (2020). PD-L1 expression by dendritic cells is a key regulator of T-cell immunity in cancer. Nature Cancer 1, 681–691. doi: 10.1038/s43018-020-0075-x
Pagès, F., Mlecnik, B., Marliot, F., Bindea, G., Ou, F.-S., Bifulco, C., et al. (2018). International validation of the consensus Immunoscore for the classification of colon cancer: a prognostic and accuracy study. Lancet 391, 2128–2139. doi: 10.1016/S0140-6736(18)30789-X
Park, J. E., Botting, R. A., Dominguez Conde, C., Popescu, D. M., Lavaert, M., Kunz, D. J., et al. (2020). A cell atlas of human thymic development defines T cell repertoire formation. Science 367:eaay3224. doi: 10.1126/science.aay3224
Pauken, K. E., Sammons, M. A., Odorizzi, P. M., Manne, S., Godec, J., Khan, O., et al. (2016). Epigenetic stability of exhausted T cells limits durability of reinvigoration by PD-1 blockade. Science 354, 1160–1165. doi: 10.1126/science.aaf2807
Peng, D., Kryczek, I., Nagarsheth, N., Zhao, L., Wei, S., Wang, W., et al. (2015). Epigenetic silencing of TH1-type chemokines shapes tumour immunity and immunotherapy. Nature 527, 249–253. doi: 10.1038/nature15520
Peterson, V. M., Zhang, K. X., Kumar, N., Wong, J., Li, L., Wilson, D. C., et al. (2017). Multiplexed quantification of proteins and transcripts in single cells. Nature Biotechnol. 35, 936–939. doi: 10.1038/nbt.3973
Petitprez, F., de Reyniès, A., Keung, E. Z., Chen, T. W.-W., Sun, C.-M., Calderaro, J., et al. (2020). B cells are associated with survival and immunotherapy response in sarcoma. Nature 577, 556–560. doi: 10.1038/s41586-019-1906-8
Philip, M., Fairchild, L., Sun, L., Horste, E. L., Camara, S., Shakiba, M., et al. (2017). Chromatin states define tumour-specific T cell dysfunction and reprogramming. Nature 545, 452–456. doi: 10.1038/nature22367
Pitt, J. M., Vetizou, M., Daillere, R., Roberti, M. P., Yamazaki, T., Routy, B., et al. (2016). Resistance mechanisms to immune-checkpoint blockade in cancer: tumor-intrinsic and -extrinsic factors. Immunity 44, 1255–1269. doi: 10.1016/j.immuni.2016.06.001
Plitas, G., Konopacki, C., Wu, K., Bos, P. D., Morrow, M., Putintseva, E. V., et al. (2016). Regulatory T cells exhibit distinct features in human breast cancer. Immunity 45, 1122–1134. doi: 10.1016/j.immuni.2016.10.032
Puram, S. V., Tirosh, I., Parikh, A. S., Patel, A. P., Yizhak, K., Gillespie, S., et al. (2017). Single-cell transcriptomic analysis of primary and metastatic tumor ecosystems in head and neck cancer. Cell 171, 1611–1624.e1624. doi: 10.1016/j.cell.2017.10.044
Robert, C., Long, G. V., Brady, B., Dutriaux, C., Maio, M., Mortier, L., et al. (2014). Nivolumab in previously untreated melanoma without BRAF mutation. N. Engl. J. Med. 372, 320–330. doi: 10.1056/NEJMoa1412082
Roberts, E. W., Broz, M. L., Binnewies, M., Headley, M. B., Nelson, A. E., Wolf, D. M., et al. (2016). Critical role for CD103(+)/CD141(+) dendritic cells bearing CCR7 for tumor antigen trafficking and priming of T cell immunity in melanoma. Cancer Cell 30, 324–336. doi: 10.1016/j.ccell.2016.06.003
Rodriques, S. G., Stickels, R. R., Goeva, A., Martin, C. A., Murray, E., Vanderburg, C. R., et al. (2019). Slide-seq: a scalable technology for measuring genome-wide expression at high spatial resolution. Science 363, 1463–1467. doi: 10.1126/science.aaw1219
Rooijers, K., Markodimitraki, C. M., Rang, F. J., de Vries, S. S., Chialastri, A., de Luca, K. L., et al. (2019). Simultaneous quantification of protein-DNA contacts and transcriptomes in single cells. Nat. Biotechnol. 37, 766–772. doi: 10.1038/s41587-019-0150-y
Rooney, M. S., Shukla, S. A., Wu, C. J., Getz, G., and Hacohen, N. (2015). Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell 160, 48–61. doi: 10.1016/j.cell.2014.12.033
Rosenberg, A. B., Roco, C. M., Muscat, R. A., Kuchina, A., Sample, P., Yao, Z., et al. (2018). Single-cell profiling of the developing mouse brain and spinal cord with split-pool barcoding. Science 360, 176–182. doi: 10.1126/science.aam8999
Rosenthal, R., Cadieux, E. L., Salgado, R., Bakir, M. A., Moore, D. A., Hiley, C. T., et al. (2019). Neoantigen-directed immune escape in lung cancer evolution. Nature 567, 479–485. doi: 10.1038/s41586-019-1032-7
Rubin, A. J., Parker, K. R., Satpathy, A. T., Qi, Y., Wu, B., Ong, A. J., et al. (2019). Coupled single-cell CRISPR screening and epigenomic profiling reveals causal gene regulatory networks. Cell 176:e317. doi: 10.1016/j.cell.2018.11.022
Sade-Feldman, M., Yizhak, K., Bjorgaard, S. L., Ray, J. P., de Boer, C. G., Jenkins, R. W., et al. (2018). Defining T cell states associated with response to checkpoint immunotherapy in melanoma. Cell 175, 998–1013.e1020. doi: 10.1016/j.cell.2018.10.038
Salmon, H., Idoyaga, J., Rahman, A., Leboeuf, M., Remark, R., Jordan, S., et al. (2016). Expansion and activation of CD103(+) dendritic cell progenitors at the tumor site enhances tumor responses to therapeutic PD-L1 and BRAF inhibition. Immunity 44, 924–938. doi: 10.1016/j.immuni.2016.03.012
Satpathy, A. T., Saligrama, N., Buenrostro, J. D., Wei, Y., Wu, B., Rubin, A. J., et al. (2018). Transcript-indexed ATAC-seq for precision immune profiling. Nat. Med. 24, 580–590. doi: 10.1038/s41591-018-0008-8
Sautès-Fridman, C., Petitprez, F., Calderaro, J., and Fridman, W. H. (2019). Tertiary lymphoid structures in the era of cancer immunotherapy. Nat. Rev. Cancer 19, 307–325. doi: 10.1038/s41568-019-0144-6
Savas, P., Virassamy, B., Ye, C., Salim, A., Mintoff, C. P., Caramia, F., et al. (2018). Single-cell profiling of breast cancer T cells reveals a tissue-resident memory subset associated with improved prognosis. Nat. Med. 24, 986–993. doi: 10.1038/s41591-018-0078-7
Scheper, W., Kelderman, S., Fanchi, L. F., Linnemann, C., Bendle, G., de Rooij, M. A. J., et al. (2019). Low and variable tumor reactivity of the intratumoral TCR repertoire in human cancers. Nat. Med. 25, 89–94. doi: 10.1038/s41591-018-0266-5
Schietinger, A., Philip, M., Krisnawan, V. E., Chiu, E. Y., Delrow, J. J., Basom, R. S., et al. (2016). Tumor-specific T cell dysfunction is a dynamic antigen-driven differentiation program initiated early during tumorigenesis. Immunity 45, 389–401. doi: 10.1016/j.immuni.2016.07.011
Siddiqui, I., Schaeuble, K., Chennupati, V., Fuertes Marraco, S. A., Calderon-Copete, S., Pais Ferreira, D., et al. (2019). Intratumoral Tcf1+PD-1+CD8+ T cells with stem-like properties promote tumor control in response to vaccination and checkpoint blockade immunotherapy. Immunity 50, 195–211.e110. doi: 10.1016/j.immuni.2018.12.021
Simoni, Y., Becht, E., Fehlings, M., Loh, C. Y., Koo, S. L., Teng, K. W. W., et al. (2018). Bystander CD8(+) T cells are abundant and phenotypically distinct in human tumour infiltrates. Nature 557, 575–579. doi: 10.1038/s41586-018-0130-2
Singhal, S., Stadanlick, J., Annunziata, M. J., Rao, A. S., Bhojnagarwala, P. S., O’Brien, S., et al. (2019). Human tumor-associated monocytes/macrophages and their regulation of T cell responses in early-stage lung cancer. Sci. Transl. Med. 11:eaat1500. doi: 10.1126/scitranslmed.aat1500
Snell, L. M., MacLeod, B. L., Law, J. C., Osokine, I., Elsaesser, H. J., Hezaveh, K., et al. (2018). CD8(+) T cell priming in established chronic viral infection preferentially directs differentiation of memory-like cells for sustained immunity. Immunity 49:e675. doi: 10.1016/j.immuni.2018.08.002
Spanjaard, B., Hu, B., Mitic, N., Olivares-Chauvet, P., Janjuha, S., Ninov, N., et al. (2018). Simultaneous lineage tracing and cell-type identification using CRISPR-Cas9-induced genetic scars. Nat. Biotechnol. 36, 469–473. doi: 10.1038/nbt.4124
Stoeckius, M., Hafemeister, C., Stephenson, W., Houck-Loomis, B., Chattopadhyay, P. K., Swerdlow, H., et al. (2017). Simultaneous epitope and transcriptome measurement in single cells. Nat. Methods 14, 865–868. doi: 10.1038/nmeth.4380
Stuart, T., Butler, A., Hoffman, P., Hafemeister, C., Papalexi, E., Mauck, W. M., et al. (2019). Comprehensive integration of single-cell data. Cell 177, 1888–1902.e1821. doi: 10.1016/j.cell.2019.05.031
Stubbington, M. J. T., Lonnberg, T., Proserpio, V., Clare, S., Speak, A. O., Dougan, G., et al. (2016). T cell fate and clonality inference from single-cell transcriptomes. Nat. Methods 13, 329–332. doi: 10.1038/nmeth.3800
Thommen, D. S., Koelzer, V. H., Herzig, P., Roller, A., Trefny, M., Dimeloe, S., et al. (2018). A transcriptionally and functionally distinct PD-1+ CD8+ T cell pool with predictive potential in non-small-cell lung cancer treated with PD-1 blockade. Nat. Med. 24, 994–1004. doi: 10.1038/s41591-018-0057-z
Tirosh, I., Izar, B., Prakadan, S. M., Wadsworth, M. H., Treacy, D., Trombetta, J. J., et al. (2016). Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 352, 189–196. doi: 10.1126/science.aad0501
Togashi, Y., Shitara, K., and Nishikawa, H. (2019). Regulatory T cells in cancer immunosuppression — implications for anticancer therapy. Nat. Rev. Clin. Oncol. 16, 356–371. doi: 10.1038/s41571-019-0175-7
Trapnell, C., Cacchiarelli, D., Grimsby, J., Pokharel, P., Li, S., Morse, M., et al. (2014). The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat. Biotechnol. 32, 381–386. doi: 10.1038/nbt.2859
Tritschler, S., Büttner, M., Fischer, D. S., Lange, M., Bergen, V., Lickert, H., et al. (2019). Concepts and limitations for learning developmental trajectories from single cell genomics. Development 146:dev170506. doi: 10.1242/dev.170506
Utzschneider, D. T., Charmoy, M., Chennupati, V., Pousse, L., Ferreira, D. P., Calderon-Copete, S., et al. (2016). T cell factor 1-expressing memory-like CD8+ T cells sustain the immune response to chronic viral infections. Immunity 45, 415–427. doi: 10.1016/j.immuni.2016.07.021
Villani, A.-C., Satija, R., Reynolds, G., Sarkizova, S., Shekhar, K., Fletcher, J., et al. (2017). Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors. Science 356:eaah4573. doi: 10.1126/science.aah4573
Wagner, J., Rapsomaniki, M. A., Chevrier, S., Anzeneder, T., Langwieder, C., Dykgers, A., et al. (2019). A single-cell atlas of the tumor and immune ecosystem of human breast cancer. Cell 177, 1330–1345.e1318. doi: 10.1016/j.cell.2019.03.005
Wang, L., Simons, D. L., Lu, X., Tu, T. Y., Solomon, S., Wang, R., et al. (2019). Connecting blood and intratumoral Treg cell activity in predicting future relapse in breast cancer. Nat. Immunol. 20, 1220–1230. doi: 10.1038/s41590-019-0429-7
Wei, S. C., Levine, J. H., Cogdill, A. P., Zhao, Y., Anang, N. A. S., Andrews, M. C., et al. (2017). Distinct cellular mechanisms underlie anti-CTLA-4 and anti-PD-1 checkpoint blockade. Cell 170, 1120–1133.e1117. doi: 10.1016/j.cell.2017.07.024
Weinreb, C., Rodriguez-Fraticelli, A., Camargo, F. D., and Klein, A. M. (2020). Lineage tracing on transcriptional landscapes links state to fate during differentiation. Science 367:eaaw3381. doi: 10.1126/science.aaw3381
Welch, J. D., Kozareva, V., Ferreira, A., Vanderburg, C., Martin, C., and Macosko, E. Z. (2019). Single-cell multi-omic integration compares and contrasts features of brain cell identity. Cell 177, 1873–1887.e1817. doi: 10.1016/j.cell.2019.05.006
Wu, T., Ji, Y., Moseman, E. A., Xu, H. C., Manglani, M., Kirby, M., et al. (2016). The TCF1-Bcl6 axis counteracts type I interferon to repress exhaustion and maintain T cell stemness. Sci. Immunol. 1:eaai8593. doi: 10.1126/sciimmunol.aai8593
Wu, T. D., Madireddi, S., de Almeida, P. E., Banchereau, R., Chen, Y.-J. J., Chitre, A. S., et al. (2020). Peripheral T cell expansion predicts tumour infiltration and clinical response. Nature 579, 274–278. doi: 10.1038/s41586-020-2056-8
Xu, J., Nuno, K., Litzenburger, U. M., Qi, Y., Corces, M. R., Majeti, R., et al. (2019). Single-cell lineage tracing by endogenous mutations enriched in transposase accessible mitochondrial DNA. eLife 8:e45105. doi: 10.7554/eLife.45105.033
Yao, C., Sun, H. W., Lacey, N. E., Ji, Y., Moseman, E. A., Shih, H. Y., et al. (2019). Single-cell RNA-seq reveals TOX as a key regulator of CD8(+) T cell persistence in chronic infection. Nat. Immunol. 20, 890–901. doi: 10.1038/s41590-019-0403-4
Yost, K. E., Satpathy, A. T., Wells, D. K., Qi, Y., Wang, C., Kageyama, R., et al. (2019). Clonal replacement of tumor-specific T cells following PD-1 blockade. Nat. Med. 25, 1251–1259. doi: 10.1038/s41591-019-0522-3
Youngblood, B., Oestreich, K. J., Ha, S. J., Duraiswamy, J., Akondy, R. S., West, E. E., et al. (2011). Chronic virus infection enforces demethylation of the locus that encodes PD-1 in antigen-specific CD8(+) T cells. Immunity 35, 400–412. doi: 10.1016/j.immuni.2011.06.015
Zemmour, D., Zilionis, R., Kiner, E., Klein, A. M., Mathis, D., and Benoist, C. (2018). Single-cell gene expression reveals a landscape of regulatory T cell phenotypes shaped by the TCR. Nat. Immunol. 19, 291–301. doi: 10.1038/s41590-018-0051-0
Zhang, L., Li, Z., Skrzypczynska, K. M., Fang, Q., Zhang, W., O’Brien, S. A., et al. (2020). Single-cell analyses inform mechanisms of myeloid-targeted therapies in colon cancer. Cell 181, 442–459.e429. doi: 10.1016/j.cell.2020.03.048
Zhang, L., Yu, X., Zheng, L., Zhang, Y., Li, Y., Fang, Q., et al. (2018). Lineage tracking reveals dynamic relationships of T cells in colorectal cancer. Nature 564, 268–272. doi: 10.1038/s41586-018-0694-x
Zhang, Q., He, Y., Luo, N., Patel, S. J., Han, Y., Gao, R., et al. (2019). Landscape and dynamics of single immune cells in hepatocellular carcinoma. Cell 179, 829–845.e820. doi: 10.1016/j.cell.2019.10.003
Zheng, C., Zheng, L., Yoo, J.-K., Guo, H., Zhang, Y., Guo, X., et al. (2017). Landscape of Infiltrating T cells in liver cancer revealed by single-cell sequencing. Cell 169, 1342–1356.e1316. doi: 10.1016/j.cell.2017.05.035
Zhu, C., Yu, M., Huang, H., Juric, I., Abnousi, A., Hu, R., et al. (2019). An ultra high-throughput method for single-cell joint analysis of open chromatin and transcriptome. Nat. Struct. Mol. Biol. 26, 1063–1070. doi: 10.1038/s41594-019-0323-x
Zilionis, R., Engblom, C., Pfirschke, C., Savova, V., Zemmour, D., Saatcioglu, H. D., et al. (2019). Single-cell transcriptomics of human and mouse lung cancers reveals conserved myeloid populations across individuals and species. Immunity 50, 1317–1334.e1310. doi: 10.1016/j.immuni.2019.03.009
Keywords: single-cell sequencing, single-cell multi-omics, tumor microenvironment, immunotherapy, tumor-specific immunity, tumor infiltrating lymphocytes (TILs), tumor infiltrating myeloid cells (TIMs)
Citation: Guo T, Li W and Cai X (2020) Applications of Single-Cell Omics to Dissect Tumor Microenvironment. Front. Genet. 11:548719. doi: 10.3389/fgene.2020.548719
Received: 03 April 2020; Accepted: 19 October 2020;
Published: 27 November 2020.
Edited by:Xianwen Ren, Peking University, China
Reviewed by:Penghui Zhou, Dana–Farber Cancer Institute, United States
Lai Wei, Zhongshan Ophthalmic Center, Sun Yat-sen University, China
Copyright © 2020 Guo, Li and Cai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.