# Medical Treatment of Lung Cancer: Can Immune Cells Predict the Response? A Systematic Review

Philippe Rochigneux, Alejandro Garcia, Brice Chanez, Anne Madroszyk,

Daniel Olive, Edward Garon

► To cite this version:

Philippe Rochigneux, Alejandro Garcia, Brice Chanez, Anne Madroszyk, Daniel Olive, et al. Medical Treatment of Lung Cancer: Can Immune Cells Predict the Response? A Systematic Review. Frontiers in Immunology, 2020, 11, 10.3389/fimmu.2020.01036. hal-03152719

HAL Id: hal-03152719

https://amu.hal.science/hal-03152719

Submitted on 25 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Medical Treatment of Lung Cancer: Can Immune Cells Predict the Response? A Systematic Review

Philippe Rochigneux¹,²,³*, Alejandro J. Garcia⁴, Brice Chanez¹, Anne Madroszyk¹, Daniel Olive² and Edward B. Garon³

¹Department of Medical Oncology, Paoli-Calmettes Institute, Marseille, France, ²Team Immunity and Cancer, Centre de Recherche en Cancérologie de Marseille (CRCM), INSERM U1068, CNRS UMR 7258, Aix-Marseille Université and Institut Paoli-Calmettes, Marseille, France, ³Division of Hematology/Oncology, Department of Medicine, David Geffen School of Medicine at the University of California, Los Angeles, CA, United States, ⁴Cytometry Core Laboratory, David Geffen School of Medicine at the University of California, Los Angeles, CA, United States

The landscape for medical treatment of lung cancer has irreversibly changed since the development of immuno-oncology (IO). Yet, while immune checkpoint blockade (ICB) revealed that T lymphocytes play a major role in lung cancer, the precise dynamic of innate and adaptive immune cells induced by anticancer treatments including chemotherapy, targeted therapy, and/or ICB is poorly understood. In lung cancer, studies evaluating specific immune cell populations as predictors of response to medical treatment are scarce, and knowledge is fragmented. Here, we review the different techniques allowing the detection of immune cells in the tumor and blood (multiplex immunohistochemistry and immunofluorescence, RNA-seq, DNA methylation pattern, mass cytometry, functional tests). In addition, we present data that consider different baseline immune cell populations as predictors of response to medical treatments of lung cancer. We also review the potential for assessing dynamic changes in cell populations during treatment as a biomarker. As powerful tools for immune cell detection and data analysis are available, clinicians and researchers could increase understanding of mechanisms of efficacy and resistance in addition to identifying new targets for IO by developing translational studies that decipher the role of different immune cell populations during lung cancer treatments.

**Keywords:** lung cancer, immune biomarker, immune monitoring, chemotherapy, immunotherapy

## INTRODUCTION

