Immunologic constant of rejection as a predictive biomarker of immune checkpoint inhibitors efficacy in non-small cell lung cancer
Journal of Translational Medicine volume 21, Article number: 637 (2023)
Anti-PD1/PDL1 immune checkpoint inhibitors (ICI) transformed the prognosis of patients with advanced non-small cell lung cancer (NSCLC). However, the response rate remains disappointing and toxicity may be life-threatening, making urgent identification of biomarkers predictive for efficacy. Immunologic Constant of Rejection signature (ICR) is a 20-gene expression signature of cytotoxic immune response with prognostic value in some solid cancers. Our objective was to assess its predictive value for benefit from anti-PD1/PDL1 in patients with advanced NSCLC.
We retrospectively profiled 44 primary tumors derived from NSCLC patients treated with ICI as single-agent in at least the second-line metastatic setting. Transcriptomic analysis was performed using the nCounter® analysis system and the PanCancer Immune Profiling Panel. We then pooled our data with clinico-biological data from four public gene expression data sets, leading to a total of 162 NSCLC patients treated with single-agent anti-PD1/PDL1. ICR was applied to all samples and correlation was searched between ICR classes and the Durable Clinical Benefit (DCB), defined as stable disease or objective response according to RECIST 1.1 for a minimum of 6 months after the start of ICI.
The DCB rate was 29%; 22% of samples were classified as ICR1, 30% ICR2, 22% ICR3, and 26% ICR4. These classes were not associated with the clinico-pathological variables, but showed enrichment from ICR1 to ICR4 in quantitative/qualitative markers of immune response. ICR2-4 class was associated with a 5.65-fold DCB rate when compared with ICR1 class. In multivariate analysis, ICR classification remained associated with DCB, independently from PDL1 expression and other predictive immune signatures. By contrast, it was not associated with disease-free survival in 556 NSCLC TCGA patients untreated with ICI.
The 20-gene ICR signature was independently associated with benefit from anti-PD1/PDL1 ICI in patients with advanced NSCLC. Validation in larger retrospective and prospective series is warranted.
Despite significant improvement in systemic treatments during the last few years, advanced non-small-cell lung cancer (NSCLC) remains the leading cause of cancer-related death worldwide . Recent advances in lung cancer biology knowledge led to molecular dismemberment of the disease, leading to the development of targeted therapies. This changed the lung cancer prognosis thanks to personalized therapeutic approach . Beyond oncogenic addictions, cancer immunology breakthrough allowed for the emergence of immune checkpoint inhibitors (ICI) . These drugs, which target the PD1-PDL1 axis, are now widely used in the first- and second-line settings for advanced NSCLC [4,5,6]. Unfortunately and in contrast with targeted therapies, only 20% of patients respond to ICI as single agent and up to 50% when combined with chemotherapy [4,5,6]. This suggests an incomplete knowledge of the molecular determinants of immune responsiveness. Furthermore, life-threatening or fatal immune-related adverse events occur in some patients. Consequently, there is still a long way to achieve immune dissection of targeted immunotherapies-mediated effects, and an urgent need to identify biomarkers able to predict anti-PD1/PDL1 ICI efficacy .
To date, evaluation of PDL1 expression on tumors cells by immunohistochemistry (IHC) using the PD-L1 IHC 22C3 pharmDx assay (CPS > 50%) is the only approved companion diagnostics for pembrolizumab, whereas the PD-L1 IHC 22C3 pharmDx and Ventana PD-L1 (SP142) have status as complementary diagnostics for nivolumab and atezolizumab. It is however an imperfect predictive biomarker which is not sufficient to face the most important existing challenge in the field of immunotherapy . The main reason given is the operator’s variability to evaluate PDL1 protein expression and its spatial heterogeneity . Furthermore, the use of different IHC tests with the different anti-PD1 and anti-PD-L1 ICI adds to these difficulties. Several other predictive biomarkers have been studied, but none of them appears to be reliable and reproducible enough to improve on the predictive significance of PDL1 expression, despite a better consideration of various biological features underlying ICI efficacy. Tumor mutational burden (TMB), thought to reflect the amount of neo-antigens on tumor cells by quantifying non-synonymous mutations in coding areas, failed to show a clear survival difference when used alone and remains difficult to standardize for routine use [10, 11]. In a same way, quantification of tumor-infiltrating lymphocytes and other circulating biomarkers appeared to be interesting prognostic biomarkers, but not predictive of ICI efficacy [12, 13]. In addition, interferon gamma (IFN-Ɣ) pathway is known to be involved in PDL1 expression and is a key feature to several anti-tumor immune expression signatures, such as the T cell-inflamed signature (TIS) . The independent predictive value of this signature was validated in a prospective series of 37 NSCLC treated with nivolumab, independently from PDL1 expression by IHC , and more recently in patients treated with pembrolizumab, in a pan-cancer clinical trial performed across 20 tumor types, including lung cancer . Finally, methylation profiles and microbiota assessment turned out to be hardly operable in daily care [17,18,19].
Thus, determinants of ICI response in NSCLC appeared to be many and various, depending on tumor and microenvironment features. Considering the relative failure of these descriptive biomarkers, more functional biomarkers such as transcriptomic signatures could become key parameters to solve this complex equation. Besides the TIS signature, other signatures including IFN-Ɣ-related genes showed also positive predictive value with atezolizumab in 115 NSCLC in the POPLAR study  and with durvalumab in 97 NSCLC . Gene expression signatures (GES) associated with enrichment in tertiary lymphoid structures (TLS) have also been generated in other cancers, including the Coppola’s 12-chemokine signature . Other prognostic and/or predictive immune signatures were reported in series ranging from 21 to 67 patients with NSCLC [23,24,25,26]. The ICR (Immunologic Constant of Rejection) signature is defined as an immune phenotype quantifying the expression of 20 genes, all involved in anti-tumor immunity. Selected genes are reflecting main anti-tumor immune pathways, such as Th1 signaling (IFNG, TBX21, CD8A/B, IL12B, STAT1 and IRF1), Th1 chemoattraction (CXCL9, CXCL10 and CCL5) and cytotoxic functions (GNLY, PRF1, GZMA, GZMB and GZMH). Interestingly, the expression of these pro-cytotoxic transcripts is associated with inhibition of suppressive mechanisms, known as ICIs (CD274, PDCD1, IDO1, CTLA4 and FOXP3) .
In early-stage breast cancer, we showed that ICR signature divided the tumors in four ICR groups reflecting a clinically and biologically relevant immune continuum, and displayed independent predictive values for metastasis-free survival (MFS) and for achievement of pathological complete response (pCR) to neo-adjuvant chemotherapy [28, 29]. ICR4 tumors strongly expressed ICR signature and showed longer MFS and higher pCR rate to chemotherapy than other ICR tumors. Similarly, we recently showed the independent prognostic value of ICR in soft tissue sarcomas . Regarding ICI, a potential positive predictive value for response was reported in melanoma .
Given these correlations in other cancers, its immune relevance, and the growing place of ICIs in advanced NSCLC treatment, we assessed the ICR signature as a potential predictive biomarker for ICI response in a cohort of 162 patients with advanced NSCLC treated with anti-PD1/PDL1 ICI. Our secondary objectives were to compare the predictive value of ICR with that of other immune signatures and to assess its prognostic value in NSCLC untreated with ICI.
Patients’ populations and gene expression profiling
Our own series included formalin-fixed paraffin-embedded (FFPE) tumor samples from the previously biopsied primary tumor of 77 consecutive patients with advanced NSCLC who had received an anti-PD1/PDL1 ICI as single agent outside clinical trials in our institution (Hôpital Nord, Marseille, France). All clinical data were collected from the patients’ electronic medical record. This non-interventional retrospective study was approved by our institutional review board under the number 2019_93. All patients had given their signed written informed consent for the use of archived material for research purpose. Tumor RNA was extracted using “Maxwell® RSC Instrument” (Promega) shortly after microtome dissection to prevent nucleic acid degradation. After extraction, RNA quantification and quality control were done using the Nano-Drop ND-1000 technology (ThermoFisher). Forty-four out of 77 tumor samples fulfilled the quality and quantity thresholds for downstream transcriptomic analysis. This latter was performed using nCounter® technology Dx analysis as recommended by Nanostring, based on microscopic imaging counting relative abundance of 770 transcripts present in the “PanCancer Immune Profiling Panel”, which included the 20 ICR genes. Briefly, the profiling required three successive steps: a first hybridization step to capture target sequences, followed by a purification step, and finally a quantification step using the Digital Analyzer processor.
In order to expand our series, we collected gene expression and clinico-pathological data from four publicly available data sets [15, 23, 32, 33] of NSCLC patients treated with a single-agent anti-PD1/PDL1 ICI and with available clinical outcome. In these published data sets, gene expression profiling had been done using Nanostring technology or RNA-Sequencing. The final pooled data set included 162 metastatic NSCLC patients clinically annotated, notably in term of clinical benefit after ICI treatment (Additional file 2: Table S1). Finally, in order to assess the eventual prognostic value of ICR outside any ICI treatment, we collected The Cancer Genome Atlas (TCGA) lung adenocarcinoma (LUAD) data set  and the TCGA lung squamous cell carcinoma (LUSC) data set  including 515 and 502 patients respectively with gene expression data.
Gene expression data analysis
Nanostring data processing and normalization were performed using the nSolver™ 4.0 analysis software. Briefly, data processing of raw counts was done with background subtraction defined by the geometric mean of the eight negative control probes. Next, normalization was done with the geometric mean algorithm using the 40 housekeeping and the six positive control probes. Processed data were then log2-transformed. Before analysis of pooled data sets, several steps of data processing were applied. The first step was the normalization of each set separately. It was done in R using Bioconductor and associated packages; we used quantile normalization for the available processed log2-transformed data. We then applied to each data set separately several multigene signatures.
First, the ICR classifier based on consensus clustering (CC) analysis of the expression levels of 20 representative immune genes (namely CCL5, CD274, CD8A, CD8B, CTLA4, CXCL9, CXCL10, FOXP3, GNLY, GZMA, GZMB, GZMH, IDO1, IFNG, IL12B, IRF1, PDCD1, PRF1, STAT1, and TBX21) as previously described . Briefly, the CC analysis was performed in R using the Bioconductor package “ConsensusClusterPlus”  setting as input parameters 5000 repetitions, 80% item resampling (pItem), a number of groups (k) fixed to 4 (in order to have all datasets stratified with the same number of classes, 4 being the optimal number of groups for the TCGA cohort , and the use of an agglomerative hierarchical clustering with ward criterion (Ward.D2) inner and complete outer linkage. Second, we applied several other transcriptional signatures related to immune response: metagenes of signatures of 28 innate and adaptative immune cell subpopulations defined by Bindea et al. , activation score of IFN-α, IFN-Ɣ, and TNF-α pathways , cytolytic activity score , and antigen processing and presentation machinery score (APMS) . We also applied the TP53 activation score , and three potential predictors of response to ICI: PDL1 (CD274) expression, T cell-inflamed signature (TIS) , and tertiary lymphoid structures (TLS) signature . These three predictors were tested as binary variables using the first quintile as cut-off, thus defining the “low” (first quintile) and “high” (four last quintiles) classes. Third, PDL1 gene expression levels (CD274) were extracted from each data set and were standardized within each data set using the NSCLC population as a reference to be comparable across data sets and to exclude bias from population heterogeneity. Of note, we verified the correlation between mRNA and protein expressions of PDL1 in lung cancer samples (cell lines and clinical tumors; Additional file 1: Fig S1).
The continuous variables were described by median and range, and the binary variables by numbers and percentage. Correlations between tumor classes and clinico-pathological or molecular variables were analyzed using the one-way analysis of variance (ANOVA) or the Fisher’s exact test when appropriate. Our primary endpoint was the Durable Clinical Benefit (DCB), a relevant criterion of ICI efficacy considering the challenge and the expected paradigm shift to overcome natural course of the disease and induce prolonged response . We defined DCB as stable disease or objective response according to RECIST 1.1 for a minimum of 6 months after the start of ICI. This definition has been mainly previously published among ICI literature and especially biomarkers literature [42,43,44]. Uni- and multivariate analyses for DCB were done using logistic regression (glm; significance estimated by specifying a binomial family for models with a logit link). The variables submitted to univariate analysis included patients’ age (continuous value), sex (male vs female), smoker status (current vs former vs non-smoker), pathological type (squamous vs non-squamous), mutational status (mutated vs non-mutated), and classifications based on ICR (ICR2-4 vs ICR1), TIS (“high” vs “low”), TLS (“high” vs “low”) and PDL1 expression (“high” vs “low”). Variables with a p-value < 0.05 in univariate analyses were tested in multivariate analyses. In this meta-analysis of five independent sets, we used the test of homogeneity in a fixed-effects model to provide evidence about whether the effect sizes are measuring a common effect size. The prognostic analysis in the TCGA set used the Disease-Free Survival (DFS) as endpoint, calculated from the date of diagnosis until the date of distant relapse. Follow-up was measured from the date of diagnosis to the date of last news for event-free patients. Survivals were calculated using the Kaplan–Meier method and curves were compared with the log-rank test. All statistical tests were two-sided at the 5% level of significance. Statistical analysis was done using the survival package (version 2.43-3) in the R software (version 3.5.1; http://www.cran.r-project.org/). We followed the reporting REcommendations for tumor MARKer prognostic studies (REMARK criteria) .
Patients’ population and ICR classification
We profiled our series of 44 tumor samples that we coupled with 118 profiled samples from four published public datasets [15, 23, 32, 33], obtaining a cohort of 162 advanced NSCLC patients treated with anti-PD1/PDL1 ICI. Baseline patients’ characteristics are summarized in Table 1. They were representative of patients treated by daily care: the median patients’ age was 60 years, 67% of informative patients were male, 89% were active or former smokers, 65% of tumors were non-squamous pathological type, and 38% displayed somatic mutations (KRAS, then EGFR, then NRAS, STK11, and ROS1). Regarding our primary endpoint, 47 patients (29%) experienced a DCB after ICI treatment and 115 (71%) did not. ICR classification of the 162 tumors defined 35 tumors (22%) as ICR1, 48 (30%) as ICR2, 36 (22%) as ICR3, and 43 (26%) as ICR4 (Fig. 1A).
ICR classification and clinico-pathological and immune variables
We first searched for correlations between the four ICR classes and different variables. There was no significant correlation with the following clinico-pathological features: patients’ age and sex, smoker status, pathological type, and mutational status (Additional file 2: Table S2). By contrast, a significant correlation existed with the achievement of DCB (p = 2.67E-02, Fisher’s exact test), with DCB rates equal to 9% in ICR1 versus 35% in ICR2, 36% in ICR3, and 33% in ICR4 (Fig. 1B).
We also found strong correlations with immunity-related features with a continuum between the four classes (from ICR1 to ICR4) for nearly all features (Fig. 1C; Additional file 2: Table S2). All but one Bindea’s signatures for immune cell subsets  showed a strong enrichment from ICR1 to ICR4, notably cytotoxic cells, T-cells, CD8 + T-cells, Th1 cells, TFH cells and activated NK CD56dim cells (p < 1.00E-15). Among T-helper cells, the Th1/Th2 ratio increased from ICR1 to ICR4. This increasing anti-tumor activation was associated with increasing immune cell subsets involved in antigen presentation, such as activated dendritic cells (aDC), DC, B-cells, and macrophages. Such enrichment in cell subsets was confirmed using more functional immune signatures, with gradual enrichment from ICR1 to ICR4 for activation scores of IFN-α, IFN-Ɣ, and TNF-α pathways (p < 1.00E-08), for the cytolytic activity score  (p < 1.00E-30) and for the APMS score  (p < 1.00E-15). Conversely, the TP53 pathway activation score  decreased from ICR1 to ICR4 (p < 1.00E-05).
ICR classification and DCB after ICI treatment
Based on the absence of difference in the DCB rate between the ICR2, 3 and 4 classes (p = 0.937, Fisher’s exact test), we pooled them into the ICR2-4 class that we then compared to the ICR1 class. The DCB rates were 9% in ICR1 (3/35 patients) versus 35% in ICR2-4 (44/127 patients), corresponding to a 5.65 Odds Ratio (OR) (95% CI 1.64–19.51; p = 6.10E-03, logit function) in ICR2-4 versus ICR1 in univariate analysis (Fig. 2). Interestingly, analysis of homogeneity in a fixed-effects model revealed homogeneity between the five data sets in term of correlation between DCB rate and ICR1/2–4 classes (p = 0.950; Additional file 1: Fig S2). Notably, in our own cohort of 44 patients, the DCB rates were 11% in ICR1 (1/9 patients) versus 40% in ICR2-4 (10/35 patients), corresponding to a 3.20 OR (95% CI 0.35–29.01). None of the clinico-pathological variables tested in univariate analysis displayed significant correlation with the achievement of DCB (Fig. 2): patients’ age and sex, smoker status, pathological type, and mutational status.
Then, we compared this predictive value of ICR for DCB to that of PDL1 expression and two immune signatures (TIS, TLS) previously reported as predictive for response to ICI. In univariate analysis (Fig. 2), the four variables showed a positive association with DCB with an OR for DCB superior to 2 in the “TIS-high”, “TLS-high”, “PDL1-high”, and ICR2-4 patients when compared with the “TIS-low”, “TLS-low”, “PDL1-low”, and ICR1 patients respectively. However, likely due to the number of patients, the correlation was not significant for TIS and PDL1. However, it was significant for TLS (p = 2.27E-02, logit function) and ICR (p = 6.10E-03, logit function) signatures. In multivariate analysis (Fig. 2), only ICR remained significant, suggesting stronger and independent predictive value.
ICR classification and correlation with survival without ICI treatment
We assessed the eventual prognostic value of ICR in an early-stage NSCLC population (TCGA dataset of 556 M0 patients including 408 with lung adenocarcinoma and 148 with lung squamous carcinoma) naive from ICI and with available post-operative disease-free survival (DFS). A total of 164 samples were classified as ICR1 (29%), 127 as ICR2 (23%), 188 as ICR3 (34%), and 77 as ICR4 (14%), with similar correlations with immune variables as those described in our population treated with ICI (Additional file 2: Table S3). With a median follow-up of 158 months, 101 (18%) displayed a DFS event and the 5 year DFS was 67% (95% CI 61–74). It was 72% (95% CI 63–82) in the ICR1 class, and 66% (95% CI 58–74) in the ICR2-4 class, suggesting no prognostic value of ICR classification (p = 0.782, log-rank test; Fig. 3).
Immune checkpoint inhibitors are considered as a revolution for patients with advanced NSCLC and became a therapeutic standard for NSCLC without oncogenic addiction in European and US Guidelines regardless of PDL1 expression [46, 47]. Nevertheless, extended efficacy with acceptable safety profile remains too scarce and identification of predictive biomarkers for ICI efficacy became a major challenge. Here, we showed the independent predictive value of the ICR signature for efficacy of single-agent anti-PD1/PDL1 ICI in the largest retrospective multicentric clinical reported cohort of patients with NSCLC.
Our approach tested ICR in an independent series of samples, thus avoiding the problem of overfitting. A total of 162 pre-therapeutic cancer samples informative for DCB after ICI and large-scale gene expression profile were available, allowing not only to test our hypothesis in uni- and multivariate analyses, but also to test many other gene signatures and modules relevant to immune response. Our series was profiled using the Nanostring® technology, already available in clinical routine for breast cancer prognostication (CE-IVD label) and more adapted for daily care, notably for small FFPE or frozen biopsies as often available in NSCLC. The DCB was chosen to explore ICR predictive value because of its clinical relevance and consistency in the immunotherapy field [43, 48].
As described in other cancers such as breast cancer , sarcomas , or colon cancer , we found an immunological continuum in NSCLC from ICR1 to ICR4 classes. There was an increasing enrichment in scores reflecting the amount of different immune cell types, notably T-cells, cytotoxic T-cells, Th1-cells, CD8 + T-cells, and antigen-presenting cells, and functional scores reflecting IFN-Ɣ pathway activation, cytolytic activity, and antigen presentation machinery. Conversely, the activation score of TP53 pathway decreased from ICR1 to ICR4, in agreement with the higher rate of inactivating TP53 mutations reported in ICR4 previously reported in breast cancer .
A strong association between ICR classes and DCB was observed: the DCB rates were 9% in ICR1 patients versus 35% in ICR2-4 patients, corresponding to a 5.65 OR. Importantly, there was homogeneity through the five pooled data sets in term of correlation between DCB rate and ICR classes. This predictive value strengthens the potential clinical role of ICR signature and the correlation between tumor immunogenicity and ICI efficacy. Interestingly, ICR was more predictive than other potential biomarkers such as PDL1 expression, TIS and TLS signatures, whereas, as expected, no tested clinico-pathological feature showed any predictive value. Although associated with the same immune and biological variables in the TCGA set than in our 162 patients, the ICR signature was not associated with survival in TCGA patients untreated with ICI, suggesting that ICR is not broadly prognostic in NSCLC and might be only informative in the presence of an ICI treatment. However, since the lymphocyte infiltration seems to have some prognostic role in NSCLC , this absence of prognostic value of ICR in our TCGA analysis may appear surprising and clearly deserves further comparative assessment of both biomarkers in the same clinical series.
All 20 genes of ICR signature are involved in tumor immune pathways and previously published data highlighted their functional implication and clinical relevance. Despite the failure of therapeutic application, IFN-Ɣ remains the most described cytokine involved in anti-tumor immune response. On one hand, IFN-Ɣ upregulation at the early phase of immune response is obviously a major step of anticancer immunity using various downstream mechanisms. On the other hand, IFN-Ɣ is also known for its regulatory role in a further step of immune response, also able to enhance tumor growth and PDL1 expression . Moreover, TBX21, known as a transcription factor of IFN-Ɣ, already showed a negative association with NSCLC prognosis by mediating tumor growth . IL12 is also involved in triggering IFN-Ɣ to switch CD4 helper differentiation into Th1 cells and was already described across expression signatures as a prognostic biomarker in NSCLC . Downstream IFN-Ɣ pathway is inducing STAT1 activation achieving various biological functions such as IRF1 expression and tumor growth regulation, but STAT1 activation is also described for enhancing tumor progression through chronic inflammation responsible for controversial therapeutic implication . ICR signature also includes genes involved in Th1 chemo-attraction through IFN-Ɣ pathway, like CXCL9 and CXCL10, regulating naïve T cells migration, differentiation and activation and inhibiting angiogenesis . Other genes included in ICR signature are mainly involved in direct cytotoxic function, such as CD8, GNLY, GZMA, GZMB and GZMH. Finally, immunosuppressive mechanisms known as checkpoint inhibitors are considered in the signature and widely involved in ICI efficacy such as CD274, PDCD1 or CTLA4 but also IDO1 and FOXP3. IDO1 is involved in tryptophan metabolism and known to enhance the PI3K/AKT signaling and also with the T regulatory cell generation in tumor microenvironment, thus associated with poor prognosis in various tumor types [56, 57]. ICR provides a surrogate for the ability of immune response to attack the cancer cells if the checkpoint is inhibited.
Our result is consistent with other studies that reported a predictive value of immune signatures for the benefit of ICI in lung cancer . These signatures, ranging from 4 to 59 genes, include genes coding for proteins involved in tumor antigenicity and T cell priming/activation, trafficking of T cells and infiltration into tumors, recognition of cancer cells by T cells, infiltration by inhibitory cells, and immune checkpoints receptors or ligands. Many of them include IFN-Ɣ-related genes, such as the 18-gene TIS associated with objective response and PFS after pembrolizumab in the KEYNOTE-028 trial , the 8-gene “T-effector and IFN-Ɣ signature” associated with overall survival (OS) after atezolizumab in 115 NSCLC patients in the POPLAR study , or the 4-gene “IFN-Ɣ signature” associated with objective response and PFS after durvalumab in 97 NSCLC patients . Other predictive immune signatures were reported in smaller NSCLC series ranging from 21 to 67 patients [23,24,25,26]. Interestingly, the ICR signature remained the sole variable significant in our multivariate analysis for DCB prediction, when confronted to the TIS and TLS signatures. The predictive value of ICR for the response to systemic anti-cancer treatments has already been reported in retrospective studies. In early-stage breast cancer, we showed that ICR signature was independently associated with achievement of pathological complete response (pCR) to neo-adjuvant chemotherapy [28, 29]. Regarding ICI, a potential positive predictive value of ICR was reported for clinical response in melanoma , and more recently for the pathological response to neo-adjuvant chemotherapy combined with pembrolizumab in breast cancer . Such predictive value should not be limited to these two cancers, and clearly further studies are warranted in other indications.
Our study displays strengths and limitations. The first limitations are related to the retrospective nature and its associated biases such as heterogeneity and missing data. For example, our population included patients treated by ICI as single agent in different metastatic lines, whereas a wide majority of patients now receive ICI in the first-line setting; the samples used for the transcriptomic analyses were obtained at various stages of the disease (early and late) leading to some heterogeneity in our population considering temporal tumor heterogeneity; relevant data were missing such as the TMB or PDL1 expression by IHC, even if we replaced it by mRNA expression that correlates with protein expression in lung cancers (Additional file 1: Fig S1). Second, we used the DCB as efficacy endpoint, rather than overall survival (OS) and progression free survival (PFS), but DCB has now been widely used among cancer immunotherapy literature, reinforcing its clinical significance. Finally, the number of samples analyzed was relatively small when compared to the frequency of NSCLC, but to our knowledge, our series is the largest series of gene expression profiling of NSCLC patients treated with ICI. The strength of our results lies in: (i) the number of 162 samples that, to our knowledge, makes our series the largest gene expression study reported so far in this setting; (ii) its originality, being the first one to describe the predictive value of ICR signature for response to ICI in NSCLC; (iii) a population of patients treated in the community setting, likely more reflective of real-life than a selected patient population within a clinical trial; (iv) the biological relevance of ICR, its predictive value in the context of ICI treatment and absence of prognostic value in absence of ICI; and (v) the small number of genes included (20 genes), which should facilitate its clinical application once validated.
We showed the predictive value of ICR in a large, composite, and multicentric cohort. Even if a validation is required in larger retrospective cohorts and in prospective trials to confirm the robustness of our results in the first-line setting and chemotherapy combination, ICR displays a promising signal of efficacy in the so far disappointing immune biomarkers field. In a near future, predicting ICI efficacy will be a major challenge with three pivotal questions: (i) when? Balanced with the new data on ICI use in early-stage NSCLC and the re-challenge in widely pretreated patients [59,60,61]; (ii) how? By assessing various techniques to develop composite predictors likely including clinical, biological, genomic, transcriptomic, proteomics, and spatial data to better consider the different drivers of immune response; and (iii) what? Considering the amount of new potential therapeutic immune targets, starting by the 20 genes of the signature, and the wide possibilities of treatment associations and sequences.
One-way analysis of variance
Antigen processing and presentation machinery score
Centre de Recherche en Cancérologie de Marseille
Durable clinical benefit
Disease free survival
Gene expression signatures
Immune checkpoint inhibitors
Immunologic constant of rejection
Lung squamous cell carcinoma
Non-small cell lung cancer
Pathological complete response
Programmed cell death protein 1
Programmed death ligand 1
REcommendations for tumor MARKer prognostic studies
The cancer genome atlas
Tumor mutational burden
T cell-inflamed signature
Tertiary lymphoid structures
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 Countries. CA Cancer J Clin. 2021;71(3):209–49.
Tan AC, Tan DSW. Targeted therapies for lung cancer patients with oncogenic driver molecular alterations. J Clin Oncol. 2022;40(6):611–25.
Miller JFAP, Sadelain M. The journey from discoveries in fundamental immunology to cancer immunotherapy. Cancer Cell. 2015;27(4):439–49.
Borghaei H, Paz-Ares L, Horn L, Spigel DR, Steins M, Ready NE, et al. Nivolumab versus docetaxel in advanced nonsquamous non–small-cell lung cancer. N Engl J Med. 2015;373(17):1627–39.
Reck M, Rodríguez-Abreu D, Robinson AG, Hui R, Csőszi 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(19):1823–33.
Gandhi L, Rodríguez-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(22):2078–92.
Mino-Kenudson M, Schalper K, Cooper W, Dacic S, Hirsch FR, Jain D, et al. Predictive biomarkers for immunotherapy in lung cancer: perspective from the international association for the study of lung cancer pathology committee. J Thorac Oncol. 2022;17(12):1335–54.
Aguiar PN, De Mello RA, Hall P, Tadokoro H, de Lima Lopes G. PD-L1 expression as a predictive biomarker in advanced non-small-cell lung cancer: updated survival data. Immunotherapy. 2017;9(6):499–506.
Kerr KM, Tsao MS, Nicholson AG, Yatabe Y, Wistuba II, Hirsch FR, et al. Programmed death-ligand 1 immunohistochemistry in lung cancer: in what state is this art? J Thorac Oncol. 2015;10(7):985–9.
Greillier L, Tomasini P, Barlesi F. The clinical utility of tumor mutational burden in non-small cell lung cancer. Transl Lung Cancer Res. 2018;7(6):639–46.
Hellmann MD, Callahan MK, Awad MM, Calvo E, Ascierto PA, Atmaca A, et al. Tumor mutational burden and efficacy of nivolumab monotherapy and in combination with ipilimumab in small-cell lung cancer. Cancer Cell. 2018;33(5):853-861.e4.
Zeng DQ, Yu YF, Ou QY, Li XY, Zhong RZ, Xie CM, et al. Prognostic and predictive value of tumor-infiltrating lymphocytes for clinical therapeutic research in patients with non-small cell lung cancer. Oncotarget. 2016;7(12):13765–81.
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(3):351–7.
Ayers M, Lunceford J, Nebozhyn M, Murphy E, Loboda A, Kaufman DR, et al. IFN-γ-related mRNA profile predicts clinical response to PD-1 blockade. J Clin Invest. 2017;127(8):2930–40.
Damotte D, Warren S, Arrondeau J, Boudou-Rouquette P, Mansuet-Lupo A, Biton J, et al. The tumor inflammation signature (TIS) is associated with anti-PD-1 treatment benefit in the CERTIM pan-cancer cohort. J Transl Med. 2019;17(1):357.
Ott PA, Bang YJ, Piha-Paul SA, Razak ARA, Bennouna J, Soria JC, et 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(4):318–27.
Dora D, Ligeti B, Kovacs T, Revisnyei P, Galffy G, Dulka E, et al. Non-small cell lung cancer patients treated with Anti-PD1 immunotherapy show distinct microbial signatures and metabolic pathways according to progression-free survival and PD-L1 status. Oncoimmunology. 2023;12(1):2204746.
Duruisseaux M, Besse B, Cadranel J, Pérol M, Mennecier B, Bigay-Game L, et al. Overall survival with crizotinib and next-generation ALK inhibitors in ALK-positive non-small-cell lung cancer (IFCT-1302 CLINALK): a French nationwide cohort retrospective study. Oncotarget. 2017;8(13):21903–17.
Routy B, Gopalakrishnan V, Daillère R, Zitvogel L, Wargo JA, Kroemer G. The gut microbiota influences anticancer immunosurveillance and general health. Nat Rev Clin Oncol. 2018;15(6):382–96.
Fehrenbacher L, Spira A, Ballinger M, Kowanetz M, Vansteenkiste J, Mazieres J, et 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(10030):1837–46.
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(16):3857–66.
Coppola D, Nebozhyn M, Khalil F, Dai H, Yeatman T, Loboda A, et al. Unique ectopic lymph node-like structures present in human primary colorectal carcinoma are identified by immune gene array profiling. Am J Pathol. 2011;179(1):37–45.
Hwang S, Kwon AY, Jeong JY, Kim S, Kang H, Park J, et 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(1):643.
Ranganath H, Jain AL, Smith JR, Ryder J, Chaudry A, Miller E, et al. Association of a novel 27-gene immuno-oncology assay with efficacy of immune checkpoint inhibitors in advanced non-small cell lung cancer. BMC Cancer. 2022;22(1):407.
Jang HJ, Lee HS, Ramos D, Park IK, Kang CH, Burt BM, et al. Transcriptome-based molecular subtyping of non-small cell lung cancer may predict response to immune checkpoint inhibitors. J Thorac Cardiovasc Surg. 2020;159(4):1598-1610.e3.
Budczies J, Kirchner M, Kluck K, Kazdal D, Glade J, Allgäuer M, et al. A gene expression signature associated with B cells predicts benefit from immune checkpoint blockade in lung adenocarcinoma. Oncoimmunology. 2021;10(1):1860586.
Galon J, Angell HK, Bedognetti D, Marincola FM. The continuum of cancer immunosurveillance: prognostic, predictive, and mechanistic signatures. Immunity. 2013;39(1):11–26.
Bertucci F, Finetti P, Simeone I, Hendrickx W, Wang E, Marincola FM, et al. The immunologic constant of rejection classification refines the prognostic value of conventional prognostic signatures in breast cancer. Br J Cancer. 2018;119(11):1383–91.
Hendrickx W, Simeone I, Anjum S, Mokrab Y, Bertucci F, Finetti P, et al. Identification of genetic determinants of breast cancer immune phenotypes by integrative genome-scale analysis. OncoImmunology. 2017;6(2):e1253654.
Bertucci F, Niziers V, de Nonneville A, Finetti P, Mescam L, Mir O, et al. Immunologic constant of rejection signature is prognostic in soft-tissue sarcoma and refines the CINSARC signature. J Immunother Cancer. 2022;10(1):e003687.
Roelands J, Hendrickx W, Zoppoli G, Mall R, Saad M, Halliwill K, et al. Oncogenic states dictate the prognostic and predictive connotations of intratumoral immune response. J Immunother Cancer. 2020. https://doi.org/10.1136/jitc-2020-000617.
Prat A, Navarro A, Paré L, Reguart N, Galván P, Pascual T, et al. Immune-related gene expression profiling after PD-1 blockade in non-small cell lung carcinoma, head and neck squamous cell carcinoma, and melanoma. Cancer Res. 2017;77(13):3540–50.
Pender A, Titmuss E, Pleasance ED, Fan KY, Pearson H, Brown SD, et al. Genome and transcriptome biomarkers of response to immune checkpoint inhibitors in advanced solid tumors. Clin Cancer Res. 2021;27(1):202–12.
Cancer Genome Atlas Research Network. Comprehensive molecular profiling of lung adenocarcinoma. Nature. 2014;511(7511):543–50.
Cancer Genome Atlas Research Network. Comprehensive genomic characterization of squamous cell lung cancers. Nature. 2012;489(7417):519–25.
Wilkerson MD, Yin X, Hoadley KA, Liu Y, Hayward MC, Cabanski CR, et al. Lung squamous cell carcinoma mRNA expression subtypes are reproducible, clinically important, and correspond to normal cell types. Clin Cancer Res. 2010;16(19):4864–75.
Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.
Gatza ML, Lucas JE, Barry WT, Kim JW, Wang Q, Crawford MD, et al. A pathway-based classification of human breast cancer. Proc Natl Acad Sci USA. 2010;107(15):6994–9.
Rooney MS, Shukla SA, Wu CJ, Getz G, Hacohen N. Molecular and genetic properties of tumors associated with local immune cytolytic activity. Cell. 2015;160(1–2):48–61.
Thompson JC, Davis C, Deshpande C, Hwang WT, Jeffries S, Huang A, et al. Gene signature of antigen processing and presentation machinery predicts response to checkpoint blockade in non-small cell lung cancer (NSCLC) and melanoma. J Immunother Cancer. 2020;8(2):e000974.
Borcoman E, Kanjanapan Y, Champiat S, Kato S, Servois V, Kurzrock R, et al. Novel patterns of response under immunotherapy. Ann Oncol. 2019;30(3):385–96.
Harutani Y, Ozawa Y, Murakami E, Sato K, Oyanagi J, Akamatsu H, et al. Pre-treatment serum protein levels predict survival of non-small cell lung cancer patients without durable clinical benefit by PD-1/L1 inhibitors. Cancer Immunol Immunother. 2022;71(9):2109–16.
Zhu G, Ren D, Lei X, Shi R, Zhu S, Zhou N, et al. Mutations associated with no durable clinical benefit to immune checkpoint blockade in non-S-cell lung cancer. Cancers. 2021;13(6):1397.
Martini DJ, Hamieh L, McKay RR, Harshman LC, Brandao R, Norton CK, et al. durable clinical benefit in metastatic renal cell carcinoma patients who discontinue PD-1/PD-L1 therapy for immune-related adverse events. Cancer Immunol Res. 2018;6(4):402–8.
McShane LM, Altman DG, Sauerbrei W, Taube SE, Gion M, Clark GM, et al. REporting recommendations for tumour MARKer prognostic studies (REMARK). Br J Cancer. 2005;93(4):387–91.
Ettinger DS, Wood DE, Aisner DL, Akerley W, Bauman JR, Bharat A, et al. Non-small cell lung cancer, version 3.2022, NCCN clinical practice guidelines in oncology. J Natl Compr Cancer Netw. 2022;20(5):497–530.
Planchard D, Popat S, Kerr K, Novello S, Smit EF, Faivre-Finn C, et al. Metastatic non-small cell lung cancer: ESMO clinical practice guidelines for diagnosis, treatment and follow-up. Ann Oncol. 2018;29(4):192–237.
Rizvi H, Sanchez-Vega F, La K, Chatila W, Jonsson P, Halpenny D, et al. Molecular determinants of response to anti-programmed cell death (PD)-1 and anti-programmed death-ligand 1 (PD-L1) blockade in patients with non-small-cell lung cancer profiled with targeted next-generation sequencing. J Clin Oncol. 2018;36(7):633–41.
Roelands J, Kuppen PJK, Ahmed EI, Mall R, Masoodi T, Singh P, et al. An integrated tumor, immune and microbiome atlas of colon cancer. Nat Med. 2023;29(5):1273–86.
Zeng DQ, Yu YF, Ou QY, Li XY, Zhong RZ, Xie CM, Hu QG. Prognostic and predictive value of tumor-infiltrating lymphocytes for clinical therapeutic research in patients with non-small cell lung cancer. Oncotarget. 2016;7(12):13765–81.
Kursunel MA, Esendagli G. The untold story of IFN-γ in cancer biology. Cytokine Growth Factor Rev. 2016;31:73–81.
Zhao S, Shen W, Yu J, Wang L. TBX21 predicts prognosis of patients and drives cancer stem cell maintenance via the TBX21-IL-4 pathway in lung adenocarcinoma. Stem Cell Res Ther. 2018;9(1):89.
Li N, Li Y, Zheng P, Zhan X. Cancer stemness-based prognostic immune-related gene signatures in lung adenocarcinoma and lung squamous cell carcinoma. Front Endocrinol. 2021;12:755805.
Zhang Y, Liu Z. STAT1 in cancer: friend or foe? Discov Med. 2017;24(130):19–29.
Tokunaga R, Zhang W, Naseem M, Puccini A, Berger MD, Soni S, et al. CXCL9, CXCL10, CXCL11/CXCR3 axis for immune activation—a target for novel cancer therapy. Cancer Treat Rev. 2018;63:40–7.
Cheong JE, Sun L. Targeting the IDO1/TDO2-KYN-AhR pathway for cancer immunotherapy—challenges and opportunities. Trends Pharmacol Sci. 2018;39(3):307–25.
Chen B, Alvarado DM, Iticovici M, Kau NS, Park H, Parikh PJ, et al. Interferon-induced IDO1 mediates radiation resistance and is a therapeutic target in colorectal cancer. Cancer Immunol Res. 2020;8(4):451–64.
Bertucci F, De Nonneville A, Finetti P, Mamessier E. Predictive power of tertiary lymphoid structure signature for neoadjuvant chemotherapy response and immunotherapy benefit in HER2-negative breast cancer. Cancer Commun. 2023;43(8):943–6.
Forde PM, Spicer J, Lu S, Provencio M, Mitsudomi T, Awad MM, et al. Neoadjuvant nivolumab plus chemotherapy in resectable lung cancer. N Engl J Med. 2022;386(21):1973–85.
Felip E, Altorki N, Zhou C, Csőszi T, Vynnychenko I, Goloborodko O, et al. Adjuvant atezolizumab after adjuvant chemotherapy in resected stage IB-IIIA non-small-cell lung cancer (IMpower010): a randomised, multicentre, open-label, phase 3 trial. Lancet. 2021;398(10308):1344–57.
Takahara Y, Tanaka T, Ishige Y, Shionoya I, Yamamura K, Sakuma T, et al. Efficacy and predictors of rechallenge with immune checkpoint inhibitors in non-small cell lung cancer. Thorac Cancer. 2022;13(4):624–30.
The study was funded by the AORC Jeunes Chercheurs Grant. Our work was supported by Institut Paoli-Calmettes and Assistance Publique-Hôpitaux de Marseille.
Ethics approval and consent to participate
This study was approved by our institutional review board under the number 2019_93. All patients had given their signed written informed consent for the use of archived material for research purpose. The in silico analysis was based upon public data from published studies in which the informed patient's consent to participate and the ethics and institutional review board were already obtained by authors.
Consent for publication
A Mogenet received consulting fees from Takeda, Viatris and Pfizer and travel fees from BMS and Pfizer. The other authors have declared no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Correlation of PDL1 mRNA (CD274) and protein expression in lung carcinoma samples. A/ Correlation in the 329 cancer cell lines (grey) and in the 66 lung cancer cell lines (blue) in the DepMap database. B/ Correlation in the 234 lung adenocarcinoma clinical samples in the TCGA series. Figure S2. Analysis of homogeneity in a fixed-effects model revealed homogeneity between the five data sets in term of correlation between DCB rate and ICR1/2-4 classes (p=0.950).
: Table S1. List of NSCLC data sets analyzed. Table S2. ICR classification and correlation with clinico-pathological and biological variables in the 162 NSCLC samples. Table S3. ICR classification and correlations with biological variables in the 556 TCGA NSCLC samples. Table S4. Normalized gene expression data for the 44 NSCLC samples included in the study.
About this article
Cite this article
Mogenet, A., Finetti, P., Denicolai, E. et al. Immunologic constant of rejection as a predictive biomarker of immune checkpoint inhibitors efficacy in non-small cell lung cancer. J Transl Med 21, 637 (2023). https://doi.org/10.1186/s12967-023-04463-2
- Lung cancer
- ICR signature
- Immune therapy
- Immune checkpoints inhibitors