Breaking the Immune Complexity of the Tumor Microenvironment Using Single-Cell Technologies

Tumors are not a simple aggregate of transformed cells but rather a complicated ecosystem containing various components, including infiltrating immune cells, tumor-related stromal cells, endothelial cells, soluble factors, and extracellular matrix proteins. Profiling the immune contexture of this intricate framework is now mandatory to develop more effective cancer therapies and precise immunotherapeutic approaches by identifying exact targets or predictive biomarkers, respectively. Conventional technologies are limited in reaching this goal because they lack high resolution. Recent developments in single-cell technologies, such as single-cell RNA transcriptomics, mass cytometry, and multiparameter immunofluorescence, have revolutionized the cancer immunology field, capturing the heterogeneity of tumor-infiltrating immune cells and the dynamic complexity of tenets that regulate cell networks in the tumor microenvironment. In this review, we describe some of the current single-cell technologies and computational techniques applied for immune-profiling the cancer landscape and discuss future directions of how integrating multi-omics data can guide a new “precision oncology” advancement.

Introduction

Cancer is one of the leading causes of death worldwide, accounting for nearly 10 million deaths in 2020 (https://www.who.int/news-room/fact-sheets/detail/cancer). Unfortunately, the pandemic COVID-19 will have consequences for cancer patients in coming years, since it has been associated with delays in diagnosis as well as interruption of therapeutic treatments and follow-up care. Hence the number of cancer victims will increase in the near future. Identifying cancer as a genetic disease characterized by a set of genomic aberrations, including in-frame insertions or deletions, missense amino acid changes, and large copy number variations, initially ingrained a “cancer cell-centric” vision in the scientific community where cancer cell-intrinsic properties exclusively drove tumorigenesis (Reddy et al., 1982; Santos et al., 1984; Lengauer et al., 1998; Futreal et al., 2004; Tomlins et al., 2005; Hanahan and Weinberg, 2011). Therefore, recognizing driver gene modifications has been a central aim of cancer research over the past 30 years, resulting in global initiatives such as The Cancer Genome Atlas (TCGA) (https://www.cancer.gov/about-nci/organization/ccg/research/structural-genomics/tcga) and the Catalog of Somatic Mutations in Cancer (COSMIC) (Tate et al., 2019) which include a broad collection of large-scale, systematic sequencing studies that constitute comprehensive catalogs of mutational abnormalities in the major tumor types. More than 100,000 somatic mutations in cancer genomes have been identified in the quarter-century since the first somatic mutation was reported (Zack et al., 2013), opening a new age for the classification and treatment protocols of some cancers, such as breast (Cancer Genome Atlas, 2012b; Ciriello et al., 2015), cervical (Cancer Genome Atlas Research Network, 2017), colorectal (Cancer Genome Atlas, 2012a), pancreatic (Consortium, 2020), gastric (Cancer Genome Atlas Research, 2014a), prostate (Cancer Genome Atlas Research, 2015), and lung carcinomas (Cancer Genome Atlas Research, 2012; Cancer Genome Atlas Research, 2014b).

In addition to the study of genetic diversity as a source of tumor cell heterogeneity, stable epigenetic changes in cancer have been received intense interest since they derive from poised or initiated chromatin states of several genes that can modulate expression of different pathways (Brown et al., 2014). Indeed, alterations of cellular activities, such as cell growth and differentiation, can be driven by epigenetic events that involve DNA methylation, histone modification, the readout of these modifications, chromatin remodeling and the effects of noncoding RNA (Feinberg et al., 2016). These alterations are temporary or yet long-lasting and they impact on tumorigenesis. For instance, aberrant DNA methylation has been, generally, associated with cancer development by inactivating gene transcription or repressing gene transcription and affecting chromatin structure (Baylin and Jones, 2011; Hanahan and Weinberg, 2011). Indeed, some gene promoters, especially key tumor suppressor genes, are unmethylated in normal tissues and highly methylated in cancer (Baylin and Jones, 2011). Interestingly, DNA methylation in cancer has generally been associated with tumor size (Feinberg et al., 2016) as well as with drug resistance and predicting response to treatment (Brown et al., 2014). Therefore, the comprehensive genetic and epigenetic analysis of cancer genomes has been for many years the most effective way to identify causative changes involved in tumorigenesis.

Despite the achievement of these extraordinary milestones in deciphering cancer cell biology, new insights demonstrate that aberrant genetic profiles of transforming cells alone are required but insufficient to nurse tumor development and progression. Indeed, cancer cells need to alter the stromal framework of the microenvironment to manipulate diverse physiological processes, such as promoting angiogenesis and vasculogenesis, to provide adequate nutrients and oxygen, and alter immune responses to avoid activating tumor-fighting elements, such as cytotoxic T lymphocytes (CTLs) (Hanahan and Weinberg, 2011). Therefore, “cancer-cell extrinsic” factors, such as local inflammation, metabolic switch, and immunity, are critical in fueling cancer growth. At this point, the key question is whether these extrinsic factors are independent of the genetic profile of cancer cells. A cornerstone study demonstrated a remarkable difference in the composition of tumor-infiltrating leukocytes in different tumor types by analyzing data on clinical outcomes and gene expression of 18,000 human tumors (Gentles et al., 2015). Interestingly, memory CD4⁺ T lymphocyte frequency correlated positively with a favorable outcome in lung cancer patients, whereas the same cell subset was associated with a worse outcome in patients affected with bladder tumors (Gentles et al., 2015), suggesting that cancer cell-intrinsic features can dictate the immune landscape of the tumor microenvironment (TME). Indeed, oncogene-driven modifications can alter tumor immunogenicity in a completely different way. Ongoing mutational processes generate either cancer neoantigens capable of activating tumor-eliminating immune cells (Balachandran et al., 2017; Luksza et al., 2017; Keenan et al., 2019) or produce immune soluble factors such as interleukins (e.g., IL-6, IL1-β), growth factors (e.g., granulocyte-macrophage colony-stimulating factor (GM-CSF)), and chemokines (e.g., CCL4), capable of differentiating immune cells into pro-tumor elements (Bronte et al., 2006; Pylayeva-Gupta et al., 2012; Calcinotto et al., 2018; Fiore et al., 2018; Vitale et al., 2019; Hofer et al., 2021). Host immune cells also influence cancer progression (Schreiber et al., 2011). The notion that tumors derived from immunodeficient hosts are more immunogenic than those derived from immunocompetent mice allowed us to hypothesize that the immune system actively shapes cancer cells, promoting the acquisition of genetic aberrations that can compromise cancer cell immunogenicity, favoring tumors escaping immune attacks. Several clinical and pre-clinical observations have validated the cancer immunoediting theory. For instance, two large studies demonstrated that tumor-infiltrating immune subsets in colorectal cancer were significant independent prognostic markers as well as microsatellite instability, long interspersed nucleotide element-1 (LINE-1) hypomethylation, and BRAF mutations (Ogino et al., 2009; Nosho et al., 2010). Moreover, the linear correlation between the density of effector memory CTLs at the site of the primary tumor and the survival of patients explicitly revealed the importance of TME immunity in cancer control (Galon et al., 2006; Fridman et al., 2017). Thus, the definition of the cancer genome and immune landscape of the TME is hierarchically equivalent and complementary to predict disease progression and therapeutic outcome.

Immune cell heterogeneity and the presence of various cell subsets complicate TME immune profiling (Binnewies et al., 2018). Moreover, the spatial distribution of immunity within the tumor mass is a critical parameter that significantly influences immune cells acquiring pro- or anti-tumor functions (Bindea et al., 2013; Sautes-Fridman et al., 2019). A multi-omic perspective considering genomics, transcriptomics, epigenomics and proteomics is necessary to unveil the immune complexity of the TME (Figures 1A–D). Unfortunately, traditional technologies used to evaluate the immune landscape of the TME, such as immunohistochemistry (IHC), flow cytometry (FC), and bulk analysis of genomic, transcriptional, and proteomic analyses, have several limitations, such as large amount of biological material requirements, fewer parameters tested simultaneously, and low analysis resolution (Peregrin et al., 1973; Bendall et al., 2012; Abel et al., 2014). Recently, high-dimensional single-cell techniques (Finotello and Eduati, 2018; Chen G. et al., 2019; Gohil et al., 2021) have revolutionized the approach to decipher the cellular diversity, cell interactions and dynamics that exists in the TME and heterogeneity across patients by resolving cell subset complexity at the single-cell level using unsupervised clustering to identify potential unknown subpopulations of cells within the populations under study (Thorsson et al., 2019; Wagner et al., 2019; Zhao J. et al., 2020; Kieffer et al., 2020; Cheng et al., 2021). Here, we review how single-cell technologies (Figure 1E) and related computational techniques have improved our knowledge about the TME and discuss future applications of these cutting-edge techniques in immune-oncology to develop more effective personalized immunotherapy.

FIGURE 1

FIGURE 1. Multi-omic perspective to study the features of the TME. (A) Genomics analysis informs about how the tumor mutational landscape influences the TME, favoring, for example, the production of cytokines (e.g., IL-6, IL1-β) and growth factors (e.g., GM-CSF) inducing the proliferation of suppressive myeloid cells and the pro-tumor differentiation of antigen-presenting cells. (B) Transcriptomic analysis enables to inspect the transcriptional machinery of the single cells of the TME, deciphering, for example, developmental trajectories, cell states and cell-cell interactions. (C) Epigenomic analysis reveals how specific switches such as histone methylation and chromatin dynamics regulate different mechanisms capable to interfere with anti-tumor immune recognition and effector functions. (D) Proteomics provides information about the state of activation of the immune cells of the TME looking at the expression of immunomodulatory proteins such as checkpoint inhibitors (e.g., PD-L1, CTLA-4). Spatial proteomics gives additional information about the localization of the cells allowing, for example, to identify cell that interact in the TME. (E) Venn diagram depicting single-cell technologies to study single (non-overlapping sets) or multimodal (overlapping sets) omics. Genomics, transcriptomics, epigenomics and proteomics are represented as blue, red, green, and purple sets, respectively.

Single-Cell Technologies to Study the TME

Single-Cell Genomics

Several studies have demonstrated the connection between the tumor mutations and the immune composition or TME. Deciphering the complexity of the mutational landscape of tumor cell clones and sub-clones that underlies intratumor heterogeneity is necessary to understand tumor patient lethal outcomes, therapeutic failures, and drug resistance (McGranahan and Swanton, 2017). Single-cell genomics technologies can make an important contribution to this goal. However, the study of tumor clonality by analyzing single-nucleotide variations (SNVs) and copy number variations (CNVs) at the single-cell level is challenging. Most methods for creating single-cell libraries rely on whole-genome amplification (WGA) to overcome the inability of sequencing technologies to capture single-cell DNA molecules in low amounts of material. Since the first method based on degenerate oligonucleotide-primed PCR (DOP-PCR) (Telenius et al., 1992), several methods have been developed to increase the coverage and uniformity of the genome to allow both SNVs and CSVs to be studied within the same experiment (Dean et al., 2001; Zong et al., 2012), given the importance of the WGA step. More recent methods based on tagmentation (Adey et al., 2010), such as direct library preparation (DLP) (Zahn et al., 2017), and linear amplification via transposon insertion (LIANTI) (Chen et al., 2017), outperformed WGA methods in terms of accuracy, time, and cost, expanding potential future application scenarios.

Finally, high-throughput methods capable of processing thousands of single cells per experiment were introduced, such as those based on microfluidics and barcoding (SiC-seq) (Lan et al., 2017), split-pool strategies (SCI-seq) (Vitak et al., 2017), high-throughput versions of linear amplification via transposon insertion (sci-L3) (Yin et al., 2019), and direct library preparation (DLP+) (Laks et al., 2019). The aim of high-throughput methods was introducing automation to increase the number of cells analyzed while maintaining accuracy of genomic information. SiC-seq exploits droplet microfluidics to encapsulate cells into microspheres in which to perform the reactions required for cell and genome processing without compromising genomic DNA. SiC-seq performs a series of steps to lyse the cells, fragment the genomes, barcode the DNA fragments and sequence them after library preparation. SCI-seq uses the strategy of transposase-based combinatorial indexing (Adey et al., 2010) to obtain barcoded libraries without using droplet microfluidics. Sci-L3 addresses the problem of genomic artifacts due to PCR amplification and low-throughput using respectively linear amplification and a three-level combinatorial indexing. DLP + takes advantage of specialized hardware and software for imaging microscopy to capture genomic information of thousands of cells per experiment. Given the complexity in obtaining single cell genomic information, commercial solutions such as those based on the 10x Genomics droplet microfluidics system have been also proposed to simplify single-cell genomics data generation (Figure 2A). Single-cell genomic technologies have been used, for example, for the deconvolution of clonal cell clusters and tracing the evolutionary trajectories of clonal breast cancer cells (Wang et al., 2014), the identification of structural and mutational events of melanoma cell line clones (Velazquez-Villarreal et al., 2020), and for an in-depth view of the intratumoral copy number alteration (CNA) heterogeneity present in breast cancer genomes (Baslan et al., 2020).

FIGURE 2

FIGURE 2. Representative single-cell technologies to study the TME. (A) 10x Genomics single-cell genomics involves two steps of encapsulation using the microfluidics system. In the first step, cells are partitioned using a cell beads polymer. The obtained cell beads are lysed to denature the genomic DNA and a second step on microfluidics chip is performed to encapsulate cell beads with barcode gel beads. After collecting single cell GEMs, amplification and barcoding of fragments is performed prior to breaking the emulsion and constructing the library for sequencing. (B) In 10x Genomics scRNA-seq, cells are encapsulated into droplets together with barcoded beads. Next, reverse transcription (RT) is performed in the collected GEMs and barcoded cDNAs are amplified for library construction and sequencing. (C) In 10x Genomics scATAC-seq, nuclei are transposed and encapsulated into droplets using the microfluidics chip. Next, the collected single nuclei GEMs are linearly amplified, and barcoded accessible DNA fragments are obtained after breaking the emulsion. Finally, DNA fragments are ready for library construction and sequencing. (D) In CyTOF, the cells are labeled using stable heavy metals, nebulized, and vaporized to form ion clouds through an argon plasma torch. Each cloud passes through a quadrupole which performs a purification step, the remaining heavy ions are quantified by a time-of-flight (TOF) mass spectrometer that determines the value of each marker. (E) In co-detection by indexing (CODEX), FFPE or FF tissues samples are stained with DNA-barcoded antibodies. Next, a multicycle reaction characterized by iteratively imaging up to three antibodies and nuclear stain, stripping and hybridizing is performed. This process is performed for all antibodies. Finally, raw images are processed and analyzed. (F) In CITE-seq, antibody-derived tags (ADTs) are used to bind the cells of interest. Next, cells are incapsulated into droplets using a microfluidics platform and after cell lysis in droplets, mRNAs and ADTs are barcoded during the RT. After amplification, cDNAs and ADTs are separated by size, converted into two independent libraries that are, finally, pooled, and sequenced.

All these aspects offer a greater understanding of the dynamics of the TME, especially in the cancer therapy response context. For example, identifying clonal and sub-clonal cell composition or the presence of specific mutated subsets would be useful for stratifying patients to understand whether they could benefit from immunotherapy. In the future, longitudinal studies coupled with other single-cell omics could help to construct a more detailed map of clone evolution kinetics and elucidate the time-dependent mechanisms underlying therapy and the emergence of more aggressive clones caused by selection pressures (McGranahan and Swanton, 2015). Furthermore, these studies could also suggest evolutionary time points with a favorable TME to perform more effective therapeutic interventions. In the future, the widespread use of single-cell genomics methods will depend on improvements of the actual technical limitations. Current technologies are not optimal to study all the different genomic aberrations such as copy number variations (CNVs), small indels, single-nucleotide variations and structural variations (SVs) at the same time (Fan et al., 2021). Advancements on the current protocols and bioinformatics solutions will be critical for a wider adoption of single-cell genomics by the scientific community.

Single-Cell Transcriptomics

Single-cell RNA-sequencing (scRNA-seq) is a key technique to explore the TME. scRNA-seq has expanded the scenarios opened by previous bulk RNA-seq technology to investigate the transcriptome of single cells inside a biological sample. Despite countless discoveries in the last decade due to the ability to sequence millions of RNA fragments from a “bulk” of cells, it is now mandatory to look at the transcription machinery of each single cell to understand the complexity of the TME. Indeed, the possibility of observing the expression of thousands of genes for each cell has allowed us to resolve cell subset heterogeneity of the TME in an unbiased way (Lambrechts et al., 2018; Zilionis et al., 2019). Furthermore, scientists can go beyond merely characterizing cell compositions and obtain information about cell state, differentiation trajectories, and cellular pathway activation (Qian et al., 2020; Zhang et al., 2020; Raghavan et al., 2021).

The TME is crucial for tumor cells to evade immune surveillance, necessitating detailed identities of cancer, immune, and non-immune cells. A central aspect in the fight against cancer requires an enhanced understanding of the role of immune cells in cancer therapies, especially in immunotherapy based on immune checkpoint inhibitors (ICIs). Hence, in recent years, scientists have tried to understand the reasons for the success or failure of immunotherapy in relation to the immune features of the TME. It is now clear that T cell infiltration into tumor tissues is a key feature of the immunotherapy response (Tang et al., 2016). However, the presence of T lymphocytes is a necessary but insufficient condition for an effective antitumor response. Indeed, the immunotherapeutic reactivity of a tumor largely depends on the functional state of infiltrated T cells and the expression of specific drug-targetable molecules, such as programmed cell death protein 1 (PD-1) and cytotoxic T lymphocyte-associated antigen 4 (CTLA-4), as well as and the lymphocyte-activation gene 3 (LAG-3), which negatively regulates T-cell proliferation and effector T-cell functions. Recently, a phase II/III, global, double-blind, randomized study named RELATIVITY-047, demonstrated that the dual inhibition of LAG-3 and PD-1, using relatlimab and nivolumab, had a synergistic effect on progression-free survival in melanoma patients highlighting the possibility to develop more effective therapy by targeting different T cell activation brakes (Tawbi et al., 2022). Furthermore, the expression of co-stimulatory molecules, such as programmed death-ligand 1 (PD-L1), by cancer cells and an inflamed TME are crucial elements for an effective immunotherapy response (Marigo et al., 2016; Duan et al., 2020; Ugel et al., 2021). In contrast, the presence of myeloid-infiltrating cells with immunosuppressive features, such as myeloid-derived suppressor cells (MDSCs) (De Sanctis et al., 2016) and a low frequency of infiltrating CD8⁺ T cells in the TME, correlate with tumors averse to immunotherapy (Zhu et al., 2017; Li J. et al., 2018).

Cancer immunology discoveries have continued in parallel with advancements in scRNA-seq technologies. Among others, methods based on droplet microfluidics and cellular barcoding, such as InDrop (Klein et al., 2015), Drop-seq (Macosko et al., 2015), and the 10x Genomics platform (Zheng et al., 2017) have gained attention in recent years. Similarly to single-cell genomics technologies, their success is due to the greater scalability offered by droplet microfluidics platforms that allow to isolate thousands of cells in a short time (Figure 2B). A drawback of these technologies is their low transcript coverage. Indeed, these methods cannot be used to study, for example, splicing and enhancer RNAs at the single-cell level (Hayashi et al., 2018) because they sequence only the 3’ end of each transcript. In contrast, full-length technologies, such as SMART-seq (Ramskold et al., 2012; Picelli et al., 2014) and MATQ-seq (Sheng et al., 2017), have a higher accuracy level in transcript detection, allowing the study of gene variants and splicing events at the cost of lower single-cell profiling throughput.

Given the complexity of the TME in terms of cellular heterogeneity, high-throughput scRNA-seq techniques have been widely used in cancer studies because they can extract thousands of cells from a biological sample. scRNA-seq has been successfully used to decipher several mechanisms by which the tumor can shape the microenvironment towards a hostile ecosystem characterized by exhausted infiltrating T cells, suppressive myeloid cell subsets, cells with pro-tumor differentiation programs, and aberrant cell-cell interaction networks (Guo et al., 2018; Elyada et al., 2019; Davidson et al., 2020; Kim et al., 2020; Marigo et al., 2020; Wu et al., 2021; De Sanctis et al., 2022). The integration of scRNA-seq data with clinical patient information has further allowed us to define cellular and molecular signatures of the TME to explain the response to immunotherapy or survival outcome in different types of cancer (Jerby-Arnon et al., 2018; Sade-Feldman et al., 2018; Ma et al., 2019; Peng et al., 2019; Di Pilato et al., 2021; Zhang et al., 2021). These results motivate future applications of scRNA-seq as a valuable tool for clinicians to perform TME screening, stratification, and identification of new druggable targets based on the integration of different patient datasets (Giladi and Amit, 2017; Binnewies et al., 2018). Data integration is of primary importance because transcriptional-level results cannot be directly translated to the functional level, requiring additional support and experimental validation at the protein level.

Recently, scRNA-seq has been combined with the innovative spatial transcriptomics (Zhang et al., 2022). This technology has revolutionized the traditional RNA-fluorescence in situ hybridization (FISH) tool able to identify target messenger RNA transcripts in tissue sections by-passing the absence of selective antibodies for unknown candidates (Facciponte et al., 2014; Cui et al., 2016). Indeed, this innovative approach allows to visualize profiles of RNA molecules in identified tissue regions, including technologies based on micro-dissected gene expression, in situ hybridization, in situ capturing and in situ sequencing technologies (Calvanese et al., 2022; Nirmal et al., 2022; Wei et al., 2022). This technology continues to aid the development of human cellular atlases of cancer, the reclassification of the immune landscape of TME, and overall, the identification of important therapeutic targets.

Single-Cell Epigenomics

The epigenome is the complete atlas of chemical modifications that can induce changes in gene expression without modifying the DNA sequence. Such modifications involve DNA, RNA, and histone proteins and can cause chromatin remodeling that can turn genes “on” or “off”. Methylation and acetylation of histones on lysine and arginine residues are the best-known epigenetic mechanisms capable of enhancing or repressing gene transcription, as exemplified by histone H3 on lysine 27 (H3K27ac) and histone H3 on lysine 9 (H3K9me) modifications. Other epigenetic mechanisms that provoke chromatin remodeling include nucleosome positioning/reorganization and DNA methylation.

Several studies have highlighted the ability of tumor cells to induce epigenetic modifications in TME-infiltrating immune cells to aid immune surveillance evasion. The epigenetic strategies implemented by tumors to avoid immune surveillance are based on the disruption of different anti-cancer immune mechanisms, such as immune recognition, signal activation, and effector functions. Common mechanisms exploited by tumors to evade immune recognition and signal triggering include the epigenetic silencing of major histocompatibility complex (MHC) genes and the inhibition of cytokines and chemokines. For example, trimethylation of H3K4 (H3K4me3) and H3K27 (H3K27me3) by polycomb repressive complex 2 (PRC2) has been linked to MHC-I repression and missed tumor recognition by CD8⁺ T cells mediated by EED and EZH2 (Burr et al., 2019). Furthermore, the lack of tumor-infiltrating lymphocytes (TILs) in several human cancers is associated with DNA methylation-induced epigenetic silencing of CCL5 (Dangaj et al., 2019). An illuminating example of how epigenetic modifications can result in the loss of immune effector function is the acetylation of H3K27 (H3K27ac). Indeed, high H3K27ac levels have been linked to a TNF-NFKB1 pathway capable of inducing CD47 upregulation and inhibiting macrophage phagocytosis of breast cancer cells (Betancur et al., 2017).

Our understanding of the cancer epigenome has evolved rapidly with the adoption of next-generation sequencing (NGS). In this context, chromatin immunoprecipitation followed by sequencing (ChIP-seq) and an assay for transposase-accessible chromatin using sequencing (ATAC-seq) have been widely used tools to study epigenetic regulation. ChIP-seq is a method for studying the interactions between proteins and DNA. It allows us to analyze chromatin states induced by histone modifications that alter gene transcription. ATAC-seq measures chromatin accessibility by directly deciphering its effects on gene transcription without detailed histone modification and chromatin state characterization. ChIP-seq and ATAC-seq have been used successfully, for example, to link epigenetic features capable of maintaining an immune cell-excluded TME and immunotherapy resistance (Benci et al., 2016; Yang et al., 2021).

The advancements of these techniques are their single-cell counterparts, scChIP-seq and scATAC-seq. Several cancer studies have shown the utility of scChIP-seq or scATAC-seq to study the epigenetic regulators responsible for tumor cell heterogeneity (Grosselin et al., 2019; LaFave et al., 2020; Taavitsainen et al., 2021). Other applications of single-cell epigenomics include adopting scATAC-seq to uncover chromatin regulators responsible for T cell exhaustion in the TME of patients treated with immunotherapy (Satpathy et al., 2019; Zhang et al., 2021). These examples motivate the use of single-cell technologies to study cancer epigenetics.

In the future, it will be pertinent to understand the link between epigenetic changes in the TME and cancer progression and how to obtain a more effective therapy response targeting epigenetic switches. One reason for the growing interest in targeting the epigenome is the possibility of identifying small molecules, such as proteolysis-targeting chimeras (PROTAC) (Sakamoto et al., 2001), that can indirectly interfere with aberrant gene expression, given the difficulties in targeting oncogenic transcription factors, such as Myc and p53 (Jones et al., 2016).

The dissemination of scChIP-seq studies depends on future technological advances. Despite recent successful attempts to improve scChIP-seq in terms of cellular sensitivity and data sparsity (Kaya-Okur et al., 2019; Wang et al., 2019), the use of this technology remains limited and largely relies on being able to increase the number of sequenced reads per cell. In contrast, scATAC-seq has a simpler and more efficient experimental protocol that requires a lower number of cells. Furthermore, introducing microfluidics approaches (Figure 2C) has dramatically improved the throughput compared to the previous technology based on combinatorial cellular indexing (sciATAC-seq) (Cusanovich et al., 2015). However, low per-cell coverage remains a weakness, even for scATAC-seq, potentially limiting the identification of significant open chromatin sites, especially in rare cell subsets (Ma and Zhang, 2020).

Single-Cell Proteomics

Despite the wide adoption of single-cell transcriptomics technologies, single-cell proteomic approaches remain the key tools for studying the functional status of TME cell populations. Indeed, in addition to identifying proteins and related isoforms, they allow us to recognize post-translational modifications that single-cell transcriptomics cannot capture.

FC is still a fundamental multiparametric technique for identifying immune cell subsets within the TME based on morphological characteristics and the expression of certain proteins. FC also allows specific cell populations to be isolated by cell sorting before analysis with other omics techniques. The main limitation of FC is the number of parameters, usually around 20, which can be intercepted simultaneously in the same experiment because of signal overlap (i.e., spillover) between the channels. This limits the number of cell-surface proteins that can be identified and, consequently, the resolution of identifying the related cell subsets present in the experiment.

A more advanced technology is Full Spectrum Flow Cytometry (FSFC), an improvement of spectral flow cytometry (SFC) (Robinson, 2004; Nolan et al., 2013). FSFC exploits high-sensitive light detectors to measure the full spectral profile of fluorophores. This technology has a higher quality and resolution than conventional FC, allowing to design multicolor panels up to 40 parameters useful for characterizing important aspects of the immune context in cancer studies (Bonilla et al., 2020).

Another recent technology for single-cell proteomics is cytometry by time-of-flight (CyTOF) (Bandura et al., 2009) (Figure 2D). Compared to FC, CyTOF allows at least 40 markers per cell to be detected in a single run and is more sensitive and less prone to errors (Bandura et al., 2009; Bendall et al., 2011). The disadvantages of CyTOF compared to FC include a lower acquisition flow rate, more critical sample preparation to avoid contamination, and the inability to perform cell sorting for populations of interest due to the final vaporization of the cells (Gadalla et al., 2019). CyTOF potentiates better resolution of TME heterogeneity, which is particularly important in the context of cancer immunology. For example, CyTOF makes finding specific TIL subsets that correlate with patient survival and response to immunotherapy, respectively, in follicular lymphoma (Yang et al., 2020) and melanoma (Subrahmanyam et al., 2018) possible. Furthermore, scientists were able to localize the expression of the T cell inhibitory molecule VISTA in CD68⁺ macrophages of human pancreatic cancer (Blando et al., 2019) and retrieve information on the composition, expansion, and activity of TILs in patients with non-small cell lung cancer (NSCLC) (Sanmamed et al., 2021). Since CyTOF dissects the cellular composition and the activation status of the immune cells that surround and infiltrate the tumor, it is an effective tool for studying the TME immune landscape. CyTOF and scRNA-seq allowed us to answer similar biological questions. However, these two technologies are not interchangeable. The main limitations of CyTOF are the resolution and bias due to prior parameter selection. These two aspects can limit the discovery of rare cell populations. An obvious advantage of CyTOF is that protein expression indicates specific functional states or activities of the cell that in scRNA-seq must be validated with other techniques. Since CyTOF is simpler than scRNA-seq and allows for better discrimination of certain immune cell subsets (Kashima et al., 2021), it represents a convenient method to monitor TME features, such as cell subsets, activation states, and immune checkpoint molecules in patient cohorts in clinical trials. It is mandatory to rapidly standardize experimental procedures, computational tools, and antibody panels to make results comparable between different institutions to accomplish this aim (Hartmann et al., 2019).

In CyTOF, protein identification strictly depends on the availability of highly specific antibodies and the quality of their interactions. Mass spectrometry-based methods, such as SCoPE-MS and the improved version SCoPE2, have been recently introduced to improve throughput and sensitivity concerning the number of proteins detected in single cells (Budnik et al., 2018; Specht et al., 2021). With SCoPE2, the authors were able to dissect cellular heterogeneity by protein expression and trace the differentiation of monocytes into macrophage-like cells in the absence of specific cytokines. This work is an important step towards future applications of single-cell technologies based on mass spectrometry for TME dissection, looking at the expression of thousands of proteins.

Single-Cell Spatial Proteomics

An important aspect to consider in the study of the TME is its spatial organization and heterogeneity. Indeed, like natural ecosystems, tumor tissues can reveal strong heterogeneity in relatively small spatial distances due to tumor cells adapting to the microenvironment or through its remodeling (Yuan, 2016). Several studies have emphasized the importance of considering spatial heterogeneity in the TME. Some of these reports have shown that the location, density, and spatial distribution of immune cells are more robust markers for predicting patient outcomes than traditional clinical parameters (Galon et al., 2006; Maley et al., 2015). These facts motivated the development and improvement of technological platforms for spatial analysis and their use in dissecting the TME.

IHC is a widely used technique for both basic research and cancer diagnosis. It is used to localize cells that express specific protein markers and study the spatial localization of cells in a tissue slide. In conventional IHC (i.e., chromogenic IHC), antibodies recognize specific antigens in the tissue and are conjugated to an enzyme to catalyze a color-producing reaction. A major problem with conventional IHC is that it only allows labeling of one marker for tissue sections. This severely limits understanding the cellular complexity of the TME. Several multiplex platforms have been introduced over the years to facilitate the analysis of cellular composition, functional states, and cell-cell interactions within the TME to address this problem. Multiplex IHC (mIHC) methods based on chromogens (Remark et al., 2016; Tsujikawa et al., 2017), fluorophores (Gerdes et al., 2013; Gorris et al., 2018; Viratham Pulsawatdi et al., 2020), metal-tagged antibodies (Angelo et al., 2014), and DNA barcodes (Goltsev et al., 2018; Manesse et al., 2020) (Figure 2E) increased the number of biomarkers to be used simultaneously by up to 50 in a single tissue section. The introduction of more effective mIHC technologies has been accompanied by the development of new software solutions capable of performing sophisticated analyses of digitalized images (e.g., segmentation and filtering) and support the work of pathologists in sample processing.

mIHC technologies have been used in different types of cancer to better understand the spatial architecture of the TME and how it can affect the response to therapy and clinical outcomes. For example, in breast cancer, mIHC was used to obtain detailed information on the spatial localization, cellular composition, and expression of regulatory proteins in the TME and recover clinically relevant characteristics (Keren et al., 2018; Jackson et al., 2020). In colorectal cancer, mIHC has been used to decipher the complex dynamic interplay between TME components (Schurch et al., 2020). In human pancreatic ductal adenocarcinoma (PDAC), mIHC was applied to assess the density and spatial distribution of myeloid and lymphoid cells in the TME and its correlation with the clinical outcome of patients (Liudahl et al., 2021).

Single-Cell Multimodal Omics

In cancer, single-cell omics aims to dissect all aspects of cellular machinery to understand its functional status and relationship with other TME cells. The use of single-cell omics technologies to study single-cell modalities, such as transcriptomics and proteomics, has increased our knowledge of cell biology in cancer without precedent. However, cells are dynamic entities whose states are characterized by a complex interplay of genomic, transcriptomic, epigenomic, and proteomic features that non-linearly contribute to the TME’s heterogeneity. The new Frontier is the simultaneous measurement of multiple modalities of the same cells to gain a better understanding of cellular and molecular mechanisms in cancer. In recent years, this has motivated the introduction of numerous single-cell approaches capable of combining two or more modalities between genomics, transcriptomics, epigenomics, and proteomics (Figure 1E). Multimodal single-cell approaches have evolved from those using tubes or microwells of plates to measure single cells (i.e., low-throughput) to those that take advantage of droplet-based technologies or combinatorial DNA barcoding strategies (i.e., high-throughput), allowing increased scalability and reduced costs per run (Zhu et al., 2020). Almost all current single-cell multimodal omics technologies extract the transcriptome of each cell. Low-throughput multimodal omics with transcriptomics paired with genomics are gDNA-mRNA sequencing (DR-seq) (Dey et al., 2015) and genome and transcriptome sequencing (G&T-seq) (Macaulay et al., 2015), while transcriptomics and epigenomics are obtained with single-cell methylome and transcriptome sequencing (scM&T) (Angermueller et al., 2016), scMT-seq (Hu et al., 2016), single-cell nucleosome, methylation, and transcription sequencing (snNMT-seq) (Clark et al., 2018), single-cell nucleosome occupancy, methylome, and RNA expression sequencing (scNOMeRe-seq) (Wang et al., 2021), and single-cell chromatin accessibility and transcriptome sequencing (scCAT-seq) (Liu et al., 2019). An interesting multimodal omics technique integrating transcriptomics, genomics, and epigenomics is single-cell triple omics sequencing (scTrio-seq) (Hou et al., 2016).

In high-throughput multimodal omics, transcriptomics is paired with epigenomics in parallel analysis of individual cells for RNA expression and DNA accessibility by sequencing (Paired-seq) (Zhu et al., 2019) and single-nucleus chromatin accessibility and mRNA expression sequencing (SNARE-seq) (Chen S. et al., 2019). Furthermore, transcriptomics is paired with epitopes in RNA expression and protein sequencing (REAP-seq) (Peterson et al., 2017) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) (Stoeckius et al., 2017) (Figure 2F). CITE-seq protocol has been extended to obtain transcriptomics and proteomics of the same cells after CRISPR-Cas9 genetic perturbations in expanded CRISPR-compatible CITE-seq (ECCITE-seq) (Mimitou et al., 2019) and Perturb-CITE-sequencing (Perturb-CITE-seq) (Frangieh et al., 2021). These techniques perform genetic perturbations through single gRNAs (sgRNA) linking them to transcriptomic and proteomic profiles of the same cells to allow demultiplexing. ECCITE-seq combines the use of CRISPR libraries and cell hashtags (Stoeckius et al., 2018) to perform genetic perturbations and pool together different experimental samples. Perturb-CITE-seq uses a method based on CROP-seq (Datlinger et al., 2017) to express sgRNAs and to link them to transcripts and surface proteins of single cells. Other interesting multimodal omics techniques that do not include transcriptomics are, for example, those that provide different layers of epigenomics (Guo et al., 2017; Pott, 2017), epigenomics paired with genomics, and CRISPR-Cas9 genetic perturbations (Rubin et al., 2019; Tedesco et al., 2021).