Lung cancer is the leading cause of cancer deaths worldwide (1), and systemic treatments (chemotherapy, targeted therapy, or immunotherapy) are required for the vast majority of patients (clinical stages Ib–IV) (2). In the last 5 years, the development of immune checkpoint blockade (ICB) has improved the outcomes of advanced non-small cell lung cancer (NSCLC) (3–5). With this new focus on immuno-oncology (IO), ~200 lung cancer immunotherapy clinical trials are active worldwide (6). The standard initial treatment for advanced NSCLC without driver mutations now includes immunotherapy [monotherapy with high [programmed death ligand 1 (PDL-1) expression or combined with platinum doublet chemotherapy at any level of PDL-1 expression] (7, 8). Despite these treatment strategies, the immune cell dynamic induced by systemic treatments

**OPEN ACCESS**

**Edited by:**

Sandra Tuyaerts,
KU Leuven, Belgium

**Reviewed by:**

Rodabe N. Amaria,
University of Texas MD Anderson
Cancer Center, United States
Catherine Sautes-Fridman,
INSERM U1138 Centre de Recherche
des Cordeliers, France

**\*Correspondence:**

Philippe Rochigneux
rochigneuxp@ipc.unicancer.fr

**Specialty section:**

This article was submitted to
Cancer Immunity and Immunotherapy,
a section of the journal
Frontiers in Immunology

**Received:** 13 December 2019

**Accepted:** 29 April 2020

**Published:** 24 June 2020

**Citation:**

Rochigneux P, Garcia AJ, Chanez B, Madroszyk A, Olive D and Garon EB (2020) *Medical Treatment of Lung Cancer: Can Immune Cells Predict the Response? A Systematic Review*. *Front. Immunol.* 11:1036. doi: 10.3389/fimmu.2020.01036
### BOX 1 | Method: search strategy.

**Database:** PubMed, Google Scholar, Scopus, Web of Science, ASCO abstracts, ESMO abstracts, WCLC abstracts

**Database key words (MeSH terms and title/abstract):** "immune cells" OR "immune biomarker" OR "immune monitoring" OR "immuno-monitoring AND "lung cancer" OR "lung neoplasm" OR "non-small cell lung cancer" OR "small cell lung cancer"

**Time limit:** last 10 years (2009–2019)

**Language:** English only

**Species:** Human only

**Inclusion criteria:** Lung cancer treated with chemotherapy (CT), targeted therapy, or immunotherapy (IT), in an adjuvant setting, locally advanced or metastatic setting. Data available for immune cell population (e.g., CD3T cells) or immune cell markers (e.g., PDL-1 expression in macrophage or dendritic cells).

(chemotherapy, targeted therapy, or ICB) is poorly understood. Additional knowledge about this immune dynamic could be important to better understand the mechanisms of both efficacy and resistance to anticancer drugs. Further, it is possible that this exploration could identify new targets for immunotherapy. Compared to the field of melanoma (9), data are limited and very fragmented among lung cancer studies. Consequently, we conducted a systematic review of lung cancer studies that look at immune cell populations (Box 1). First, we will describe techniques that allow immune cells detection, focusing on recent discoveries. Second, we will review the immune predictors of response to medical treatments in lung cancer, studied at baseline or during treatment.

# CLASSICAL AND NEW TECHNICS OF IMMUNE MONITORING

## New Developments in Immunohistochemistry and Immunofluorescence

Since its invention by Albert Coons in the 1940s, immunohistochemistry (IHC) has been the gold standard in studying immune cell infiltration inside the tumor and the tumor immune microenvironment (TIME) (10) (Figure 1). Briefly, the classical IHC is the staining of formalin-fixed,

paraffin-embedded (FFPE) tissues with antibodies linked to an enzyme or fluorescent dye (14). This approach allows preservation of the tumor and TIME architecture, allowing spatial resolution. However, the main inconvenience of classical IHC is that a maximum of 2 antibodies are stained on one slide, necessitating a lot of tumor material and complicating the study of multiple cell populations. “Multiplex IHC” is a new technique allowing the visualization and quantification of specific immune cell populations by using multiple markers to identify different subsets [e.g., subset of dendritic cells (DCs) or CD8⁺ T cells] (15, 16). Multiplex IHC allows consecutive staining on a single slide (up to 10 antibodies) by using multiple rounds of staining and destaining (15). Several commercial solutions are available, mainly based on proprietary fluorescent probes and frozen materials (16). A nonproprietary assay called MICSSS (multiplexed immunohistochemical consecutive staining on single slide) has been developed based on chromogens and virtual color assignment (17). For example, in NSCLC patients (n = 75), MICSSS allowed staining of T cells, regulatory T cells (TReg), B cells, DCs, macrophages, and neutrophils together on a single slide, allowing for meaningful co-localizations (17). To improve the semiquantitative aspect of IHC, several algorithms have been developed to automatically analyze the slides of FFPE tumors (e.g., AQUA®), both for immunofluorescence (IF) and IHC, with an increasing use in translational research (18, 19). In the context of immune checkpoint inhibitors, quantitative IF and AQUA® were used to determine a “dormant” tumor-infiltrating lymphocyte (TIL) signature (elevated TILs with low activation and proliferation) associated with survival benefit (20). Recently, the performance of several biomarkers of anti-PDL-1 was studied in a meta-analysis of 45 studies; multiplex IHC/IF was associated with improved performance over PDL-1 IHC, tumor mutational burden, or gene expression (21). Finally, a new technique called “imaging mass cytometry” couples the principles of IHC and mass cytometry; tissue sections are stained with antibodies (up to 40) linked to rare metal isotopes, and an ultraviolet laser ablates the material spot by spot, which is then sent to the CyTOF (cytometry by time-of-flight) mass detector (see below cytometry section) (22).

## Immune Gene Expression Profiling

As cell types have distinct transcriptional profiles, it is possible to define immune cell populations by gene expression analysis. Traditional techniques, such as microarrays or NanoString nCounter® systems already, allow this identification, but the development of RNA-seq has transformed gene expression analyses into a powerful tool to identify cell populations. The main advantage of RNA-seq is to determine the primary sequence and relative abundance of each RNA molecule without previous knowledge of the sequence (using retrotranscription to cDNA and next-generation sequencing) (23). Furthermore, with microfluidics and barcodes (24), samples can be analyzed on a single-cell basis (single cell (sc) RNA-seq), allowing precise characterization of cell types in samples with cellular heterogeneity, such as lung (25–27). Sc RNA-seq also has the potential to define novel cell subtypes in blood (28) or solid tissues (29) and to follow cell differentiation with

**Abbreviations:** CyTOF, cytometry by time-of-flight; DCs, dendritic cells; ELISPOT, enzyme-linked immunospot; FFPE, formalin-fixed, paraffin-embedded; HR, hazard ratio; ICB, immune checkpoint blockade; IF, immunofluorescence; IFN, interferon; IHC, immunohistochemistry; IO, immuno-oncology; LAG3, lymphocyte-activation gene 3; MDSC, myeloid-derived suppressor cells; Mo-MDSC, monocytic myeloid-derived suppressor cells; NGS, next-generation sequencing; NLR, neutrophil-to-lymphocyte ratio; NSCLC, non-small cell lung cancer; OS, overall survival; OR, odds ratio; PBMC, peripheral blood mononuclear cell; PD1, programmed death 1; PDL-1, programmed death ligand 1; PFS, progression-free survival; PR, partial response; RNA-seq, RNA sequencing; SD, stable disease; STING, stimulator of interferon genes; TIGIT, T cell immunoreceptor with Ig and ITIM domains; TIL, tumor-infiltrating lymphocytes; TIME, tumor immune microenvironment; TKI, tyrosine kinase inhibitors; TLS, tertiary lymphoid structures; TMB, tumor mutational burden; TReg, regulatory T cells.
<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | New techniques for immune monitoring, in blood or tumor microenvironment of lung cancer. Colors in the text (bellow each image): green is for the main advantage; red is for the main disadvantage. Illustrations are adapted from the following references: (A) Rakaee M et al. (11); (B) Stern et al. (12); (C) Papalexi et al. (13); (D) Epicentral.com. Ab, antibody; IHC, immunohistochemistry; NGS, next-generation sequencing; scRNA-seq, single-cell RNA sequencing.</figcaption></figure>
RNA velocity (30). Of course, RNA-seq requires important computational statistical analysis, but machine learning methods such as Cibersort, XCell, and MetaNeighbor have been recently developed to simplify the characterization of cell composition from transcriptome data (31–33).

## Methylation Patterns

Epigenetic modifications, particularly DNA methylation, are crucial biological processes, allowing for the expression of specific cellular phenotypes from a common genetic background (34). DNA methylation is the addition of a methyl group to the C5 carbon residue of cytosines by DNA methyltransferases. Interestingly, these methylation patterns are cell type specific, and several studies describe that the methylome distinguishes cell lineages with high sensitivity and specificity (35–37). Technically, methylome analysis begins with a bisulfite conversion, as sodium bisulfite converts cytosines into uracils, whereas methylcytosines remain unmodified (38, 39). Subsequent amplification gives rise to two polymerase chain reaction products that are sequenced. From there, the sequences are aligned to a reference sequence, which can prove to be challenging (40, 41). As the tools for DNA methylation mapping are improving, and the required amount of DNA is decreasing, DNA methylome can now predict cell compositions in plasma (42). Consequently, plasmatic immuno-monitoring studies could be published in the future.

## Cytometry (Flow and Mass)

Invented in the 1960's, flow cytometry is a technique that studies the properties of a single cell in a liquid environment using fluorophore-linked antibodies (43). Lasers excite the fluorophore-linked antibodies at a certain spectrum, and the detectors record the emission spectrum (44). This signal is proportional to the expression of intracellular or extracellular cell markers, allowing for the identification of cell types-based markers of interest (45). Flow cytometry has quickly become a routine technique not only in malignant hematology and infectious disease studies, but also in drug development and drug monitoring (45–47). Successive technical improvements have broadened the number of excitation lasers (up to 10; BD InFlux) and the available fluorophores (notably with tandem and brilliant violet) allowing classical 18+ antibody panels. However, the multiplication of antibodies in a limited wavelength range (350–550 nm) leads to spectral overlap of fluorophores, requiring complex and time-consuming compensation setup before analyzing the cytometry data (48).

Around 2010, DVS Sciences Company and the Nolan Lab at Stanford University developed mass cytometry (CyTOF), a technique using antibodies tagged with rare earth metal isotopes (lanthanide series of the periodic table) (49). After staining with antibodies, the cells are nebulized and ionized with an argon inductively coupled plasma, and the ratio charge/mass is used to get a specific time of flight (TOF) proportional to the marker of interest (50). Thanks to the ability to distinguish isotopes clearly, up to 40 parameters can be studied in a single cell simultaneously, which is very useful for precious samples (51). Comparative analysis with flow cytometry showed that mass cytometry had a strong accuracy and reproducibility (52). The development of

CyTOF and its very large amount of biological data initiated a new era of high dimensional analyses (53). Several algorithms were created (Spade, viSNE, Citrus) to automatically cluster cell populations and provide a global map of cell signaling responses to interventions (54–57). Use of CyTOF led to several interesting results in predictors of response to ICB in melanoma (58, 59). Although the cost of the instrument is high, and the amount of the data generated requires some statistical expertise, the value of mass cytometry as a tool for immune monitoring is quickly increasing.

## Functional Tests

While understanding the cellular composition of TIME is important, knowing the functional status of the infiltrating immune cells tremendously enhances the value of the data. Researchers can study the proliferation of lymphocytes [a prognostic marker in several tumor types (60, 61)], thanks to the incorporation of ³H thymidine after stimulation with mitogens (62, 63). Cytotoxicity properties of CD8⁺ T cells or natural killer (NK) cells are also frequently measured with radioactive chromium (⁵¹Cr) assays in which target cells are loaded *in vitro* with radioactive chromium, and lysis is determined by measuring chromium released by dying cells in the supernatant (64, 65). Developed in 1968, this technique is still the gold standard for evaluation of cytotoxicity but requires handling radioactivity and autologous tumor targets. Interestingly, nonradioactive tests for cytotoxicity are emerging, based on flow cytometry (66) or microscopy of fluorescent target cells (67). Another classical way to measure T cell and NK cell cytotoxicity is enzyme-linked immunospot (ELISPOT), a sensitive immunoassay that measures the frequency of cytokine-secreting cells at the single-cell level (68, 69). As cytotoxic effector cells, on the one hand, will induce killing by different mechanisms (perforin, granzyme, Trail) but at the same time produce cytokines important for immune and inflammatory functions [such as interferon (IFN) and tumor necrosis factor (TNF)], IFN-γ ELISPOT is often used as a surrogate marker for cytotoxic properties of effector T cells. Other ELISPOT analyses also include granzyme B or perforin, two secreted proteins involved in perforation and caspase activation of the target cells (70). Lastly, measuring cytokine production with ELISPOT or intracellular flow cytometry allows assessment of the activation of CD4⁺ T cells (IL-2 production) or differentiation of subpopulations, such as CD4⁺T<sub>H1</sub> (IFN-γ, IL-2, IL-6, IL-12, IL-21, and TNF-α), T<sub>H2</sub> (IL-4, IL-5, and IL-13), and T<sub>H17</sub> (IL-17) cells.

## BASELINE PREDICTORS OF RESPONSE

### Immune Predictors of Chemotherapy and Targeted Therapy Response

The majority of studied biomarkers for efficacy of platinum doublets are not immune biomarkers but tumor biomarkers, studied at a genomic (71), transcriptomic (72), or protein level (73) (Figure 2). However, a few studies have raised the question of TILs in the context of chemotherapy. For example, in an analysis of 1,586 resected lung cancers treated with platinum-based adjuvant chemotherapy, an intense infiltration
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Baseline immune predictors of response to lung cancer medical treatments, in blood and tumor microenvironment. "Dormant TIL" is defined as CD3<sup>high</sup> granzyme B<sup>low</sup> Ki-67<sup>low</sup>. Abbreviations: CD4 THD, highly differentiated CD4+ T cells (CD27<sup>-</sup> CD28<sup>low/-</sup>); Mo-MDSC, monocytic myeloid-derived suppressor cells; PDL-1, programmed death ligand 1; TIME, tumor immune microenvironment; T<sub>Reg</sub>, regulatory T cells.</figcaption></figure>

of TIL in IHC was a positive prognostic marker, but had no predictive value for efficacy of the platinum doublets (74). In a metastatic setting, of 159 patients analyzed for TILs based on IHC, none of the T-cell subsets alone (CD8, CD4, T<sub>Reg</sub>) were associated with tumor response, but a low T<sub>Reg</sub>/CD8<sup>+</sup> ratio was associated with more tumor response to platinum doublets in multivariate analysis [odds ratio = 4.17, 95% confidence interval (CI) = 1.02–13.37, *p* = 0.029] (75). Finally, even if less studied, B cells can organize in the stroma with T cells and DCs into tertiary lymphoid structures (TLSs) that are ectopic lymphoid organs at site of inflammation (76). In 122 NSCLCs treated with neoadjuvant chemotherapy and surgery, density of follicular B cells and DCs organized in TLS was associated with an improved disease-specific survival (DSS) after 50 months’ follow-up (median DSS for B cells <sup>Hi</sup>/DCs <sup>Hi</sup> >60 vs. 21 months for B cells <sup>Low</sup>/DCs <sup>Low</sup>, *p* = 0.007) (77).

Among immune cells biomarkers in the blood, one of the most studied is monocytic myeloid-derived suppressor cells (Mo-MDSCs: CD33<sup>+</sup>, HLADR<sup>-</sup>, CD11b<sup>+</sup>, CD14<sup>+</sup>), an immature myeloid cell population inhibiting proliferation and cytotoxicity of T cells. In a cohort of 24 stage IV patients treated with chemotherapy (cisplatin/pemetrexed), progressors had higher rates of a subset of Mo-MDSCs, CD11b<sup>+</sup>, CD14<sup>+</sup>, and S100A9<sup>+</sup> (damage-associated molecular pattern molecules). Patients with a decreased frequency of these cells under the median had significantly longer progression-free survival (PFS) (9.2 vs. 3 months, *p* < 0.001) (78). In another study of 60 patients with advanced NSCLC treated with chemotherapy, progressive disease was associated with more baseline Mo-MDSCs (HLA-DR<sup>-</sup>/low, CD14<sup>+</sup>), and a baseline number of Mo-MDSCs under the median

was also associated with longer median PFS (9 vs. 3 months, *p* < 0.001) (79). This effect of Mo-MDSCs on PFS was consistent both in frequency (%) and absolute number (cells/µL). Another interesting blood biomarker is the pretreatment neutrophil-to-lymphocyte ratio (NLR). In 182 stage IV patients treated with platinum doublets (80), a high NLR (>2.63) was associated with worse PFS [hazard ratio (HR) = 1.81, *p* = 0.018] and overall survival (OS) (HR = 1.76, *p* = 0.02) in multivariate analysis, suggesting an independent detrimental effect of inflammation in response to chemotherapy. Next to NLR, a 2019 abstract suggests that high prechemotherapy absolute lymphocyte count is associated with favorable outcome in stage IB–III NSCLC patients who received adjuvant chemotherapy after surgical resection (81).

In the field of targeted therapy (e.g., EGFR and ALK inhibitors), data about baseline immune biomarkers are very scarce. Similar to chemotherapy, a retrospective analysis of 152 stage III/IV patients treated with EGFR tyrosine kinase inhibitors (TKIs) also found that a high NLR (>2.11) was an independent prognostic factor for longer OS (HR = 1.07, *p* = 0.03) (82). Additionally, in a prospective cohort of 33 patients treated with EGFR-TKI (mainly erlotinib) (83), flow cytometry analyses on peripheral blood mononuclear cells reported that a high baseline PDL-1<sup>+</sup> CD3<sup>+</sup> T cells predicted shorter OS in multivariate analysis (HR = 3.52, 95% CI = 1.09–11.4, *p* = 0.036).

## Immune Predictors of Checkpoint Inhibitors Response

In the TIME, several baseline biomarkers are associated with response to checkpoint inhibitors, which we can schematically
resume in 2 categories: (1) immune infiltrate or immune exclusion and (2) immune function or immune dysfunction of TIL infiltrate.

First, the pattern of immune infiltration (immune infiltrate vs. immune exclusion) seems crucial, as pejorative outcomes were described for lung tumors without immune CD8+ infiltration (immunological ignorance) or with immune CD8+ infiltrate outside of the tumor (excluded infiltrate) (84, 85). However, compared to melanoma, few data about immune exclusion are available in the ICB area in lung cancer (86). Recently, in 39 NSCLC patients treated with diverse ICB, the level of T-cell infiltration (CD3+) assessed by quantitative multiplex IF was 2.4-fold higher in patients with durable clinical benefit (n = 16) (20).

Second, data suggest that the presence of infiltration *per se* is not predictive, but depends of the precise nature of TIL infiltration and TIL functional state. Indeed, in the previous study, the highest response to ICB was observed in a specific subgroup called “dormant TIL signature” (CD3 high granzyme B low Ki-67 low, n = 7) (20). Interestingly, this signature was independent of PDL-1 expression in tumor cells or tumor mutational burden. Similarly, in the phase 1 study of atezolizumab in NSCLC patients (n = 53), best responses were observed in a specific subgroup of tumor-infiltrating immune cells with PDL-1 expression (macrophages, DCs, and lymphoid cells) (84). The presence of this immune infiltrate correlated with an improved objective response rate (ORR) (83% of ORR in the subgroup of immune PDL-1 expression >10%, compared to 14%-20% of ORR otherwise). Moreover, this association with atezolizumab response was stronger for PDL-1 expression in tumor-infiltrating immune cells than for PDL-1 expression in tumor cells. These data suggest that PDL-1 has a major role in pretreatment immunosuppression that can be reversed by that checkpoint inhibitor. Gene expression profile can also assess the functional state of TILs: in a prospective setting, IMpower150 study of PDL-1 inhibitor atezolizumab used the biomarker CD8+ Teff gene signature, defined by PDL-1, CXCL9, and IFN-γ mRNA expression from baseline tumor tissue (87). The benefit of the combination treatment atezolizumab, bevacizumab, carboplatin, paclitaxel (ABCP) vs. bevacizumab, carboplatin, paclitaxel (BCP) was particularly important in the Teff-high subgroup (median PFS ABCP = 11.3 months vs. median PFS BCP = 6.8 months) compared to the general population (median PFS ABCP = 8.3 months vs. median PFS BCP = 6.8 months). If confirmed, Teff mRNA signature is a promising surrogate of exhausted peritumoral T cells that can be reversed by ICB (85). Interestingly, this transcriptomic IFN-γ signature was independent of PDL-1 tumor expression (IHC). Increasing IFN-γ response is a current goal of pharmacological development to boost ICB, notably with STING (stimulator of interferon genes) agonist (88). Finally, there are limited data for alternative checkpoints such as LAG-3 expression on T cells, associated with shorter OS with PD-1 inhibitors (89). In conclusion, an approach combining both T cell exclusion and dysfunction [as the gene expression TIDE computational method (85)] may be particularly relevant to predict ICB response.

In the blood, contrary to the field of melanoma where
several baseline immune predictors of response are described

[e.g., baseline CD14⁺ CD16⁻ HLA-DR<sup>hi</sup> monocytes (58) or baseline Ki67⁺ PD-1⁺ CD8⁺ T cells (59)], only one large study is available in NSCLC, including 466 patients treated with diverse PD-1/PDL-1 inhibitors. The authors studied the impact of baseline Lung Immune Prognostic Index (LIPI), combining derived neutrophils/(leukocytes minus neutrophils) ratio (dNLR) and Lactate dehydrogenase (LDH) (90). Poor baseline LIPI, combining dNLR greater than 3 and LDH greater than upper limit of normal, was correlated with worse outcomes for ICB treatment in patients with NSCLC, but not with chemotherapy. Median OS was 4.8 (95% CI, 3.6-7.7) vs. 10.0 (95% CI, 7.3-12.6) vs. 16.5 (95% CI, 11.4-34.0) months for the poor, intermediate, and good LIPI groups, respectively. Similarly to chemotherapy (80, 82), these data raise the important question of the detrimental effect of baseline inflammation for ICB treatment in NSCLC. Recently, our group presented mass cytometry analysis about baseline predictors of pembrolizumab efficacy in NSCLC on KEYNOTE-001 using machine-learning algorithm (91). Three predictors of response were identified (classical monocytes perforin⁺ granzyme⁺/central memory CD4⁺ T cells ICOS⁺ CD28⁺ PD1⁺/41BB⁺ and perforin⁺ effector CD8⁺ T cells), and prospective validation is ongoing. Interestingly, the positive impact of classical monocyte in baseline has been previously published in melanoma (58).

Recently, the importance of myeloid cells is emerging in lung cancer; in 32 NSCLC patients treated with ICB (pembrolizumab, nivolumab, atezolizumab), a high proportion of myeloid cells expressing PDL-1 (PDL-1⁺ CD11b⁺ myeloid cells) was associated with objective response (92). Additionally, a functional CD4 immunity also seems important; among 51 NSCLC patients treated with anti-PDL-1, those with an increased proportion of highly differentiated CD4 (T<sub>HD</sub>: CD27<sup>-</sup> CD28<sup>low/-</sup>) had improved PFS/OS (93). The majority of these CD4⁺ T<sub>HD</sub> cells corresponded to nonsenescent, nonexhausted memory CD4 cells.

Finally, some data are emerging about the worrying phenomenon of hyperprogression; in 263 NSCLC patients treated with PD-1/PDL-1 inhibitors, a lower frequency of effector/memory CD8+ T cells (CCR7<sup>-</sup> CD45RA<sup>-</sup>) and a higher frequency of severely exhausted populations (TIGIT<sup>+</sup> T cells among PD-1<sup>+</sup> CD8<sup>+</sup> T cells) were associated with hyperprogression (n = 55/263) and inferior survival rate (94).

DYNAMIC PREDICTORS OF RESPONSE

Baseline predictive biomarkers may be sufficient to identify patients benefiting from medical treatments (Figure 3). However, dynamic biomarkers have theoretical added value. Immune system is a highly dynamic system with many switches, thresholds, and feedforward and feedback loops (95). Consequently, immunity is very sensitive to initial conditions, and minuscule differences may go undetected. Moreover, immune system has been described as a complex system with randomness and stochastic variations (e.g., lymphocyte fate
decision) (96), where small differences can lead to massive consequences downstream (minor bacterial exposure leading to a septic shock, targeting one single molecule leading to cytokines storm). If dynamic biomarkers are just emerging, they are promising tools to identify secondary mechanisms of resistance, especially if the selected timepoints are chosen wisely (at the time of progression or response). Finally, understanding these resistance mechanisms could help develop new IO combination therapies (97).

## Immune Predictors of Chemotherapy and Targeted Therapy Response

Similar to baseline markers, one of the most studied dynamic immune cells blood biomarker is MDSCs (CD33⁺ CD11b⁺). Liu et al. (98) reported variation of granulocytic MDSCs (CD33⁺ CD11b⁺ CD14⁻ CD15⁺) in advanced NSCLC patients treated with chemotherapy; nonprogressors (partial response or stable disease, n = 41) decreased their Gr-MDSC frequency compared

to baseline ($p < 0.0001$), contrary to progressors (n = 37), who kept similar frequency. The authors also described a negative correlation between Gr-MDSC frequency and CD8⁺ T cells in blood, consistent with the known immunosuppressive effect of MDSCs on T cells. Although less convincing statistically, another study of 46 patients with unresectable NSCLC, treated with platinum-based chemotherapy, reported that disease progression was associated with significantly higher levels of MDSC subpopulations (CD15⁺ and CD14⁺) compared to patients with disease control (100).

Except limited data about cytokines variations (101), the only immune variations described during EGFR inhibition are those described in the previously cited study about PDL-1⁺ CD3⁺ T cells (83). A value over the median after 1 week of EGFR inhibitors is associated with a decreased OS (HR = 6.49, 95% CI = 1.9–21.8, $p = 0.002$). These findings may reflect an immune resistance mechanism occurring in the PD-1/PDL-1 pathway after initiation of the targeted therapy (83).

## Immune Predictors of Checkpoint Inhibitors Response

Unfortunately, despite several studies describing the evolution of systemic immunity during checkpoint inhibitors in melanoma [e.g., increase of central memory CD4⁺ T cells (102)], there are few data available for lung cancer patients. In the previously cited phase 1 study of atezolizumab, an increase of CD8⁺ HLA-DR⁺ Ki67⁺ T cells in blood was seen at C2D1 (second infusion), but this variation was not correlated with atezolizumab response (84). Similarly, after 6 weeks of treatment, durvalumab significantly increased tumor gene expression of T cell chemotactic chemokine CXCL9, the checkpoint molecule LAG3, and IFN-γ, but the clinical effects of these immune variations was unclear (103).

However, some positive results are emerging; in an American Society of Clinical Oncology 2019 abstract, in 88 NSCLC patients who received anti-PD1 therapy in a single institution, NLR and relative lymphocyte count (RLC) were recorded at baseline and during the treatment (104). Using median RLC

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Dynamic changes in immune cell populations during chemotherapy or immune checkpoint blockade in lung cancer (data available in blood only). Abbreviations: CEMo, chemotherapy; ICB, immune checkpoint blockade; MDSC, myeloid-derived suppressor cells; PD-1, programmed death 1; PDL-1, programmed death ligand 1; RLC, relative lymphocyte count.</figcaption></figure>

<table><caption>**TABLE 1** | Immune cell populations predicting efficacy in lung cancer medical treatments (including cells function and references).</caption><thead><tr><th>Immune biomarker</th><th>Definition</th><th>Location</th><th>Outcome</th><th>References</th></tr></thead><tbody><tr><td>Treg/ CD8<sup>+</sup> ratio</td><td>CD4<sup>+</sup> FoxP3<sup>+</sup></td><td>Tumor</td><td>↑ Cemo response rate</td><td>(75)</td></tr><tr><td>Tertiary Lymphoid Structures</td><td>B cells, DCs, CD4<sup>+</sup> CD8<sup>+</sup> T cells</td><td>Tumor</td><td>↑ Cemo PFS</td><td>(77)</td></tr><tr><td>"T<sup>eff</sup>" lymphocytes signature</td><td>PD-L1, CXCL9, and IFNγ mRNA</td><td>Tumor</td><td>↑ ICB PFS/OS (signature)</td><td>(87)</td></tr><tr><td>"Dormant" lymphocytes</td><td>CD3<sup>+</sup> Granzyme B<sup>-</sup> Ki67<sup>-</sup></td><td>Tumor</td><td>↑ ICB response rate</td><td>(20)</td></tr><tr><td>CD4<sup>+</sup> highly differentiated</td><td>CD27<sup>-</sup> CD28<sup>low/-</sup></td><td>Blood</td><td>↑ ICB response rate</td><td>(93)</td></tr><tr><td>PDL1<sup>+</sup> immune cells</td><td>PDL1<sup>+</sup> macro, DCs, lymphocytes</td><td>Tum/Blood</td><td>↑ ICB response rate</td><td>(83, 84, 92)</td></tr><tr><td>"Reactivated" lymphocytes</td><td>CD8<sup>+</sup> PD1<sup>+</sup> Ki67<sup>+</sup></td><td>Blood</td><td>↑ ICB response rate</td><td>(99)</td></tr><tr><td>Neutrophil lymphocyte ratio (NLR)</td><td>Neutrophil/Lymphocyte</td><td>Blood</td><td>↓ Cemo PFS/OS (high NLR)</td><td>(80)</td></tr><tr><td>Lung immune prognostic index (LIPI)</td><td>dNLR >3; LDH > ULN</td><td>Blood</td><td>↓ ICB PFS/OS (high LIPI)</td><td>(90)</td></tr><tr><td>Myeloid derived suppressor cells</td><td>CD33<sup>+</sup> HLADR<sup>-</sup> CD11b<sup>+</sup></td><td>Blood</td><td>↓ Cemo/ICB response rate</td><td>(78, 98)</td></tr></tbody></table>

Data are presented for published studies (not for abstract only). Immune populations with a positive or detrimental outcome are highlighted in green or red, respectively. CXCL9, chemokine (C-X-C motif) ligand 9; DCs, dendritic cells; IFN, interferon; Cemo, chemotherapy; ICB, immune checkpoint blockade; NLR, neutrophil-to-lymphocyte ratio; PD-1, programmed death 1; PDL-1, programmed death ligand 1; PFS, progression-free survival; OS, overall survival, w, weeks.
at 4 weeks as a threshold, patients with high RLC at 4 weeks had significantly favorable survival (log-rank *p* < 0.0001). For patients with acquired resistance to therapy, RLC increased early during treatment followed by a decrease at the time of progression. Moreover, it is known that PD-1 inhibitors can rescue exhausted T cells (59). Interestingly, in a cohort of 29 NSCLC (responders: *n* = 10) treated with three different anti-PD-(L)1 agents (pembrolizumab, nivolumab, and atezolizumab), 80% of NSCLC patients with partial clinical responses presented early proliferative CD8+ T cells, which were both PD-1+ (exhausted) and Ki67+ (in replication) (99). Patients with an early proliferation of these PD-1+ Ki67+ CD8+ T cells in the blood within 4 weeks of treatment initiation had higher response rates (99). In another study in NSCLC patients (*n* = 13), the same team described that PD-1+ CD8+ T cells activated by PD-1 therapy were mostly CD28+, suggesting a potential role for this costimulation molecule in ICB response (105). Altogether, these studies describe a positive outcome of an early PD-1+ CD8+ T-cell blood response unleashed by blockade of the PD-1 pathway.

## CONCLUDING REMARKS

To be concise, we voluntarily limited this review to medical treatments of lung cancer and immune cells, even if some interesting immunological data are emerging about neoantigens or T Cell Receptor (TCR) repertoire (106). Inside this field, baseline data showed that a low number of Mo-MDSCs, a low NLR, and a high number of TLS improve platinum-based chemotherapy outcomes. Dynamic data showed that the decrease of MDSCs and the increase of PD1+ CD8+ T cells improve chemotherapy or immune checkpoint inhibitors outcomes (Table 1). However, the knowledge about immune

dynamics induced by medical treatments is still scarce, and none of these markers have sufficient level of evidence to be used in clinical practice. Academic research needs to lead in designing translational studies capable of rigorously evaluating this kind of data. Deciphering immune dynamics induced by medical treatments would aid in understanding the mechanisms of resistance and could also identify new targets for immunotherapy. Powerful and robust tools such as multiplex IHC/IF, RNA-seq, methylome, or mass cytometry are now available, together with software for data analysis. In the future, multi-omics approach will help integrate the different data about immune cells biomarkers (107). In that perspective, systems cancer immunology may soon help guide clinical decision-making (108). Lung cancer physicians and researchers should seize this opportunity to pursue immune monitoring together with drug development in IO.

## AUTHOR CONTRIBUTIONS

PR, DO, and EG designed the review. AG helped with the methodological and technical aspects. BC and AM helped with designing the figures and improving medical relevance. PR and EG wrote the paper with input from all authors. All authors contributed to the article and approved the submitted version.

## FUNDING

PR was supported by Fondation de France UNICANCER. EG was in part supported by NIH grant (R01 CA208403). DO was supported by research grant of FRM (Equipe FRM 2018 DEQ20180339209).

## REFERENCES

1. Research Topics-Lung Cancer. Available online at: http://www.iarc.fr/en/cancertopics/lungtopics.php (accessed August 9, 2018).

2. NCCN Guidelines-Non-Small Cell Lung Cancer. Available online at: https://www.nccn.org/professionals/physician_gls/pdf/nscl.pdf (accessed February 2020).

3. Garon EB, Rizvi NA, Hui R, Leighl N, Balmanoukian AS, Eder JP, et al. Pembrolizumab for the treatment of non-small-cell lung cancer. *N Engl J Med.* (2015) 372:2018-28. doi: 10.1056/NEJMoa1501824

4. Garon EB. The race for combined checkpoint inhibition in NSCLC. *Lancet Oncol.* (2016) 17:259-60. doi: 10.1016/S1470-2045(15)00580-X

5. Ellis PM, Vella ET, Ung YC. Immune checkpoint inhibitors for patients with advanced non-small-cell lung cancer: a systematic review. *Clin Lung Cancer.* (2017) 18:444-59.e1. doi: 10.1016/j.cllc.2017.02.001

6. ClinicalTrials. Search of: Immunotherapy | Recruiting, Active, Not Recruiting Studies | Lung Cancer - List Results. Available online at: https://clinicaltrials.gov (accessed August 9, 2018).

7. Gandhi L, Rodriguez-Abreu D, Gadgeel S, Esteban E, Felip E, De Angelis F, et al. Pembrolizumab plus chemotherapy in metastatic non-small-cell lung cancer. *N Engl J Med.* (2018) 378:2078-92. doi: 10.1056/NEJMoa1801005

8. Reck M, Rodriguez-Abreu D, Robinson AG, Hui R, Csoszi T, Fülöp A, et al. Pembrolizumab versus chemotherapy for PD-L1-positive non-small-cell lung cancer. *N Engl J Med.* (2016) 375:1823-33. doi: 10.1056/NEJMoa1606774

9. Hegde PS, Karanikas V, Evers S. The where, the when, and the how of immune monitoring for cancer immunotherapies in the era of checkpoint inhibition. *Clin Cancer Res.* (2016) 22:1865-74. doi: 10.1158/1078-0432.CCR-15-1507

10. Coons AH, Creech HJ, Jones RN. Immunological properties of an antibody containing a fluorescent group. *Proc Soc Exp Biol Med.* (1941) 47:200-2. doi: 10.3181/00379727-47-13084P

11. Rakaee M, Busund L-TR, Jamaly S, Paulsen E-E, Richardsen E, Andersen S, et al. Prognostic value of macrophage phenotypes in resectable non-small cell lung cancer assessed by multiplex immunohistochemistry. *Neoplasia.* (2019) 21:282-93. doi: 10.1016/j neo.2019.01.005

12. Stern L, McGuire H, Avdic S, Rizzetto S, Fazekas de St Groth B, Luciani F, et al. Mass cytometry for the assessment of immune reconstitution after hematopoietic stem cell transplantation. *Front Immunol.* (2018) 9:1672. doi: 10.3389/fimmu.2018.01672

13. Papalexi E, Satija R. Single-cell RNA sequencing to explore immune cell heterogeneity. *Nat Rev Immunol.* (2018) 18:35-45. doi: 10.1038/nri.2017.76

14. Mason DY, Sammons R. Alkaline phosphatase and peroxidase for double immunoenzymatic labelling of cellular constituents. *J Clin Pathol.* (1978) 31:454-60. doi: 10.1136/jcp.31.5.454

15. Gerdes MJ, Sevinsky CJ, Sood A, Adak S, Bello MO, Bordwell A, et al. Highly multiplexed single-cell analysis of formalin-fixed, paraffin-embedded cancer tissue. *Proc Natl Acad Sci USA.* (2013) 110:11982-87. doi: 10.1073/pnas.1300136110
16. Dixon AR, Bathany C, Tsuei M, White J, Barald KF, Takayama S. Recent developments in multiplexing techniques for immunohistochemistry. *Expert Rev Mol Diagn.* (2015) 15:1171-86. doi: 10.1586/14737159.2015.1069182

17. Remark R, Merghoub T, Grabe N, Litjens G, Damotte D, Wolchok JD, et al. In-depth tissue profiling using multiplexed immunohistochemical consecutive staining on single slide. *Sci Immunol.* (2016) 1:aaf6925. doi: 10.1126/sciimmunol.aaf6925

18. Gustavson MD, Bourke-Martin B, Reilly DM, Cregger M, Williams C, Tedeschi G, et al. Development of an unsupervised pixel-based clustering algorithm for compartmentalization of immunohistochemical expression using Automated QUantitative Analysis. *Appl Immunohistochem Mol Morphol AIMM.* (2009) 17:329-37. doi: 10.1097/PAI.0b013e318195ecaa

19. Bejnordi BE, Veta M, van Diest PJ, van Ginneken B, Karssemeijer N, Litjens G, et al. Diagnostic assessment of deep learning algorithms for detection of lymph node metastases in women with breast cancer. *JAMA.* (2017) 318:2199-210. doi: 10.1001/jama.2017.14585

20. Gettinger SN, Choi J, Mani N, Sanmamed MF, Datar I, Sowell R, et al. A dormant TIL phenotype defines non-small cell lung carcinomas sensitive to immune checkpoint blockers. *Nat Commun.* (2018) 9:3196. doi: 10.1038/s41467-018-05032-8

21. Lu S, Stein JE, Rimm DL, Wang DW, Bell JM, Johnson DB, et al. Comparison of biomarker modalities for predicting response to PD-1/PD-L1 checkpoint blockade: a systematic review and meta-analysis. *JAMA Oncol.* (2019) 5:1195-204. doi: 10.1001/jamaoncol.2019.1549

22. Giesen C, Wang HAO, Schapiro D, Zivanovic N, Jacobs A, Hattendorf B, et al. Highly multiplexed imaging of tumor tissues with subcellular resolution by mass cytometry. *Nat Methods.* (2014) 11:417-22. doi: 10.1038/nmeth.2869

23. Hrdlickova R, Toloue M, Tian B. RNA-Seq methods for transcriptome analysis. *Wiley Interdiscip Rev RNA.* (2017) 8. doi: 10.1002/wrna.1364

24. Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, et al. Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets. *Cell.* (2015) 161:1202-14. doi: 10.1016/j.cell.2015.05.002

25. Treutlein B, Brownfield DG, Wu AR, Neff NF, Mantalas GL, Espinoza FH, et al. Reconstructing lineage hierarchies of the distal lung epithelium using single-cell RNA-seq. *Nature.* (2014) 509:371-5. doi: 10.1038/nature13173

26. Muraro MJ, Dharmadhikari G, Grün D, Groen N, Dielen T, Jansen E, et al. A single-cell transcriptome atlas of the human pancreas. *Cell Syst.* (2016) 3:385-94.e3. doi: 10.1016/j.cels.2016.09.002

27. Habib N, Li Y, Heidenreich M, Swiech L, Avraham-Davidi I, Trombetta JJ, et al. Div-Seq: Single-nucleus RNA-Seq reveals dynamics of rare adult newborn neurons. *Science.* (2016) 353:925-8. doi: 10.1126/science.aad7038

28. Villani A-C, Satija R, Reynolds G, Sarkizova S, Shekhar K, Fletcher J, et al. Single-cell RNA-seq reveals new types of human blood dendritic cells, monocytes, and progenitors. *Science.* (2017) 356:eaah4573. doi: 10.1126/science.aah4573

29. La Manno G, Gyllborg D, Codeluppi S, Nishimura K, Salto C, Zeisel A, et al. Molecular diversity of midbrain development in mouse, human, and stem cells. *Cell.* (2016) 167:566-80.e19. doi: 10.1016/j.cell.2016.09.027

30. Manno GL, Soldatov R, Zeisel A, Braun E, Hochgerner H, Petukhov V, et al. RNA velocity of single cells. *Nature.* (2018) 560:494-8. doi: 10.1038/s41586-018-0414-6

31. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. *Nat Methods.* (2015) 12:453-7. doi: 10.1038/nmeth.3337

32. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. *Genome Biol.* (2017) 18:220. doi: 10.1186/s13059-017-1349-1

33. Crow M, Paul A, Ballouz S, Huang ZJ, Gillis J. Characterizing the replicability of cell types defined by single cell RNA-sequencing data using MetaNeighbor. *Nat Commun.* (2018) 9:884. doi: 10.1038/s41467-018-03282-0

34. Esteller M. Epigenetics in Cancer. *N Engl J Med.* (2008) 358:1148-59. doi: 10.1056/NEJMra072067

35. Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, et al. DNA methylation arrays as surrogate measures of cell mixture distribution. *BMC Bioinformatics.* (2012) 13:86. doi: 10.1186/1471-2105-13-86

36. Schuyler RP, Merkel A, Raineri E, Altucci L, Vellenga E, Martens JHA, et al. Distinct Trends of DNA methylation patterning

in the innate and adaptive immune systems. *Cell Rep.* (2016) 17:2101-11. doi: 10.1016/j.celrep.2016.10.054

37. Delacher M, Imbusch CD, Weichenhan D, Breiling A, Hotz-Wagenblatt A, Träger U, et al. Genome-wide DNA-methylation landscape defines specialization of regulatory T cells in tissues. *Nat Immunol.* (2017) 18:1160-72. doi: 10.1038/nli.3799

38. Krueger F, Kreck B, Franke A, Andrews SR. DNA methylome analysis using short bisulfite sequencing data. *Nat Methods.* (2012) 9:145-51. doi: 10.1038/nmeth.1828

39. Feng S, Rubbi L, Jacobsen SE, Pellegrini M. Determining DNA methylation profiles using sequencing. In: *High-Throughput Next Generation Sequencing Methods in Molecular Biology*. Totowa, NJ: Humana Press. p. 223-238. doi: 10.1007/978-1-61779-089-8_16

40. Ji H, Ehrlich LIR, Seita J, Murakami P, Doi A, Lindau P, et al. Comprehensive methylome map of lineage commitment from haematopoietic progenitors. *Nature.* (2010) 467:338-42. doi: 10.1038/nature09367

41. Farlik M, Sheffield NC, Nuzzo A, Datlinger P, Schönegger A, Klughammer J, et al. Single-cell DNA methylome sequencing and bioinformatic inference of epigenomic cell-state dynamics. *Cell Rep.* (2015) 10:1386-97. doi: 10.1016/j.celrep.2015.02.001

42. Moss J, Magenheim J, Neiman D, Zemmour H, Loyfer N, Korach A, et al. Comprehensive human cell-type methylation atlas reveals origins of circulating cell-free DNA in health and disease. *Nat Commun.* (2018) 9:5068. doi: 10.1038/s41467-018-07466-6

43. Herzenberg LA, Parks D, Sahaf B, Perez O, Roederer M, Herzenberg LA. The history and future of the fluorescence activated cell sorting and flow cytometry: a view from stanford. *Clin Chem.* (2002) 48:1819-27. doi: 10.1093/clinchem/48.10.1819

44. Cossarizza A, Chang H-D, Radbruch A, Akdis M, Andrä I, Annunziato F, et al. Guidelines for the use of flow cytometry and cell sorting in immunological studies. *Eur J Immunol.* (2017) 47:1584-797. doi: 10.1002/eji.201646632

45. Woo J, Baumann A, Arguello V. Recent advancements of flow cytometry: new applications in hematology and oncology. *Expert Rev Mol Diagn.* (2014) 14:67-81. doi: 10.1586/14737159.2014.862153

46. Krutzik PO, Crane JM, Clutter MR, Nolan GP. High-content single-cell drug screening with phosphospecific flow cytometry. *Nat Chem Biol.* (2008) 4:132-42. doi: 10.1038/nchembio.2007.59

47. Thibult M-L, Mamessier E, Gertner-Dardenne J, Pastor S, Just-Landi S, Xerri L, et al. PD-1 is a novel regulator of human B-cell activation. *Int Immunol.* (2013) 25:129-37. doi: 10.1093/intimm/dxs098

48. Szalóki G, Goda K. Compensation in multicolor flow cytometry. *Cytometry A.* (2015) 87:982-5. doi: 10.1002/cyto.a.22736

49. Bendall SC, Simonds EF, Qiu P, Amir ED, Krutzik PO, Finck R, et al. Single-cell mass cytometry of differential immune and drug responses across a human hematopoietic continuum. *Science.* (2011) 332:687-96. doi: 10.1126/science.1198704

50. Bendall SC, Nolan GP, Roederer M, Chattopadhyay PK. A deep profiler's guide to cytometry. *Trends Immunol.* (2012) 33:323-32. doi: 10.1016/j.it.2012.02.010

51. Spitzer MH, Nolan GP. Mass cytometry: single cells, many features. *Cell.* (2016) 165:780-91. doi: 10.1016/j.cell.2016.04.019

52. Gadalla R, Noamani B, MacLeod BL, Dickson RI, Guo M, Xu W, et al. Validation of CyTOF against flow cytometry for immunological studies and monitoring of human cancer clinical trials. *Front Oncol.* (2019) 9:415. doi: 10.3389/fonc.2019.00415

53. Newell EW, Sigal N, Bendall SC, Nolan GP, Davis MM. Cytometry by time-of-flight shows combinatorial cytokine expression and virus-specific cell niches within a continuum of CD8+ T cell phenotypes. *Immunity.* (2012) 36:142-52. doi: 10.1016/j.immuni.2012.01.002

54. Qiu P, Simonds EF, Bendall SC, Gibbs Jr KD, Bruggner RV, Linderman MD, et al. Extracting a cellular hierarchy from high-dimensional cytometry data with SPADE. *Nat Biotechnol.* (2011) 29:886-91. doi: 10.1038/nbt.1991

55. Amir ED, Davis KL, Tadmor MD, Simonds EF, Levine JH, Bendall SC, et al. viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. *Nat Biotechnol.* (2013) 31:545-52. doi: 10.1038/nbt.2594
56. Bruggner RV, Bodenmiller B, Dill DL, Tibshirani RJ, Nolan GP. Automated identification of stratifying signatures in cellular subpopulations. *Proc Natl Acad Sci USA.* (2014) 111:E2770-E7. doi: 10.1073/pnas.1408792111

57. Gondois-Rey F, Granjeaud S, Rouillier P, Rioualen C, Bidaut G, Olive D. Multi-parametric cytometry from a complex cellular sample: improvements and limits of manual versus computational-based interactive analyses. *Cytom Part J Int Soc Anal Cytol.* (2016) 89:480-90. doi: 10.1002/cyto.a.22850

58. Krieg C, Nowicka M, Guglietta S, Schindler S, Hartmann FJ, Weber LM, et al. High-dimensional single-cell analysis predicts response to anti-PD-1 immunotherapy. *Nat Med.* (2018) 24:144-53. doi: 10.1038/nm.4466

59. Huang AC, Postow MA, Orlowski RJ, Mick R, Bengsch B, Manne S, et al. T-cell invigoration to tumour burden ratio associated with anti-PD-1 response. *Nature.* (2017) 545:60-5. doi: 10.1038/nature22079

60. Nakano O, Sato M, Naito Y, Suzuki K, Orikasa S, Aizawa M, et al. Proliferative Activity of Intratumoral CD8+ T-lymphocytes as a prognostic factor in human renal cell carcinoma: clinicopathologic demonstration of antitumor immunity. *Cancer Res.* (2001) 61:5132-6.

61. Filaci G, Fenoglio D, Fravega M, Ansaldo G, Borgonovo G, Traverso P, et al. CD8+CD28- T regulatory lymphocytes inhibiting t cell proliferative and cytotoxic functions infiltrate human cancers. *J Immunol.* (2007) 179:4323-34. doi: 10.4049/jimmunol.179.7.4323

62. Taylor JH, Woods PS, Hughes WL. The organization and duplication of chromosomes as revealed by autoradiographic studies using tritium-labeled thymidine. *Proc Natl Acad Sci USA.* (1957) 43:122-8.

63. Wilson DB, Blyth JL, Nowell PC. Quantitative studies on the mixed lymphocyte interaction in rats. *J Exp Med.* (1968) 128:1157-81.

64. Brunner KT, Mauel J, Cerottini JC, Chapuis B. Quantitative assay of the lytic action of immune lymphoid cells on 51-Cr-labelled allogeneic target cells in vitro; inhibition by isoantibody and by drugs. *Immunology.* (1968) 14:181-96.

65. Pasero C, Gravis G, Guerin M, Granjeaud S, Thomassin-Piana J, Rocchi P, et al. Inherent and tumor-driven immune tolerance in the prostate microenvironment impairs natural killer cell antitumor activity. *Cancer Res.* (2016) 76:2153-65. doi: 10.1158/0008-5472.CAN-15-1965

66. Zaritskaya L, Shurin MR, Sayers TJ, Malyguine AM. New flow cytometric assays for monitoring cell-mediated cytotoxicity. *Expert Rev Vaccines.* (2010) 9:601-16. doi: 10.1586/erv.10.49

67. Fassy J, Tsalkitzi K, Goncalves-Maia M, Braud VM. A real-time cytotoxicity assay as an alternative to the standard chromium-51 release assay for measurement of human NK and T cell cytotoxic activity. *Curr Protoc Immunol.* (2017) 118:7.42.1-12. doi: 10.1002/cpim.28

68. Malyguine AM, Strobl S, Dunham K, Shurin MR, Sayers TJ. ELISPOT assay for monitoring cytotoxic T lymphocytes (CTL) activity in cancer vaccine clinical trials. *Cells.* (2012) 1:111-126. doi: 10.3390/cells1020111

69. Barabas S, Spindler T, Kiener R, Tonar C, Lugner T, Batzilla J, et al. An optimized IFN-γ ELISpot assay for the sensitive and standardized monitoring of CMV protein-reactive effector cells of cell-mediated immunity. *BMC Immunol.* (2017) 18:14. doi: 10.1186/s12865-017-0195-y

70. Ranieri E, Popescu I, Gigante M. CTL ELISPOT assay. *Methods Mol Biol Clifton NJ.* (2014) 1186:75-86. doi: 10.1007/978-1-4939-1158-5_6

71. Zhu X, Lin MCM, Fan W, Tian L, Wang J, Ng SS, et al. An intronic polymorphism in grp78 improves chemotherapeutic prediction in non-small cell lung cancer. *Chest.* (2012) 141:1466-72. doi: 10.1378/chest.11-0469

72. Zhang K, Keymeulen S, Nelson R, Tong TR, Yuan Y-C, Yun X, et al. Overexpression of flap endonuclease 1 correlates with enhanced proliferation and poor prognosis of non-small-cell lung cancer. *Am J Pathol.* (2018) 188:242-51. doi: 10.1016/j.ajpath.2017.09.017

73. Wang X, Zhu X, Zhang H, Fan X, Xue X, Chen Y, et al. ERCC1_202 Is A prognostic biomarker in advanced stage non-small cell lung cancer patients treated with platinum-based chemotherapy. *J Cancer.* (2017) 8:2846-53. doi: 10.7150/jca.19897

74. Brambilla E, Le Teuff G, Marguet S, Lantuejoul S, Dunant A, Graziano S, et al. Prognostic effect of tumor lymphocytic infiltration in resectable non-small-cell lung cancer. *J Clin Oncol.* (2016) 34:1223-30. doi: 10.1200/JCO.2015.63.0970

75. Liu H, Zhang T, Ye J, Li H, Huang J, Li X, et al. Tumor-infiltrating lymphocytes predict response to chemotherapy in patients with advance

non-small cell lung cancer. *Cancer Immunol Immunother.* (2012) 61:1849-56. doi: 10.1007/s00262-012-1231-7

76. Sautès-Fridman C, Petitprez F, Calderaro J, Fridman WH. Tertiary lymphoid structures in the era of cancer immunotherapy. *Nat Rev Cancer.* (2019) 19:307-25. doi: 10.1038/s41568-019-0144-6

77. Germain C, Gnjatic S, Tamzalit F, Knockaert S, Remark R, Goc J, et al. Presence of B cells in tertiary lymphoid structures is associated with a protective immunity in patients with lung cancer. *Am J Respir Crit Care Med.* (2014) 189:832-44. doi: 10.1164/rccm.201309-1611OC

78. Feng P-H, Lee K-Y, Chang Y-L, Chan Y-F, Kuo L-W, Lin T-Y, et al. CD14+S100A9+ monocytic myeloid-derived suppressor cells and their clinical relevance in non-small cell lung cancer. *Am J Respir Crit Care Med.* (2012) 186:1025. doi: 10.1164/rccm.201204-0636OC

79. Huang A, Zhang B, Wang B, Zhang F, Fan K-X, Guo Y-J. Increased CD14(+)-HLA-DR (-/low) myeloid-derived suppressor cells correlate with extrathoracic metastasis and poor response to chemotherapy in non-small cell lung cancer patients. *Cancer Immunol Immunother CII.* (2013) 62:1439-451. doi: 10.1007/s00262-013-1450-6

80. Yao Y, Yuan D, Liu H, Gu X, Song Y. Pretreatment neutrophil to lymphocyte ratio is associated with response to therapy and prognosis of advanced non-small cell lung cancer patients treated with first-line platinum-based chemotherapy. *Cancer Immunol Immunother CII.* (2013) 62:471-9. doi: 10.1007/s00262-012-1347-9

81. Lee HW, Kang SY, Han JH, Yoo BM. Not preoperative but prechemotherapy lymphocyte count predicts outcome of non-small cell lung cancer patients treated with adjuvant chemotherapy. *J Glob Oncol.* (2019) 5:76-76. doi: 10.1200/JGO.2019.5.suppl.76

82. Minami S, Ogata Y, Ihara S, Yamamoto S, Komuta K. Neutrophil-to-lymphocyte ratio predicts overall survival of advanced non-small cell lung cancer harboring mutant epidermal growth factor receptor. *World J Oncol.* (2017) 8:180-7. doi: 10.14740/wjon1069w

83. Meniawy TM, Lake RA, McDonnell AM, Millward MJ, Nowak AK. PD-L1 on peripheral blood T lymphocytes is prognostic in patients with non-small cell lung cancer (NSCLC) treated with EGFR inhibitors. *Lung Cancer.* (2016) 93:9-16. doi: 10.1016/j.lungcan.2015.12.006

84. Herbst RS, Soria J-C, Kowanetz M, Fine GD, Hamid O, Gordon MS, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. *Nature.* (2014) 515:563-7. doi: 10.1038/nature14011

85. Jiang P, Gu S, Pan D, Fu J, Sahu A, Hu X, et al. Signatures of T cell dysfunction and exclusion predict cancer immunotherapy response. *Nat Med.* (2018) 1. doi: 10.1038/s41591-018-0136-1

86. Bonaventura P, Shekarian T, Alcazer V, Valladeau-Guilemond J, Valsesia-Wittmann S, Amigorena S, et al. Cold tumors: a therapeutic challenge for immunotherapy. *Front Immunol.* (2019) 10:168. doi: 10.3389/fimmu.2019.00168

87. Socinski MA, Jotte RM, Cappuzzo F, Orlandi F, Stroyakovskiy D, Nogami N, et al. Atezolizumab for first-line treatment of metastatic nonsquamous NSCLC. *N Engl J Med.* (2018) 378:2288-301. doi: 10.1056/NEJMoa1716948

88. Ramanjulu JM, Pesiridis GS, Yang J, Concha N, Singhaus R, Zhang S-Y, et al. Design of amidobenzimidazole STING receptor agonists with systemic activity. *Nature.* (2018) 564:439-43. doi: 10.1038/s41586-018-0705-y

89. Datar I, Sanmamed M, Choi J, Wang J, S Henick B, Badri T, et al. In patients with advanced non-small cell lung cancer (NSCLC) LAG-3 is expressed on activated TILs and predicts resistance to PD-1 axis blockers. *Ann Oncol.* (2017) 28: doi: 10.1093/annonc/mdx710.006

90. Mezquita L, Auclin E, Ferrara R, Charrier M, Remon J, Planchard D, et al. Association of the lung immune prognostic index with immune checkpoint inhibitor outcomes in patients with advanced non-small cell lung cancer. *JAMA Oncol.* (2018) 4:351-7. doi: 10.1001/jamaoncol.2017.4771

91. Rochigneux P, Lisberg A, Garcia A, Chretien AS, Fattori S, Madroszyk A, et al. P1.04-33 deep phenotyping of immune populations reveals baseline predictors of pembrolizumab efficacy in NSCLC on KEYNOTE-001. *J Thorac Oncol.* (2019) 14:S452-S3. doi: 10.1016/j.jtho.2019.08.936

92. Bocanegra A, Fernandez-Hinojal G, Zuazo-Ibarra M, Arasanz H, Garcia-Granda MJ, Hernandez C, et al. PD-L1 expression in systemic immune cell populations as a potential predictive biomarker of responses to
PD-L1/PD-1 blockade therapy in lung cancer. *Int J Mol Sci.* (2019) 20:1631. doi: 10.3390/ijms20071631

93. Zuazo M, Arasanz H, Fernández-Hinojal G, García-Granda MJ, Gato M, Bocanegra A, et al. Functional systemic CD4 immunity is required for clinical responses to PD-L1/PD-1 blockade therapy. *EMBO Mol Med.* (2019) 11:e10293. doi: 10.15252/emmm.201910293

94. Kim CG, Kim KH, Pyo K-H, Xin C-F, Hong MH, Ahn B-C, et al. Hyperprogressive disease during PD-1/PD-L1 blockade in patients with non-small-cell lung cancer. *Ann Oncol.* (2019) 30:1104–113. doi: 10.1093/annonc/mdz123

95. Lesterhuis WJ, Bosco A, Millward MJ, Small M, Nowak AK, Lake RA. Dynamic versus static biomarkers in cancer immune checkpoint blockade: unravelling complexity. *Nat Rev Drug Discov.* (2017) 16:264–72. doi: 10.1038/nrd.2016.233

96. Hodgkin PD. A probabilistic view of immunology: Drawing parallels with physics. *Immunol Cell Biol.* (2007) 85:295–9. doi: 10.1038/sj.icb.7100061

97. Yu JX, Hodge JP, Oliva C, Neftelinov ST, Hubbard-Lucey VM, Tang J. Trends in clinical development for PD-1/PD-L1 inhibitors. *Nat Rev Drug Discov.* (2019) 19:163–4. doi: 10.1038/d41573-019-00182-w

98. Liu C-Y, Wang Y-M, Wang C-L, Feng P-H, Ko H-W, Liu Y-H, et al. Population alterations of l-arginase- and inducible nitric oxide synthase-expressed CD11b+/CD14-/CD15+/CD33+ myeloid-derived suppressor cells and CD8+ T lymphocytes in patients with advanced-stage non-small cell lung cancer. *J Cancer Res Clin Oncol.* (2010) 136:35–45. doi: 10.1007/s00432-009-0634-0

99. Kamphorst AO, Pillai RN, Yang S, Nasti TH, Akondy RS, Wieland A, et al. Proliferation of PD-1+ CD8 T cells in peripheral blood after PD-1–targeted therapy in lung cancer patients. *Proc Natl Acad Sci USA.* (2017) 114:4993–8. doi: 10.1073/pnas.1705327114

100. Koinis F, Vetsika EK, Aggouraki D, Skalidaki E, Koutoulaki A, Gkioulmvasani M, et al. Effect of first-line treatment on myeloid-derived suppressor cells’ subpopulations in the peripheral blood of patients with non-small cell lung cancer. *J Thorac Oncol Off Publ Int Assoc Study Lung Cancer.* (2016) 11:1263–72. doi: 10.1016/j.jtho.2016.04.026

101. Lisberg A, Horton J, Goldman JW, Wolf B, Tucker DA, Carroll J, et al. P3.02b-042 reduction in peripheral blood cytokine levels observed in egfr mutant (egfrm) patients treated with erlotinib: topic: EGFR biomarkers. *J Thorac Oncol.* (2017) 12:S1212–S3. doi: 10.1016/j.jtho.2016.11.1709

102. Takeuchi Y, Tanemura A, Tada Y, Katayama I, Kumanogoh A, Nishikawa H. Clinical response to PD-1 blockade correlates with a sub-fraction of peripheral central memory CD4+ T cells in patients with malignant melanoma. *Int Immunol.* (2018) 30:13–22. doi: 10.1093/intimm/dxx073

103. Higgs BW, Morehouse CA, Streicher K, Brohawn PZ, Pilataxi F, Gupta A, et al. Interferon gamma messenger RNA signature in tumor biopsies predicts outcomes in patients with non-small cell lung carcinoma or urothelial cancer treated with durvalumab. *Clin Cancer Res.* (2018) 24:3857–66. doi: 10.1158/1078-0432.CCR-17-3451

104. Rosner S, Forde PM, Naidoo J, Marrone K, Reuss JE, Feliciano JL, et al. Early shifts in immune cell subsets to predict response to immune checkpoint blockade in non-small cell lung cancer (NSCLC). *J Clin Oncol.* (2019) 37:105. doi: 10.1200/JCO.2019.37.8_suppl.105

105. Kamphorst AO, Wieland A, Nasti T, Yang S, Zhang R, Barber DL, et al. Rescue of exhausted CD8T cells by PD-1-targeted therapies is CD28-dependent. *Science.* (2017) doi: 10.1126/science. aaf0683

106. Anagnostou V, Forde PM, White JR, Niknafs N, Hraban C, Naidoo J, et al. Dynamics of tumor and immune responses during immune checkpoint blockade in non-small cell lung cancer. *Cancer Res.* (2019) 79:1214–25. doi: 10.1158/0008-5472.CAN-18-1127

107. Lee JS, Ruppin E. Multiomics prediction of response rates to therapies to inhibit programmed cell death 1 and programmed cell death 1 ligand 1. *JAMA Oncol.* (2019)5:1614–18. doi: 10.1001/jamaoncol. 2019.2311

108. Greenplate AR, Johnson DB, Ferrell PB, Irish JM. Systems immune monitoring in cancer therapy. *Eur J Cancer Oxf Engl 1990.* (2016) 61:77–84. doi: 10.1016/j.ejca.2016.03.085

**Conflict of Interest:** EG reports funds to his institution from AstraZeneca, Boehringer Ingelheim, Bristol Myers Squibb, Eli Lilly, Genentech, Mirati, Merck, and Novartis. DO was co-founder of Imcheck Therapeutics and Emergence Therapeutics and has research funds from Imcheck Therapeutics, Alderaan, BMS, Talix, Cellectis and Emergence Therapeutics.

The remaining 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.

Copyright © 2020 Rochigneux, Garcia, Chanez, Madroszyk, Olive and Garon. 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.