Abstract
Immunotherapy with checkpoint blockers (ICBs), aimed at unleashing the immune response toward tumor cells, has shown a great improvement in overall patient survival compared to standard therapy, but only in a subset of patients. While a number of recent studies have significantly improved our understanding of mechanisms playing an important role in the tumor microenvironment (TME), we still have an incomplete view of how the TME works as a whole. This hampers our ability to effectively predict the large heterogeneity of patients' response to ICBs. Systems approaches could overcome this limitation by adopting a holistic perspective to analyze the complexity of tumors. In this Mini Review, we focus on how an integrative view of the increasingly available multi-omics experimental data and computational approaches enables the definition of new systems-based predictive biomarkers. In particular, we will focus on three facets of the TME toward the definition of new systems biomarkers. First, we will review how different types of immune cells influence the efficacy of ICBs, not only in terms of their quantification, but also considering their localization and functional state. Second, we will focus on how different cells in the TME interact, analyzing how inter- and intra-cellular networks play an important role in shaping the immune response and are responsible for resistance to immunotherapy. Finally, we will describe the potential of looking at these networks as dynamic systems and how mathematical models can be used to study the rewiring of the complex interactions taking place in the TME.
A Change in the Landscape of Biomarkers Discovery
Tumor cells are able to activate several mechanisms to evade the immune response by disguising themselves as “self” cells. Binding to inhibitory checkpoint molecules (i.e., immune checkpoints) they can block antitumor activities of the immune system. Immunotherapy with immune checkpoint blockers (ICBs) uses antibodies to target immune checkpoints, such as PD1, PD-L1, and CTLA-4, unleashing the immune response. In clinical trials, ICB therapy has been shown to achieve durable therapeutic response and to increase patient survival in different cancer types, although still a small number of ICBs are FDA-approved (1, 2). Even if clinically approved, ICB therapy is effective for a small subset of patients. Given the potential immunological toxicity (3, 4) and the elevated costs (>US$100,000 per patient per year) (5) associated with ICBs, it is of paramount importance to be able to predict which patients will likely respond to the therapy, in order to administer the optimal treatment based on biomarkers.
The investigation of mechanisms supporting immune resistance has provided a great opportunity for biomarker discovery of patient response to ICBs (Figure 1). Two biomarkers have been clinically approved for PD-1/PD-L1 blockade therapy: the first is immunohistochemistry (IHC) staining of PD-L1 in non-small-cell lung cancer (NSCLC), melanoma, renal cell carcinoma (RCC), urothelial cancer, and triple-negative breast cancer (TNBC) (6); and the second is high microsatellite instability/defective mismatch repair (MSI-H/dMMR) regardless of tumor type (7, 8). Other emerging predictive biomarkers such as tumor mutational burden (TMB) (9, 10), signatures of a T cell inflamed tumor microenvironment (TME) either alone (10) or in combination (11), and neoantigen load (12–14) are still undergoing clinical trials. In addition, T cell receptor (TCR) diversity has been used as a biomarker to monitor the clonal expansion of T cells in breast cancer, glioma, cervical cancer, and leukemia/lymphoma (15–18). Further efforts both to exploit the utility of these biomarkers and to search for additional ones are still ongoing. For a complete review of these biomarkers and in which tumors they work, we refer to Havel et al. (19).
Figure 1
Despite being promising, these biomarkers also present some limitations. For instance, IHC enables measuring PD-L1 expressed on tumor cells, however the expression of this biomarker fluctuates over time and varies between different tumor sites. This variability undermines the ability to evaluate PD-1/PD-L1 therapies effectiveness based on IHC, as reviewed in Topalian et al. (20) and Camidge et al. (21). Another example is TMB, which is known to correlate imperfectly with clinical response (12, 13, 22). Neoantigen burden should partially overcome this issue, however most computational tools fail to estimate true neoantigens (19, 20, 23), and additional features should be considered to better determine neoantigen immunogenicity as reviewed in Finotello et al. (24).
Above-mentioned examples shed light upon the conceptual problem of looking only at individual components of the TME. While the characterization of different parts playing a role in the interaction between tumor and immune system has been essential to elucidate the most important actionable mechanisms, further research is required to define biomarkers harnessing a more coordinated joint action of these mechanisms. Predictive biomarkers for immunotherapy with ICBs have been extensively reviewed previously (19, 20, 23, 25). In this Mini Review we focus on how a holistic profiling of the TME can provide new opportunities for identifying systems-based biomarkers built on existing synergies between the different individual components of the TME. Such a shift toward multifaceted strategies has been favored by increasingly available multi-omics data from bulk populations, individual cells, and imaging technologies (26), that can be integrated using computational approaches. In the following sections we will describe how biomarkers can be derived by considering three increasing levels of complexity. The first is the cellular component, focusing on the immune contexture of tumors, such as immune cells quantification, functionality, and localization. The second is the network of communication between and within cells of the TME. Finally, we will elaborate on how mathematical models can be used to take the dynamic nature of these networks into account.
The Role of the Immune Contexture on ICB Efficacy
It is well-known that different types of immune cells can play a different role in the response to ICBs (27). For example, while the presence of CD8+ T cells within the TME is a good biomarker of ICBs efficacy, a high abundance of regulatory T (Treg) cells is generally associated with poor prognosis. Different tools have been developed to quantify tumor-infiltrating immune cells from bulk (RNA-seq) and single-cell (scRNA-seq) RNA sequencing measurements, as extensively reviewed in Finotello and Eduati (26) and Finotello and Trajanoski (28).
Apart from quantification of immune cells, their spatial localization also plays a pivotal role in the response to immunotherapy (29). For instance, CD8+ T cells not only need to be present, but also to be infiltrated (hot tumor) for the ICB therapy to work. In fact, pure quantification of CD8+ T cells is not always associated with favorable prognosis (30). Imaging techniques can be used to explore the spatial patterns of immune infiltration. A notable example of a biomarker assessing through IHC, both the abundance and the location (tumor center and invasive margin) of two lymphocyte populations (CD3+ and CD8+ T cells) is the immunoscore (31), that was shown to accurately predict patient survival in colorectal cancer patients (32). More recently, spatial information of CD8+ T cells from IHC was integrated with transcriptomics data to study the effect of lymphocyte infiltration in patients with TNBC, providing predictive biomarkers of ICBs response (33). Automatic approaches for image analysis could reveal useful in the future for high-throughput identification of spatial biomarkers. A first attempt in this direction was the development of tumor infiltrating lymphocytes maps by using deep learning on images from the cancer genome atlas (TCGA) (34).
Another important factor that affects patients' response to ICBs is the functional state of the different immune cells (35). Dysfunctional states of T cells can be characterized from bulk and single-cell RNA-seq (36–38) and epigenetic profiling (39–41). ICBs aim at rescuing dysfunctional T cells, therefore the investigation of their functional state can inform on ICBs therapy success and limitations (36–39, 41). Depending on the type of stimulatory signal, macrophages (42, 43), and B cells (44, 45) can develop into functional subsets that have either positive or negative effects on tumors. Another example are dendritic cells (DCs), that normally control cancer antigen presentation, priming and activation of T cell responses, however the TME can compromise their ability to stimulate the immune response (46, 47). Certain computational tools for cell-type quantification can also unmask the phenotypic state of cell subpopulations in the TME by inferring the transcriptomics profiles of individual cells (48, 49). A promising research direction for biomarkers discovery is also given by new technologies that allow generation of omics data from tissue slides preserving cell spatial identity (50, 51). These approaches would result in combined localization and characterization of the cells in the TME.
Analysis on the immune infiltrate quantification, functionality, and localization can help both to explain the diversity of the tumor immune milieu and develop informative biomarkers for ICBs (27, 52, 53). Pointing in this direction, different efforts have recently explored the use of bulk transcriptomics data to derive more complex immune-related scores to assess the likelihood of a patient to respond to ICBs (38, 54–63).
Intra- and Inter-Cellular Networks Orchestrate the Immune Response
The functional state of cells in the TME is defined by a complex system of communication between molecules within the cells (intra-cellular networks) and among different cells (inter-cellular networks). Looking at intra- and inter-cellular networks can provide a more holistic perspective of the TME and inform a new class of biomarkers for immunotherapy and its potential combination with other targeted therapies (64).
Intra-cellular signaling pathways play a part in shaping the interaction with the immune system [(65, 66); Figure 2]. Abnormalities in tumor-intrinsic signaling, involving oncogenes and tumor suppressor genes, have been associated with mechanisms of inherent immune resistance (67). Examples are PTEN loss (68) or EGFR gain of function (69), both causing PI3K-Akt pathway activation and leading to over-expression of PD-L1 and consequent immunoresistance. Due to the complexity of signaling pathways, with numerous cross-talks and feedback loops, the adoption of individual oncogenic drivers as biomarkers is not expected to be effective in most cases (20). In fact, PD-L1 signal is directly regulated by numerous oncogenic pathways such as Ras, mTOR, EGFR, MEK, ERK, and MAPK (70). Besides pathways regulating immune checkpoints, other signaling cross-talks control the immune response from different perspectives, like inactivation of TP53 or activation of β-catenin pathway, both reducing chemokine production by tumor cells and thereby reducing recruitment of immune cells into the TME (71, 72).
Figure 2
In addition, cancer cells receive signals from other cells in the TME through ligand-receptor interactions. These inter-cellular communications lead to changes in the phenotype of the regulated cells thus playing an important role in both progression and prognosis of cancer (73, 74). An example is the response elicited on cancer cells by two cytokines (TNF and IFNγ) produced by activated T cells. These cytokines induce PD-L1 expression through JAK-STAT and NF-kB signaling, inducting acquired resistance to the immune response (75, 76). Another study identified a relationship between high expression of NOS2 and prolonged IFN signaling in tumors resistant to PD-1 blockade (77).
While collections of intra- (78) and inter-cellular (79) interactions can be derived from literature and databases, additional data are required to characterize the networks for each patient or group of patients. Transcriptomics and proteomics data can provide the basis to study intra- and inter-cellular signaling networks. Imaging data can also be integrated to improve our understanding on spatial localization of interacting cells. Computational methods have been developed to infer integrated inter- and intra-cellular networks from bulk (80, 81) and single-cell (81, 82) RNA-seq data. These tools could be exploited to derive biomarkers for immunotherapy by studying the functional effect of cell-cell communication. In a recent study, a curated database of ligand-receptor interactions (79) was integrated with gene expression data to deconvolute the transcriptional profile of cancer and stromal cells and infer cross-talks in the TME (83). Interestingly, the authors show that for different cancer types, PD-L1 expression is higher on cancer or stromal cells which nicely correlates with the general responsiveness to immunotherapy. Further research is required to assess if this holds also for individual patients, making it potentially a more effective biomarker than bulk PD-L1 expression. In another recent publication (84), researchers performed an extensive literature curation to derive a comprehensive signaling network of innate immune response in cancer, including cell type-specific signaling in macrophages, DCs, myeloid-derived suppressor cells, and natural killer cells. Such network was then integrated with scRNAseq data from macrophages and natural killer cells in melanoma to study the heterogeneity of innate immune cell types and could potentially be used to predict patient survival and response to immunotherapies. Finally, Worzfeld et al. combined parallel bulk transcriptomics and proteomics data on tumor cell spheroids, tumor-associated T cells and macrophages to derive inter-cellular signaling networks in the ovarian cancer microenvironment (85). Such networks included several immune checkpoint regulators and appeared to have potential clinical relevance. Overall, these studies have demonstrated the enormous benefit that holistic approaches combining complex multicellular networks can bring into the immuno-oncology field, and we expect that in the forthcoming future more research efforts will be spent in this direction. The recent developments of 3D cell culture models resembling the TME, are expected to be a powerful tool for further in vitro and ex vivo investigation of intra-cellular communication, and to study their effect on the response to ICBs (86).
The Potential of Looking at the Dynamicity and Plasticity of the TME
It is well-known that the cellular functional state changes dynamically in response to environmental changes and perturbations such as drug treatment (87, 88), calling for identification of the dynamic properties of the networks. The ideal data for dynamic functional characterization of the system's response are obtained upon perturbation (89). Functional screening of the effect of cancer drugs has been so far focused on cancer cell lines. While cell lines are a debatable model system, they proved to be a valuable tool to explore novel biomarkers of drug response (90, 91). High-throughput drug screening studies are now also being increasingly performed on organoids (92) or other 3D experimental models (86), which are more physiological human cancer models of the TME. These efforts open new ways for pre-clinical investigation of the effect of immunotherapy. Finally, more recent technologies allow screening also of patient biopsies without need for culturing steps (93–95) paving the way for functional characterization of ex vivo tumor samples potentially improving personalized cancer treatment.
To capture the functional context of the immune response, statistical, and mathematical approaches are developing into more compendious methods that integrate multi-omics data and prior knowledge on network structure (Figure 2). While mathematical models do not fall into the standard definition of biomarkers, they can provide predictions of response to immunotherapy. Additionally they can be used to define dynamic biomarkers based on properties of the modeled system, as opposed to static biomarkers that only consider the initial conditions of the system (88).
Dynamic mathematical models can be used to study intra-cellular networks of the different cell types populating the TME (96). To characterize these networks at the patient-specific level, models of signaling pathways in cancer cells have been trained from perturbation experiments (97, 98), gene expression data (99), or integrating multi-omics data (100). The resulting parameters corresponding to these personalized models can be relevant biomarkers of clinical outcome (99–101). Mathematical models have also been used to study intra-cellular signaling in T cells. This includes the investigation of how PD-1 leads to deactivation of the T cell receptor signaling (102) or mechanistic understanding of T cell exhaustion (103). PD-1 is one of the main targets of ICB, and exhausted T cells have a higher number of targetable checkpoint proteins like PD-1 and CTLA-4, therefore the investigation of these aspects could be relevant to identify possible biomarkers.
More studies are now focusing on mathematical models incorporating inter-cellular interactions to better capture the complexity of the TME. Agent-based models can be used to simulate the interactions between cells in the tumor microenvironment seen as a 2D or a 3D grid (104). Each cell is seen as an agent that can perform different tasks with a certain probability (e.g., cells can non-proliferate, divide, or die). Since the immune response can be seen as a probabilistic outcome of a complex system (88), agent-based models are an adequate mathematical approximation to capture this stochasticity. These models can be refined using a multitude of data types and used to simulate the effect of immunotherapy (105, 106), providing a variety of possible outcomes given the same initial conditions that can be interpreted as probability of success. It has been shown that tumor-bearing inbred mice, which have only minimal differences, can respond differently to immunotherapy (88), therefore having models that can incorporate stochasticity provides an interesting approximation of the in vivo situation. Another approach to model cell-cell communication is by using response-time modeling (107), where cells are modeled as a black-box that can receive inputs (e.g., cytokines) from other cells, process them, and change state (e.g., immune cells can switch between inactive and active) accordingly with a certain probability. Recently, Grandclaudon et al. combined perturbation data with a multivariate quantitative model to study context dependent interactions between DCs and helper T cells (108). A different approach based on quantitative systems pharmacology was recently used to simulate the effect of ICB therapy in metastatic breast cancer patients using a four compartments (central, peripheral, tumor-draining lymph node, and tumor) model (109).
Additionally, combining mathematical models with longitudinal data, i.e., data collected at different time points, can be used to investigate the evolutionary dynamics of treatment response. This aspect is particularly relevant, especially to be able to distinguish at an early stage real tumor progression (patient should be assigned to a different treatment) from what is called pseudoprogression, i.e., temporary progression followed by a response to the treatment (patient should be kept on ICB). The latter behavior has been described using a model of immune activation incorporating the dynamics of antigen presentation (110). Based on a system of three ordinary differential equations to describe the interaction between tumor cells, Treg cells, and cytotoxic T cells, this model could explain why, in response to ICBs, the tumor can worsen before starting regressing. Other multi-cellular models have been used to derive in silico patients to test different possible dynamics of treatment response (111, 112), that could be compared with longitudinal measurements of tumor load from PET/CT imaging (112). Longitudinal data are often limited to non-invasive imaging and, in a few cases, to transcriptomics, IHC, TCR, and genome sequencing data (113, 114) for a limited number of time points due to invasiveness of biopsies. Computational modeling of longitudinal data is still at its infancy, but we envision that in the future more mechanistic dynamic models will be able to exploit this type of data for definition of dynamic biomarkers.
Conclusions and Future Perspectives
Current limitations in identifying predictive biomarkers for ICB therapy are partially due to overlooking the complexity of the TME. Following the advancements in technologies to measure multi-omics data, measurements of bulk populations, individual cells, and spatial information have paved the way to a more comprehensive view of the TME. Recent efforts are focused on searching for signatures of response to ICBs that consider quantification, localization, and functionality of different immune cells in the TME, showing improved predictive power with respect to simpler biomarkers (115). However, they still miss an integrative strategy that takes a view of the whole TME, rather than examining each factor in isolation. In this respect, mechanistic models incorporating existing biological basis, e.g., on intra- and inter-cellular pathways, can accompany both therapy and biomarker development in immuno-oncology (116).
There is compelling evidence that the interplay of the immune system, tumors, organs, and external environment, harmonizes antitumor immune responses (117). Therefore, we envision that novel systems medicine approaches entailing mathematical models can gradually build up a profile of the TME, both in the lab and, more importantly, in the clinic. To this end, building patient specific models have become of increasing importance, especially when based on data that can be measured in clinical settings. Moreover, systems approaches can especially be useful to provide rationale for alternative personalized treatments such as combinatorial therapy.
Statements
Author contributions
ÓL-S and FE wrote and edited the article. All authors contributed to the article and approved the submitted version.
Funding
Costs related to this publication are covered by the Computational Biology group of the Department of Biomedical Engineering of Eindhoven University of Technology.
Acknowledgments
The authors would like to thank the members of the CBio group for their helpful proofreading suggestions.
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.
- CTLA-4
cytotoxic T lymphocyte antigen 4
- DC
dendritic cell
- ICB
immune checkpoint blocker
- IFNγ
interferon gamma
- IHC
immunohistochemistry
- MMR
DNA mismatch repair
- MSI-H
high microsatellite instability
- NOS2
nitric oxide synthase 2
- NSCLC
non-small-cell lung cancer
- PD-L1
programmed cell death-ligand 1
- PD-1
programmed cell death protein 1
- RCC
renal cell carcinoma
- RNA-seq
RNA sequencing
- scRNA-seq
single-cell RNA sequencing
- TCGA
the cancer genome atlas
- TCR
T cell receptor
- TMB
tumor mutational burden
- TME
tumor microenvironment
- TNBC
triple-negative breast cancer
- TNF
tumor necrosis factor
- Treg
regulatory T cell.
Abbreviations
References
1.
HargadonKMJohnsonCEWilliamsCJ. Immune checkpoint blockade therapy for cancer: an overview of FDA-approved immune checkpoint inhibitors. Int Immunopharmacol. (2018) 62:29–39. 10.1016/j.intimp.2018.06.001
2.
SunCMezzadraRSchumacherTN. Regulation and function of the PD-L1 checkpoint. Immunity. (2018) 48:434–52. 10.1016/j.immuni.2018.03.014
3.
BoutrosCTarhiniARoutierELambotteOLadurieFLCarbonnelFet al. Safety profiles of anti-CTLA-4 and anti-PD-1 antibodies alone and in combination. Nat Rev Clin Oncol. (2016) 13:473–86. 10.1038/nrclinonc.2016.58
4.
PostowMASidlowRHellmannMD. Immune-related adverse events associated with immune checkpoint blockade. New Engl J Med. (2018) 378:158–68. 10.1056/NEJMra1703481
5.
SchmidtC. The benefits of immunotherapy combinations. Nature. (2017) 552:S67–9. 10.1038/d41586-017-08702-7
6.
AroraSVelichinskiiRLeshRWAliUKubiakMBansalPet al. Existing and emerging biomarkers for immune checkpoint immunotherapy in solid tumors. Adv Ther. (2019) 36:2638–78. 10.1007/s12325-019-01051-z
7.
LeDTUramJNWangHBartlettBKemberlingHEyringAet al. PD-1 blockade in tumors with mismatch repair deficiency. N Engl J Med. (2015) 372:2509–20. 10.1056/NEJMoa1500596
8.
LeDTDurhamJNSmithKNWangHBartlettBRAulakhLKet al. Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade. Science. (2017) 357:409–13. 10.1126/science.aan6733
9.
LegrandFAGandaraDRMariathasanSPowlesTHeXZhangWet al. Association of high tissue TMB and atezolizumab efficacy across multiple tumor types. J Clin Oncol. (2018) 36(Suppl. 15):12000. 10.1200/JCO.2018.36.15_suppl.12000
10.
OttPABangY-JPiha-PaulSAAbdul RazakARBennounaJSoriaJ-Cet al. T-cell–inflamed gene-expression profile, programmed death ligand 1 expression, and tumor mutational burden predict efficacy in patients treated with pembrolizumab across 20 cancers: KEYNOTE-028. J Clin Oncol. (2019) 37:318–27. 10.1200/JCO.2018.78.2276
11.
CristescuRMoggRAyersMAlbrightAMurphyEYearleyJet al. Pan-tumor genomic biomarkers for PD-1 checkpoint blockade-based immunotherapy. Science. (2018) 362:eaar3593. 10.1126/science.aar3593
12.
RizviNAHellmannMDSnyderAKvistborgPMakarovVHavelJJet al. Mutational landscape determines sensitivity to PD-1 blockade in non–small cell lung cancer. Science. (2015) 348:124–8. 10.1126/science.aaa1348
13.
Van AllenEMMiaoDSchillingBShuklaSABlankCZimmerLet al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science. (2015) 350:207–11. 10.1126/science.aad0095
14.
LuTWangSXuLZhouQSinglaNGaoJet al. Tumor neoantigenicity assessment with CSiN score incorporates clonality and immunogenicity to predict immunotherapy outcomes. Sci Immunol. (2020) 5:aaz3199. 10.1126/sciimmunol.aaz3199
15.
PageDBYuanJRedmondDWenYHDurackJCEmersonRet al. Deep sequencing of T-cell receptor DNA as a biomarker of clonally expanded TILs in breast cancer after immunotherapy. Cancer Immunol Res. (2016) 4:835–44. 10.1158/2326-6066.CIR-16-0013
16.
SimsJSGrinshpunBFengYUngTHNeiraJASamanamudJLet al. Diversity and divergence of the glioma-infiltrating T-cell receptor repertoire. Proc Natl Acad Sci USA. (2016) 113:E3529–37. 10.1073/pnas.1601012113
17.
CuiJHLinKRYuanSHJinYBChenXPSuXKet al. TCR Repertoire as a novel indicator for immune monitoring and prognosis assessment of patients with cervical cancer. Front Immunol. (2018) 9:2729. 10.3389/fimmu.2018.02729
18.
FarmanbarAKnellerRFirouziS. RNA sequencing identifies clonal structure of T-cell repertoires in patients with adult T-cell leukemia/lymphoma. NPJ Genom Med. (2019) 4:10. 10.1038/s41525-019-0084-9
19.
HavelJJChowellDChanTA. The evolving landscape of biomarkers for checkpoint inhibitor immunotherapy. Nat Rev Cancer. (2019) 19:133–50. 10.1038/s41568-019-0116-x
20.
TopalianSLTaubeJMAndersRAPardollDM. Mechanism-driven biomarkers to guide immune checkpoint blockade in cancer therapy. Nat Rev Cancer. (2016) 16:275–87. 10.1038/nrc.2016.36
21.
CamidgeDRDoebeleRCKerrKM. Comparing and contrasting predictive biomarkers for immunotherapy and targeted therapy of NSCLC. Nat Rev Clin Oncol. (2019) 16:341–55. 10.1038/s41571-019-0173-9
22.
SnyderAMakarovVMerghoubTYuanJZaretskyJMDesrichardAet al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. (2014) 371:2189–99. 10.1056/NEJMoa1406498
23.
KeenanTEBurkeKPVan AllenEM. Genomic correlates of response to immune checkpoint blockade. Nat Med. (2019) 25:389–402. 10.1038/s41591-019-0382-x
24.
FinotelloFRiederDHacklHTrajanoskiZ. Next-generation computational tools for interrogating cancer immunity. Nat Rev Genet. (2019) 20:724–46. 10.1038/s41576-019-0166-7
25.
DuffyMJCrownJ. Biomarkers for predicting response to immunotherapy with immune checkpoint inhibitors in cancer patients. Clin Chem. (2019) 65:1228–38. 10.1373/clinchem.2019.303644
26.
FinotelloFEduatiF. Multi-omics profiling of the tumor microenvironment: paving the way to precision immuno-oncology. Front Oncol. (2018) 8:430. 10.3389/fonc.2018.00430
27.
FridmanWHZitvogelLSautès-FridmanCKroemerG. The immune contexture in cancer prognosis and treatment. Nat Rev Clin Oncol. (2017) 14:717–34. 10.1038/nrclinonc.2017.101
28.
FinotelloFTrajanoskiZ. Quantifying tumor-infiltrating immune cells from transcriptomics data. Cancer Immunol Immunother. (2018) 67:1031–40. 10.1007/s00262-018-2150-z
29.
GalonJBruniD. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies. Nat Rev Drug Discov. (2019) 18:197–218. 10.1038/s41573-018-0007-y
30.
GiraldoNABechtEPagèsFSklirisGVerkarreVVanoYet al. Orchestration and prognostic significance of immune checkpoints in the microenvironment of primary and metastatic renal cell cancer. Clin Cancer Res. (2015) 21:3031–40. 10.1158/1078-0432.CCR-14-2926
31.
GalonJCostesASanchez-CaboFKirilovskyAMlecnikBLagorce-PagèsCet al. Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science. (2006) 313:1960–64. 10.1126/science.1129139
32.
PagèsFMlecnikBMarliotFBindeaGOuF-SBifulcoCet al. International validation of the consensus immunoscore for the classification of colon cancer: a prognostic and accuracy study. Lancet. (2018) 391:2128–39. 10.1016/S0140-6736(18)30789-X
33.
GruossoTGigouxMManemVSKBertosNZuoDPerlitchIet al. Spatially distinct tumor immune microenvironments stratify triple-negative breast cancers. J Clin Invest. (2019) 129:1785–1800. 10.1172/JCI96313
34.
SaltzJGuptaRHouLKurcTSinghPNguyenVet al. Spatial organization and molecular correlation of tumor-infiltrating lymphocytes using deep learning on pathology images. Cell Rep. (2018) 23:181–93.e7. 10.1016/j.celrep.2018.03.086
35.
ThorssonVGibbsDLBrownSDWolfDBortoneDSOu YangT-Het al. The immune landscape of cancer. Immunity. (2019) 51:411–12. 10.1016/j.immuni.2019.08.004
36.
Sade-FeldmanMYizhakKBjorgaardSLRayJPde BoerCGJenkinsRWet al. Defining T cell states associated with response to checkpoint immunotherapy in melanoma. Cell. (2018) 175:998–1013.e20. 10.1016/j.cell.2018.10.038
37.
ZhengCZhengLYooJ-KGuoHZhangYGuoXet al. Landscape of infiltrating T cells in liver cancer revealed by single-cell sequencing. Cell. (2017) 169:1342–56.e16. 10.1016/j.cell.2017.05.035
38.
JiangPGuSPanDFuJSahuAHuXet al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. Nat Med. (2018) 24:1550–58. 10.1038/s41591-018-0136-1
39.
PaukenKESammonsMAOdorizziPMManneSGodecJKhanOet al. Epigenetic stability of exhausted T cells limits durability of reinvigoration by PD-1 blockade. Science. (2016) 354:1160–65. 10.1126/science.aaf2807
40.
SenDRKaminskiJBarnitzRAKurachiMGerdemannUYatesKBet al. The epigenetic landscape of T cell exhaustion. Science. (2016) 354:1165–9. 10.1126/science.aae0491
41.
PhilipMFairchildLSunLHorsteELCamaraSShakibaMet al. Chromatin states define tumour-specific T cell dysfunction and reprogramming. Nature. (2017) 545:452–6. 10.1038/nature22367
42.
QuarantaVRainerCNielsenSRRaymantMLAhmedMSEngleDDet al. Macrophage-derived granulin drives resistance to immune checkpoint inhibition in metastatic pancreatic cancer. Cancer Res. (2018) 78:4253–69. 10.1158/0008-5472.CAN-17-3876
43.
KlugFPrakashHHuberPESeibelTBenderNHalamaNet al. Low-dose irradiation programs macrophage differentiation to an iNOS /M1 phenotype that orchestrates effective T cell immunotherapy. Cancer Cell. (2013) 24:589–602. 10.1016/j.ccr.2013.09.014
44.
LindnerSDahlkeKSontheimerKHagnMKaltenmeierCBarthTFEet al. Interleukin 21-induced granzyme B-expressing B cells infiltrate tumors and regulate T cells. Cancer Res. (2013) 73:2468–79. 10.1158/0008-5472.CAN-12-3450
45.
SelitskySRMoseLESmithCCChaiSHoadleyKADittmerDPet al. Prognostic value of B cells in cutaneous melanoma. Genome Med. (2019) 11:36. 10.1186/s13073-019-0647-5
46.
TesoneAJRutkowskiMRBrencicovaESvoronosNPerales-PuchaltAStephenTLet al. Satb1 overexpression drives tumor-promoting activities in cancer-associated dendritic cells. Cell Rep. (2016) 14:1774–86. 10.1016/j.celrep.2016.01.056
47.
FuertesMBKachaAKKlineJWooS-RKranzDMMurphyKMet al. Host type I IFN signals are required for antitumor CD8 T cell responses through CD8α dendritic cells. J Exp Med. (2011) 208:2005–16. 10.1084/jem.20101159
48.
NewmanAMSteenCBLiuCLGentlesAJChaudhuriAASchererFet al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. (2019) 37:773–82. 10.1038/s41587-019-0114-2
49.
ZaitsevKBambouskovaMSwainAArtyomovMN. Complete deconvolution of cellular mixtures based on linearity of transcriptional signatures. Nat Commun. (2019) 10:2209. 10.1038/s41467-019-09990-5
50.
GiesenCWangHAOSchapiroDZivanovicNJacobsAHattendorfBet al. Highly multiplexed imaging of tumor tissues with subcellular resolution by mass cytometry. Nat Methods. (2014) 11:417–22. 10.1038/nmeth.2869
51.
AngeloMBendallSCFinckRHaleMBHitzmanCBorowskyADet al. Multiplexed ion beam imaging of human breast tumors. Nat Med. (2014) 20:436–42. 10.1038/nm.3488
52.
GalluzziLChanTAKroemerGWolchokJDLópez-SotoA. The hallmarks of successful anticancer immunotherapy. Sci Transl Med. (2018) 10:eaat7807. 10.1126/scitranslmed.aat7807
53.
BinnewiesMRobertsEWKerstenKChanVFearonDFMeradMet al. Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat Med. (2018) 24:541–50. 10.1038/s41591-018-0014-x
54.
AuslanderNZhangGLeeJSFrederickDTMiaoBMollTet al. Publisher Correction: Robust prediction of response to immune checkpoint blockade therapy in metastatic melanoma. Nat Med. (2018) 24:1942. 10.1038/s41591-018-0247-8
55.
DavoliTUnoHWootenECElledgeSJ. Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy. Science. (2017) 355:eaaf8399. 10.1126/science.aaf8399
56.
FehrenbacherLSpiraABallingerMKowanetzMVansteenkisteJMazieresJet al. Atezolizumab versus docetaxel for patients with previously treated non-small-cell lung cancer (POPLAR): a multicentre, open-label, phase 2 randomised controlled trial. Lancet. (2016) 387:1837–46. 10.1016/S0140-6736(16)00587-0
57.
HuangACOrlowskiRJXuXMickRGeorgeSMYanPKet al. A single dose of neoadjuvant PD-1 blockade predicts clinical outcomes in resectable melanoma. Nat Med. (2019) 25:454–61. 10.1038/s41591-019-0357-y
58.
MessinaJLFenstermacherDAEschrichSQuXBerglundAELloydMCet al. 12-Chemokine gene signature identifies lymph node-like structures in melanoma: potential for patient selection for immunotherapy?Sci Rep. (2012) 2:765. 10.1038/srep00765
59.
RohWChenP-LReubenASpencerCNPrietoPAMillerJPet al. Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci Transl Med. (2017) 9:eaah3560. 10.1126/scitranslmed.aah3560
60.
RooneyMSShuklaSAWuCJGetzGHacohenN. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. (2015) 160:48–61. 10.1016/j.cell.2014.12.033
61.
AyersMLuncefordJNebozhynMMurphyELobodaAKaufmanDRet al. IFN-γ-related mRNA profile predicts clinical response to PD-1 blockade. J Clin Invest. (2017) 127:2930–40. 10.1172/JCI91190
62.
OckC-YHwangJ-EKeamBKimS-BShimJ-JJangH-Jet al. Genomic landscape associated with potential response to anti-CTLA-4 treatment in cancers. Nat Commun. (2017) 8:1050. 10.1038/s41467-017-01018-0
63.
CharoentongPFinotelloFAngelovaMMayerCEfremovaMRiederDet al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. (2017) 18:248–62. 10.1016/j.celrep.2016.12.019
64.
SzetoGLFinleySD. Integrative approaches to cancer immunotherapy. Trends Cancer Res. (2019) 5:400–10. 10.1016/j.trecan.2019.05.010
65.
WellensteinMDde VisserKE. Cancer-cell-intrinsic mechanisms shaping the tumor immune landscape. Immunity. (2018) 48:399–416. 10.1016/j.immuni.2018.03.004
66.
SprangerSGajewskiTF. Impact of oncogenic pathways on evasion of antitumour immune responses. Nat Rev Cancer. (2018) 18:139–47. 10.1038/nrc.2017.117
67.
SharmaPHu-LieskovanSWargoJARibasA. Primary, adaptive, and acquired resistance to cancer immunotherapy. Cell. (2017) 168:707–23. 10.1016/j.cell.2017.01.017
68.
PengWChenJQLiuCMaluSCreasyCTetzlaffMTet al. Loss of PTEN promotes resistance to T cell-mediated immunotherapy. Cancer Discov. (2016) 6:202–16. 10.1158/1538-7445.AM2016-4363
69.
AkbayEAKoyamaSCarreteroJAltabefATchaichaJHChristensenCLet al. Activation of the PD-1 pathway contributes to immune escape in EGFR-driven lung tumors. Cancer Discov. (2013) 3:1355–63. 10.1158/2159-8290.CD-13-0310
70.
EscorsDGato-CañasMZuazoMArasanzHGarcía-GrandaMJVeraRet al. The intracellular signalosome of PD-L1 in cancer cells. Signal Transduct Target Ther. (2018) 3:26. 10.1038/s41392-018-0022-9
71.
QuigleyDSilwal-PanditLDannenfelserRLangerødAVollanHKMVaskeCet al. Lymphocyte invasion in IC10/basal-like breast tumors is associated with wild-type TP53. Mol Cancer Res. (2015) 13:493–501. 10.1158/1541-7786.MCR-14-0387
72.
SprangerSBaoRGajewskiTF. Melanoma-intrinsic β-catenin signalling prevents anti-tumour immunity. Nature. (2015) 523:231–5. 10.1038/nature14404
73.
MamanSWitzIP. A history of exploring cancer in context. Nat Rev Cancer. (2018) 18:359–76. 10.1038/s41568-018-0006-7
74.
Altan-BonnetGMukherjeeR. Cytokine-mediated communication: a quantitative appraisal of immune complexity. Nat Rev Immunol. (2019) 19:205–17. 10.1038/s41577-019-0131-x
75.
ZaretskyJMGarcia-DiazAShinDSEscuin-OrdinasHHugoWHu-LieskovanSet al. Mutations associated with acquired resistance to PD-1 blockade in melanoma. N Engl J Med. (2016) 375:819–29. 10.1056/NEJMoa1604958
76.
BertrandFMontfortAMarcheteauEImbertCGilhodesJFilleronTet al. TNFα blockade overcomes resistance to anti-PD-1 in experimental melanoma. Nat Commun. (2017) 8:2256. 10.1038/s41467-017-02358-7
77.
JacquelotNYamazakiTRobertiMPDuongCPMAndrewsMCVerlingueLet al. Sustained type I interferon signaling as a mechanism of resistance to PD-1 blockade. Cell Res. (2019) 29:846–61. 10.1038/s41422-019-0224-x
78.
TüreiDKorcsmárosTSaez-RodriguezJ. OmniPath: guidelines and gateway for literature-curated signaling pathway resources. Nat Methods. (2016) 13:966–67. 10.1038/nmeth.4077
79.
RamilowskiJAGoldbergTHarshbargerJKloppmannELizioMSatagopamVPet al. A draft network of ligand-receptor-mediated multicellular signalling in human. Nat Commun. (2015) 6:7866. 10.1038/ncomms8866
80.
ChoiHShengJGaoDLiFDurransARyuSet al. Transcriptome analysis of individual stromal cell populations identifies stroma-tumor crosstalk in mouse lung cancer model. Cell Rep. (2015) 10:1187–201. 10.1016/j.celrep.2015.01.040
81.
BrowaeysRSaelensWSaeysY. NicheNet: modeling intercellular communication by linking ligands to target genes. Nat Methods. (2020) 17:159–62. 10.1038/s41592-019-0667-5
82.
WangSKarikomiMMacLeanALNieQ. Cell lineage and communication network inference via optimization for single-cell transcriptomics. Nucleic Acids Res. (2019) 47:e66. 10.1093/nar/gky882
83.
GhoshdastiderUNaeiniMMRohatgiNRevkovE. Data-driven inference of crosstalk in the tumor microenvironment. BioRxiv. (2019). 10.1101/835512
84.
KondratovaMCzerwinskaUSompairacNAmigorenaSDSoumelisVBarillotEet al. A multiscale signalling network map of innate immune response in cancer reveals cell heterogeneity signatures. Nat Commun. (2019) 10:4808. 10.1038/s41467-019-12270-x
85.
WorzfeldTFinkernagelFReinartzSKonzerAAdhikaryTNistAet al. Proteotranscriptomics reveal signaling networks in the ovarian cancer microenvironment. Mol Cell Proteomics. (2018) 17:270–89. 10.1074/mcp.RA117.000400
86.
ModugnoFDDi ModugnoFColosiCTronoPAntonacciGRuoccoGet al. 3D models in the new era of immune oncology: focus on T cells, CAF and ECM. J Exp Clin Cancer Res. (2019) 38:117. 10.1186/s13046-019-1086-2
87.
BarabásiALOltvaiZN. Network biology: understanding the cell's functional organization. Nat Rev Genet. (2004) 5:101–13. 10.1038/nrg1272
88.
LesterhuisWJBoscoAMillwardMJSmallMNowakAKLakeRA. Dynamic versus static biomarkers in cancer immune checkpoint blockade: unravelling complexity. Nat Rev Drug Discov. (2017) 16:264–72. 10.1038/nrd.2016.233
89.
LetaiA. Functional precision cancer medicine-moving beyond pure genomics. Nat Med. (2017) 23:1028–35. 10.1038/nm.4389
90.
BarretinaJCaponigroGStranskyNVenkatesanKMargolinAAKimSet al. The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. (2012) 483:603–7. 10.1038/nature11003
91.
IorioFKnijnenburgTAVisDJBignellGRMendenMPSchubertMet al. A landscape of pharmacogenomic interactions in cancer. Cell. (2016) 166:740–54. 10.1016/j.cell.2016.06.017
92.
Bar-EphraimYEKretzschmarKCleversH. Organoids in immunological research. Nat Rev Immunol. (2020) 20:279–93. 10.1038/s41577-019-0248-y
93.
WongAHHLiHJiaYMakPIMartinsRPdaSLiuYet al. Drug screening of cancer cell lines and human primary tumors using droplet microfluidics. Sci Rep. (2017) 7:9109. 10.1038/s41598-017-08831-z
94.
EduatiFUtharalaRMadhavanDNeumannUPLongerichTCramerTet al. A microfluidics platform for combinatorial drug screening on cancer biopsies. Nat Commun. (2018) 9:2434. 10.1038/s41467-018-04919-w
95.
MonteroJSarosiekKADeAngeloJDMaertensORyanJErcanDet al. Drug-induced death signaling strategy rapidly predicts cancer response to chemotherapy. Cell. (2015) 160:977–89. 10.1016/j.cell.2015.01.042
96.
RohrsJAWangPFinleySD. Understanding the dynamics of T-cell activation in health and disease through the lens of computational modeling. JCO Clin Cancer Inform. (2019) 3:1–8. 10.1200/CCI.18.00057
97.
TognettiMGaborAYangMCappellettiVWindhagerJCharmpiKet al. Deciphering the signaling network landscape of breast cancer improves drug sensitivity prediction. biorxiv. (2020). 10.1101/2020.01.21.907691
98.
EduatiFDoldàn-MartelliVKlingerBCokelaerTSieberAKogeraFet al. Drug resistance mechanisms in colorectal cancer dissected with cell type-specific dynamic logic models. Cancer Res. (2017) 77:3364–75. 10.1158/0008-5472.CAN-17-0078
99.
FeyDHalaszMDreidaxDKennedySPHastingsJFRauchNet al. Signaling pathway models as biomarkers: patient-specific simulations of JNK activity predict the survival of neuroblastoma patients. Sci Signal. (2015) 8:ra130. 10.1126/scisignal.aab0990
100.
BéalJMontagudATraynardPBarillotECalzoneL. Personalization of logical models with multi-omics data allows clinical stratification of patients. Front Physiol. (2018) 9:1965. 10.3389/fphys.2018.01965
101.
EduatiFJaaksPWapplerJCramerTMertenCAGarnettMJet al. Patient-specific logic models of signaling pathways from screenings on cancer biopsies to prioritize personalized combination therapies. Mol Syst Biol. (2020) 16:e8664. 10.15252/msb.20188664
102.
ArulrajTBarikD. Mathematical modeling identifies Lck as a potential mediator for PD-1 induced inhibition of early TCR signaling. PLoS ONE. (2018) 13:e0206232. 10.1371/journal.pone.0206232
103.
BolouriHYoungMBeilkeJJohnsonRFoxBHuangLet al. Integrative network modeling reveals mechanisms underlying T cell exhaustion. Sci Rep. (2020) 10:1915. 10.1038/s41598-020-58600-8
104.
NortonK-AGongCJamalianSPopelAS. Multiscale agent-based and hybrid modeling of the tumor immune microenvironment. Processes. (2019) 7:37. 10.3390/pr7010037
105.
KatherJNPoleszczukJSuarez-CarmonaMKrisamJCharoentongPValousNAet al. Modeling of immunotherapy and stroma-targeting therapies in human colorectal cancer. Cancer Res. (2017) 77:6442–52. 10.1158/0008-5472.CAN-17-2006
106.
KatherJNCharoentongPSuarez-CarmonaMHerpelEKluppFUlrichAet al. High-throughput screening of combinatorial immunotherapies with patient-specific models of metastatic colorectal cancer. Cancer Res. (2018) 78:5155–63. 10.1158/0008-5472.CAN-18-1126
107.
ThurleyKWuLFAltschulerSJ. Modeling cell-to-cell communication networks using response-time distributions. Cell Syst. (2018) 6:355–67.e5. 10.1016/j.cels.2018.01.016
108.
GrandclaudonMPerrot-DockèsMTrichotCKarpfLAbouzidOChauvinCet al. A quantitative multivariate model of human dendritic cell-T helper cell communication. Cell. (2019) 179:432–47.e21. 10.1016/j.cell.2019.09.012
109.
WangHMilbergOBartelinkIHViciniPWangBNarwalRet al. In silico simulation of a clinical trial with anti-CTLA-4 and anti-PD-L1 immunotherapies in metastatic breast cancer using a systems pharmacology model. R Soc Open Sci. (2019) 6:190366. 10.1098/rsos.190366
110.
SontagED. A dynamic model of immune responses to antigen presentation predicts different regions of tumor or pathogen elimination. Cell Syst. (2017) 4:231–41.e11. 10.1016/j.cels.2016.12.003
111.
SorribesICBasuABradyREnriquez-NavasPMFengXKatherJNet al. Harnessing patient-specific response dynamics to optimize evolutionary therapies for metastatic clear cell renal cell carcinoma – Learning to adapt. biorxiv. (2019). 10.1101/563130
112.
PerlsteinDShlagmanOKoganYHalevi-TobiasKYakobsonALazarevIet al. Personal response to immune checkpoint inhibitors of patients with advanced melanoma explained by a computational model of cellular immunity, tumor growth, and drug. PLoS ONE. (2019) 14:e0226869. 10.1371/journal.pone.0226869
113.
RiazNHavelJJMakarovVDesrichardAUrbaWJSimsJSet al. Tumor and microenvironment evolution during immunotherapy with nivolumab. Cell. (2017) 171:934–49.e16. 10.1016/j.cell.2017.09.028
114.
ChenPLRohWReubenACooperZASpencerCNPrietoPAet al. Analysis of immune signatures in longitudinal tumor samples yields insight into biomarkers of response and mechanisms of resistance to immune checkpoint blockade. Cancer Discov. (2016) 6:827–37. 10.1158/2159-8290.CD-15-1545
115.
HwangSKwonAYJeongJYKimSKangHParkJet al. Immune gene signatures for predicting durable clinical benefit of anti-PD-1 immunotherapy in patients with non-small cell lung cancer. Sci Rep. (2020) 10:643. 10.1038/s41598-019-57218-9
116.
PeskovKAzarovIChuLVoronovaVKosinskyYHelmlingerG. Quantitative mechanistic modeling in support of pharmacological therapeutics development in immuno-oncology. Front Immunol. (2019) 10:924. 10.3389/fimmu.2019.00924
117.
EisensteinM. Making cancer immunotherapy a surer bet. Nature. (2017) 552:S72–S73. 10.1038/d41586-017-08704-5
Summary
Keywords
tumor microenvironment, precision immuno-oncology, multi-omics profiling, systems biology, predictive biomarkers, cancer signaling networks, immune checkpoint blockers
Citation
Lapuente-Santana Ó and Eduati F (2020) Toward Systems Biomarkers of Response to Immune Checkpoint Blockers. Front. Oncol. 10:1027. doi: 10.3389/fonc.2020.01027
Received
30 January 2020
Accepted
22 May 2020
Published
24 June 2020
Volume
10 - 2020
Edited by
Enrica Calura, University of Padova, Italy
Reviewed by
Howard Donninger, University of Louisville, United States; Mirjana Efremova, Wellcome Sanger Institute (WT), United Kingdom
Updates
Copyright
© 2020 Lapuente-Santana and Eduati.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Federica Eduati f.eduati@tue.nl
This article was submitted to Cancer Genetics, a section of the journal Frontiers in Oncology
Disclaimer
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.