Low-throughput multimodal omics were successfully adopted in cancer studies to better understand heterogeneity, complexity, and the evolution of cancer cells by integrating genomic, epigenomic, and transcriptomic features of the same cells (Macaulay et al., 2015; Hou et al., 2016; Bian et al., 2018; Zhu et al., 2021). However, the possibility of obtaining only tens or hundreds of cells provided by low-throughput multimodal omics and their costs have limited their application in TME studies.

Furthermore, given the importance of immunomodulatory proteins (e.g., PD-L1, CTLA-4) in response to cancer immunotherapy, high-throughput multimodal technologies combining transcriptomics and proteomics have been applied in TME studies. In a recent report, CITE-seq successfully discovered new macrophage populations expressing PD-L1 and PD-L2 surface proteins linked to survival in breast cancer (Wu et al., 2021). In other studies, ECCITE-seq and Perturb-CITE-seq were used to define new clinically relevant resistance mechanisms to ICIs in human cancer cell lines and melanoma by exploiting CRISPR-Cas9 screens with multimodal single-cell readouts (Frangieh et al., 2021; Papalexi et al., 2021). Furthermore, the possibility of having both transcript and protein expression from the same cells also helped scientists increase the robustness of their scRNA-seq results. Indeed, they were able to validate the expression of known protein markers and identify novel proteins expressed by immune subsets in different human cancers (Leader et al., 2021; Pombo Antunes et al., 2021). The progress in single cell technologies will guarantee in the near future the possibility to explore at single cell level other omics such as glycomics, lipidomics, metabolomics and microbiomics able to highlight essential cell functions and biological proprieties of tissue components generating a more detailed map of the immune landscape of TME.

