Next Article in Journal
The Current Evidence of Intensity-Modulated Radiotherapy for Hepatocellular Carcinoma: A Systematic Review and Meta-Analysis
Next Article in Special Issue
Pan-Cancer Profiling of Intron Retention and Its Clinical Significance in Diagnosis and Prognosis
Previous Article in Journal
Intranasal Delivery of Oncolytic Adenovirus XVir-N-31 via Optimized Shuttle Cells Significantly Extends Survival of Glioblastoma-Bearing Mice
Previous Article in Special Issue
Time-Dependent ROC Curve Analysis for Assessing the Capability of Radiation-Induced CD8 T-Lymphocyte Apoptosis to Predict Late Toxicities after Adjuvant Radiotherapy of Breast Cancer Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors †

1
Departments of Cutaneous Oncology and Immunology, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
2
Department of Machine Learning, H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
3
Department of Internal Medicine, Division of Medical Oncology, Comprehensive Cancer Center, The Ohio State University, Columbus, OH 43210, USA
4
Center for Public Health Genomics, School of Medicine, University of Virginia, Charlottesville, VA 22903, USA
5
Division of Surgical Oncology, Department of Surgery, School of Medicine, University of Colorado, Aurora, CO 80045, USA
6
USC Norris Comprehensive Cancer Center, Los Angeles, CA 90033, USA
7
Department of Neurosurgery, School of Medicine, University of Utah, Salt Lake City, UT 84132, USA
8
Huntsman Cancer Institute, Salt Lake City, UT 84132, USA
9
Oklahoma University Health Stephenson Cancer Center, Oklahoma City, OK 73104, USA
10
Department of Medicine, Roswell Park Comprehensive Cancer Center, Buffalo, NY 14263, USA
11
University of Kentucky Markey Cancer Center, Lexington, KY 40536, USA
12
Clinical & Life Sciences Department, Aster Insights, Hudson, FL 34667, USA
13
H. Lee Moffitt Cancer Center and Research Institute, Tampa, FL 33612, USA
14
Aster Insights, Hudson, FL 34667, USA
15
Department of Internal Medicine, Carver College of Medicine, University of Iowa Health Care, Iowa City, IA 52242, USA
*
Authors to whom correspondence should be addressed.
This work was presented at the Society for Immunotherapy of Cancer 37th Annual Meeting (SITC 2022), Boston, MA, USA, 8–12 November 2022; Trial Registration: NCT03977402.
These authors contributed equally to this work.
Cancers 2023, 15(20), 4913; https://doi.org/10.3390/cancers15204913
Submission received: 24 August 2023 / Revised: 14 September 2023 / Accepted: 29 September 2023 / Published: 10 October 2023

Abstract

:

Simple Summary

This study utilized real-world data from patients with advanced malignancies who underwent immune checkpoint inhibitor (ICI) treatment. We used transcriptomic data to derive an immunoscore based on CD3+ and CD8+ T cell densities using CIBERSORTx and the LM22 gene signature matrix. The imputed immunoscore effectively predicted overall survival (OS); patients with an intermediate–high immunoscore achieved better OS than those with a low immunoscore. Therefore, the T cell immunoscore represents a promising signature for estimating OS with ICIs and can be used as a reference for future machine learning-based biomarker development.

Abstract

Background: We aimed to determine the prognostic value of an immunoscore reflecting CD3+ and CD8+ T cell density estimated from real-world transcriptomic data of a patient cohort with advanced malignancies treated with immune checkpoint inhibitors (ICIs) in an effort to validate a reference for future machine learning-based biomarker development. Methods: Transcriptomic data was collected under the Total Cancer Care Protocol (NCT03977402) Avatar® project. The real-world immunoscore for each patient was calculated based on the estimated densities of tumor CD3+ and CD8+ T cells utilizing CIBERSORTx and the LM22 gene signature matrix. Then, the immunoscore association with overall survival (OS) was estimated using Cox regression and analyzed using Kaplan–Meier curves. The OS predictions were assessed using Harrell’s concordance index (C-index). The Youden index was used to identify the optimal cut-off point. Statistical significance was assessed using the log-rank test. Results: Our study encompassed 522 patients with four cancer types. The median duration to death was 10.5 months for the 275 participants who encountered an event. For the entire cohort, the results demonstrated that transcriptomics-based immunoscore could significantly predict patients at risk of death (p-value < 0.001). Notably, patients with an intermediate–high immunoscore achieved better OS than those with a low immunoscore. In subgroup analysis, the prediction of OS was significant for melanoma and head and neck cancer patients but did not reach significance in the non-small cell lung cancer or renal cell carcinoma cohorts. Conclusions: Calculating CD3+ and CD8+ T cell immunoscore using real-world transcriptomic data represents a promising signature for estimating OS with ICIs and can be used as a reference for future machine learning-based biomarker development.

1. Introduction