Computational Approaches to Analyze Single-Cell Data of the TME

With the widespread use of high-dimensional biological datasets, the scientific community designed specific computational techniques capable of extracting knowledge from complex multi-omic data. Compared to its bulk counterpart, single-cell data analysis is particularly challenging because of the high dimensionality given by the number of cells and markers and the presence of peculiar technical and biological factors that are important to keep in mind. For example, an underestimated problem in single-cell analysis is computational power. As the number of cells increases, all the data analysis steps become more computationally intensive, and some of these steps may require more scalable computational methods and architectures. Machine learning (ML) techniques, especially deep learning (DL), are an emerging class computational methods in single-cell data analysis for their capability to manage complex datasets and the possibility to be implemented into high-parallel architectures (e.g., GPU). Additionally, single-cell data is noisier and more subject to batch effects compared to the bulk counterpart due, for example, to more critical experimental procedures such as single cell isolation, the technical variability in the number of reads and cells sequenced in each sample, and the biological variability caused by heterogeneity in cell composition. TME studies particularly exemplify the latter aspect, in which the presence of cancer and immune and non-immune cells contributes to the complexity of the experiment. Potential confounding factors must be removed or included in statistical models during data analysis. Typical computational steps performed during single-cell data analysis include 1) preprocessing and harmonization steps in which multiple datasets or modalities are combined to perform an integrated analysis after the removal of outliers or low-quality cells, 2) applying dimensionality reducing techniques useful for visualization, 3) clustering, 4) cell annotation, and 5) cellular and molecular functional analysis.

Pre-Processing and Harmonization of Single-Cell Data

Pre-processing is essential before the downstream analysis of single-cell data. This phase encompasses many computational steps ranging from raw to processed data through various types and file formats. Sequencing- and mass spectrometry-based single-cell data involve several steps before quantifying the features-by-cell matrices. Sequencing-based methods have common procedures, including processing raw FASTQ files containing the reads and alignments to the reference genome or transcriptome. Next, for both sequencing-based and mass spectrometry-based approaches, feature detection and quantification (e.g., exons, peaks, and peptides) are common stages to obtain the final features for analysis. Prior to matrix processing and analysis, preliminary quality control is required to assess, for example, the quality of the reads, the percentage of valid barcodes, and reads mapped to the genome. The software suite of specific platforms (e.g., 10x Genomics Cell Ranger) often provides this information. Next, the matrices are reduced by removing low-quality and outlier cells. For example, cells with an unexpectedly high or low number of features detected or with poor quantification are removed from further analyses. In droplet-based single-cell technologies, a filtering step removes cells with a hybrid transcriptome (e.g., doublets), that is, two or more cells incapsulated in the same droplet, using specialized software (McGinnis et al., 2019; Wolock et al., 2019; DePasquale et al., 2020). A specific scRNA-seq pre-processing step involves removing cells with a high percentage of mitochondrial and ribosomal genes expressed because they are usually considered low quality. Thresholds on mitochondrial and ribosomal expression must be chosen carefully, especially when dissecting the immune complexity of the TME, to avoid removing certain cell subsets (Zilionis et al., 2019; Osorio and Cai, 2021; Subramanian et al., 2021).

In FC and CyTOF, pre-processing steps and downstream analyses are typically performed using flow cytometry standard (FCS) files. Typical FC pre-processing steps include data compensation and transformation (e.g., biexponential, generalized Box-Cox) to correct the channel spillover and the effects of outliers and distorted distributions, respectively. Like droplet-based single-cell technologies, an important FC pre-processing step involves removing doublets. This happens when the cytometer cannot discriminate between 2 cells because they pass too closely through the trigger laser. In FC, single cells are differentiated at the beginning of the gating strategy using the 2D plot with Forward Side Channel-Aria (FSC-A) and Forward Side Channel-Height (FSC-H). Cells that do not display a linear correlation of these two parameters are marked as doublets and are excluded from the analysis. In CyTOF, the discrimination of doublets is more complex because the cell size parameters used in FC are not available. Here, DNA intercalators and event length are used to obtain single cells (Gadalla et al., 2019). In spatial proteomics, the pre-processing step involves some adjustment (e.g., illumination and contrast correction) to the images before cell identification through the segmentation process. After segmentation, the cells are classified and quantified for subsequent statistical analysis.

As part of the pre-processing phase, normalization is a fundamental procedure in single-cell data analysis to make all the cells comparable. In sequencing-based single-cell data, the variability of reads sequenced per cell and data sparsity have been carefully considered in computational pipelines to avoid technical effects confounding biological heterogeneity. In scRNA-seq, this has motivated the design of different normalization methods, for example, those based on the estimation of size factors to correct gene expression (Lun et al., 2016; Wolf et al., 2018; Stuart et al., 2019) and others based on linear regression (Bacher et al., 2017; Yip et al., 2017; Hafemeister and Satija, 2019). Text-mining techniques are used to normalize scATAC-seq datasets (Cusanovich et al., 2015; Cusanovich et al., 2018; Fang et al., 2021). Normalization is also an important step in single-cell proteomics data to remove the technical variability due, for example, to differences in instrument performance over acquisition time, particularly in mass cytometry (Rybakowska et al., 2020). Bead-based normalization was introduced to correct these technical artifacts (Finck et al., 2013). This method uses information obtained through standard calibration beads that track changes in the signal over the acquisition time to adjust the marker values. However, technical and biological differences may arise from the disparate technical and biological aspects of data generation in single-cell experiments. A key step in studying the TME through single-cell technologies is the study of changes in the cell subset composition of different individuals. Harmonization techniques have been proposed for single-cell data to facilitate the comparison of multiple samples. These techniques aim to overcome noise due to the variability among cells, individuals, species, and protocols trying to maintain the true biological signals. An important advantage of data harmonization is the possibility to aggregate multiple samples into a single dataset making faster identifying shared or sample-specific cell subsets. This approach makes the analysis simpler and less error prone because it avoids processing each dataset individually or merging them without considering potential biases. A drawback of harmonization is the possibility to lose true biological signals due to the “correction” procedure. In single-cell proteomics, methods have been introduced to correct unwanted variability in FC and mass cytometry, such as aligning the marker intensity distributions across samples or performing cell-type-specific normalization using shared controls across multiple batches (Hahne et al., 2009; Finak et al., 2014; Van Gassen et al., 2020). Several methods ranging from more conventional methods to highly scalable and fast ML approaches have been proposed to align and make cell subsets comparable to gene expression and chromatin accessibility of different datasets, considering diverse technical and biological variation sources (Luecken et al., 2022). The advent of single-cell multimodal omics has posed additional computational challenges in extrapolating useful information from the different layers measured for each cell. This has motivated the introduction of harmonization methods for performing a joint analysis, such as clustering, to exploit the power given by all available cell modalities (Argelaguet et al., 2020; Wang et al., 2020; Gayoso et al., 2021; Hao et al., 2021; Singh et al., 2021; Zuo and Chen, 2021). It is worth to remark that DL approaches (e.g., deep generative models) represent a significant part of harmonization methods for single-cell datasets. Their power is given by their capacity to learn complex mechanisms of biological systems from multiple biological datasets and modalities (Li Y. et al., 2018). This makes them important tools to manage the increasing complexity of single-cell data.

Dimensionality Reduction, Clustering, and Cell Annotation of Single-Cell Data

Single-cell datasets are intrinsically high dimensional. A common step in single-cell data analysis is projecting the data into a low-dimensional space using dimensionality reducing techniques to dissect the complexity and understand the TME’s cellular composition, cell states, and trajectories. The dimensions calculated by these algorithms are useful for visualizing the cell subsets, usually using 2D scatter plots, and as inputs of other computational techniques. Methods such as principal component analysis (PCA) and singular value decomposition (SVD) are commonly used to decompose high-dimensional datasets into several important axes of variation. Owing to their linearity, these methods cannot accurately represent the complex structure of single-cell data. However, the dimensions extracted by these approaches are commonly injected into other dimensionality reducing methods, such as t-distributed stochastic neighbor embedding (t-SNE) and uniform manifold approximation and projection (UMAP), to improve their accuracy (Kobak and Berens, 2019; Baek and Lee, 2020; Kobak and Linderman, 2021). t-SNE and UMAP are currently the most commonly used techniques for visualizing single-cell datasets intuitively. With these methods, each cell is represented as a point in a 2D scatter plot, where the proximity to other cells indicates a similarity in expression profiles. These algorithms can position cells in a biologically meaningful way allowing the correct interpretation of the data. A vital property of these algorithms is to preserve the local and global structure of the data, namely, the distances between points (e.g., cells) within the same cluster and between different clusters. The local topology ensures that the cells in a cluster are homogeneous and represent, for example, the same immune cell subset. In contrast, the global structure of data can provide important biological insights into cell subset relationships, such as monocyte-to-macrophage or epithelial-to-malignant cell transitions, and it is the most debated feature of t-SNE and UMAP. Although previous reports have shown UMAP to be better than t-SNE in preserving the organization of cell clusters (Becht et al., 2018), recent studies suggest that more effective t-SNE parameterizations make it as good as UMAP for conserving global data geometry (Kobak and Berens, 2019; Kobak and Linderman, 2021). Considering the importance of studying the state transition and differentiation trajectories of single cells, other dimensionality reducing approaches have been proposed to better preserve the global features of the data and, consequently, provide biologically relevant information (Weinreb et al., 2018; Moon et al., 2019). After dimensionality reduction, clustering and cell annotation are usually the next steps in single-cell data analysis. Clustering algorithms are mainly used to define homogenous cell subsets to be annotated using specific molecular measures, such as gene expression, proteins, or peak counts. In single-cell multimodal omics data, all the provided molecular measures can be combined to perform an integrated analysis, as the previous section introduced.

Among the large toolboxes of clustering techniques made available by computational sciences, several methods are more commonly used in the single-cell data analysis context. Graph-based clustering has become very popular in single-cell data analysis owing to its minimal required assumptions compared to other well-known techniques, such as k-means. These include methods based on clique detection (Xu and Su, 2015), spectral clustering (Ng et al., 2002), and community detection algorithms (Blondel et al., 2008; Waltman and Van Eck, 2013; Traag et al., 2019). The latter has been highly appreciated in recent years because they can scale the number of cells and be implemented in popular R and Python packages, such as Seurat (Stuart et al., 2019), Monocle (Qiu et al., 2017), and SCANPY (Wolf et al., 2018). In recent years, ML techniques based on self-organizing maps (SOMs) (Kohonen, 1990) have also been introduced for both mass cytometry (Van Gassen et al., 2015) and scRNA-seq data (Camp et al., 2017). These tools usually provide useful visualization features for an intuitive interpretation of expression similarity among cell clusters.

Cluster analysis is usually used for cell annotation following an iterative process in which each cluster is mapped to a biologically relevant cell type by observing the expression of multiple markers. Several clustering algorithms allow us to set a level of granularity (i.e., resolution) based on the level of detail at which the cell subsets are dissected. The use of clustering to annotate cells has several drawbacks. For example, the number of clusters may be overestimated, underestimated, or not reproducible. Furthermore, manual annotation may be error-prone or not correspond to similar annotated cells in the literature. Alternatively, several classification techniques have been introduced in recent years to perform automatic cell annotation (Zhao X. et al., 2020) based on well-annotated reference datasets containing, for example, fluorescence-activated cell sorting (FACS)-sorted cell populations. These methods typically use bulk or single-cell references from large general-purpose (Consortium, 2012; Mabbott et al., 2013; Martens and Stunnenberg, 2013), immune-specific (Heng et al., 2008), or tumor-specific atlases (Rozenblatt-Rosen et al., 2020) to infer the cell types present in dataset by correlating their expression profiles. Most of these tools make feature selections before performing cell classifications. However, the unbiased nature of these approaches makes them critical for improving the reproducibility and consistency among single-cell studies. With the spread of large-scale disease datasets, an important goal will be to integrate these datasets into harmonized and batch-corrected references to be queried efficiently. Scalable architectures based on ML will be valuable tools for integrating human references to study the TME (Lotfollahi et al., 2022).

Analysis of Cell Function and Differentiation of Single-Cell Data

Analyzing changes in the cellular composition of the TME is necessary for understanding the various mechanisms of cancer. However, a deep molecular characterization is crucial for identifying the key drivers of functional cellular changes. Differential expression (DE) analysis is the main statistical technique for detecting functional perturbations caused by changes in gene or protein expression, chromatin accessibility, and genomic aberrations. Although general-purpose statistical tests or DE methods for bulk datasets have been widely used for DE analysis of different single-cell omics datasets, several specialized techniques have been adopted over the years to deal with the heterogeneity and sparsity of scRNA-seq (Kharchenko et al., 2014; Finak et al., 2015). However, DE methods based on pseudo-bulk aggregation of biological replicates have recently gained attention for their capacity to extrapolate more robust results than general-purpose and specialized single-cell DE methods (Squair et al., 2021). After differential expression analysis, a typical task in the study of the TME is to extract a list of biological processes linked to the molecular changes induced by cancer or therapy in different cell subsets. Methods for gene set analysis that have been widely used in recent years for bulk data, such as over-representation, are commonly used in single-cell data. These approaches take into account a list of differentially expressed molecules and gene sets from Gene Ontology (GO) (Gene Ontology, 2015), the Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2016), Reactome (Fabregat et al., 2018), and the Molecular Signature Database (MSigDB) (Liberzon et al., 2011) to identify altered molecular pathways.

In addition to molecular changes concerning pathway regulation, characterizing molecular states between different cells is a common step in single-cell data analysis. While clustering and cell annotation techniques provide information about the cell types present in the experiment, they supply no information on the relationships between cell types and clusters. Recently, several methods have been proposed to reconstruct differentiation trajectory maps by ordering cells based on expression pattern similarity (Saelens et al., 2019) or transcriptional dynamics (La Manno et al., 2018; Bergen et al., 2020). These techniques have been successfully applied to study the developmental trajectories of different cell types and TMEs (Zhang et al., 2019; Chen et al., 2020; Qian et al., 2020; He et al., 2021; Liu et al., 2021). With the improvement in single-cell omics technology performance and their integration ability, we need to address the problem of cell number scalability and the generalizability of these methods to different omics (Stassen et al., 2021).

Pro- and anti-tumor mechanisms are governed by changes driven by complex molecular and cellular interactions within the TME (Anderson and Simon, 2020). Various methods have been introduced in single-cell data analyses to study the intra- and inter-cellular interactome through the inference of gene regulatory networks (GRNs) and ligand-receptor pairs to understand the biological processes underlying these mechanisms. Several methods for GRN inference have been proposed to decipher intracellular networks, including tools originally designed for bulk transcriptomics (Huynh-Thu and Sanguinetti, 2015; Moerman et al., 2019) and techniques specifically designed for single-cell transcriptomics that exploit additional information, such as pseudo-temporal ordering (Matsumoto et al., 2017; Specht and Li, 2017; Deshpande et al., 2021) and information about transcription factors and their targets (Aibar et al., 2017). The first class of methods tries to learn the gene regulatory structure without prior information, and the second class uses pseudo-temporal information to better explain gene regulation during cell differentiation and development. A systematic evaluation of these techniques was recently published (Pratapa et al., 2020). In addition to the intracellular regulation state of a cell, a crucial aspect in understanding the TME involves exploring the cell-to-cell interactome induced by cancer cells (Whiteside, 2008). In spatial technologies, this task can be accomplished by observing cellular co-localization by inspecting appropriate cell markers in histological regions of interest (ROIs) of tissue sections. Without spatial information, cellular interactions must be inferred from the ligand and respective cognate receptor expression levels.

Several computational methods based on different mathematical models have been proposed to identify cell-to-cell interactions. Among others, methods based on the permutation of expression have been widely used (Armingol et al., 2021). These methods typically calculate the communication score of a list of ligand-receptor pairs obtained from curated databases and evaluate the significance of the interactions through cluster label permutation and statistical tests. Ligand-receptor interactions inferred from single-cell transcriptomics may provide interesting hypotheses that need to be further validated using other technologies.

Conclusion and Future Perspectives

The shift from bulk to single-cell sequencing has allowed us to move forward from a general molecular signature in which the contribution of each cell is averaged to the complete molecular fingerprint of each sequenced cell. This is particularly important when complex samples characterized by heterogeneous cell compositions, such as tumor tissues, are analyzed. Accordingly, the employment of single-cell technologies has radically improved the understanding of the cancer framework both quantitatively and qualitatively. Indeed, a single-cell platform can resolve the plasticity of tumor cells and decode tumor phenotypes (invasiveness, stemness, proliferation, and apoptosis), revealing the composition of the TME and the differentiation of immune and stromal cells towards anti- or pro-tumor phenotypes. In addition, single-cell sequencing can track the evolutionary trajectories of neoplastic clones in primary tumors, with results that challenge the original vision of gradual neoplastic evolution (Gao et al., 2016) and improve the understanding of the metastatic spreading process by profiling circulating tumor cells or metastatic lesions (Leung et al., 2017). Similarly, molecular tumor fingerprinting can predict the response to target therapy (Tirosh et al., 2016; Rambow et al., 2018). Although single-cell technologies were initially developed for research purposes and contributed significantly to dissecting cancer evolution mechanisms, clinical settings will soon use them. Investigation at the single-cell level can improve early tumor detection, prognostic biomarker identification, and patient risk stratification, thus supporting a more tumor-tailored therapy. Similarly, dissecting the complexity of the TME can help design the best immunotherapy approach, reverting local immune suppression or empowering the fitness and killing abilities of tumor-infiltrating effector cells. Finally, single-cell sequencing platforms can be employed as potent diagnostic tools for non-invasive monitoring of tumor evolution and patient relapse by profiling circulating tumor cells. This approach can detect, for instance, the appearance of clones resistant to targeted therapy, promptly driving the clinical decision towards an alternative therapeutic solution. However, some limitations need to be overcome to make this technology available for mainstream clinical purposes. First, its usage requires a more complex team, including surgeons, oncologists, pathologists, and researchers working in a fast and coordinated manner and the development of robust tissue processing protocols for primary tumors. Moreover, sample processing requires loss of tissue architecture, whereas spatial single-cell technologies combining molecular and histological information do not guarantee the same resolution of single-cell sequencing on suspension cells. Secondly, scRNA-seq data depend on their intrinsic noisy since eukaryotic transcription does not occur at a persistent basal rate but it takes place in pulses (Chubb et al., 2006). Therefore, a failure to uncover a transcript of a specific gene in a cell at a single time point is an ambiguous result since it can be considered as a result of either permanent gene inactivation or timely limitation of gene transcription detection where the gene is active but the transcript in the time window of the sampling is not present. To avoid possible serious faults, the interpretation of scRNA-seq results should directed on pathways analysis and gene-set enrichment rather than single gene expression. Finally, batch effects can take places when aggregating multiple samples. Several methods of correcting batch effects have been optimized but their use must be balanced against the risk of eclipsing true biological differences. Nonetheless, as recently exemplified by the introduction of bulk NGS, single-cell technology use will soon be extended to patients to support cancer diagnosis and treatment in the near future.

Author Contributions

All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.

Funding

This work was supported by the PRIN programs of the Italian Ministry of Education, University and Research (MIUR, CUP: B38D19000140006, CUP: B39J22001200001) and Fondazione Associazione Italiana per la Ricerca sul Cancro (AIRC, Project: 21509).

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

Figures 1A–D were adapted from the templates “Cancer Cell Mutations Affect Tumor Microenvironment”, “Epigenetic Deregulation in Cancer”, “Cancer Evolution with Treatment”, “ATAC Sequencing”, “Immune Mechanisms Affected by COVID-19″, “TI-1 Antigens Stimulate B cells in Conjuntion with TLR Signal”, “Icon Pack–Immunology” by BioRender.com (2022), “CITE-seq (Cellular Indexing of Transcriptomes & Epitopes by Sequencing) Workflow”. Retrieved from "https://app.biorender.com/biorender-templates"

References

Abel, E. J., Bauman, T. M., Weiker, M., Shi, F., Downs, T. M., Jarrard, D. F., et al. (2014). Analysis and Validation of Tissue Biomarkers for Renal Cell Carcinoma Using Automated High-Throughput Evaluation of Protein Expression. Hum. Pathol. 45 (5), 1092–1099. doi:10.1016/j.humpath.2014.01.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Adey, A., Morrison, H. G., AsanXun, A., Xun, X., Kitzman, J. O., Turner, E. H., et al. (2010). Rapid, Low-Input, Low-Bias Construction of Shotgun Fragment Libraries by High-Density In Vitro Transposition. Genome Biol. 11 (12), R119. doi:10.1186/gb-2010-11-12-r119

PubMed Abstract | CrossRef Full Text | Google Scholar

Aibar, S., González-Blas, C. B., Moerman, T., Huynh-Thu, V. A., Imrichova, H., Hulselmans, G., et al. (2017). SCENIC: Single-Cell Regulatory Network Inference and Clustering. Nat. Methods 14 (11), 1083–1086. doi:10.1038/nmeth.4463

PubMed Abstract | CrossRef Full Text | Google Scholar

Anderson, N. M., and Simon, M. C. (2020). The Tumor Microenvironment. Curr. Biol. 30 (16), R921–R925. doi:10.1016/j.cub.2020.06.081

PubMed Abstract | CrossRef Full Text | Google Scholar

Angelo, M., Bendall, S. C., Finck, R., Hale, M. B., Hitzman, C., Borowsky, A. D., et al. (2014). Multiplexed Ion Beam Imaging of Human Breast Tumors. Nat. Med. 20 (4), 436–442. doi:10.1038/nm.3488

PubMed Abstract | CrossRef Full Text | Google Scholar

Angermueller, C., Clark, S. J., Lee, H. J., Macaulay, I. C., Teng, M. J., Hu, T. X., et al. (2016). Parallel Single-Cell Sequencing Links Transcriptional and Epigenetic Heterogeneity. Nat. Methods 13 (3), 229–232. doi:10.1038/nmeth.3728

PubMed Abstract | CrossRef Full Text | Google Scholar

Argelaguet, R., Arnol, D., Bredikhin, D., Deloro, Y., Velten, B., Marioni, J. C., et al. (2020). MOFA+: a Statistical Framework for Comprehensive Integration of Multi-Modal Single-Cell Data. Genome Biol. 21 (1), 111. doi:10.1186/s13059-020-02015-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Armingol, E., Officer, A., Harismendy, O., and Lewis, N. E. (2021). Deciphering Cell-Cell Interactions and Communication from Gene Expression. Nat. Rev. Genet. 22 (2), 71–88. doi:10.1038/s41576-020-00292-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Bacher, R., Chu, L.-F., Leng, N., Gasch, A. P., Thomson, J. A., Stewart, R. M., et al. (2017). SCnorm: Robust Normalization of Single-Cell RNA-Seq Data. Nat. Methods 14 (6), 584–586. doi:10.1038/nmeth.4263

PubMed Abstract | CrossRef Full Text | Google Scholar

Baek, S., and Lee, I. (2020). Single-cell ATAC Sequencing Analysis: From Data Preprocessing to Hypothesis Generation. Comput. Struct. Biotechnol. J. 18, 1429–1439. doi:10.1016/j.csbj.2020.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Balachandran, V. P., Luksza, M., Łuksza, M., Zhao, J. N., Makarov, V., Moral, J. A., et al. (2017). Identification of Unique Neoantigen Qualities in Long-Term Survivors of Pancreatic Cancer. Nature 551 (7681), 512–516. doi:10.1038/nature24462

PubMed Abstract | CrossRef Full Text | Google Scholar

Bandura, D. R., Baranov, V. I., Ornatsky, O. I., Antonov, A., Kinach, R., Lou, X., et al. (2009). Mass Cytometry: Technique for Real Time Single Cell Multitarget Immunoassay Based on Inductively Coupled Plasma Time-Of-Flight Mass Spectrometry. Anal. Chem. 81 (16), 6813–6822. doi:10.1021/ac901049w

PubMed Abstract | CrossRef Full Text | Google Scholar

Baslan, T., Kendall, J., Volyanskyy, K., McNamara, K., Cox, H., D'Italia, S., et al. (2020). Novel Insights into Breast Cancer Copy Number Genetic Heterogeneity Revealed by Single-Cell Genome Sequencing. Elife 9. doi:10.7554/eLife.51480

CrossRef Full Text | Google Scholar

Baylin, S. B., and Jones, P. A. (2011). A Decade of Exploring the Cancer Epigenome - Biological and Translational Implications. Nat. Rev. Cancer 11 (10), 726–734. doi:10.1038/nrc3130

PubMed Abstract | CrossRef Full Text | Google Scholar

Becht, E., McInnes, L., Healy, J., Dutertre, C.-A., Kwok, I. W. H., Ng, L. G., et al. (2018). Dimensionality Reduction for Visualizing Single-Cell Data Using UMAP. Nat. Biotechnol. 37, 38–44. doi:10.1038/nbt.4314

PubMed Abstract | CrossRef Full Text | Google Scholar

Benci, J. L., Xu, B., Qiu, Y., Wu, T. J., Dada, H., Twyman-Saint Victor, C., et al. (2016). Tumor Interferon Signaling Regulates a Multigenic Resistance Program to Immune Checkpoint Blockade. Cell 167 (6), 1540–1554. doi:10.1016/j.cell.2016.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Bendall, S. C., Nolan, G. P., Roederer, M., and Chattopadhyay, P. K. (2012). A Deep Profiler's Guide to Cytometry. Trends Immunol. 33 (7), 323–332. doi:10.1016/j.it.2012.02.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Bendall, S. C., Simonds, E. F., Qiu, P., Amir, E.-a. D., Krutzik, P. O., Finck, R., et al. (2011). Single-cell Mass Cytometry of Differential Immune and Drug Responses across a Human Hematopoietic Continuum. Science 332 (6030), 687–696. doi:10.1126/science.1198704

PubMed Abstract | CrossRef Full Text | Google Scholar

Bergen, V., Lange, M., Peidli, S., Wolf, F. A., and Theis, F. J. (2020). Generalizing RNA Velocity to Transient Cell States through Dynamical Modeling. Nat. Biotechnol. 38 (12), 1408–1414. doi:10.1038/s41587-020-0591-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Betancur, P. A., Abraham, B. J., Yiu, Y. Y., Willingham, S. B., Khameneh, F., Zarnegar, M., et al. (2017). A CD47-Associated Super-enhancer Links Pro-inflammatory Signalling to CD47 Upregulation in Breast Cancer. Nat. Commun. 8, 14802. doi:10.1038/ncomms14802

PubMed Abstract | CrossRef Full Text | Google Scholar

Bian, S., Hou, Y., Zhou, X., Li, X., Yong, J., Wang, Y., et al. (2018). Single-cell Multiomics Sequencing and Analyses of Human Colorectal Cancer. Science 362 (6418), 1060–1063. doi:10.1126/science.aao3791

PubMed Abstract | CrossRef Full Text | Google Scholar

Bindea, G., Mlecnik, B., Tosolini, M., Kirilovsky, A., Waldner, M., Obenauf, A. C., et al. (2013). Spatiotemporal Dynamics of Intratumoral Immune Cells Reveal the Immune Landscape in Human Cancer. Immunity 39 (4), 782–795. doi:10.1016/j.immuni.2013.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (5), 541–550. doi:10.1038/s41591-018-0014-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Blando, J., Sharma, A., Higa, M. G., Zhao, H., Vence, L., Yadav, S. S., et al. (2019). Comparison of Immune Infiltrates in Melanoma and Pancreatic Cancer Highlights VISTA as a Potential Target in Pancreatic Cancer. Proc. Natl. Acad. Sci. U.S.A. 116 (5), 1692–1697. doi:10.1073/pnas.1811067116

PubMed Abstract | CrossRef Full Text | Google Scholar

Blondel, V. D., Guillaume, J.-L., Lambiotte, R., and Lefebvre, E. (2008). Fast Unfolding of Communities in Large Networks. J. Stat. Mech. 2008 (10), P10008–P15468. doi:10.1088/1742-5468/2008/10/P10008

CrossRef Full Text | Google Scholar

Bonilla, D. L., Reinin, G., and Chua, E. (2020). Full Spectrum Flow Cytometry as a Powerful Technology for Cancer Immunotherapy Research. Front. Mol. Biosci. 7, 612801. doi:10.3389/fmolb.2020.612801

PubMed Abstract | CrossRef Full Text | Google Scholar

Bronte, V., Cingarlini, S., Marigo, I., De Santo, C., Gallina, G., Dolcetti, L., et al. (2006). Leukocyte Infiltration in Cancer Creates an Unfavorable Environment for Antitumor Immune Responses: a Novel Target for Therapeutic Intervention. Immunol. Investig. 35 (3-4), 327–357. doi:10.1080/08820130600754994

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, R., Curry, E., Magnani, L., Wilhelm-Benartzi, C. S., and Borley, J. (2014). Poised Epigenetic States and Acquired Drug Resistance in Cancer. Nat. Rev. Cancer 14 (11), 747–753. doi:10.1038/nrc3819

PubMed Abstract | CrossRef Full Text | Google Scholar

Budnik, B., Levy, E., Harmange, G., and Slavov, N. (2018). SCoPE-MS: Mass Spectrometry of Single Mammalian Cells Quantifies Proteome Heterogeneity during Cell Differentiation. Genome Biol. 19 (1), 161. doi:10.1186/s13059-018-1547-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Burr, M. L., Sparbier, C. E., Chan, K. L., Chan, Y.-C., Kersbergen, A., Lam, E. Y. N., et al. (2019). An Evolutionarily Conserved Function of Polycomb Silences the MHC Class I Antigen Presentation Pathway and Enables Immune Evasion in Cancer. Cancer Cell 36 (4), 385–401. doi:10.1016/j.ccell.2019.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Calcinotto, A., Spataro, C., Zagato, E., Di Mitri, D., Gil, V., Crespo, M., et al. (2018). IL-23 Secreted by Myeloid Cells Drives Castration-Resistant Prostate Cancer. Nature 559 (7714), 363–369. doi:10.1038/s41586-018-0266-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Calvanese, V., Capellera-Garcia, S., Ma, F., Fares, I., Liebscher, S., Ng, E. S., et al. (2022). Mapping Human Haematopoietic Stem Cells from Haemogenic Endothelium to Birth. Nature 604, 534–540. doi:10.1038/s41586-022-04571-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Camp, J. G., Sekine, K., Gerber, T., Loeffler-Wirth, H., Binder, H., Gac, M., et al. (2017). Multilineage Communication Regulates Human Liver Bud Development from Pluripotency. Nature 546 (7659), 533–538. doi:10.1038/nature22796

PubMed Abstract | CrossRef Full Text | Google Scholar

Cancer Genome Atlas, N. (2012a). Comprehensive Molecular Characterization of Human Colon and Rectal Cancer. Nature 487 (7407), 330–337. doi:10.1038/nature11252

CrossRef Full Text | Google Scholar

Cancer Genome Atlas, N. (2012b). Comprehensive Molecular Portraits of Human Breast Tumours. Nature 490 (7418), 61–70. doi:10.1038/nature11412

CrossRef Full Text | Google Scholar

Cancer Genome Atlas Research, N. (2012). Comprehensive Genomic Characterization of Squamous Cell Lung Cancers. Nature 489 (7417), 519–525. doi:10.1038/nature11404

PubMed Abstract | CrossRef Full Text | Google Scholar

Cancer Genome Atlas Research, N. (2014a). Comprehensive Molecular Characterization of Gastric Adenocarcinoma. Nature 513 (7517), 202–209. doi:10.1038/nature13480

CrossRef Full Text | Google Scholar

Cancer Genome Atlas Research, N. (2014b). Comprehensive Molecular Profiling of Lung Adenocarcinoma. Nature 511 (7511), 543–550. doi:10.1038/nature13385

CrossRef Full Text | Google Scholar

Cancer Genome Atlas Research, N. (2015). The Molecular Taxonomy of Primary Prostate Cancer. Cell 163 (4), 1011–1025. doi:10.1016/j.cell.2015.10.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Cancer Genome Atlas Research, Network (2017). Analytical BiologicalIntegrated Genomic and Molecular Characterization of Cervical Cancer. Nature 543 (7645), 378–384. doi:10.1038/nature21386

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Xing, D., Tan, L., Li, H., Zhou, G., Huang, L., et al. (2017). Single-cell Whole-Genome Analyses by Linear Amplification via Transposon Insertion (LIANTI). Science 356 (6334), 189–194. doi:10.1126/science.aak9787

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, G., Ning, B., and Shi, T. (2019a). Single-Cell RNA-Seq Technologies and Related Computational Data Analysis. Front. Genet. 10, 317. doi:10.3389/fgene.2019.00317

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Lake, B. B., and Zhang, K. (2019b). High-throughput Sequencing of the Transcriptome and Chromatin Accessibility in the Same Cell. Nat. Biotechnol. 37 (12), 1452–1457. doi:10.1038/s41587-019-0290-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Z., Zhou, L., Liu, L., Hou, Y., Xiong, M., Yang, Y., et al. (2020). Single-cell RNA Sequencing Highlights the Role of Inflammatory Cancer-Associated Fibroblasts in Bladder Urothelial Carcinoma. Nat. Commun. 11 (1), 5077. doi:10.1038/s41467-020-18916-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, S., Li, Z., Gao, R., Xing, B., Gao, Y., Yang, Y., et al. (2021). A Pan-Cancer Single-Cell Transcriptional Atlas of Tumor Infiltrating Myeloid Cells. Cell 184 (3), 792–809. doi:10.1016/j.cell.2021.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Chubb, J. R., Trcek, T., Shenoy, S. M., and Singer, R. H. (2006). Transcriptional Pulsing of a Developmental Gene. Curr. Biol. 16 (10), 1018–1025. doi:10.1016/j.cub.2006.03.092

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciriello, G., Gatza, M. L., Beck, A. H., Wilkerson, M. D., Rhie, S. K., Pastore, A., et al. (2015). Comprehensive Molecular Portraits of Invasive Lobular Breast Cancer. Cell 163 (2), 506–519. doi:10.1016/j.cell.2015.09.033

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (1), 781. doi:10.1038/s41467-018-03149-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Consortium, E. P. (2012). An Integrated Encyclopedia of DNA Elements in the Human Genome. Nature 489 (7414), 57–74. doi:10.1038/nature11247

PubMed Abstract | CrossRef Full Text | Google Scholar

Consortium, I. T. P.-C. A. o. W. G. (2020). Pan-cancer Analysis of Whole Genomes. Nature 578 (7793), 82–93. doi:10.1038/s41586-020-1969-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, C., Shu, W., and Li, P. (2016). Fluorescence In Situ Hybridization: Cell-Based Genetic Diagnostic and Research Applications. Front. Cell Dev. Biol. 4, 89. doi:10.3389/fcell.2016.00089

PubMed Abstract | CrossRef Full Text | Google Scholar

Cusanovich, D. A., Daza, R., Adey, A., Pliner, H. A., Christiansen, L., Gunderson, K. L., et al. (2015). Multiplex Single-Cell Profiling of Chromatin Accessibility by Combinatorial Cellular Indexing. Science 348 (6237), 910–914. doi:10.1126/science.aab1601

PubMed Abstract | CrossRef Full Text | Google Scholar

Cusanovich, D. A., Hill, A. J., Aghamirzaie, D., Daza, R. M., Pliner, H. A., Berletch, J. B., et al. (2018). A Single-Cell Atlas of In Vivo Mammalian Chromatin Accessibility. Cell 174 (5), 1309–1324. e1318. doi:10.1016/j.cell.2018.06.052

PubMed Abstract | CrossRef Full Text | Google Scholar

Dangaj, D., Bruand, M., Grimm, A. J., Ronet, C., Barras, D., Duttagupta, P. A., et al. (2019). Cooperation between Constitutive and Inducible Chemokines Enables T Cell Engraftment and Immune Attack in Solid Tumors. Cancer Cell 35 (6), 885–900. doi:10.1016/j.ccell.2019.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Datlinger, P., Rendeiro, A. F., Schmidl, C., Krausgruber, T., Traxler, P., Klughammer, J., et al. (2017). Pooled CRISPR Screening with Single-Cell Transcriptome Readout. Nat. Methods 14 (3), 297–301. doi:10.1038/nmeth.4177

PubMed Abstract | CrossRef Full Text | Google Scholar

Davidson, S., Efremova, M., Riedel, A., Mahata, B., Pramanik, J., Huuhtanen, J., et al. (2020). Single-Cell RNA Sequencing Reveals a Dynamic Stromal Niche that Supports Tumor Growth. Cell Rep. 31 (7), 107628. doi:10.1016/j.celrep.2020.107628

PubMed Abstract | CrossRef Full Text | Google Scholar

De Sanctis, F., Bronte, V., and Ugel, S. (2016). Tumor-Induced Myeloid-Derived Suppressor Cells. Microbiol. Spectr. 4 (3). doi:10.1128/microbiolspec.MCHD-0016-2015

CrossRef Full Text | Google Scholar

De Sanctis, F., Lamolinara, A., Boschi, F., Musiu, C., Caligola, S., Trovato, R., et al. (2022). Interrupting the Nitrosative Stress Fuels Tumor-specific Cytotoxic T Lymphocytes in Pancreatic Cancer. J. Immunother. Cancer 10 (1), e003549. doi:10.1136/jitc-2021-003549

PubMed Abstract | CrossRef Full Text | Google Scholar

Dean, F. B., Nelson, J. R., Giesler, T. L., and Lasken, R. S. (2001). Rapid Amplification of Plasmid and Phage DNA Using Phi29 DNA Polymerase and Multiply-Primed Rolling Circle Amplification. Genome Res. 11 (6), 1095–1099. doi:10.1101/gr.180501

PubMed Abstract | CrossRef Full Text | Google Scholar

DePasquale, E. A. K., Schnell, D., Chetal, K., and Salomonis, N. (2020). Protocol for Identification and Removal of Doublets with DoubletDecon. Star. Protoc. 1 (2), 100085. doi:10.1016/j.xpro.2020.100085

PubMed Abstract | CrossRef Full Text | Google Scholar

Deshpande, A., Chu, L.-F., Stewart, R., and Gitter, A. (2021). Network Inference with Granger Causality Ensembles on Single-Cell Transcriptomic Data. Cell Rep. 38, 534834. doi:10.1016/j.celrep.2022.110333

CrossRef Full Text | Google Scholar

Dey, S. S., Kester, L., Spanjaard, B., Bienko, M., and van Oudenaarden, A. (2015). Integrated Genome and Transcriptome Sequencing of the Same Cell. Nat. Biotechnol. 33 (3), 285–289. doi:10.1038/nbt.3129

PubMed Abstract | CrossRef Full Text | Google Scholar

Di Pilato, M., Kfuri-Rubens, R., Pruessmann, J. N., Ozga, A. J., Messemaker, M., Cadilha, B. L., et al. (2021). CXCR6 Positions Cytotoxic T Cells to Receive Critical Survival Signals in the Tumor Microenvironment. Cell 184 (17), 4512–4530. e4522. doi:10.1016/j.cell.2021.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, Q., Zhang, H., Zheng, J., and Zhang, L. (2020). Turning Cold into Hot: Firing up the Tumor Microenvironment. Trends Cancer 6 (7), 605–618. doi:10.1016/j.trecan.2020.02.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Elyada, E., Bolisetty, M., Laise, P., Flynn, W. F., Courtois, E. T., Burkhart, R. A., et al. (2019). Cross-Species Single-Cell Analysis of Pancreatic Ductal Adenocarcinoma Reveals Antigen-Presenting Cancer-Associated Fibroblasts. Cancer Discov. 9 (8), 1102–1123. doi:10.1158/2159-8290.CD-19-0094

PubMed Abstract | CrossRef Full Text | Google Scholar

Fabregat, A., Jupe, S., Matthews, L., Sidiropoulos, K., Gillespie, M., Garapati, P., et al. (2018). The Reactome Pathway Knowledgebase. Nucleic Acids Res. 46 (D1), D649–D655. doi:10.1093/nar/gkx1132

PubMed Abstract | CrossRef Full Text | Google Scholar

Facciponte, J. G., Ugel, S., De Sanctis, F., Li, C., Wang, L., Nair, G., et al. (2014). Tumor Endothelial Marker 1-specific DNA Vaccination Targets Tumor Vasculature. J. Clin. Invest. 124 (4), 1497–1511. doi:10.1172/JCI67382

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, X., Yang, C., Li, W., Bai, X., Zhou, X., Xie, H., et al. (2021). SMOOTH-seq: Single-Cell Genome Sequencing of Human Cells on a Third-Generation Sequencing Platform. Genome Biol. 22 (1), 195. doi:10.1186/s13059-021-02406-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, R., Preissl, S., Li, Y., Hou, X., Lucero, J., Wang, X., et al. (2021). Comprehensive Analysis of Single Cell ATAC-Seq Data with SnapATAC. Nat. Commun. 12 (1), 1337. doi:10.1038/s41467-021-21583-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Feinberg, A. P., Koldobskiy, M. A., and Göndör, A. (2016). Epigenetic Modulators, Modifiers and Mediators in Cancer Aetiology and Progression. Nat. Rev. Genet. 17 (5), 284–299. doi:10.1038/nrg.2016.13

PubMed Abstract | CrossRef Full Text | Google Scholar

Finak, G., Frelinger, J., Jiang, W., Newell, E. W., Ramey, J., Davis, M. M., et al. (2014). OpenCyto: an Open Source Infrastructure for Scalable, Robust, Reproducible, and Automated, End-To-End Flow Cytometry Data Analysis. PLoS Comput. Biol. 10 (8), e1003806. doi:10.1371/journal.pcbi.1003806

PubMed Abstract | CrossRef Full Text | Google Scholar

Finak, G., McDavid, A., Yajima, M., Deng, J., Gersuk, V., Shalek, A. K., et al. (2015). MAST: a Flexible Statistical Framework for Assessing Transcriptional Changes and Characterizing Heterogeneity in Single-Cell RNA Sequencing Data. Genome Biol. 16, 278. doi:10.1186/s13059-015-0844-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Finck, R., Simonds, E. F., Jager, A., Krishnaswamy, S., Sachs, K., Fantl, W., et al. (2013). Normalization of Mass Cytometry Data with Bead Standards. Cytometry 83A (5), 483–494. doi:10.1002/cyto.a.22271

PubMed Abstract | CrossRef Full Text | Google Scholar

Finotello, F., and Eduati, F. (2018). Multi-Omics Profiling of the Tumor Microenvironment: Paving the Way to Precision Immuno-Oncology. Front. Oncol. 8, 430. doi:10.3389/fonc.2018.00430

PubMed Abstract | CrossRef Full Text | Google Scholar

Fiore, A., Ugel, S., De Sanctis, F., Sandri, S., Fracasso, G., Trovato, R., et al. (2018). Induction of Immunosuppressive Functions and NF-Κb by FLIP in Monocytes. Nat. Commun. 9 (1), 5193. doi:10.1038/s41467-018-07654-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Frangieh, C. J., Melms, J. C., Thakore, P. I., Geiger-Schuller, K. R., Ho, P., Luoma, A. M., et al. (2021). Multimodal Pooled Perturb-CITE-Seq Screens in Patient Models Define Mechanisms of Cancer Immune Evasion. Nat. Genet. 53 (3), 332–341. doi:10.1038/s41588-021-00779-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Fridman, W. H., Zitvogel, L., Sautès–Fridman, C., and Kroemer, G. (2017). The Immune Contexture in Cancer Prognosis and Treatment. Nat. Rev. Clin. Oncol. 14 (12), 717–734. doi:10.1038/nrclinonc.2017.101

PubMed Abstract | CrossRef Full Text | Google Scholar

Futreal, P. A., Coin, L., Marshall, M., Down, T., Hubbard, T., Wooster, R., et al. (2004). A Census of Human Cancer Genes. Nat. Rev. Cancer 4 (3), 177–183. doi:10.1038/nrc1299

PubMed Abstract | CrossRef Full Text | Google Scholar

Gadalla, R., Noamani, B., MacLeod, B. L., Dickson, R. J., Guo, M., Xu, W., et al. (2019). Validation of CyTOF against Flow Cytometry for Immunological Studies and Monitoring of Human Cancer Clinical Trials. Front. Oncol. 9, 415. doi:10.3389/fonc.2019.00415

PubMed Abstract | CrossRef Full Text | Google Scholar