Cancer is one of the world’s most complicated biological systems [1]. Notably, there are numerous conventional treatment modalities, including surgery, radiotherapy, molecularly targeted therapy, and chemotherapy, that have been developed [2]. Because of these advances, many patients can live longer and/or with a better quality of life. Other patients, however, do not respond to these treatments, so additional strategies are necessary [3]. One alternative approach is to induce an immune-mediated antitumor response. In this context, the discovery of immune checkpoints like cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) and programmed death/ligand 1 (PD-1/PD-L1) has enabled the development of immune checkpoint inhibitors (ICIs), antibodies that target these checkpoints for the treatment of cancer [4]. Nevertheless, while ICIs can induce long-term remissions, even in patients with solid metastatic tumors, the majority of patients with other cancers do not achieve durable clinical responses [5,6,7]. One important reason for these low survival rates is the lack of prognostic biomarkers required to select patients for ICI monotherapies or combinations [8,9].
The development of ICIs clearly demonstrated that the immunogenicity of the tumor microenvironment (TME) plays an essential role in terms of the likelihood of response [10]. The TME can be infiltrated by various immune cells of potential immune-activating or suppressive effects with variable prognostic significance [11,12]. Indeed, several studies have indicated that helper T lymphocytes, cytotoxic T lymphocytes, and B lymphocytes play a crucial role as the prognostic markers of cancer patients [11,12,13,14,15,16,17,18]. The first evidence of tumors infiltrating lymphocytes (TILs) affecting survival was reported in 1921 [12]. Today, almost a century later, the prognostic value of TILs has become relevant to the staging of colon adenocarcinoma and with growing evidence to support prognostic value in melanoma [19].
The immunoscore, as validated by Galon et al. [20], estimates CD3+ and CD8+ T cell density in the core of the tumor (CT) and the invasive margin (IM). Evidence supports that a higher density of TILs (i.e., a high immunoscore) is associated with a better prognosis [11,21,22,23,24,25,26,27,28,29,30]. Hence, it is essential to note that immunoscore may reflect the immunogenicity of the TME and the tumor’s susceptibility to immunotherapy [31,32]. In fact, immunoscore has been widely evaluated in colorectal cancer patients to determine its effects on both overall survival (OS) and disease-free survival (DFS). For instance, a meta-analysis published recently found that a low immunoscore significantly correlated with poor OS (HR = 1.74, 95% CI: 1.43–2.13) and DFS (HR = 1.82, 95% CI: 1.64–2.03) [33]. Whereas a search of the literature revealed that few studies have investigated the role of immunoscore as a predictor for outcomes in other types of cancers, including esophageal cancer [34], bladder cancer [35], and non-small cell lung cancer (NSCLC), or in a cross-cancer capacity as a pan-cancer prognostic signature [36]. Furthermore, research into the intricate networks of cell–cell interactions in the TME and how this may affect several types of cancer responses to immunotherapies, such as ICIs, is a rapidly expanding area of research [36]. It relies on a limited repertoire of phenotypic markers and may result in cell loss or damage in immunohistochemistry or fluorescence-activated cell sorting [37]. To overcome these barriers and to take advantage of RNA sequencing data, system biology methods like transcriptome deconvolution were developed to estimate the relative densities of different cell types [36]. In this regard, evaluating TIL status in the tumor will be extremely helpful in investigating tumor immune cell states and susceptibility to immunotherapy. For this purpose, the specific objective of this study is to determine the prognostic value of immunoscore based on real-world transcriptomics data from TILs in patients with advanced malignancies treated with ICIs in an effort to explore its utility as a reference or control for machine learning-based biomarker development.

2. Methods and Materials

2.1. Patients and Datasets

For this study, real-world clinical and transcriptomic data retrospectively collected under the Total Cancer Care Protocol (NCT03977402) and Avatar® project within the Oncology Research Information Exchange Network (ORIEN) of 18 collaborating cancer centers was utilized. An IRB-approved informed consent was obtained from all subjects at their participating institutions. The patients had to be 18 years old and above with cancer and have been treated with ICI. As shown in Figure 1, the construction process of the real-world immunoscore calculation is divided into different steps.
Briefly, this study involves the collection of normal and tumor tissue, blood and/or fluid samples with DNA in the form of frozen and/or formalin-fixed, paraffin-embedded (FFPE) tissue that was obtained prior to ICI treatment initiation from consenting subjects. Also, an additional tumor sample was obtained at the time of planned diagnostic biopsies or from previously collected or stored tumor tissue (if available). All patients’ related data, such as survey data, medical records data, cancer registry data, and other related data, were collected from the time they first joined the study until the time of conducting the current study.

2.2. RNA-Sequencing and Data Processing

RNA sequencing and data processing were completed as previously published https://www.asterinsights.com/white-paper/renal-cell-carcinoma-rwd-data/ (accessed on 28 September 2023). Subsequently, the RNA expression profiles were identified in the ORIEN database by downloading a series of matrix files that contain transcript per million (TPM) at the gene level. To avoid the value of zero during the normalization of TPM, a value of 1 was added to TPM, then Log10 transformation was used to normalize data value (TPM to (log2(TPM+1)); finally, the results were exponentiated to be a linear scale. Afterward, using linear regression, CIBERSORTx selected genes from the input matrix based on the LM22 signature matrix to deconvolve a given mixture. The input matrix of reference gene expression signatures was made using the standard annotation file. Finally, the CIBERSORTx algorithm runs in Python with 20 five-fold cross-validation permutations using the LM22 signature. The cut-off for statistical significance was set at p < 0.05. The input matrix of reference gene expression signatures was created using the standard annotation file.

2.3. Immunoscore Imputation

A limitation of our real-world dataset was that the direct calculation of the immunoscore was impractical, as our data did not include direct measurements of the CD3+ and CD8+ densities. However, it was shown that the gene expression levels can be used to impute cell type abundance. For this purpose, CIBERSORTx was employed using the leukocyte signature matrix (LM22), a gene signature matrix to impute the abundance of member cell types in the mixed cell population using gene expression data. Mainly, CIBERSORTx is an analytical tool that imputes gene expression profiles based on an input gene signature matrix. LM22 is a signature matrix containing 547 genes capable of accurately distinguishing 22 mature populations of human hematopoietic cells. This includes seven T cell types: naive and memory B, plasma, natural killer (NK), and myeloid subsets. Even though LM22 was made and tested with data from gene expression microarrays, it can also be used with data from RNA sequences [38]. Upon the imputation of CD3 and CD8 densities, patients’ immunoscore was calculated as categorical variables with two states of “low” and “intermediate and high”, following the procedure previously proposed by Galon et al. [1].

2.4. Construction of the Immunoscore

Following the percentile cut-off that is presented in Figure 2, the immunoscore was calculated as the average of the percentile of CD3CT and CD8CT densities of the patients among the training population. In order to avoid over-fitting, the immunoscore was computed using 20 iterations of five-fold validation such that in each iteration, four folds are employed to define the percentile functions, which are used to calculate the immunoscore of the patients in the remaining fold. Instantly, using this approach, the 20th percentile was calculated for each patient, with the ultimate immunoscore being computed based on their average.
Regarding the immunoscore categories, we considered two different cut-off points. The first cut-off point was the 25th percentile, as proposed by Galon et al. [1], and also shown in Figure S1 in the Supplementary Materials. For the second cut-off point, we employed Youden’s J statistic to define the optimal cut-off point for our data. To do so, for each iteration of our cross-validations, a cut-off point was computed using Youden’s J statistic, with the final cut-off point being the average.