Galon, J., Costes, A., Sanchez-Cabo, F., Kirilovsky, A., Mlecnik, B., Lagorce-Pagès, C., et al. (2006). Type, Density, and Location of Immune Cells within Human Colorectal Tumors Predict Clinical Outcome. Science 313 (5795), 1960–1964. doi:10.1126/science.1129139

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, R., Davis, A., McDonald, T. O., Sei, E., Shi, X., Wang, Y., et al. (2016). Punctuated Copy Number Evolution and Clonal Stasis in Triple-Negative Breast Cancer. Nat. Genet. 48 (10), 1119–1130. doi:10.1038/ng.3641

PubMed Abstract | CrossRef Full Text | Google Scholar

Gayoso, A., Steier, Z., Lopez, R., Regier, J., Nazor, K. L., Streets, A., et al. (2021). Joint Probabilistic Modeling of Single-Cell Multi-Omic Data with totalVI. Nat. Methods 18 (3), 272–282. doi:10.1038/s41592-020-01050-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gene Ontology, C. (2015). Gene Ontology Consortium: Going Forward. Nucleic Acids Res. 43, D1049–D1056. doi:10.1093/nar/gku1179

PubMed Abstract | CrossRef Full Text | Google Scholar

Gentles, A. J., Newman, A. M., Liu, C. L., Bratman, S. V., Feng, W., Kim, D., et al. (2015). The Prognostic Landscape of Genes and Infiltrating Immune Cells across Human Cancers. Nat. Med. 21 (8), 938–945. doi:10.1038/nm.3909

PubMed Abstract | CrossRef Full Text | Google Scholar

Gerdes, M. J., Sevinsky, C. J., Sood, A., Adak, S., Bello, M. O., Bordwell, A., et al. (2013). Highly Multiplexed Single-Cell Analysis of Formalin-Fixed, Paraffin-Embedded Cancer Tissue. Proc. Natl. Acad. Sci. U.S.A. 110 (29), 11982–11987. doi:10.1073/pnas.1300136110

PubMed Abstract | CrossRef Full Text | Google Scholar

Giladi, A., and Amit, I. (2017). Immunology, One Cell at a Time. Nature 547 (7661), 27–29. doi:10.1038/547027a

PubMed Abstract | CrossRef Full Text | Google Scholar

Gohil, S. H., Iorgulescu, J. B., Braun, D. A., Keskin, D. B., and Livak, K. J. (2021). Applying High-Dimensional Single-Cell Technologies to the Analysis of Cancer Immunotherapy. Nat. Rev. Clin. Oncol. 18 (4), 244–256. doi:10.1038/s41571-020-00449-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Goltsev, Y., Samusik, N., Kennedy-Darling, J., Bhate, S., Hale, M., Vazquez, G., et al. (2018). Deep Profiling of Mouse Splenic Architecture with CODEX Multiplexed Imaging. Cell 174 (4), 968–981. e915. doi:10.1016/j.cell.2018.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorris, M. A. J., Halilovic, A., Rabold, K., van Duffelen, A., Wickramasinghe, I. N., Verweij, D., et al. (2018). Eight-Color Multiplex Immunohistochemistry for Simultaneous Detection of Multiple Immune Checkpoint Molecules within the Tumor Microenvironment. J. I. 200 (1), 347–354. doi:10.4049/jimmunol.1701262

PubMed Abstract | CrossRef Full Text | Google Scholar

Grosselin, K., Durand, A., Marsolier, J., Poitou, A., Marangoni, E., Nemati, F., et al. (2019). High-throughput Single-Cell ChIP-Seq Identifies Heterogeneity of Chromatin States in Breast Cancer. Nat. Genet. 51 (6), 1060–1066. doi:10.1038/s41588-019-0424-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, F., Li, L., Li, J., Wu, X., Hu, B., Zhu, P., et al. (2017). Single-cell Multi-Omics Sequencing of Mouse Early Embryos and Embryonic Stem Cells. Cell Res. 27 (8), 967–988. doi:10.1038/cr.2017.82

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (7), 978–985. doi:10.1038/s41591-018-0045-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hafemeister, C., and Satija, R. (2019). Normalization and Variance Stabilization of Single-Cell RNA-Seq Data Using Regularized Negative Binomial Regression. Genome Biol. 20 (1), 296. doi:10.1186/s13059-019-1874-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hahne, F., LeMeur, N., Brinkman, R. R., Ellis, B., Haaland, P., Sarkar, D., et al. (2009). flowCore: a Bioconductor Package for High Throughput Flow Cytometry. BMC Bioinforma. 10, 106. doi:10.1186/1471-2105-10-106

PubMed Abstract | CrossRef Full Text | Google Scholar

Hanahan, D., and Weinberg, R. A. (2011). Hallmarks of Cancer: the Next Generation. Cell 144 (5), 646–674. doi:10.1016/j.cell.2011.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Hao, Y., Hao, S., Andersen-Nissen, E., Mauck, W. M., Zheng, S., Butler, A., et al. (2021). Integrated Analysis of Multimodal Single-Cell Data. Cell 184 (13), 3573–3587. doi:10.1016/j.cell.2021.04.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Hartmann, F. J., Babdor, J., Gherardini, P. F., Amir, E.-A. D., Jones, K., Sahaf, B., et al. (2019). Comprehensive Immune Monitoring of Clinical Trials to Advance Human Immunotherapy. Cell Rep. 28 (3), 819–831. e814. doi:10.1016/j.celrep.2019.06.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayashi, T., Ozaki, H., Sasagawa, Y., Umeda, M., Danno, H., and Nikaido, I. (2018). Single-cell Full-Length Total RNA Sequencing Uncovers Dynamics of Recursive Splicing and Enhancer RNAs. Nat. Commun. 9 (1), 619. doi:10.1038/s41467-018-02866-0

PubMed Abstract | CrossRef Full Text | Google Scholar

He, D., Wang, D., Lu, P., Yang, N., Xue, Z., Zhu, X., et al. (2021). Single-cell RNA Sequencing Reveals Heterogeneous Tumor and Immune Cell Populations in Early-Stage Lung Adenocarcinomas Harboring EGFR Mutations. Oncogene 40 (2), 355–368. doi:10.1038/s41388-020-01528-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Heng, T. S. P., Painter, M. W., Painter, M. W., Elpek, K., Lukacs-Kornek, V., Mauermann, N., et al. (2008). The Immunological Genome Project: Networks of Gene Expression in Immune Cells. Nat. Immunol. 9 (10), 1091–1094. doi:10.1038/ni1008-1091

PubMed Abstract | CrossRef Full Text | Google Scholar

Hofer, F., Di Sario, G., Musiu, C., Sartoris, S., De Sanctis, F., and Ugel, S. (2021). A Complex Metabolic Network Confers Immunosuppressive Functions to Myeloid-Derived Suppressor Cells (MDSCs) within the Tumour Microenvironment. Cells 10 (10), 2700. doi:10.3390/cells10102700

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (3), 304–319. doi:10.1038/cr.2016.23

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Y., Huang, K., An, Q., Du, G., Hu, G., Xue, J., et al. (2016). Simultaneous Profiling of Transcriptome and DNA Methylome from a Single Cell. Genome Biol. 17, 88. doi:10.1186/s13059-016-0950-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Huynh-Thu, V. A., and Sanguinetti, G. (2015). Combining Tree-Based and Dynamical Systems for the Inference of Gene Regulatory Networks. Bioinformatics 31 (10), 1614–1622. doi:10.1093/bioinformatics/btu863

PubMed Abstract | CrossRef Full Text | Google Scholar

Jackson, H. W., Fischer, J. R., Zanotelli, V. R. T., Ali, H. R., Mechera, R., Soysal, S. D., et al. (2020). The Single-Cell Pathology Landscape of Breast Cancer. Nature 578 (7796), 615–620. doi:10.1038/s41586-019-1876-x

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (4), 984–997. doi:10.1016/j.cell.2018.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, P. A., Issa, J.-P. J., and Baylin, S. (2016). Targeting the Cancer Epigenome for Therapy. Nat. Rev. Genet. 17 (10), 630–641. doi:10.1038/nrg.2016.93

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanehisa, M., Sato, Y., Kawashima, M., Furumichi, M., and Tanabe, M. (2016). KEGG as a Reference Resource for Gene and Protein Annotation. Nucleic Acids Res. 44 (D1), D457–D462. doi:10.1093/nar/gkv1070

PubMed Abstract | CrossRef Full Text | Google Scholar

Kashima, Y., Togashi, Y., Fukuoka, S., Kamada, T., Irie, T., Suzuki, A., et al. (2021). Potentiality of Multiple Modalities for Single-Cell Analyses to Evaluate the Tumor Microenvironment in Clinical Specimens. Sci. Rep. 11 (1), 341. doi:10.1038/s41598-020-79385-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaya-Okur, H. S., Wu, S. J., Codomo, C. A., Pledger, E. S., Bryson, T. D., Henikoff, J. G., et al. (2019). CUT&Tag for Efficient Epigenomic Profiling of Small Samples and Single Cells. Nat. Commun. 10 (1), 1930. doi:10.1038/s41467-019-09982-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Keenan, T. E., Burke, K. P., and Van Allen, E. M. (2019). Genomic Correlates of Response to Immune Checkpoint Blockade. Nat. Med. 25 (3), 389–402. doi:10.1038/s41591-019-0382-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Keren, L., Bosse, M., Marquez, D., Angoshtari, R., Jain, S., Varma, S., et al. (2018). A Structured Tumor-Immune Microenvironment in Triple Negative Breast Cancer Revealed by Multiplexed Ion Beam Imaging. Cell 174 (6), 1373–1387. doi:10.1016/j.cell.2018.08.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Kharchenko, P. V., Silberstein, L., and Scadden, D. T. (2014). Bayesian Approach to Single-Cell Differential Expression Analysis. Nat. Methods 11 (7), 740–742. doi:10.1038/nmeth.2967

PubMed Abstract | CrossRef Full Text | Google Scholar

Kieffer, Y., Hocine, H. R., Gentric, G., Pelon, F., Bernard, C., Bourachot, B., et al. (2020). Single-Cell Analysis Reveals Fibroblast Clusters Linked to Immunotherapy Resistance in Cancer. Cancer Discov. 10 (9), 1330–1351. doi:10.1158/2159-8290.CD-19-1384

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, N., Kim, H. K., Lee, K., Hong, Y., Cho, J. H., Choi, J. W., et al. (2020). Single-cell RNA Sequencing Demonstrates the Molecular and Cellular Reprogramming of Metastatic Lung Adenocarcinoma. Nat. Commun. 11 (1), 2285. doi:10.1038/s41467-020-16164-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Klein, A. M., Mazutis, L., Akartuna, I., Tallapragada, N., Veres, A., Li, V., et al. (2015). Droplet Barcoding for Single-Cell Transcriptomics Applied to Embryonic Stem Cells. Cell 161 (5), 1187–1201. doi:10.1016/j.cell.2015.04.044

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobak, D., and Berens, P. (2019). The Art of Using T-SNE for Single-Cell Transcriptomics. Nat. Commun. 10 (1), 5416. doi:10.1038/s41467-019-13056-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kobak, D., and Linderman, G. C. (2021). Initialization Is Critical for Preserving Global Data Structure in Both T-SNE and UMAP. Nat. Biotechnol. 39 (2), 156–157. doi:10.1038/s41587-020-00809-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Kohonen, T. (1990). The Self-Organizing Map. Proc. IEEE 78 (9), 1464–1480. doi:10.1109/5.58325

CrossRef Full Text | Google Scholar

La Manno, G., Soldatov, R., Zeisel, A., Braun, E., Hochgerner, H., Petukhov, V., et al. (2018). RNA Velocity of Single Cells. Nature 560 (7719), 494–498. doi:10.1038/s41586-018-0414-6

PubMed Abstract | CrossRef Full Text | Google Scholar

LaFave, L. M., Kartha, V. K., Ma, S., Meli, K., Del Priore, I., Lareau, C., et al. (2020). Epigenomic State Transitions Characterize Tumor Progression in Mouse Lung Adenocarcinoma. Cancer Cell 38 (2), 212–228. e213. doi:10.1016/j.ccell.2020.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Laks, E., McPherson, A., Zahn, H., Lai, D., Steif, A., Brimhall, J., et al. (2019). Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing. Cell 179 (5), 1207–1221. e1222. doi:10.1016/j.cell.2019.10.026

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (8), 1277–1289. doi:10.1038/s41591-018-0096-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Lan, F., Demaree, B., Ahmed, N., and Abate, A. R. (2017). Single-cell Genome Sequencing at Ultra-high-throughput with Microfluidic Droplet Barcoding. Nat. Biotechnol. 35 (7), 640–646. doi:10.1038/nbt.3880

PubMed Abstract | CrossRef Full Text | Google Scholar

Leader, A. M., Grout, J. A., Maier, B. B., Nabet, B. Y., Park, M. D., Tabachnikova, A., et al. (2021). Single-cell Analysis of Human Non-small Cell Lung Cancer Lesions Refines Tumor Classification and Patient Stratification. Cancer Cell 39 (12), 1594–1609. e1512. doi:10.1016/j.ccell.2021.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Lengauer, C., Kinzler, K. W., and Vogelstein, B. (1998). Genetic Instabilities in Human Cancers. Nature 396 (6712), 643–649. doi:10.1038/25292

PubMed Abstract | CrossRef Full Text | Google Scholar

Leung, M. L., Davis, A., Gao, R., Casasent, A., Wang, Y., Sei, E., et al. (2017). Single-cell DNA Sequencing Reveals a Late-Dissemination Model in Metastatic Colorectal Cancer. Genome Res. 27 (8), 1287–1299. doi:10.1101/gr.209973.116

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Byrne, K. T., Yan, F., Yamazoe, T., Chen, Z., Baslan, T., et al. (2018a). Tumor Cell-Intrinsic Factors Underlie Heterogeneity of Immune Cell Infiltration and Response to Immunotherapy. Immunity 49 (1), 178–193. e177. doi:10.1016/j.immuni.2018.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Wu, F.-X., and Ngom, A. (2018b). A Review on Machine Learning Principles for Multi-View Biological Data Integration. Brief. Bioinform 19 (2), bbw113–340. doi:10.1093/bib/bbw113

CrossRef Full Text | Google Scholar

Liberzon, A., Subramanian, A., Pinchback, R., Thorvaldsdottir, H., Tamayo, P., and Mesirov, J. P. (2011). Molecular Signatures Database (MSigDB) 3.0. Bioinformatics 27 (12), 1739–1740. doi:10.1093/bioinformatics/btr260

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, B., Hu, X., Feng, K., Gao, R., Xue, Z., Zhang, S., et al. (2021). Temporal Single-Cell Tracing Reveals Clonal Revival and Expansion of Precursor Exhausted T Cells during Anti-PD-1 Therapy in Lung Cancer. Nat. Cancer 3, 108–121. doi:10.1038/s43018-021-00292-8

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (1), 470. doi:10.1038/s41467-018-08205-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Liudahl, S. M., Betts, C. B., Sivagnanam, S., Morales-Oyarvide, V., da Silva, A., Yuan, C., et al. (2021). Leukocyte Heterogeneity in Pancreatic Ductal Adenocarcinoma: Phenotypic and Spatial Features Associated with Clinical Outcome. Cancer Discov. 11 (8), 2014–2031. doi:10.1158/2159-8290.CD-20-0841

PubMed Abstract | CrossRef Full Text | Google Scholar

L. Lun, A. T., Bach, K., and Marioni, J. C. (2016). Pooling across Cells to Normalize Single-Cell RNA Sequencing Data with Many Zero Counts. Genome Biol. 17, 75. doi:10.1186/s13059-016-0947-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lotfollahi, M., Naghipourfar, M., Luecken, M. D., Khajavi, M., Büttner, M., Wagenstetter, M., et al. (2022). Mapping Single-Cell Data to Reference Atlases by Transfer Learning. Nat. Biotechnol. 40 (1), 121–130. doi:10.1038/s41587-021-01001-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Luecken, M. D., Büttner, M., Chaichoompu, K., Danese, A., Interlandi, M., Mueller, M. F., et al. (2022). Benchmarking Atlas-Level Data Integration in Single-Cell Genomics. Nat. Methods 19 (1), 41–50. doi:10.1038/s41592-021-01336-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Łuksza, M., Riaz, N., Makarov, V., Balachandran, V. P., Hellmann, M. D., Solovyov, A., et al. (2017). A Neoantigen Fitness Model Predicts Tumour Response to Checkpoint Blockade Immunotherapy. Nature 551 (7681), 517–520. doi:10.1038/nature24473

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, L., Hernandez, M. O., Zhao, Y., Mehta, M., Tran, B., Kelly, M., et al. (2019). Tumor Cell Biodiversity Drives Microenvironmental Reprogramming in Liver Cancer. Cancer Cell 36 (4), 418–430. e416. doi:10.1016/j.ccell.2019.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, S., and Zhang, Y. (2020). Profiling Chromatin Regulatory Landscape: Insights into the Development of ChIP-Seq and ATAC-Seq. Mol. Biomed. 1 (1), 9. doi:10.1186/s43556-020-00009-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Mabbott, N. A., Baillie, J., Brown, H., Freeman, T. C., and Hume, D. A. (2013). An Expression Atlas of Human Primary Cells: Inference of Gene Function from Coexpression Networks. BMC Genomics 14, 632. doi:10.1186/1471-2164-14-632

PubMed Abstract | CrossRef Full Text | Google Scholar

Macaulay, I. C., Haerty, W., Kumar, P., Li, Y. I., Hu, T. X., Teng, M. J., et al. (2015). G&T-seq: Parallel Sequencing of Single-Cell Genomes and Transcriptomes. Nat. Methods 12 (6), 519–522. doi:10.1038/nmeth.3370

PubMed Abstract | CrossRef Full Text | Google Scholar

Macosko, E. Z., Basu, A., Satija, R., Nemesh, J., Shekhar, K., Goldman, M., et al. (2015). Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell 161 (5), 1202–1214. doi:10.1016/j.cell.2015.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Maley, C. C., Koelble, K., Natrajan, R., Aktipis, A., and Yuan, Y. (2015). An Ecological Measure of Immune-Cancer Colocalization as a Prognostic Factor for Breast Cancer. Breast Cancer Res. 17 (1), 131. doi:10.1186/s13058-015-0638-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Manesse, M., Patel, K. K., Bobrow, M., and Downing, S. R. (2020). The InSituPlex Staining Method for Multiplexed Immunofluorescence Cell Phenotyping and Spatial Profiling of Tumor FFPE Samples. Methods Mol. Biol. 2055, 585–592. doi:10.1007/978-1-4939-9773-2_26

PubMed Abstract | CrossRef Full Text | Google Scholar

Marigo, I., Trovato, R., Hofer, F., Ingangi, V., Desantis, G., Leone, K., et al. (2020). Disabled Homolog 2 Controls Prometastatic Activity of Tumor-Associated Macrophages. Cancer Discov. 10 (11), 1758–1773. doi:10.1158/2159-8290.CD-20-0036

PubMed Abstract | CrossRef Full Text | Google Scholar

Marigo, I., Zilio, S., Desantis, G., Mlecnik, B., Agnellini, A. H. R., Ugel, S., et al. (2016). T Cell Cancer Therapy Requires CD40-Cd40l Activation of Tumor Necrosis Factor and Inducible Nitric-Oxide-Synthase-Producing Dendritic Cells. Cancer Cell 30 (4), 651. doi:10.1016/j.ccell.2016.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Martens, J. H. A., and Stunnenberg, H. G. (2013). BLUEPRINT: Mapping Human Blood Cell Epigenomes. Haematologica 98 (10), 1487–1489. doi:10.3324/haematol.2013.094243

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsumoto, H., Kiryu, H., Furusawa, C., Ko, M. S. H., Ko, S. B. H., Gouda, N., et al. (2017). SCODE: an Efficient Regulatory Network Inference Algorithm from Single-Cell RNA-Seq during Differentiation. Bioinformatics 33 (15), 2314–2321. doi:10.1093/bioinformatics/btx194

PubMed Abstract | CrossRef Full Text | Google Scholar

McGinnis, C. S., Murrow, L. M., and Gartner, Z. J. (2019). DoubletFinder: Doublet Detection in Single-Cell RNA Sequencing Data Using Artificial Nearest Neighbors. Cell Syst. 8 (4), 329–337. doi:10.1016/j.cels.2019.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

McGranahan, N., and Swanton, C. (2015). Biological and Therapeutic Impact of Intratumor Heterogeneity in Cancer Evolution. Cancer Cell 27 (1), 15–26. doi:10.1016/j.ccell.2014.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

McGranahan, N., and Swanton, C. (2017). Clonal Heterogeneity and Tumor Evolution: Past, Present, and the Future. Cell 168 (4), 613–628. doi:10.1016/j.cell.2017.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (5), 409–412. doi:10.1038/s41592-019-0392-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Moerman, T., Aibar Santos, S., Bravo González-Blas, C., Simm, J., Moreau, Y., Aerts, J., et al. (2019). GRNBoost2 and Arboreto: Efficient and Scalable Inference of Gene Regulatory Networks. Bioinformatics 35 (12), 2159–2161. doi:10.1093/bioinformatics/bty916

PubMed Abstract | CrossRef Full Text | Google Scholar

Moon, K. R., van Dijk, D., Wang, Z., Gigante, S., Burkhardt, D. B., Chen, W. S., et al. (2019). Visualizing Structure and Transitions in High-Dimensional Biological Data. Nat. Biotechnol. 37 (12), 1482–1492. doi:10.1038/s41587-019-0336-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ng, A. Y., Jordan, M. I., and Weiss, Y. (2002). On Spectral Clustering: Analysis and an Algorithm, 849–856.

Google Scholar

Nirmal, A. J., Maliga, Z., Vallius, T., Quattrochi, B., Chen, A. A., Jacobson, C. A., et al. (2022). The Spatial Landscape of Progression and Immunoediting in Primary Melanoma at Single Cell Resolution. Cancer Discov. doi:10.1158/2159-8290.CD-21-1357

CrossRef Full Text | Google Scholar

Nolan, J. P., Condello, D., Duggan, E., Naivar, M., and Novo, D. (2013). Visible and Near Infrared Fluorescence Spectral Flow Cytometry. Cytometry 83A (3), 253–264. doi:10.1002/cyto.a.22241

PubMed Abstract | CrossRef Full Text | Google Scholar

Nosho, K., Baba, Y., Tanaka, N., Shima, K., Hayashi, M., Meyerhardt, J. A., et al. (2010). Tumour-infiltrating T-Cell Subsets, Molecular Changes in Colorectal Cancer, and Prognosis: Cohort Study and Literature Review. J. Pathol. 222 (4), 350–366. doi:10.1002/path.2774

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogino, S., Nosho, K., Irahara, N., Meyerhardt, J. A., Baba, Y., Shima, K., et al. (2009). Lymphocytic Reaction to Colorectal Cancer Is Associated with Longer Survival, Independent of Lymph Node Count, Microsatellite Instability, and CpG Island Methylator Phenotype. Clin. Cancer Res. 15 (20), 6412–6420. doi:10.1158/1078-0432.CCR-09-1438

PubMed Abstract | CrossRef Full Text | Google Scholar

Osorio, D., and Cai, J. J. (2021). Systematic Determination of the Mitochondrial Proportion in Human and Mice Tissues for Single-Cell RNA-Sequencing Data Quality Control. Bioinformatics 37 (7), 963–967. doi:10.1093/bioinformatics/btaa751

PubMed Abstract | CrossRef Full Text | Google Scholar

Papalexi, E., Mimitou, E. P., Butler, A. W., Foster, S., Bracken, B., Mauck, W. M., et al. (2021). Characterizing the Molecular Regulation of Inhibitory Immune Checkpoints with Multimodal Single-Cell Screens. Nat. Genet. 53 (3), 322–331. doi:10.1038/s41588-021-00778-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, J., Sun, B.-F., Chen, C.-Y., Zhou, J.-Y., Chen, Y.-S., Chen, H., et al. (2019). Single-cell RNA-Seq Highlights Intra-tumoral Heterogeneity and Malignant Progression in Pancreatic Ductal Adenocarcinoma. Cell Res. 29 (9), 725–738. doi:10.1038/s41422-019-0195-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Peregrin, J., Svĕrák, J., and Doskocil, Z. (1973). Photostimulator for the Study of Visual Functions. Sb. Ved. P. R. Lek. Fak. Karlovy Univerzity Hradci Kralove Suppl. 16 (3), 203–207.

Google Scholar

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. Nat. Biotechnol. 35 (10), 936–939. doi:10.1038/nbt.3973

PubMed Abstract | CrossRef Full Text | Google Scholar

Picelli, S., Faridani, O. R., Björklund, Å. K., Winberg, G., Sagasser, S., and Sandberg, R. (2014). Full-length RNA-Seq from Single Cells Using Smart-Seq2. Nat. Protoc. 9 (1), 171–181. doi:10.1038/nprot.2014.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Pombo Antunes, A. R., Scheyltjens, I., Lodi, F., Messiaen, J., Antoranz, A., Duerinck, J., et al. (2021). Single-cell Profiling of Myeloid Cells in Glioblastoma across Species and Disease Stage Reveals Macrophage Competition and Specialization. Nat. Neurosci. 24 (4), 595–610. doi:10.1038/s41593-020-00789-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Pott, S. (2017). Simultaneous Measurement of Chromatin Accessibility, DNA Methylation, and Nucleosome Phasing in Single Cells. Elife 6. doi:10.7554/eLife.23203

PubMed Abstract | CrossRef Full Text | Google Scholar

Pratapa, A., Jalihal, A. P., Law, J. N., Bharadwaj, A., and Murali, T. M. (2020). Benchmarking Algorithms for Gene Regulatory Network Inference from Single-Cell Transcriptomic Data. Nat. Methods 17 (2), 147–154. doi:10.1038/s41592-019-0690-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Pylayeva-Gupta, Y., Lee, K. E., Hajdu, C. H., Miller, G., and Bar-Sagi, D. (2012). Oncogenic Kras-Induced GM-CSF Production Promotes the Development of Pancreatic Neoplasia. Cancer Cell 21 (6), 836–847. doi:10.1016/j.ccr.2012.04.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, J., Olbrecht, S., Boeckx, B., Vos, H., Laoui, D., Etlioglu, E., et al. (2020). A Pan-Cancer Blueprint of the Heterogeneous Tumor Microenvironment Revealed by Single-Cell Profiling. Cell Res. 30 (9), 745–762. doi:10.1038/s41422-020-0355-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiu, X., Mao, Q., Tang, Y., Wang, L., Chawla, R., Pliner, H. A., et al. (2017). Reversed Graph Embedding Resolves Complex Single-Cell Trajectories. Nat. Methods 14 (10), 979–982. doi:10.1038/nmeth.4402

PubMed Abstract | CrossRef Full Text | Google Scholar

Raghavan, S., Winter, P. S., Navia, A. W., Williams, H. L., DenAdel, A., Lowder, K. E., et al. (2021). Microenvironment Drives Cell State, Plasticity, and Drug Response in Pancreatic Cancer. Cell 184 (25), 6119–6137. e6126. doi:10.1016/j.cell.2021.11.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Rambow, F., Rogiers, A., Marin-Bejar, O., Aibar, S., Femel, J., Dewaele, M., et al. (2018). Toward Minimal Residual Disease-Directed Therapy in Melanoma. Cell 174 (4), 843–855. doi:10.1016/j.cell.2018.06.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramsköld, D., Luo, S., Wang, Y.-C., Li, R., Deng, Q., Faridani, O. R., et al. (2012). Full-length mRNA-Seq from Single-Cell Levels of RNA and Individual Circulating Tumor Cells. Nat. Biotechnol. 30 (8), 777–782. doi:10.1038/nbt.2282

PubMed Abstract | CrossRef Full Text | Google Scholar

Reddy, E. P., Reynolds, R. K., Santos, E., and Barbacid, M. (1982). A Point Mutation Is Responsible for the Acquisition of Transforming Properties by the T24 Human Bladder Carcinoma Oncogene. Nature 300 (5888), 149–152. doi:10.1038/300149a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Remark, R., Merghoub, T., Grabe, N., Litjens, G., Damotte, D., Wolchok, J. D., et al. (2016). In-depth Tissue Profiling Using Multiplexed Immunohistochemical Consecutive Staining on Single Slide. Sci. Immunol. 1 (1), aaf6925. doi:10.1126/sciimmunol.aaf6925

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, J. P. (2004). Multispectral Cytometry: the Next Generation. Biophot. Int. 11, 36–41.

Google Scholar

Rozenblatt-Rosen, O., Regev, A., Oberdoerffer, P., Nawy, T., Hupalowska, A., Rood, J. E., et al. (2020). The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution. Cell 181 (2), 236–249. doi:10.1016/j.cell.2020.03.053

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (1-2), 361–376. doi:10.1016/j.cell.2018.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Rybakowska, P., Alarcón-Riquelme, M. E., and Marañón, C. (2020). Key Steps and Methods in the Experimental Design and Data Analysis of Highly Multi-Parametric Flow and Mass Cytometry. Comput. Struct. Biotechnol. J. 18, 874–886. doi:10.1016/j.csbj.2020.03.024

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (4), 998–1013. doi:10.1016/j.cell.2018.10.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Saelens, W., Cannoodt, R., Todorov, H., and Saeys, Y. (2019). A Comparison of Single-Cell Trajectory Inference Methods. Nat. Biotechnol. 37 (5), 547–554. doi:10.1038/s41587-019-0071-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakamoto, K. M., Kim, K. B., Kumagai, A., Mercurio, F., Crews, C. M., and Deshaies, R. J. (2001). Protacs: Chimeric Molecules that Target Proteins to the Skp1-Cullin-F Box Complex for Ubiquitination and Degradation. Proc. Natl. Acad. Sci. U.S.A. 98 (15), 8554–8559. doi:10.1073/pnas.141230798

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanmamed, M. F., Nie, X., Desai, S. S., Villaroel-Espindola, F., Badri, T., Zhao, D., et al. (2021). A Burned-Out CD8+ T-Cell Subset Expands in the Tumor Microenvironment and Curbs Cancer Immunotherapy. Cancer Discov. 11 (7), 1700–1715. doi:10.1158/2159-8290.CD-20-0962

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, E., Martin-Zanca, D., Reddy, E. P., Pierotti, M. A., Della Porta, G., and Barbacid, M. (1984). Malignant Activation of a K- Ras Oncogene in Lung Carcinoma but Not in Normal Tissue of the Same Patient. Science 223 (4637), 661–664. doi:10.1126/science.6695174

PubMed Abstract | CrossRef Full Text | Google Scholar

Satpathy, A. T., Granja, J. M., Yost, K. E., Qi, Y., Meschi, F., McDermott, G. P., et al. (2019). Massively Parallel Single-Cell Chromatin Landscapes of Human Immune Cell Development and Intratumoral T Cell Exhaustion. Nat. Biotechnol. 37 (8), 925–936. doi:10.1038/s41587-019-0206-z

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (6), 307–325. doi:10.1038/s41568-019-0144-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Schreiber, R. D., Old, L. J., and Smyth, M. J. (2011). Cancer Immunoediting: Integrating Immunity's Roles in Cancer Suppression and Promotion. Science 331 (6024), 1565–1570. doi:10.1126/science.1203486

PubMed Abstract | CrossRef Full Text | Google Scholar

Schürch, C. M., Bhate, S. S., Barlow, G. L., Phillips, D. J., Noti, L., Zlobec, I., et al. (2020). Coordinated Cellular Neighborhoods Orchestrate Antitumoral Immunity at the Colorectal Cancer Invasive Front. Cell 182 (5), 1341–1359. doi:10.1016/j.cell.2020.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheng, K., Cao, W., Niu, Y., Deng, Q., and Zong, C. (2017). Effective Detection of Variation in Single-Cell Transcriptomes Using MATQ-Seq. Nat. Methods 14 (3), 267–270. doi:10.1038/nmeth.4145

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, R., Hie, B. L., Narayan, A., and Berger, B. (2021). Schema: Metric Learning Enables Interpretable Synthesis of Heterogeneous Single-Cell Modalities. Genome Biol. 22 (1), 131. doi:10.1186/s13059-021-02313-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Specht, A. T., and Li, J. (2017). LEAP: Constructing Gene Co-expression Networks for Single-Cell RNA-Sequencing Data Using Pseudotime Ordering. Bioinformatics 33 (5), btw729–766. doi:10.1093/bioinformatics/btw729

PubMed Abstract | CrossRef Full Text | Google Scholar

Specht, H., Emmott, E., Petelski, A. A., Huffman, R. G., Perlman, D. H., Serra, M., et al. (2021). Single-cell Proteomic and Transcriptomic Analysis of Macrophage Heterogeneity Using SCoPE2. Genome Biol. 22 (1), 50. doi:10.1186/s13059-021-02267-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Squair, J. W., Gautier, M., Kathe, C., Anderson, M. A., James, N. D., Hutson, T. H., et al. (2021). Confronting False Discoveries in Single-Cell Differential Expression. Nat. Commun. 12 (1), 5692. doi:10.1038/s41467-021-25960-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Stassen, S. V., Yip, G. G. K., Wong, K. K. Y., Ho, J. W. K., and Tsia, K. K. (2021). Generalized and Scalable Trajectory Inference in Single-Cell Omics Data with VIA. Nat. Commun. 12 (1), 5528. doi:10.1038/s41467-021-25773-3

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (9), 865–868. doi:10.1038/nmeth.4380

PubMed Abstract | CrossRef Full Text | Google Scholar

Stoeckius, M., Zheng, S., Houck-Loomis, B., Hao, S., Yeung, B. Z., Mauck, W. M., et al. (2018). Cell Hashing with Barcoded Antibodies Enables Multiplexing and Doublet Detection for Single Cell Genomics. Genome Biol. 19 (1), 224. doi:10.1186/s13059-018-1603-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Stuart, T., Butler, A., Hoffman, P., Hafemeister, C., Papalexi, E., Mauck, W. M., et al. (2019). Comprehensive Integration of Single-Cell Data. Cell 177 (7), 1888–1902. doi:10.1016/j.cell.2019.05.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Subrahmanyam, P. B., Dong, Z., Gusenleitner, D., Giobbie-Hurder, A., Severgnini, M., Zhou, J., et al. (2018). Distinct Predictive Biomarker Candidates for Response to Anti-CTLA-4 and Anti-PD-1 Immunotherapy in Melanoma Patients. J. Immunother. cancer 6 (1), 18. doi:10.1186/s40425-018-0328-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Subramanian, A., Alperovich, M., Yang, Y., and Li, B. (2021). Biology-inspired Data-Driven Quality Control for Scientific Discovery in Single-Cell Transcriptomics. bioRxiv. doi:10.1101/2021.10.27.466176

CrossRef Full Text | Google Scholar

Taavitsainen, S., Engedal, N., Cao, S., Handle, F., Erickson, A., Prekovic, S., et al. (2021). Single-cell ATAC and RNA Sequencing Reveal Pre-existing and Persistent Cells Associated with Prostate Cancer Relapse. Nat. Commun. 12 (1), 5307. doi:10.1038/s41467-021-25624-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, H., Wang, Y., Chlewicki, L. K., Zhang, Y., Guo, J., Liang, W., et al. (2016). Facilitating T Cell Infiltration in Tumor Microenvironment Overcomes Resistance to PD-L1 Blockade. Cancer Cell 29 (3), 285–296. doi:10.1016/j.ccell.2016.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Tate, J. G., Bamford, S., Jubb, H. C., Sondka, Z., Beare, D. M., Bindal, N., et al. (2019). COSMIC: the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res. 47 (D1), D941–D947. doi:10.1093/nar/gky1015

PubMed Abstract | CrossRef Full Text | Google Scholar

Tawbi, H. A., Schadendorf, D., Lipson, E. J., Ascierto, P. A., Matamala, L., Castillo Gutiérrez, E., et al. (2022). Relatlimab and Nivolumab versus Nivolumab in Untreated Advanced Melanoma. N. Engl. J. Med. 386 (1), 24–34. doi:10.1056/NEJMoa2109970

PubMed Abstract | CrossRef Full Text | Google Scholar

Tedesco, M., Giannese, F., Lazarević, D., Giansanti, V., Rosano, D., Monzani, S., et al. (2021). Chromatin Velocity Reveals Epigenetic Dynamics by Single-Cell Profiling of Heterochromatin and Euchromatin. Nat. Biotechnol. 40, 235–244. doi:10.1038/s41587-021-01031-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Telenius, H. k., Carter, N. P., Bebb, C. E., Nordenskjo¨ld, M., Ponder, B. A. J., and Tunnacliffe, A. (1992). Degenerate Oligonucleotide-Primed PCR: General Amplification of Target DNA by a Single Degenerate Primer. Genomics 13 (3), 718–725. doi:10.1016/0888-7543(92)90147-k

PubMed Abstract | CrossRef Full Text | Google Scholar

Thorsson, V., Gibbs, D. L., Brown, S. D., Wolf, D., Bortone, D. S., Ou Yang, T. H., et al. (2019). The Immune Landscape of Cancer. Immunity 51 (2), 411–412. doi:10.1016/j.immuni.2019.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (6282), 189–196. doi:10.1126/science.aad0501

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomlins, S. A., Rhodes, D. R., Perner, S., Dhanasekaran, S. M., Mehra, R., Sun, X.-W., et al. (2005). Recurrent Fusion of TMPRSS2 and ETS Transcription Factor Genes in Prostate Cancer. Science 310 (5748), 644–648. doi:10.1126/science.1117679

PubMed Abstract | CrossRef Full Text | Google Scholar

Traag, V. A., Waltman, L., and van Eck, N. J. (2019). From Louvain to Leiden: Guaranteeing Well-Connected Communities. Sci. Rep. 9 (1), 5233. doi:10.1038/s41598-019-41695-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Tsujikawa, T., Kumar, S., Borkar, R. N., Azimi, V., Thibault, G., Chang, Y. H., et al. (2017). Quantitative Multiplex Immunohistochemistry Reveals Myeloid-Inflamed Tumor-Immune Complexity Associated with Poor Prognosis. Cell Rep. 19 (1), 203–217. doi:10.1016/j.celrep.2017.03.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Ugel, S., Canè, S., De Sanctis, F., and Bronte, V. (2021). Monocytes in the Tumor Microenvironment. Annu. Rev. Pathol. Mech. Dis. 16, 93–122. doi:10.1146/annurev-pathmechdis-012418-013058

CrossRef Full Text | Google Scholar

Van Gassen, S., Callebaut, B., Van Helden, M. J., Lambrecht, B. N., Demeester, P., Dhaene, T., et al. (2015). FlowSOM: Using Self-Organizing Maps for Visualization and Interpretation of Cytometry Data. Cytometry 87 (7), 636–645. doi:10.1002/cyto.a.22625

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Gassen, S., Gaudilliere, B., Angst, M. S., Saeys, Y., and Aghaeepour, N. (2020). CytoNorm: A Normalization Algorithm for Cytometry Data. Cytometry 97 (3), 268–278. doi:10.1002/cyto.a.23904

PubMed Abstract | CrossRef Full Text | Google Scholar

Velazquez-Villarreal, E. I., Maheshwari, S., Sorenson, J., Fiddes, I. T., Kumar, V., Yin, Y., et al. (2020). Single-cell Sequencing of Genomic DNA Resolves Sub-clonal Heterogeneity in a Melanoma Cell Line. Commun. Biol. 3 (1), 318. doi:10.1038/s42003-020-1044-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Viratham Pulsawatdi, A., Craig, S. G., Bingham, V., McCombe, K., Humphries, M. P., Senevirathne, S., et al. (2020). A Robust Multiplex Immunofluorescence and Digital Pathology Workflow for the Characterisation of the Tumour Immune Microenvironment. Mol. Oncol. 14 (10), 2384–2402. doi:10.1002/1878-0261.12764

PubMed Abstract | CrossRef Full Text | Google Scholar

Vitak, S. A., Torkenczy, K. A., Rosenkrantz, J. L., Fields, A. J., Christiansen, L., Wong, M. H., et al. (2017). Sequencing Thousands of Single-Cell Genomes with Combinatorial Indexing. Nat. Methods 14 (3), 302–308. doi:10.1038/nmeth.4154

PubMed Abstract | CrossRef Full Text | Google Scholar

Vitale, I., Manic, G., Coussens, L. M., Kroemer, G., and Galluzzi, L. (2019). Macrophages and Metabolism in the Tumor Microenvironment. Cell Metab. 30 (1), 36–50. doi:10.1016/j.cmet.2019.06.001

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (5), 1330–1345. e1318. doi:10.1016/j.cell.2019.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Waltman, L., and Van Eck, N. J. (2013). A Smart Local Moving Algorithm for Large-Scale Modularity-Based Community Detection. Eur. Phys. J. B 86 (11), 11434–146036. doi:10.1140/epjb/e2013-40829-0

CrossRef Full Text | Google Scholar

Wang, Q., Xiong, H., Ai, S., Yu, X., Liu, Y., Zhang, J., et al. (2019). CoBATCH for High-Throughput Single-Cell Epigenomic Profiling. Mol. Cell 76 (1), 206–216. doi:10.1016/j.molcel.2019.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Sun, Z., Zhang, Y., Xu, Z., Xin, H., Huang, H., et al. (2020). BREM-SC: a Bayesian Random Effects Mixture Model for Joint Clustering Single Cell Multi-Omics Data. Nucleic Acids Res. 48 (11), 5814–5824. doi:10.1093/nar/gkaa314

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Waters, J., Leung, M. L., Unruh, A., Roh, W., Shi, X., et al. (2014). Clonal Evolution in Breast Cancer Revealed by Single Nucleus Genome Sequencing. Nature 512 (7513), 155–160. doi:10.1038/nature13600

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Yuan, P., Yan, Z., Yang, M., Huo, Y., Nie, Y., et al. (2021). Single-cell Multiomics Sequencing Reveals the Functional Regulatory Landscape of Early Embryos. Nat. Commun. 12 (1), 1247. doi:10.1038/s41467-021-21409-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, R., He, S., Bai, S., Sei, E., Hu, M., Thompson, A., et al. (2022). Spatial Charting of Single-Cell Transcriptomes in Tissues. Nat. Biotechnol. doi:10.1038/s41587-022-01233-1

CrossRef Full Text | Google Scholar

Weinreb, C., Wolock, S., and Klein, A. M. (2018). SPRING: a Kinetic Interface for Visualizing High Dimensional Single-Cell Expression Data. Bioinformatics 34 (7), 1246–1248. doi:10.1093/bioinformatics/btx792

PubMed Abstract | CrossRef Full Text | Google Scholar

Whiteside, T. L. (2008). The Tumor Microenvironment and its Role in Promoting Tumor Growth. Oncogene 27 (45), 5904–5912. doi:10.1038/onc.2008.271

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolf, F. A., Angerer, P., and Theis, F. J. (2018). SCANPY: Large-Scale Single-Cell Gene Expression Data Analysis. Genome Biol. 19 (1), 15. doi:10.1186/s13059-017-1382-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolock, S. L., Lopez, R., and Klein, A. M. (2019). Scrublet: Computational Identification of Cell Doublets in Single-Cell Transcriptomic Data. Cell Syst. 8 (4), 281–291. doi:10.1016/j.cels.2018.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, S. Z., Al-Eryani, G., Roden, D. L., Junankar, S., Harvey, K., Andersson, A., et al. (2021). A Single-Cell and Spatially Resolved Atlas of Human Breast Cancers. Nat. Genet. 53 (9), 1334–1347. doi:10.1038/s41588-021-00911-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, C., and Su, Z. (2015). Identification of Cell Types from Single-Cell Transcriptomes Using a Novel Clustering Method. Bioinformatics 31 (12), 1974–1980. doi:10.1093/bioinformatics/btv088

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, W., Feng, Y., Zhou, J., Cheung, O. K.-W., Cao, J., Wang, J., et al. (2021). A Selective HDAC8 Inhibitor Potentiates Antitumor Immunity and Efficacy of Immune Checkpoint Blockade in Hepatocellular Carcinoma. Sci. Transl. Med. 13 (588). doi:10.1126/scitranslmed.aaz6804

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z.-Z., Kim, H. J., Wu, H., Jalali, S., Tang, X., Krull, J. E., et al. (2020). TIGIT Expression Is Associated with T-Cell Suppression and Exhaustion and Predicts Clinical Outcome and Anti-PD-1 Response in Follicular Lymphoma. Clin. Cancer Res. 26 (19), 5217–5231. doi:10.1158/1078-0432.CCR-20-0558

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, Y., Jiang, Y., Lam, K.-W. G., Berletch, J. B., Disteche, C. M., Noble, W. S., et al. (2019). High-Throughput Single-Cell Sequencing with Linear Amplification. Mol. Cell 76 (4), 676–690. e610. doi:10.1016/j.molcel.2019.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Yip, S. H., Wang, P., Kocher, J.-P. A., Sham, P. C., and Wang, J. (2017). Linnorm: Improved Statistical Analysis for Single Cell RNA-Seq Expression Data. Nucleic Acids Res. 45 (22), e179. doi:10.1093/nar/gkx828

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, Y. (2016). Spatial Heterogeneity in the Tumor Microenvironment. Cold Spring Harb. Perspect. Med. 6 (8), a026583. doi:10.1101/cshperspect.a026583

PubMed Abstract | CrossRef Full Text | Google Scholar

Zack, T. I., Schumacher, S. E., Carter, S. L., Cherniack, A. D., Saksena, G., Tabak, B., et al. (2013). Pan-cancer Patterns of Somatic Copy Number Alteration. Nat. Genet. 45 (10), 1134–1140. doi:10.1038/ng.2760

PubMed Abstract | CrossRef Full Text | Google Scholar

Zahn, H., Steif, A., Laks, E., Eirew, P., VanInsberghe, M., Shah, S. P., et al. (2017). Scalable Whole-Genome Single-Cell Library Preparation without Preamplification. Nat. Methods 14 (2), 167–173. doi:10.1038/nmeth.4140

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Chen, D., Song, D., Liu, X., Zhang, Y., Xu, X., et al. (2022). Clinical and Translational Values of Spatial Transcriptomics. Sig Transduct. Target Ther. 7 (1), 111. doi:10.1038/s41392-022-00960-w

CrossRef Full Text | Google Scholar

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 (2), 442–459. doi:10.1016/j.cell.2020.03.048

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (4), 829–845. doi:10.1016/j.cell.2019.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Chen, H., Mo, H., Hu, X., Gao, R., Zhao, Y., et al. (2021). Single-cell Analyses Reveal Key Immune Cell Subsets Associated with Response to PD-L1 Blockade in Triple-Negative Breast Cancer. Cancer Cell 39 (12), 1578–1593. e1578. doi:10.1016/j.ccell.2021.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Zhang, S., Liu, Y., He, X., Qu, M., Xu, G., et al. (2020a). Single-cell RNA Sequencing Reveals the Heterogeneity of Liver-Resident Immune Cells in Human. Cell Discov. 6, 22. doi:10.1038/s41421-020-0157-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, X., Wu, S., Fang, N., Sun, X., and Fan, J. (2020b). Evaluation of Single-Cell Classifiers for Single-Cell RNA Sequencing Data Sets. Brief. Bioinform 21 (5), 1581–1595. doi:10.1093/bib/bbz096

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, G. X. Y., Terry, J. M., Belgrader, P., Ryvkin, P., Bent, Z. W., Wilson, R., et al. (2017). Massively Parallel Digital Transcriptional Profiling of Single Cells. Nat. Commun. 8, 14049. doi:10.1038/ncomms14049

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, C., Preissl, S., and Ren, B. (2020). Single-cell Multimodal Omics: the Power of Many. Nat. Methods 17 (1), 11–14. doi:10.1038/s41592-019-0691-5

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (11), 1063–1070. doi:10.1038/s41594-019-0323-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J., Powis de Tenbossche, C. G., Cané, S., Colau, D., van Baren, N., Lurquin, C., et al. (2017). Resistance to Cancer Immunotherapy Mediated by Apoptosis of Tumor-Infiltrating Lymphocytes. Nat. Commun. 8 (1), 1404. doi:10.1038/s41467-017-00784-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, Z., Wang, W., Lin, F., Jordan, T., Li, G., Silverman, S., et al. (2021). Genome Profiles of Pathologist-Defined Cell Clusters by Multiregional LCM and G&T-seq in One Triple-Negative Breast Cancer Patient. Cell Rep. Med. 2 (10), 100404. doi:10.1016/j.xcrm.2021.100404

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (5), 1317–1334. doi:10.1016/j.immuni.2019.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zong, C., Lu, S., Chapman, A. R., and Xie, X. S. (2012). Genome-wide Detection of Single-Nucleotide and Copy-Number Variations of a Single Human Cell. Science 338 (6114), 1622–1626. doi:10.1126/science.1229164

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, C., and Chen, L. (2021). Deep-joint-learning Analysis Model of Single Cell Transcriptome and Open Chromatin Accessibility Data. Brief. Bioinform 22 (4). doi:10.1093/bib/bbaa287

CrossRef Full Text | Google Scholar