2.5. Statistical Analysis

Kaplan–Meier curves with 95% confidence interval (CI) survival plots as shaded areas were used to visualize the survival curves of patients with different immunoscores. In the survival analysis, the starting time point was the first exposure to ICI drugs, while death was considered an event and alive or lost follow-up was considered right-censored. A log-rank test was also employed to assess the significance of the difference in the survival distributions of the identified immunoscore categories. The analysis was performed first considering our entire dataset and second reflecting each cancer category individually. Statistical analyses were conducted using Python 3.8 libraries, Delaware, United States.

3. Results

3.1. Baseline Characteristics

Based on our methodology, data from 522 patients were included in our analysis. Of those patients, just over half (275 (53.0%)) experienced a study event, with 10.5 months as the median time to death. It is apparent in Table 1 that there is an imperceptible difference in the distribution of the four distinct malignant tumor types among our cohort, where each type represents just above or below a quarter of the total cohort, namely renal cell carcinoma (RCC), representing 28.5% of the overall cohort, followed by NSCLC (24.5%) and melanoma (23.9%), then 23% for head and neck cancer. In addition, data from this table indicates that the majority of the patients (80%) were treated with either nivolumab or pembrolizumab (42.0% and 38%, respectively), while a minority of patients (13.2%) were treated with the combination of ipilimumab and nivolumab. Moreover, only a smaller minority of participants (6.0%) received either ipilimumab, avelumab or cemiplimab.

3.2. Entire Cohort

The first set of investigations examined the immunoscore analysis on our entire cohort using the 25th percentile cut-off. This identified 122 (23.4%) patients with a low immunoscore. However, a superior result was achieved when employing Youden’s J statistic, which identified the 43.5th percentile as the optimal cut-off point. This resulted in approximately doubling the number of patients (222 (42.5%)) with a low immunoscore.
Figure 2A,B compare the Kaplan–Meier curves of the stratified patients based on the 25th and 43.5th percentile cut-off points, correspondingly. Both figures illustrate that those patients with an intermediate–high immunoscore (immunoscore = 1) achieved a significantly better survival time than patients with a low immunoscore (immunoscore = 0) (p < 0.001 for both cut-off points); however, based on Figure 2B, the 43.5th percentile resulted in a visually better stratification. Table 2 summarizes Harrel’s average c-index and the associated 95% confidence interval (CI).

3.3. Cancer Categories

Next, we considered patients with the same histology as independent cohorts. We performed immunoscore sub-group analyses on the four cancer categories included in our dataset. The results obtained from the survival analysis and displayed in Kaplan–Meier curves of four different cancers histologies stratified using a 25th percentile cut-off point are presented in Figure 3. Remarkably, in Figure 3A,B, it is clear that patients with melanoma and head and neck tumors with an intermediate–high immunoscore (immunoscore = 1) had a significantly superior OS time than patients with a low immunoscore (immunoscore = 0) (p = 0.009 and p = 0.04, respectively). On the contrary, as shown in Figure 3C,D for OS time between an intermediate–high immunoscore (immunoscore = 1) and a low immunoscore (immunoscore = 0) in patients with NSCLC and RCC (p = 0.77 and p = 0.17, correspondingly), while we observe survival curve separation, the differences were not statistically significant.
Using Youden’s J statistic, different cut-off points were identified for the four cancer categories. Specifically, the optimal cut-off point of patients with melanoma was the 49.13th percentile; for head and neck cancer, it was the 43.7th; for NSCLC, it was the 57.75th; and for RCC, it was the 36.9th. Figure 4 presents the survival analyses of the different cancer histologies stratified based on the optimized cut-off points. Looking at Figure 4A,B, the statistically improved separation between the two immunoscore score categories for patients with melanoma and head and neck with log-rank test p-values of 0.01 and <0.001, respectively, are apparent. This could be interpreted as patients with melanoma or head and neck cancer, who achieved an intermediate–high immunoscore of 1, had a higher OS time than patients in the same cancer category who had a low immunoscore of 0. In contradiction to the above results, which indicated a significant separation of patients according to their immunoscore, Figure 4C,D demonstrated that adjusting the cut-off point of immunoscore in patients with NSCLC and RCC cancer did not result in a significant improvement in OS time between patients with an intermediate–high immunoscore versus patients with a low immunoscore (p = 0.37 and p = 0.25, respectively). Table 3 shows Harrel’s average C-index and the associated 95% CI together, and these results provide crucial insights into the role of immunoscore in predicting OS in various cancer types.

4. Discussion

Although patients with select solid metastatic tumors have been shown to achieve long-term remissions with ICIs, most tumor types do not experience clinical benefits [5,6,7]. These low response rates may be partly attributed to the lack of prognostic biomarkers for identifying patients who have the capacity to benefit from therapy with ICIs [8,9]. The T cell immunoscore has been investigated as a prognostic biomarker that may also inform the status of an immune response to a tumor [30,31]. Due to the multiple parameters required for immunophenotyping, immunohistochemistry is limited in the number of immune groups that can be highlighted [39]. Also, flow cytometry measures a limited number of markers and requires strict technical methods, and therefore, closely related cell types may be missed [40].
In biomarker development with regard to ICI therapy, there is a need for a relatively simple biomarker signature of reasonable prognostic value that can be used as a reference, particularly in the rapidly advancing field of artificial intelligence and machine learning. This becomes more important when utilizing real-world data that are becoming increasingly available to researchers interested in ICI biomarker research. Clearly, bulk tissue gene expression profiling does not rely on surface markers or have cellular dissociation errors. Although immune-enriched gene expression signatures have prognostic value, connecting these signatures to specific TIL characteristics can be challenging [41]. Mathematically splitting bulk tumor gene expression patterns into individual cell types can inform this problem [20]. In addition to this, in reviewing the literature, little to no real-world data were found concerning the analysis of the association between immunoscore and ICIs outcomes in variant types of tumors, applying RNA expression data to calculate CD3+ and CD8+. Therefore, we attempted to validate the value of the CD3+ and CD8+ T cell-imputed immunoscore as a prognostic biomarker in cancer patients treated with ICIs and to support its utility as a “reference” or “control” in machine learning (ML)/artificial intelligence-based studies for biomarker development in the context of ICI and other immunotherapies when using transcriptomic data. So, a unique method was employed to accomplish this objective, utilizing transcriptomic data and CIBERSORTx, a computational approach that implements ML to determine imputed immunoscore and suggests immunoscore cut-off percentiles for different histologic cancer types [42].
One interesting finding is that immunoscore based on transcriptomics was able to significantly distinguish between the OS among patients with different cancer types as a total cohort using the cut-off point estimated by Galon et al. [1] or Youden’s J statistic, with better performance from the latter. The results showed that patients with intermediate–high immunoscore had better OS than patients with low immunoscore. A possible explanation for this might be that patients with a higher immunoscore have a higher density of TIL known to have antitumor activity [43]. Stratification was visually improved when Youden’s J statistic was applied. This result could cast a new light on the validity of our method of using transcriptomic data for analyzing immunoscore based on gene expression data from tumor tissue; CIBERSORTx might be a useful tool for determining the cellular composition of TIL. This includes both innate and adaptive compartments [42].
It is worth discussing the interesting facts revealed by the analysis results of each malignancy in our cohort. The analysis of immunoscore in patients with melanoma revealed a significant difference between patients with low versus intermediate–high immunoscore in both cut-offs: the 25th percentile and the 49.13th percentile cut-off point was estimated using Youden’s J statistic. This also aligns with previous reports that showed a significant association between higher CD8+, CD45+, and CD3+ cell counts and better OS (p = 0.001, p = 0.004, and p = 0.009, respectively) [44]. Even with these studies, it is becoming harder and harder to define immunoscore in melanoma based on pathological tumor-node-metastasis (TNM) staging because it is built on complex immune reactions intratumorally. Alternatively, metastatic lymph node tissue could be the best source for evaluating the immunoscore because, in many cases, the metastatic lymph nodes from a lymphectomy are the only available tissue [15]. A similar pattern of results was obtained by examining OS in patients with head and neck cancer based on their immunoscore with the two cut-off points used. The observed result is in agreement with the outcomes reported by Furgiuele et al. [45], who evaluated the prognostic role of TIL in developing an immunoscore in head and neck cancer patients. According to that study, the CD8+ density was an independent prognostic marker for recurrence-free survival (RFS) and OS. Furthermore, patients with high CD8+, CD68+, and FoxP3 T cell density had better OS.
This study has not demonstrated a significant correlation between the generated T cell immunoscore and OS in patients with NSCLC or RCC when both cut-off percentiles were applied. The lack of significance in patients with NSCLC is contrary to that of X.-T. Li et al. (2021) [36], who found that patients with NSCLC could benefit from an immunoscore as a powerful, independent, and significant prognostic indicator. This could be related to the sample sizes in the different cohorts. Similarly, while our results in patients with RCC did not reach statistical significance, other studies support a significant correlation [46]. Researchers found that patients with a high immunoscore had prolonged disease-free survival (DFS), progression-free survival (PFS), and OS (HR 2.652, 2.848, and 2.933, respectively; all p < 0.001) compared to those with a lower score. These results were consistent in the sub-group analysis in patients with different Fuhrman grades and pathological TNM stages.
A strength of the current research is that it is based on real-world data and investigated the prognostic value of immunoscore in multiple types of malignancies from transcriptomic data using ML for the first time. However, this study is not without limitations. Importantly, as RNA sequencing was performed on the specimens collected from the tumor, the immunoscore reflects the immune infiltration in a part of the tumor that could be central or peripheral, depending on the tumor biopsy obtained. This is unlike the original proposition of Galon, where the CD3+ and CD8+ T cell densities were best estimated via immunohistochemistry within the invasive margin [33,34,35,36]. In addition, there is no independent validation of the sample in this study because of its retrospective design. Future studies with larger cohorts are required to confirm the findings obtained in this study.

5. Conclusions

In conclusion, the present research aimed to examine the prognostic value of immunoscore using real-world clinical and transcriptomics data in patients with various advanced malignancies treated with ICIs. The findings clearly indicate that the T cell immunoscore based on real-world transcriptomic data is able to predict patients at risk of death at a significant level and can be utilized as a reference for developing biomarkers based on machine learning.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15204913/s1, Figure S1. Shows immunoscore calculation using two categories of “low” and “intermediate–high” using the imputed densities of CD3+ and CD8+ T cells.

Author Contributions

Conceptualization, A.A.T., I.M.E.N., P.G.S. and I.E.; data curation, all authors; formal analysis, A.A.T., I.M.E.N., P.G.S. and I.E.; investigation, all authors; methodology, all authors; project administration, all authors; resources, all authors; financial support, A.A.T. and I.M.E.N.; supervision, A.A.T. and I.M.E.N.; writing—original draft, all authors; writing—review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Funding

The original clinical trial within which this project was nested was supported by an ORIEN FOUNDATION NOVA (New Oncologic Visionary Award) Grant and Community Foundation of Tampa Bay: 69-21295-01-01. In addition, it is partly supported by the National Institute of Health (NIH) grant R01-CA233487 and its supplement (CA233487-05S1).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of the Total Cancer Care Protocol (IRB# Pro00014441).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the original clinical trial within which this study was nested.

Data Availability Statement

The data sets used and analyzed for the current study are available from the corresponding author on request.

Conflicts of Interest

I.E., P.G.S., M.L.C., J.C., James Chen., A.P.I., W.S.D., G.J.W., J.C., John Carpten and P.H. declare no conflicts of interest. M.M. declares contracted research grants with Merck, Taiho, and the National Comprehensive Cancer Network. I.P. declares research funding: NIH/NCI P30CA016056 Cancer Center Support Grant Program Leader; 1.2 calendar months consulting/advisory role: Nouscom, Iovance, Nektar, and Regeneron; and stock ownership: Ideaya, Inc. H.C. declares advisory board/consultant: Best Doctors/Teladoc, Orbus Therapeutics, Bristol Meyers Squibb, Regeneron, Novocure, and PPD/Chimerix; research funding (site PI/institutional contract): Orbus, GCAR, Array BioPharma, Karyopharm Therapeutics, Nuvation Bio, Bayer, Bristol Meyer Squib, and Sumitomo Dainippon Pharma. To date, S.M.A. is a member of ASCO COI. I.M.E.N. is on the scientific advisory board of Endectra, L.L., acts as deputy editor for the Journal of Medical Physics, and receives funding from NIH and DoD. Martin McCarter declares receiving research support from Merck, Taiho, and NCCN. A.A.T. declares contracted research grants with the institution from Bristol Myers Squib, Genentech-Roche, Regeneron, Sanofi-Genzyme, Nektar, Clinigen, Merck, Acrotech, Pfizer, Checkmate, OncoSec, Scholar Rock, InflaRx GmbH, and Agenus; personal consultant/advisory board fees (less than USD 10k per year) from Bristol Myers Squibb, Merck, Easai, Instil Bio, Clinigin, Regeneron, Sanofi-Genzyme, Novartis, Partner Therapeutics, Genentech/Roche, BioNTech, Concert AI, and AstraZeneca outside the submitted work.

References

  1. Galon, J.; Lanzi, A. Immunoscore and its introduction in clinical practice. Q. J. Nucl. Med. Mol. Imaging 2020, 64, 152–161. [Google Scholar] [CrossRef]
  2. Yadav, D.; Kwak, M.; Chauhan, P.S.; Puranik, N.; Lee, P.C.W.; Jin, J.O. Cancer immunotherapy by immune checkpoint blockade and its advanced application using bio-nanomaterials. Semin. Cancer Biol. 2022, 86 Pt 2, 909–922. [Google Scholar] [CrossRef] [PubMed]
  3. Vale, C.L. Neoadjuvant chemotherapy in invasive bladder cancer: Update of a systematic review and meta-analysis of individual patient data advanced bladder cancer (ABC) meta-analysis collaboration. Eur. Urol. 2005, 48, 202–205; discussion 5-6. [Google Scholar] [CrossRef]
  4. Bondhopadhyay, B.; Sisodiya, S.; Chikara, A.; Khan, A.; Tanwar, P.; Afroze, D.; Singh, N.; Agrawal, U.; Mehrotra, R.; Hussain, S. Cancer immunotherapy: A promising dawn in cancer research. Am. J. Blood Res. 2020, 10, 375–385. [Google Scholar] [PubMed]
  5. Fares, C.M.; Van Allen, E.M.; Drake, C.G.; Allison, J.P.; Hu-Lieskovan, S. Mechanisms of Resistance to Immune Checkpoint Blockade: Why Does Checkpoint Inhibitor Immunotherapy Not Work for All Patients? Am. Soc. Clin. Oncol. Educ. Book 2019, 39, 147–164. [Google Scholar] [CrossRef] [PubMed]
  6. Haslam, A.; Prasad, V. Estimation of the Percentage of US Patients With Cancer Who Are Eligible for and Respond to Checkpoint Inhibitor Immunotherapy Drugs. JAMA Netw. Open 2019, 2, e192535. [Google Scholar] [CrossRef]
  7. Kalbasi, A.; Ribas, A. Tumour-intrinsic resistance to immune checkpoint blockade. Nat. Rev. Immunol. 2020, 20, 25–39. [Google Scholar] [CrossRef] [PubMed]
  8. Reck, M.; Rodriguez-Abreu, D.; Robinson, A.G.; Hui, R.; Csoszi, T.; Fulop, A.; Gottfried, M.; Peled, N.; Tafreshi, A.; Cuffe, S.; et al. Pembrolizumab versus Chemotherapy for PD-L1-Positive Non-Small-Cell Lung Cancer. N. Engl. J. Med. 2016, 375, 1823–1833. [Google Scholar] [CrossRef] [PubMed]
  9. Horn, L.; Spigel, D.R.; Vokes, E.E.; Holgado, E.; Ready, N.; Steins, M.; Poddubskaya, E.; Borghaei, H.; Felip, E.; Paz-Ares, L.; et al. Nivolumab Versus Docetaxel in Previously Treated Patients With Advanced Non-Small-Cell Lung Cancer: Two-Year Outcomes From Two Randomized, Open-Label, Phase III Trials (CheckMate 017 and CheckMate 057). J. Clin. Oncol. 2017, 35, 3924–3933. [Google Scholar] [CrossRef]
  10. Hinshaw, D.C.; Shevde, L.A. The Tumor Microenvironment Innately Modulates Cancer Progression. Cancer Res. 2019, 79, 4557–4566. [Google Scholar] [CrossRef] [PubMed]
  11. Bruni, D.; Angell, H.K.; Galon, J. The immune contexture and Immunoscore in cancer prognosis and therapeutic efficacy. Nat. Rev. Cancer 2020, 20, 662–680. [Google Scholar] [CrossRef] [PubMed]
  12. Galon, J.; Bruni, D. Tumor Immunology and Tumor Evolution: Intertwined Histories. Immunity 2020, 52, 55–81. [Google Scholar] [CrossRef] [PubMed]
  13. Ascierto, P.A.; Capone, M.; Urba, W.J.; Bifulco, C.B.; Botti, G.; Lugli, A.; Marincola, F.M.; Ciliberto, G.; Galon, J.; Fox, B.A. The additional facet of immunoscore: Immunoprofiling as a possible predictive tool for cancer treatment. J. Transl. Med. 2013, 11, 54. [Google Scholar] [CrossRef] [PubMed]
  14. Galon, J.; Bruni, D. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies. Nat. Rev. Drug Discov. 2019, 18, 197–218. [Google Scholar] [CrossRef]
  15. Galon, J.; Fox, B.A.; Bifulco, C.B.; Masucci, G.; Rau, T.; Botti, G.; Marincola, F.M.; Ciliberto, G.; Pages, F.; Ascierto, P.A.; et al. Immunoscore and Immunoprofiling in cancer: An update from the melanoma and immunotherapy bridge 2015. J. Transl. Med. 2016, 14, 273. [Google Scholar] [CrossRef]
  16. Kirilovsky, A.; Marliot, F.; El Sissy, C.; Haicheur, N.; Galon, J.; Pages, F. Rational bases for the use of the Immunoscore in routine clinical settings as a prognostic and predictive biomarker in cancer patients. Int. Immunol. 2016, 28, 373–382. [Google Scholar] [CrossRef]
  17. Pages, F.; Galon, J.; Fridman, W.H. The essential role of the in situ immune reaction in human colorectal cancer. J. Leukoc. Biol. 2008, 84, 981–987. [Google Scholar] [CrossRef]
  18. Fridman, W.H.; Dieu-Nosjean, M.C.; Pages, F.; Cremer, I.; Damotte, D.; Sautes-Fridman, C.; Galon, J. The immune microenvironment of human tumors: General significance and clinical impact. Cancer Microenviron. 2013, 6, 117–122. [Google Scholar] [CrossRef]
  19. Marliot, F.; Pages, F.; Galon, J. Usefulness and robustness of Immunoscore for personalized management of cancer patients. Oncoimmunology 2020, 9, 1832324. [Google Scholar] [CrossRef]
  20. Galon, J.; Costes, A.; Sanchez-Cabo, F.; Kirilovsky, A.; Mlecnik, B.; Lagorce-Pages, C.; Tosolini, M.; Camus, M.; Berger, A.; Wind, P.; et al. Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science 2006, 313, 1960–1964. [Google Scholar] [CrossRef]
  21. Galon, J.; Pages, F.; Marincola, F.M.; Thurin, M.; Trinchieri, G.; Fox, B.A.; Gajewski, T.F.; Ascierto, P.A. The immune score as a new possible approach for the classification of cancer. J. Transl. Med. 2012, 10, 1. [Google Scholar] [CrossRef]
  22. Anitei, M.G.; Zeitoun, G.; Mlecnik, B.; Marliot, F.; Haicheur, N.; Todosi, A.M.; Kirilovsky, A.; Lagorce, C.; Bindea, G.; Ferariu, D.; et al. Prognostic and predictive values of the immunoscore in patients with rectal cancer. Clin. Cancer Res. 2014, 20, 1891–1899. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, Z.Y.; Shen, W.; Yue, J.Q.; Yao, W.Y.; Liu, S.L.; Jin, Y.P.; Dong, P.; Ma, F.; Wu, X.S.; Gong, W. Combining Immunoscore with Clinicopathologic Features in Cholangiocarcinoma: An Influential Prognostic Nomogram. Onco Targets Ther. 2020, 13, 11359–11376. [Google Scholar] [CrossRef] [PubMed]
  24. Galon, J.; Mlecnik, B.; Bindea, G.; Angell, H.K.; Berger, A.; Lagorce, C.; Lugli, A.; Zlobec, I.; Hartmann, A.; Bifulco, C.; et al. Towards the introduction of the ‘Immunoscore’ in the classification of malignant tumours. J. Pathol. 2014, 232, 199–209. [Google Scholar] [CrossRef] [PubMed]
  25. Pages, F.; Mlecnik, B.; Marliot, F.; Bindea, G.; Ou, F.S.; Bifulco, C.; Lugli, A.; Zlobec, I.; Rau, T.T.; Berger, M.D.; et al. International validation of the consensus Immunoscore for the classification of colon cancer: A prognostic and accuracy study. Lancet 2018, 391, 2128–2139. [Google Scholar] [CrossRef]
  26. Jiang, Y.; Zhang, Q.; Hu, Y.; Li, T.; Yu, J.; Zhao, L.; Ye, G.; Deng, H.; Mou, T.; Cai, S.; et al. ImmunoScore Signature: A Prognostic and Predictive Tool in Gastric Cancer. Ann. Surg. 2018, 267, 504–513. [Google Scholar] [CrossRef] [PubMed]
  27. Li, X.D.; Huang, C.W.; Liu, Z.F.; Jiang, L.J.; Chen, J.W.; Xie, D.; Zhou, F.J.; Lu, H.M.; Liu, Z.W. Prognostic Role of the Immunoscore for Patients with Urothelial Carcinoma of the Bladder Who Underwent Radical Cystectomy. Ann. Surg. Oncol. 2019, 26, 4148–4156. [Google Scholar] [CrossRef] [PubMed]
  28. Angell, H.K.; Bruni, D.; Barrett, J.C.; Herbst, R.; Galon, J. The Immunoscore: Colon Cancer and Beyond. Clin. Cancer Res. 2020, 26, 332–339. [Google Scholar] [CrossRef]
  29. Zhang, X.; Yang, J.; Du, L.; Zhou, Y.; Li, K. The prognostic value of Immunoscore in patients with cancer: A pooled analysis of 10,328 patients. Int. J. Biol. Markers 2020, 35, 3–13. [Google Scholar] [CrossRef]
  30. Nie, R.C.; Yuan, S.Q.; Wang, Y.; Chen, Y.B.; Cai, Y.Y.; Chen, S.; Li, S.M.; Zhou, J.; Chen, G.M.; Luo, T.Q.; et al. Robust immunoscore model to predict the response to anti-PD1 therapy in melanoma. Aging 2019, 11, 11576–11590. [Google Scholar] [CrossRef]
  31. Stoll, G.; Enot, D.; Mlecnik, B.; Galon, J.; Zitvogel, L.; Kroemer, G. Immune-related gene signatures predict the outcome of neoadjuvant chemotherapy. Oncoimmunology 2014, 3, e27884. [Google Scholar] [CrossRef] [PubMed]
  32. Vacchelli, E.; Galluzzi, L.; Fridman, W.H.; Galon, J.; Sautes-Fridman, C.; Tartour, E.; Kroemer, G. Trial watch: Chemotherapy with immunogenic cell death inducers. Oncoimmunology 2012, 1, 179–188. [Google Scholar] [CrossRef] [PubMed]
  33. Sun, G.; Dong, X.; Tang, X.; Qu, H.; Zhang, H.; Zhao, E. The prognostic value of immunoscore in patients with colorectal cancer: A systematic review and meta-analysis. Cancer Med. 2019, 8, 182–189. [Google Scholar] [CrossRef] [PubMed]
  34. Noma, T.; Makino, T.; Ohshima, K.; Sugimura, K.; Miyata, H.; Honma, K.; Yamashita, K.; Saito, T.; Tanaka, K.; Yamamoto, K.; et al. Immunoscore Signatures in Surgical Specimens and Tumor-Infiltrating Lymphocytes in Pretreatment Biopsy Predict Treatment Efficacy and Survival in Esophageal Cancer. Ann. Surg. 2023, 277, e528–e537. [Google Scholar] [CrossRef]
  35. Nassif, E.F.; Thibault, C.; Oudard, S.; Galon, J. Precision immunity: Immunoscore and neoadjuvant treatment in bladder cancer. Oncoimmunology 2021, 10, 1888488. [Google Scholar] [CrossRef]
  36. Li, X.T.; Zhang, J.T.; Yan, H.H.; Su, J.; Cheng, M.L.; Sun, Q.H.; Zhong, W.Z.; Wu, Y.L.; Zhang, D.X.; Hou, D.J. Gene co-expression modules integrated with immunoscore predicts survival of non-small cell lung cancer. Cancer Treat. Res. Commun. 2021, 26, 100297. [Google Scholar] [CrossRef]
  37. Shen-Orr, S.S.; Gaujoux, R. Computational deconvolution: Extracting cell type-specific information from heterogeneous samples. Curr. Opin. Immunol. 2013, 25, 571–578. [Google Scholar] [CrossRef]
  38. Chen, B.; Khodadoust, M.S.; Liu, C.L.; Newman, A.M.; Alizadeh, A.A. Profiling Tumor Infiltrating Immune Cells with CIBERSORT. Methods Mol. Biol. 2018, 1711, 243–259. [Google Scholar]
  39. Tung, J.W.; Heydari, K.; Tirouvanziam, R.; Sahaf, B.; Parks, D.R.; Herzenberg, L.A.; Herzenberg, L.A. Modern flow cytometry: A practical approach. Clin. Lab. Med. 2007, 27, 453–468. [Google Scholar] [CrossRef]
  40. Abbas, A.R.; Wolslegel, K.; Seshasayee, D.; Modrusan, Z.; Clark, H.F. Deconvolution of blood microarray data identifies cellular activation patterns in systemic lupus erythematosus. PLoS ONE. 2009, 4, e6098. [Google Scholar] [CrossRef]
  41. Curtis, C.; Shah, S.P.; Chin, S.F.; Turashvili, G.; Rueda, O.M.; Dunning, M.J.; Speed, D.; Lynch, A.G.; Samarajiwa, S.; Yuan, Y.; et al. The genomic and transcriptomic architecture of 2000 breast tumours reveals novel subgroups. Nature 2012, 486, 346–352. [Google Scholar] [CrossRef] [PubMed]
  42. Newman, A.M.; Liu, C.L.; Green, M.R.; Gentles, A.J.; Feng, W.; Xu, Y.; Hoang, C.D.; Diehn, M.; Alizadeh, A.A. Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 2015, 12, 453–457. [Google Scholar] [CrossRef] [PubMed]
  43. Weeks, J.C.; Cook, E.F.; O’Day, S.J.; Peterson, L.M.; Wenger, N.; Reding, D.; Harrell, F.E.; Kussin, P.; Dawson, N.V.; Connors, A.F., Jr.; et al. Relationship between cancer patients’ predictions of prognosis and their treatment preferences. JAMA 1998, 279, 1709–1714. [Google Scholar] [CrossRef] [PubMed]
  44. Erdag, G.; Schaefer, J.T.; Smolkin, M.E.; Deacon, D.H.; Shea, S.M.; Dengel, L.T.; Patterson, J.W.; Slingluff, C.L., Jr. Immunotype and immunohistologic characteristics of tumor-infiltrating immune cells are associated with clinical outcome in metastatic melanoma. Cancer Res. 2012, 72, 1070–1080. [Google Scholar] [CrossRef] [PubMed]
  45. Furgiuele, S.; Descamps, G.; Lechien, J.R.; Dequanter, D.; Journe, F.; Saussez, S. Immunoscore Combining CD8, FoxP3, and CD68-Positive Cells Density and Distribution Predicts the Prognosis of Head and Neck Cancer Patients. Cells 2022, 11, 2050. [Google Scholar] [CrossRef] [PubMed]
  46. Selvi, I.; Demirci, U.; Bozdogan, N.; Basar, H. The prognostic effect of immunoscore in patients with clear cell renal cell carcinoma: Preliminary results. Int. Urol. Nephrol. 2020, 52, 21–34. [Google Scholar] [CrossRef]
Figure 1. The construction process of the real-world immunoscore calculation. This figure demonstrates the different steps for calculating immunoscore based on real-world transcription data using CIBERSORTx as a learning tool.
Figure 1. The construction process of the real-world immunoscore calculation. This figure demonstrates the different steps for calculating immunoscore based on real-world transcription data using CIBERSORTx as a learning tool.
Cancers 15 04913 g001
Figure 2. Kaplan–Meier plot of the entire cohort (n = 522) stratified based on (A) 25th percentile cut-off. (B) Kaplan–Meier plot of the entire cohort stratified based on (A) 43.5th percentile cut-off. p-value < 0.05 is considered statically significant.
Figure 2. Kaplan–Meier plot of the entire cohort (n = 522) stratified based on (A) 25th percentile cut-off. (B) Kaplan–Meier plot of the entire cohort stratified based on (A) 43.5th percentile cut-off. p-value < 0.05 is considered statically significant.
Cancers 15 04913 g002
Figure 3. Kaplan–Meier plot of (A) patients with melanoma (n = 125) stratified based on the 25th percentile cut-off. (B) Patients with head and neck cancer (n = 120) stratified based on 25th percentile cut-offs. (C) Patients with non-small cell lung cancer (n = 128) stratified based on the 25th percentile cut-off. (D) Patients with renal cell carcinoma (n = 149) stratified based on the 25th percentile cut-off. p-value < 0.05 is considered statically significant.
Figure 3. Kaplan–Meier plot of (A) patients with melanoma (n = 125) stratified based on the 25th percentile cut-off. (B) Patients with head and neck cancer (n = 120) stratified based on 25th percentile cut-offs. (C) Patients with non-small cell lung cancer (n = 128) stratified based on the 25th percentile cut-off. (D) Patients with renal cell carcinoma (n = 149) stratified based on the 25th percentile cut-off. p-value < 0.05 is considered statically significant.
Cancers 15 04913 g003
Figure 4. Kaplan–Meier plot of (A) patients with melanoma (n = 125) stratified based on a 49.13th percentile cut-off. (B) Patients with head and neck cancer (n = 120) stratified based on a 43.7th percentile cut-off. (C) Patients with non-small cell lung cancer (n = 128) stratified based on a 57.75th percentile cut-off. (D) Patients with renal cell carcinoma (n = 149) stratified based on a 36.9th percentile cut-off. p-value <0.05 is considered statically significant.
Figure 4. Kaplan–Meier plot of (A) patients with melanoma (n = 125) stratified based on a 49.13th percentile cut-off. (B) Patients with head and neck cancer (n = 120) stratified based on a 43.7th percentile cut-off. (C) Patients with non-small cell lung cancer (n = 128) stratified based on a 57.75th percentile cut-off. (D) Patients with renal cell carcinoma (n = 149) stratified based on a 36.9th percentile cut-off. p-value <0.05 is considered statically significant.
Cancers 15 04913 g004
Table 1. Baseline disease and immunotherapy characteristics.
Table 1. Baseline disease and immunotherapy characteristics.
VariablePatients (N = 522)
Age (in years)
 Median (range)63 (19–90)
Sex, n (%)360 (69)
 Male162 (31)
 Female
Race, n (%)496 (95)
 White16 (3)
 Black10 (2)
 Other
ECOG performance status at diagnosis, n (%)
 099 (19)
 184 (16)
 210 (2)
 Unknown329 (63)
Cancer type, n (%)
 Renal cell carcinoma149 (28.5%)
 Non-small cell lung cancer128 (24.5%)
 Melanoma125 (23.9%)
 Head and neck cancer120 (23.0%)
Prior systemic therapy, n (%)
 1 Prior line198 (38)
 2+ Prior line324 (62)
First Immune checkpoint inhibitors, n (%)
 Nivolumab219 (42.0%)
 Pembrolizumab202 (38.7%)
 Ipilimumab + nivolumab69 (13.2%)
 Ipilimumab30 (5.6%)
 Avelumab1 (0.2%)
 Cemiplimab1 (0.2%)
ECOG: Eastern Cooperative Oncology Group. 0: Fully active; no performance restrictions; 1: strenuous physical activity restricted; fully ambulatory and able to carry out light work; 2: capable of all self-care but unable to carry out any work activities, and up and about for >50% of waking hours; 3: capable of only limited self-care, and confined to bed or chair for >50% of waking hours; 4: completely disabled, cannot carry out any self-care, and totally confined to bed or chair.
Table 2. Average C-index and associated 95% Cl and log-rank p-values for the entire cohort.
Table 2. Average C-index and associated 95% Cl and log-rank p-values for the entire cohort.
Percentile Cut-OffAvg. C-Index (95% CI)Log-Rank Test p-Value *
Cut-off = 25th percentile0.5402 (0.5345, 0.5459)<0.001
Cut-off = 43.5th percentile0.5528 (0.5466, 0.5591)<0.001
* p-value < 0.05 was tested for significance using a log-rank test.
Table 3. Average C-index and associated 95% Cl and log-rank p-values per each cancer category.
Table 3. Average C-index and associated 95% Cl and log-rank p-values per each cancer category.
Cancer CategoryPercentile Cut-OffAvg. C-Index (95% CI)Log-Rank Test p-Value *
Head and neck43.7th0.55 (0.54, 0.56)0.04
Renal Cell Carcinoma36.9th0.56 (0.55, 0.58)0.17
Non-small cell lung 57.75th0.47 (0.46, 0.48)0.77
Melanoma49.13th0.58 (0.57, 0.59)0.009
* p-value < 0.05 was tested using a log-rank test.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Eljilany, I.; Saghand, P.G.; Chen, J.; Ratan, A.; McCarter, M.; Carpten, J.; Colman, H.; Ikeguchi, A.P.; Puzanov, I.; Arnold, S.; et al. The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors. Cancers 2023, 15, 4913. https://doi.org/10.3390/cancers15204913

AMA Style

Eljilany I, Saghand PG, Chen J, Ratan A, McCarter M, Carpten J, Colman H, Ikeguchi AP, Puzanov I, Arnold S, et al. The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors. Cancers. 2023; 15(20):4913. https://doi.org/10.3390/cancers15204913

Chicago/Turabian Style

Eljilany, Islam, Payman Ghasemi Saghand, James Chen, Aakrosh Ratan, Martin McCarter, John Carpten, Howard Colman, Alexandra P. Ikeguchi, Igor Puzanov, Susanne Arnold, and et al. 2023. "The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors" Cancers 15, no. 20: 4913. https://doi.org/10.3390/cancers15204913

APA Style

Eljilany, I., Saghand, P. G., Chen, J., Ratan, A., McCarter, M., Carpten, J., Colman, H., Ikeguchi, A. P., Puzanov, I., Arnold, S., Churchman, M., Hwu, P., Conejo-Garcia, J., Dalton, W. S., Weiner, G. J., El Naqa, I. M., & Tarhini, A. A. (2023). The T Cell Immunoscore as a Reference for Biomarker Development Utilizing Real-World Data from Patients with Advanced Malignancies Treated with Immune Checkpoint Inhibitors. Cancers, 15(20), 4913. https://doi.org/10.3390/cancers15204913

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop