1. Introduction
In 2020, gastric cancer was the sixth-most-frequent cancer with 1,089,103 new patients and it was the third-leading cause of cancer death with 768,793 deaths worldwide [1]. In Japan, it is the second-most common cancer with 129,475 new cases in 2017 and the third-leading cause of cancer-related deaths with 42,931 fatalities in 2019 [2]. Efforts to combat gastric cancer have led to the exploration of immune-checkpoint inhibitors (ICIs) in various clinical trials. Notably, the phase Ib study (KEYNOTE-012) demonstrated a response rate of 22.2% in programmed cell death ligand 1 (PD-L1)-positive patients treated with pembrolizumab [3]. Similarly, the phase Ib study (JAVELIN) revealed a response rate of 9.0% for avelumab, an anti-PD-L1 antibody, in patients with advanced gastric or gastroesophageal-junction cancer [4]. Furthermore, the phase III trial (ATTRACTION-2) demonstrated the efficacy of nivolumab, an anti-programmed death receptor 1 (PD-1) antibody, in improving OS by 1.2 months (5.3 months vs. 4.1 months, hazard ratio (HR) 0.63 (95% CI 0.50–0.78), p < 0.001) compared to placebo in patients with unresectable advanced or recurrent gastric cancer refractory to standard therapy [5]. Subsequently, based on these encouraging results, nivolumab gained approval in Japan on 22 September 2017, for treating unresectable advanced, or recurrent gastric cancer that has progressed after chemotherapy. Notably, nivolumab was designated as the standard of care for third-line treatment in the Japanese gastric cancer treatment guidelines 2018 (5th edition) [6].
A limited response rate characterizes the therapeutic efficacy of nivolumab in gastric cancer, but those patients who do respond often achieve long-term survival. Conversely, patients who are unlikely to respond are compelled to switch to alternative chemotherapy strategies, such as trifluridine/tipiracil, to achieve stable disease [7]. Identifying biomarkers that can predict both initial and long-term treatment responses is essential to improve treatment outcomes. Additionally, understanding the mechanisms of treatment resistance is crucial for selecting appropriate combination therapies. Therefore, there is a clear need to establish reliable biomarkers and investigate treatment-resistance mechanisms to develop more effective therapies.
In our recent studies [8,9], we developed an “immunogram” using RNA sequencing (RNA-Seq) data to assess the intratumoral immune response characteristics in individual patients. Building upon this, we applied the immunogram to a cohort of 29 gastric cancer surgery cases, successfully classifying gastric cancer into four distinct groups [10]. In addition to utilizing RNA-Seq data, we employed flow cytometry and liquid factor analysis to comprehensively analyze the intratumoral immune response. This multi-layered approach revealed that our classification did not align with conventional clinical classifications. However, it provided valuable insights into various aspects of the intratumoral immune response, including cancer antigen counts, epithelial–mesenchymal transition, genetic mutations, tumor microenvironment, tumor-infiltrating T lymphocyte (TIL) dysfunction, and infiltration exclusion. As a result, our approach proved to be a valuable tool for evaluating the intratumoral immune response.
In this prospective study, we used the aforementioned methods to compare the intratumoral immune response in gastric cancer patients before nivolumab therapy and during early-on treatment with the drug. Our primary objectives were to uncover the immunological mechanisms contributing to treatment resistance and its underlying factors and to identify biomarkers that can predict the efficacy of nivolumab treatment for gastric cancer. By investigating these aspects, we aimed to enhance our understanding of nivolumab’s effectiveness and shed light on potential strategies to overcome treatment resistance in gastric cancer.
2. Results
2.1. Study Design and Baseline Characteristic
This study was a single-arm, prospective interventional study at the University of Tokyo Hospital (Figure 1). The patients in the study received nivolumab monotherapy as their third- or later-line therapy at either 3 mg/kg (until October 2018) or 240 mg/body (from November 2018). Endoscopic biopsies were performed before initiating treatment (pre-treatment) and after completing three courses of nivolumab therapy (early-on treatment). The purpose of these biopsies was to analyze the intratumoral immune response at an early stage of nivolumab treatment in patients with unresectable advanced or recurrent gastric cancer, with the aim of identifying predictive biomarkers for nivolumab efficacy. The primary endpoint of the clinical trial was the durable clinical benefit rate (DCBR), with DCB defined as disease control for more than 4 months. The secondary endpoints included the ORR, disease control rate (DCR), PFS, OS, and immunological response.
From 1 October 2017, a total of 19 patients (10 males, median age 67 years) were enrolled, all of whom had received nivolumab as the third- or later-line of therapy. The recruitment for the study was supposed to reach 50 patients, but it was halted before the planned completion on 31 March 2023. This was due to the challenges in recruiting participants for endoscopy studies during the COVID-19 pandemic. Additionally, the standard therapy for eligible patients changed, making it even more difficult to recruit new patients. Clinical data were censored on 31 March 2023; the median follow-up time was 5.1 months (0.8–64.6). All patients underwent pre-biopsy, and in 18 cases, sufficient tumor tissue was obtained for analysis (hereafter referred to as immunological-analysis cases). Fourteen patients could complete three courses of nivolumab therapy and post-biopsy could be performed; sufficient tumor tissue for analysis was obtained in all cases (Table 1 and Table S1).
2.2. Efficacy of Nivolumab Monotherapy
The Kaplan–Meier survival analysis revealed that the median OS and PFS for all patients were 5.5 and 2.3 months, respectively (Figure 2). The median OS and PFS for patients included in the immunological analysis were 5.5 and 2.4 months, respectively (Figure S1).
The DCBR of all patients was 15.8% (95% confidence interval (CI), 3.4–39.6) (n = 3). Among the 18 patients in the immunological analysis, the DCBR was 16.7% (n = 3) (Table 2). ORR, DCR, PFS, and OS are shown in Table 2.
2.3. Immunogram Analysis
To assess the anti-tumor immune response elicited by anti-PD-1 monotherapy, we extracted RNA from biopsy specimens and conducted RNA-Seq analysis. We focused on ten gene sets associated with anti-tumor immunity, proliferation, and metabolism (Table S2). Subsequently, we performed single-sample gene-set enrichment analysis (ssGSEA) on the tumor samples and converted the scores of the ten gene sets into immunogram scores following the methodology outlined in the Section 4. The resulting immunogram scores were plotted on a radar chart to generate individual immunograms (Figure 3), which exhibited considerable variation among patients, indicating that the immune responses within the tumor and tumor microenvironment (TME) were distinct for each case. Notably, in certain patients, the orange line representing the early-on-treatment tumor expanded around the blue line representing the pre-treatment tumor, effectively enclosing it. This observation suggests that anti-tumor immune responses were activated by anti-PD-1 monotherapy in these particular patients (Figure 3). However, immunogram scores were low and did not show improvement with the treatment in the majority of patients. The immunogram score changes before and during early-on treatment for each axis of the immunogram exhibited no significant alterations in most of the axes. However, there was a significant increase in the immunogram score of Tregs in all immunological cases and non-DCB cases. No significant change was observed in DCB cases (Figure S2).
2.4. Biomarker for Anti-PD-1 Monotherapy
We investigated the molecular biomarkers for DCB in anti-PD-1 monotherapy using pre-treatment and early-on treatment samples. Immunogram scores for each axis were compared between patients with and without DCB (DCB and non-DCB). There was no statistically significant difference between these two groups (Figures S3 and S4).
As depicted in Figure S2, both T cell and Treg scores exhibited an increase in many patients following the initiation of nivolumab therapy. The rise in Treg scores was more pronounced than that of T cell scores. Moreover, the escalation of Treg scores in non-DCB patients surpassed those in DCB patients. Consequently, the ratio of T cell/Treg immunogram scores in pre-treatment samples was higher in patients with DCB than in non-DCB patients (Figure 4a). The T cell/Treg signature ratio in early-on treatment samples was also not significant, but had a higher trend in patients with DCB than in non-DCB patients (Figure 4b).
Furthermore, patients with a T cell/Treg signature ratio in immunogram score ≥ 1.2 exhibited longer OS than patients with a ratio < 1.2 (Figure 4c). Additionally, a comparison of T cell/Treg signature ratios in immunogram scores of the PRJEB25780 study in Korea [11] revealed higher scores in responders (Figure 4d).
The ssGSEA analysis was performed on patients with DCB and non-DCB using 140 gene sets from the Human Molecular Signatures Database (MSigDB) (Table S3). Four activated processes were enriched in patients with DCB, while two were enriched in non-DCB patients (Table S4 and Figure 5).
These six gene sets were then applied to responders and non-responders in the PRJEB25780 study (Table S5). Genes down-regulated in response to ultraviolet (UV) radiation (HALLMRK_UV_RESPONSE_DN) was significantly lower in responders than in non-responders (Figure 6).
2.5. TCR Repertoire Analysis
In addition to transcriptome analysis, TCR β genes were amplified from cDNA, and the TCR repertoire was evaluated based on the complementarity-determining region 3 (CDR3) sequences. The Shannon–Weaver index (SWI), which represents the diversity of the TCR repertoire, was compared (Table 3 and Figure S5). SWIs of DCB patients were found to be larger than those of non-DCB patients in pre-treatment tumors, while SWIs of DCB patients were smaller than those of non-DCB patients in early-on-treatment tumors. In patients with DCB, the SWI of early-on-treatment tumors was smaller than that of pre-treatment tumors, whereas in patients without DCB, the SWI of early-on-treatment tumors was larger than that of pre-treatment samples (Table 3 and Figure S5). These results indicate that anti-PD-1 monotherapy increased T cell clonality in DCB patients, while T cell diversity was increased in non-DCB patients.
The clonal replacement of TILs by ICIs has been reported [12,13,14,15]. Therefore, we compared the frequencies of individual T cell clones in tumors before and during treatment using the CDR3 sequences. After anti-PD-1 monotherapy, TCR clones were categorized into three groups: lost, shared, and emerged clones (Figure 7). Next, we compared the number of clonotypes between DCB and non-DCB patients. There were no significant differences in the numbers of lost, shared, or emerged clonotypes between DCB and non-DCB patients.
2.6. Immunological Classification of Gastric Cancer
Previously, we conducted an examination of surgically resected gastric cancer in the BKT study and proposed an immunological classification of gastric cancer based on immunogram scores. Surgically resected gastric cancer demonstrated four immune signatures representing the main subtypes: Hot1, Hot2, Intermediate, and Cold. To create a straightforward approach for immunological subtyping of gastric cancer without resorting to cluster analysis, a decision tree was formulated. The distinction between Immune-Hot and Immune-Cold tumors was based on the sum of immunogram scores for innate immunity (IGS1), priming and activation (IGS2), T cells (IGS3), IFN-γ response (IGS4), inhibitory molecules (IGS5), and Tregs (IGS6), with a threshold of <18.21 or >18.21. Intermediate tumors were subsequently identified by evaluating the recognition of tumor cells (IGS8 in this study), with a threshold of <3.78 or >3.78. Lastly, Hot1 and Hot2 categories were determined based on glycolysis (IGS10 in this study), with a threshold of <2.11 or >2.11 [10].
Consequently, we applied the decision tree for immunological subtypes of gastric cancer to 18 pre-treatment biopsy samples. Out of the 18 pre-treatment tumors in this study, 15 were classified as Cold tumors, 3 were classified as Intermediate subtype, and none were classified as Hot1 or Hot2 tumors. These results indicate that the tumors in patients who received nivolumab monotherapy as their third- or later-line therapy were more advanced and immunologically suppressed compared to the surgically resectable tumors.
3. Discussion
In this study, we analyzed the tumor immune response in a group that achieved DCB with nivolumab monotherapy as a third-line treatment for gastric cancer, compared to a group that did not achieve DCB. DCB was observed in 15.8% of the intent-to-treat population. The secondary endpoints, including an ORR of 15.8%, median OS of 5.5 months, and median PFS of 2.3 months, demonstrated comparable results to the ATTRACTION-2 trial [5]. Despite a slightly unfavorable patient background with Eastern Cooperative Oncology Group (ECOG) performance status (PS) 2 of 37% in the real-world setting, the efficacy of nivolumab was confirmed.
No significant changes were observed in individual immunograms before or early-on in treatment, and no consistent changes were identified among cases that achieved DCB. Furthermore, TCR repertoire analysis revealed that newly emerged clonotypes were detected in the early-on-treatment tumors. However, clonal replacement was not associated with the efficacy of the treatment. As mentioned earlier (Figure 8), at the stage of third- or later-line therapy, the tumors had already transitioned into Cold tumors, and it was believed that it was no longer a stage where the immune response could be altered simply by inhibiting PD-1. In a considerable number of cases with early disease progression, collecting early-on-treatment tissue samples was not feasible. If these cases were included in the analysis, the mechanism of resistance might have been observed between the periods before treatment and early on in the treatment.
To investigate potential biomarkers predicting treatment efficacy, we conducted an analysis of gene expression profiles from pre-treatment tumor transcriptomes in two separate groups: the DCB group (cases from this study) and the responder group (cases from the PRJEB25780 study). In terms of the signature of tumor-infiltrating immune cells, we observed that a high T cell/Treg ratio and a low signature of genes down-regulated in response to UV radiation were associated with successful DCB and treatment response. Despite the limited number of cases analyzed, we observed that clinical benefits were accompanied by an immune response.
The efficacy of anti-PD-1 therapy in exerting an anti-tumor effect is believed to rely on the balance between reactivating effector T cells and enhancing their proliferation, while concurrently suppressing PD-1+ Treg cells [16]. A previous investigation, employing flow cytometry to analyze TIL extracted from biopsy specimens of malignant melanoma (n = 12), non-small cell lung cancer (n = 27), and gastric cancer (n = 48) before anti-PD-1 monotherapy, demonstrated that the ratio of PD-1+CD8+ T cells to PD-1+ Treg cells in the tumor microenvironment could predict the clinical efficacy of anti-PD-1 monotherapy [17]. Consistent with these earlier findings, our study, utilizing RNA-seq on tumor tissue, suggested that the pre-treatment T cell/Treg signature ratio could also predict the clinical efficacy of anti-PD-1 monotherapy. Although our emphasis was on the intratumoral immune environment, it has been reported that soluble PD-1 and PD-L1 also play a role in modulating anti-tumor immunity and the efficacy of anti-PD-1 therapy [18,19,20]. In addition, the influence of genetic alterations on the tumor microenvironment is also implicated. The inclusion of these factors will enhance our understanding of the anti-tumor immune response and the mechanisms of resistance to anti-PD-1 therapy.
Compared to surgical cases, the majority of cases in this study were characterized by Cold tumors (Figure 8). All cases were in the third-line or later treatments, and it is possible that they developed into Cold tumors as a result of tumor progression or the influence of past chemotherapy. In gastric cancer, it may be effective to use immune-checkpoint inhibitors at a more frontline setting before the tumor transitions from a Hot to Cold state. In fact, the results of two phase III (CheckMate 649 and ATTRACTION-4) trials that evaluated the combination of chemotherapy with nivolumab in first-line treatment have established it as the standard treatment [21,22]. Moreover, for such Cold tumors with exhausted T cells, chimeric antigen receptor T cell (CAR-T) therapy, in which activated T cells themselves are infused, may be useful. Interim analysis results of a phase 1 study of CAR-T therapy redirecting claudin (CLDN) 18.2, which is a protein that tightly joins gastric epithelial cells and is expressed by around 60% of gastric cancers, in patients with CLDN18.2-positive gastrointestinal cancers, showed promising efficacy with an acceptable safety profile especially in gastric cancer [23].
As a limitation of this study, due to the impact of the COVID-19 pandemic, we were only able to accumulate a significantly smaller number of cases than originally planned (50 cases). As a result, there is a possibility that the statistical power of the analysis is insufficient. Additionally, a considerable number of cases were unable to have early-on-treatment tissue samples collected primarily due to disease progression, potentially limiting the ability for adequate before-and-after comparisons. Now that the COVID-19 situation has improved, we have initiated a prospective study to evaluate the tumor immune response before and after nivolumab-containing treatment in primary therapy for gastric cancer.
4. Materials and Methods
4.1. Patient Selection
Patients with unresectable advanced or recurrent gastric or esophagogastric-junction adenocarcinoma were eligible for inclusion in the study if their tumor tissues could be accessed through endoscopy at the University of Tokyo Hospital. Other criteria for eligibility included being 20 years of age or older, having an ECOG PS of 0–2, and having an evaluable lesion according to RECIST version 1.1 (measurable lesions were not required). Additionally, patients needed to be at least 42 days post-failure of standard chemotherapy and at least 14 days after their last chemotherapy infusion. Adequate organ and marrow functions were determined through laboratory tests, including assessments of neutrophil count, hemoglobin (Hb) levels, platelet count (PLT), total bilirubin (T-Bil), aspartate aminotransferase (AST), alanine aminotransferase (ALT), and serum creatinine (Cre).
The study had certain exclusion criteria, which were as follows: patients with synchronous or metachronous double cancers, except for intramucosal tumors that had been curatively resected through local therapy within the past 5 years; patients with active infections requiring systemic therapy; patients with active autoimmune diseases or a history of chronic or recurring autoimmune diseases; patients with a history of interstitial pneumonia, pulmonary fibrosis, or irradiation pneumonitis; patients with active diverticulitis or inflammatory bowel disease; patients with poorly controlled diabetes mellitus or thyroid diseases; patients with unstable angina within the past 3 weeks or a history of acute myocardial infarction within the past 3 months; patients with severe psychological illness; pregnant or lactating women or women of childbearing potential; patients within 4 weeks of receiving a live vaccination or 2 weeks of receiving an inactivated vaccination; and patients deemed unfit to participate in the study by the investigator.
4.2. Clinical Sample Processing and RNA Extraction
The tissues were collected immediately after endoscopic biopsy, stored in RNAlater Stabilization Solution (Thermo Fisher Scientific K.K., Tokyo, Japan) and cryopreserved until use. Total RNA samples from tissues were extracted using the AllPrep DNA/RNA/miRNA Universal Kits (Qiagen, Hilden, Germany), following the manufacturer’s instructions. The extracted RNAs were then assessed for quality and quantity. For next-generation sequencing (NGS), RNA samples meeting the following criteria were selected: a concentration of ≥20.0 ng/μL, a total amount of ≥0.4 μg, and an RNA integrity number (RIN) of ≥7.0, as assessed using the Agilent 2200 TapeStation (Agilent Technologies, Santa Clara, CA, USA).
4.3. RNA-Sequence (RNA-Seq)
For RNA-Seq library preparation, we used the NEBNext® UltraTM RNA Library Prep Kit for Illumina® (Agilent Technologies), following the manufacturer’s protocols. The prepared libraries were sequenced as 150 bp paired-end reads on the NovaSeq platform (Illumina, San Diego, CA, USA) at VERITAS (Danvers, MA, USA). On average, each sample yielded approximately 35.1 million reads of 150 base pairs in length. The obtained reads were aligned to the reference genome (GRCh38/hg38) using STAR (v.2.5.2b). Expression values were calculated as fragments per kilobase of exon per million fragments mapped (FPKM) using HTSeq (v.0.6.1) and the R programming language (version 3.4.3):
4.4. Computational Methods to Analyze RNA-Seq Data
We calculated a ssGSEA score using R version 3.6.2 with the GSVA package version 1.38.2. To depict the immunological status of the tumor in each patient, we constructed an immunogram based on RNA-Seq data [9]. Incorporated gene sets are listed in Table S2. Similar ssGSEA was applied to the Cancer Genome Atlas (TCGA) mRNA data of 375 gastric cancer patients. We obtained the mean (M) and standard deviation (SD) of the ssGSEA score of these 375 gastric cancer patients for each gene set. The score for each axis of the immunogram in each patient was calculated as the (immunogram score) = 3 + 1.5 × (ssGSEA score − M)/SD. This formula was applied for all axes of the immunogram of a patient.
4.5. TCR Repertoire Analysis
TCR genes were amplified using adaptor ligation-mediated PCR [24]. High-throughput sequencing was performed using the Illumina Miseq paired-end platform (2 × 300 bp) (Illumina, San Diego, CA). Assignment of TRBV and TRBJ segments in TCR genes was performed based on the international ImMunoGeneTics information system® (IMGT) database (
4.6. Hierarchical Clustering
We utilized an unsupervised hierarchical clustering algorithm for the transcriptome analysis data. This analysis used R version 3.6.2 with the pheatmap package version 1.0.12. To generate the hierarchical clustering, we calculated the squared Euclidean distance between the samples. This distance measure quantifies the dissimilarity between samples based on their transcriptome profiles. We then applied an agglomerative algorithm with Ward’s method, which iteratively merges clusters to minimize the within-cluster variance.
4.7. Statistical Analysis
PFS and OS were estimated using the Kaplan–Meier method and the log-rank test. Data were censored on 31 March 2023. Patients who were lost to follow-up were censored at the date of last contact or follow-up. PFS was calculated from the date of study enrollment to the date of disease progression or death from any cause. OS was calculated from the date of study enrollment to the date of death from any cause. Patients who were alive on 31 March 2023, were censored for OS analysis. Tumor response was evaluated according to the Response Evaluation Criteria in Solid Tumors, version 1.1 [25], based on computed tomography (CT) findings. The best overall response was assessed as complete response (CR), partial response (PR), stable disease (SD), non-CR/non-PD, or progressive disease (PD). Patients with clinically progressed disease status were defined as PD without undergoing a CT scan in this study. Disease control was defined as no PD.
The Kaplan–Meier method with a log-rank test and Cox regression analysis were performed to evaluate relapse-free survival and OS. All statistical analyses were performed with EZR (Saitama Medical Center, Jichi Medical University, Saitama, Japan), which is a graphical user interface for R (The R Foundation for Statistical Computing, Vienna, Austria, version 3.2.1). More precisely, it is a modified version of R commander designed to add statistical functions frequently used in biostatistics [26].
5. Conclusions
In conclusion, this study underscores the importance of the tumor immune response in influencing the effectiveness of nivolumab for gastric cancer. Given the challenging scenario where the tumor-infiltrating immune response may be compromised by previous treatments in third-line therapy, initiating treatment at an earlier stage is deemed advisable to enhance immunotherapeutic outcomes. Nevertheless, for the optimal utilization of immune-checkpoint inhibitors in gastric cancer treatment, it is crucial to conduct additional validation studies involving larger cohorts. Ongoing research in this direction will contribute to the development of more effective and personalized therapeutic strategies.
Conceptualization, Y.S. (Yasuyoshi Sato), Y.S. (Yasuyuki Seto) and K.K.; methodology, H.Y.; software, T.H.; investigation Y.K., K.N., A.K., N.S., R.I., Y.O., K.Y., S.A. and S.N.; formal analysis T.K.; data curation, Y.S. (Yasuyoshi Sato) and Y.K.; writing—original draft preparation, Y.S. (Yasuyoshi Sato) and K.K.; writing—review and editing, H.Y., Y.K., K.N. and Y.S. (Yasuyuki Seto); visualization, Y.S. (Yasuyoshi Sato) and Y.K.; supervision, Y.S. (Yasuyuki Seto); project administration, K.K.; funding acquisition, Y.S. (Yasuyuki Seto) and K.K. All authors have read and agreed to the published version of the manuscript.
The research protocol received approval from the Ethical Committee of the University of Tokyo Hospital (P2017008). The study is registered with the University Hospital Medical Information Network Clinical Trials Registry (UMIN-CTR) under the registration number UMIN000033110. This study was conducted in compliance with the principles outlined in the Declaration of Helsinki and the Good Clinical Practice guidelines established by the International Conference on Harmonization of Technical Requirements for Registration of Pharmaceuticals for Human Use.
Informed consent was obtained from all subjects involved in the study.
Data are deposited on the Japanese Genotype–Phenotype Archive (Accession no. JGAS000639) [
The authors thank Yaeko Furuhashi and Mikiko Shibuya for excellent technical assistance. The authors also thank Kazuhiko Yamada (Department of Surgery, Center Hospital of the National Center for Global Health and Medicine) for referral of patients.
Sato reports personal fees from ONO Pharmaceutical Co., Ltd.; Bristol-Myers Squibb Company; MSD KK; TAIHO Pharmaceutical Co., Ltd.; and Daiichi Sankyo Co., Ltd. outside the submitted work. Kakimi reports grants from TAKARA BIO Inc. outside the submitted work. The Department of Immunotherapeutics, The University of Tokyo Hospital, is an endowed department by TAKARA BIO Inc. The other authors have no competing interests to disclose. Author Tetsuro Hisayoshi is employed by the company cBioinformatics, Inc. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflicts of interest.
DCB | Durable clinical benefit |
TCR | T cell receptor |
ORR | Objective response rate |
PFS | Progression-free survival |
OS | Overall survival |
Treg | Regulatory T cell |
ICI | Immune-checkpoint inhibitor |
PD-L1 | Programmed cell death ligand 1 |
PD-1 | Programmed death receptor 1 |
HR | Hazard ratio |
RNA-seq | RNA sequencing |
TIL | Tumor-infiltrating T lymphocyte |
DCBR | Durable clinical benefit rate |
DCR | Disease control rate |
HER2 | Human epidermal growth factor receptor 2 |
ssGSEA | Single-sample gene-set enrichment analysis |
TME | Tumor microenvironment |
MSigDB | Molecular Signatures Database |
UV | Ultraviolet |
cDNA | Complementary DNA |
CDR3 | Complementarity-determining region 3 |
SWI | Shannon–Weaver index |
ECOG | Eastern Cooperative Oncology Group |
PS | Performance status |
CAR-T | Chimeric antigen receptor T cell |
CLDN | Claudin |
Hb | Hemoglobin |
PLT | Platelet count |
T-Bil | Total bilirubin |
AST | Aspartate aminotransferase |
ALT | Alanine aminotransferase |
Cre | Creatinine |
NGS | Next-generation sequencing |
FPKM | Fragments per kilobase of exon per million fragments mapped |
TCGA | The Cancer Genome Atlas |
M | Mean |
SD | Standard deviation |
CT | Computed tomography |
CR | Complete response |
PR | Partial response |
SD | Stable disease |
PD | Progressive disease |
Footnotes
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.
Figure 2. The Kaplan–Meier survival analysis. Survival curves of OS (a) and PFS (b) of all patients (n = 19).
Figure 3. Immunograms for cancer–immunity interactions in 18 patients with gastric cancer. The immunograms illustrated the differences between pre-treatment and early-on-treatment tumors using blue and orange lines, respectively. Red squares indicate the DCB cases. The best overall response for each case is listed in the lower-right corner.
Figure 4. T cell/Treg signature ratios. T cell/Treg signature ratios in immunogram scores in cases with and without DCB (DCB vs. non-DCB) (a) before nivolumab treatment, and (b) during early-on nivolumab treatment. (c) The Kaplan–Meier survival curves of OS stratified by T cell/Treg signature ratio before nivolumab treatment (≥1.2; <1.2). (d) Pre-treatment T cell/Treg signature ratios in immunogram scores before ICI treatment and in 45 earlier-ICI-treated cases in the PRJEB25780 study.
Figure 5. ssGSEA score of 6 selected gene sets from the Human MSigDB. ssGSEA scores were compared between DCB and non-DCB patients. (a) ImSig_CLASSICAL_M1_vs._ALTERNATIVE_M2_MACROPHAGE. (b) ImSig_NONSUPPRESSIVE_TCELL_vs._ACTIVATED_TREG_DN. (c) ImSig_PD1HIGH_vs._PD1LOW_CD8_TCELL_DN. (d) LM22_B cells naïve. (e) HALLMRK _UV_RESPONSE_DN. (f) ImSig_TH1_vs._TH2_CD4_TCELL_UP.
Figure 6. ssGSEA score of 6 selected gene sets from the Human MSigDB using PRJEB25780 data set. ssGSEA scores were compared between responder and non-responder patients. (a) ImSig_CLASSICAL_M1_vs._ALTERNATIVE_M2_MACROPHAGE. (b) ImSig_NONSUPPRESSIVE_TCELL_vs._ACTIVATED_TREG_DN. (c) ImSig_PD1HIGH_vs._PD1LOW_CD8_TCELL_DN. (d) LM22_B cells naïve. (e) HALLMRK _UV_RESPONSE_DN. (f) ImSig_TH1_vs._TH2_CD4_TCELL_UP.
Figure 7. Pre- and early-on treatment (Tx) T cell clonotypes by TCR repertoire analysis in the tumor. (a) Pre- and early-on-treatment T cell clonotypes by TCR repertoire analysis of OG002 is shown. Numbers of emerged clonotypes (b), lost clonotypes (c), and shared clonotypes (d) were compared between DCB and non-DCB patients.
Figure 8. A decision tree for the immunological subtypes of gastric cancer was implemented on 18 pre-treatment biopsy samples. The cutoff value for identifying the Cold subtype was set at a sum of IGS1 to IGS6 of <18.21. Similarly, IGS8 < 3.78 was established as the cutoff value for the Intermediate subtype, and IGS10 < 2.11 was utilized to distinguish between Hot1 and Hot2 tumors.
Patients’ characteristics.
Characteristic | All (n = 19) | Immunological-Analysis Cases (n = 18) |
---|---|---|
Sex, n (%) | ||
Male | 10 (53) | 9 (50) |
Age, years | ||
Median (range) | 67 (49–87) | 66 (49–87) |
ECOG PS, n (%) | ||
0 | 4 (21) | 4 (22) |
1 | 8 (42) | 8 (44) |
2 | 7 (37) | 6 (33) |
Histology, n (%) | ||
Intestinal | 5 (26) | 4 (22) |
Diffuse | 10 (53) | 10 (56) |
Mixed | 2 (11) | 2 (11) |
Special type | 0 (0) | 0 (0) |
Unknown | 2 (11) | 2 (11) |
HER2, n (%) | ||
Positive | 2 (11) | 2 (11) |
Negative | 13 (68) | 12 (67) |
Unknown | 4 (21) | 4 (22) |
Number of previous chemotherapies, n (%) | ||
0 | 0 (0) | 0 (0) |
1 | 0 (0) | 0 (0) |
2 | 14 (74) | 13 (72) |
3 | 3 (16) | 3 (17) |
≥4 | 2 (11) | 2 (11) |
HER2, human epidermal growth factor receptor 2.
Efficacy outcomes.
All (n = 19) | Immunological-Analysis Cases (n = 18) | |
---|---|---|
Primary endpoint | ||
DCB, n (% (95% CI)) | 3 (15.8 (3.4–39.6)) | 3 (16.7) |
Co-primary endpoints | ||
Objective response, n (%) | 3 (15.8) | 3 (16.7) |
Disease control, n (%) | 4 (21.0) | 4 (22.2) |
PFS | ||
Median (95% CI), months | 2.3 (1.3–2.8) | 2.4 (1.4–2.8) |
1-year rate (% (95% CI)) | 5.3 (0.4–21.4) | 5.6 (0.4–22.4) |
OS | ||
Median (95% CI), months | 5.5 (1.6–7.6) | 5.5 (2.5–8.7) |
1-year rate (95% CI), % | 17.1 (4.3–37.3) | 18.1% (4.5–39.0) |
Shannon–Weaver index by TCR repertoire analysis in tumors.
Pre-Treatment | Early-On Treatment | p Value | ||
---|---|---|---|---|
DCB | OG002 | 5.752 | 5.329 | |
OG010 | 5.630 | 5.320 | ||
OG017 | 3.908 | 3.674 | ||
Average ± SD | 5.097 ± 1.032 | 4.774 ± 0.953 | 0.028 | |
Non-DCB | OG003 | 5.143 | 5.482 | |
OG007 | 4.181 | 5.564 | ||
OG008 | 3.034 | 4.562 | ||
OG009 | 3.702 | 3.320 | ||
OG011 | 5.055 | 6.008 | ||
OG014 | 4.899 | 5.618 | ||
OG015 | 6.060 | 5.371 | ||
OG016 | 4.314 | 4.216 | ||
OG018 | 3.844 | 4.005 | ||
OG019 | 4.451 | 4.577 | ||
Average ± SD | 4.468 ± 0.860 | 4.872 ± 0.864 | 0.115 | |
All | Average ± SD | 4.613 ± 0.899 | 4.850 ± 0.845 | 0.254 |
Supplementary Materials
The supporting information can be downloaded at:
References
1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA. Cancer J. Clin.; 2021; 71, pp. 209-249. [DOI: https://dx.doi.org/10.3322/caac.21660]
2. Cancer Information Service. Cancer Statistics in Japan 2021; Foundation for Promotion of Cancer Research: Tokyo, Japan, 2021.
3. Muro, K.; Chung, H.C.; Shankaran, V.; Geva, R.; Catenacci, D.; Gupta, S.; Eder, J.P.; Golan, T.; Le, D.T.; Burtness, B. et al. Pembrolizumab for patients with PD-L1-positive advanced gastric cancer (KEYNOTE-012): A multicentre, open-label, phase 1b trial. Lancet Oncol.; 2016; 17, pp. 717-726. [DOI: https://dx.doi.org/10.1016/S1470-2045(16)00175-3]
4. Chung, H.C.; Arkenau, H.T.; Lee, J.; Rha, S.Y.; Oh, D.Y.; Wyrwicz, L.; Kang, Y.K.; Lee, K.W.; Infante, J.R.; Lee, S.S. et al. Avelumab (anti-PD-L1) as first-line switch-maintenance or second-line therapy in patients with advanced gastric or gastroesophageal junction cancer: Phase 1b results from the JAVELIN Solid Tumor trial. J. Immunother. Cancer; 2019; 7, 30. [DOI: https://dx.doi.org/10.1186/s40425-019-0508-1]
5. Kang, Y.K.; Boku, N.; Satoh, T.; Ryu, M.H.; Chao, Y.; Kato, K.; Chung, H.C.; Chen, J.S.; Muro, K.; Kang, W.K. et al. Nivolumab in patients with advanced gastric or gastro-oesophageal junction cancer refractory to, or intolerant of, at least two previous chemotherapy regimens (ONO-4538-12, ATTRACTION-2): A randomised, double-blind, placebo-controlled, phase 3 trial. Lancet; 2017; 390, pp. 2461-2471. [DOI: https://dx.doi.org/10.1016/S0140-6736(17)31827-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28993052]
6. Japanese Gastric Cancer Association Japanese gastric cancer treatment guidelines 2018 (5th edition). Gastric Cancer; 2020; 24, pp. 1-21. [DOI: https://dx.doi.org/10.1007/s10120-020-01042-y]
7. Shitara, K.; Doi, T.; Dvorkin, M.; Mansoor, W.; Arkenau, H.T.; Prokharau, A.; Alsina, M.; Ghidini, M.; Faustino, C.; Gorbunova, V. et al. Trifluridine/tipiracil versus placebo in patients with heavily pretreated metastatic gastric cancer (TAGS): A randomised, double-blind, placebo-controlled, phase 3 trial. Lancet Oncol.; 2018; 19, pp. 1437-1448. [DOI: https://dx.doi.org/10.1016/S1470-2045(18)30739-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30355453]
8. Karasaki, T.; Nagayama, K.; Kuwano, H.; Nitadori, J.I.; Sato, M.; Anraku, M.; Hosoi, A.; Matsushita, H.; Morishita, Y.; Kashiwabara, K. et al. An Immunogram for the Cancer-Immunity Cycle: Towards Personalized Immunotherapy of Lung Cancer. J. Thorac. Oncol.; 2017; 12, pp. 791-803. [DOI: https://dx.doi.org/10.1016/j.jtho.2017.01.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28088513]
9. Kobayashi, Y.; Kushihara, Y.; Saito, N.; Yamaguchi, S.; Kakimi, K. A novel scoring method based on RNA-Seq immunograms describing individual cancer-immunity interactions. Cancer Sci.; 2020; 111, pp. 4031-4040. [DOI: https://dx.doi.org/10.1111/cas.14621]
10. Sato, Y.; Wada, I.; Odaira, K.; Hosoi, A.; Kobayashi, Y.; Nagaoka, K.; Karasaki, T.; Matsushita, H.; Yagi, K.; Yamashita, H. et al. Integrative immunogenomic analysis of gastric cancer dictates novel immunological classification and the functional status of tumor-infiltrating cells. Clin. Transl. Immunol.; 2020; 9, e1194. [DOI: https://dx.doi.org/10.1002/cti2.1194]
11. Kim, S.T.; Cristescu, R.; Bass, A.J.; Kim, K.; Odegaard, J.I.; Kim, K.; Liu, X.Q.; Sher, X.; Jung, H.; Lee, M. et al. Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancer. Nat. Med.; 2018; 24, pp. 1449-1458. [DOI: https://dx.doi.org/10.1038/s41591-018-0101-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30013197]
12. Yost, K.E.; Satpathy, A.T.; Wells, D.K.; Qi, Y.; Wang, C.; Kageyama, R.; McNamara, K.L.; Granja, J.M.; Sarin, K.Y.; Brown, R.A. et al. Clonal replacement of tumor-specific T cells following PD-1 blockade. Nat. Med.; 2019; 25, pp. 1251-1259. [DOI: https://dx.doi.org/10.1038/s41591-019-0522-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31359002]
13. Wu, T.D.; Madireddi, S.; de Almeida, P.E.; Banchereau, R.; Chen, Y.J.J.; Chitre, A.S.; Chiang, E.Y.; Iftikhar, H.; O’Gorman, W.E.; Au-Yeung, A. et al. Peripheral T cell expansion predicts tumour infiltration and clinical response. Nature; 2020; 579, pp. 274-278. [DOI: https://dx.doi.org/10.1038/s41586-020-2056-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32103181]
14. Luoma, A.M.; Suo, S.; Wang, Y.; Gunasti, L.; Porter, C.B.M.; Nabilsi, N.; Tadros, J.; Ferretti, A.P.; Liao, S.; Gurer, C. et al. Tissue-resident memory and circulating T cells are early responders to pre-surgical cancer immunotherapy. Cell; 2022; 185, pp. 2918-2935.e29. [DOI: https://dx.doi.org/10.1016/j.cell.2022.06.018] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35803260]
15. Nagasaki, J.; Inozume, T.; Sax, N.; Ariyasu, R.; Ishikawa, M.; Yamashita, K.; Kawazu, M.; Ueno, T.; Irie, T.; Tanji, E. et al. PD-1 blockade therapy promotes infiltration of tumor-attacking exhausted T cell clonotypes. Cell Rep.; 2022; 38, 110331. [DOI: https://dx.doi.org/10.1016/j.celrep.2022.110331] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35108529]
16. Kamada, T.; Togashi, Y.; Tay, C.; Ha, D.; Sasaki, A.; Nakamura, Y.; Sato, E.; Fukuoka, S.; Tada, Y.; Tanaka, A. et al. PD-1 + regulatory T cells amplified by PD-1 blockade promote hyperprogression of cancer. Proc. Natl. Acad. Sci. USA; 2019; 116, 201822001. [DOI: https://dx.doi.org/10.1073/pnas.1822001116] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31028147]
17. Kumagai, S.; Togashi, Y.; Kamada, T.; Sugiyama, E.; Nishinakamura, H.; Takeuchi, Y.; Vitaly, K.; Itahashi, K.; Maeda, Y.; Matsui, S. et al. The PD-1 expression balance between effector and regulatory T cells predicts the clinical efficacy of PD-1 blockade therapies. Nat. Immunol.; 2020; 21, pp. 1346-1358. [DOI: https://dx.doi.org/10.1038/s41590-020-0769-3]
18. Serratì, S.; Guida, M.; Di Fonte, R.; De Summa, S.; Strippoli, S.; Iacobazzi, R.M.; Quarta, A.; De Risi, I.; Guida, G.; Paradiso, A. et al. Circulating extracellular vesicles expressing PD1 and PD-L1 predict response and mediate resistance to checkpoint inhibitors immunotherapy in metastatic melanoma. Mol. Cancer; 2022; 21, 20. [DOI: https://dx.doi.org/10.1186/s12943-021-01490-9]
19. Ohkuma, R.; Ieguchi, K.; Watanabe, M.; Takayanagi, D.; Goshima, T.; Onoue, R.; Hamada, K.; Kubota, Y.; Horiike, A.; Ishiguro, T. et al. Increased plasma soluble pd-1 concentration correlates with disease progression in patients with cancer treated with anti-pd-1 antibodies. Biomedicines; 2021; 9, 1929. [DOI: https://dx.doi.org/10.3390/biomedicines9121929]
20. Sagawa, R.; Sakata, S.; Gong, B.; Seto, Y.; Takemoto, A.; Takagi, S.; Ninomiya, H.; Yanagitani, N.; Nakao, M.; Mun, M. et al. Soluble PD-L1 works as a decoy in lung cancer immunotherapy via alternative polyadenylation. JCI Insight.; 2022; 7, e153323. [DOI: https://dx.doi.org/10.1172/jci.insight.153323]
21. Janjigian, Y.Y.; Shitara, K.; Moehler, M.; Garrido, M.; Salman, P.; Shen, L.; Wyrwicz, L.; Yamaguchi, K.; Skoczylas, T.; Campos Bragagnoli, A. et al. First-line nivolumab plus chemotherapy versus chemotherapy alone for advanced gastric, gastro-oesophageal junction, and oesophageal adenocarcinoma (CheckMate 649): A randomised, open-label, phase 3 trial. Lancet; 2021; 398, pp. 27-40. [DOI: https://dx.doi.org/10.1016/S0140-6736(21)00797-2]
22. Kang, Y.K.; Chen, L.T.; Ryu, M.H.; Oh, D.Y.; Oh, S.C.; Chung, H.C.; Lee, K.W.; Omori, T.; Shitara, K.; Sakuramoto, S. et al. Nivolumab plus chemotherapy versus placebo plus chemotherapy in patients with HER2-negative, untreated, unresectable advanced or recurrent gastric or gastro-oesophageal junction cancer (ATTRACTION-4): A randomised, multicentre, double-blind, placebo-contr. Lancet Oncol.; 2022; 23, pp. 234-247. [DOI: https://dx.doi.org/10.1016/S1470-2045(21)00692-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35030335]
23. Qi, C.; Gong, J.; Li, J.; Liu, D.; Qin, Y.; Ge, S.; Zhang, M.; Peng, Z.; Zhou, J.; Cao, Y. et al. Claudin18.2-specific CAR T cells in gastrointestinal cancers: Phase 1 trial interim results. Nat. Med.; 2022; 28, pp. 1189-1198. [DOI: https://dx.doi.org/10.1038/s41591-022-01800-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35534566]
24. Kitaura, K.; Shini, T.; Matsutani, T.; Suzuki, R. A new high-throughput sequencing method for determining diversity and similarity of T cell receptor (TCR) α and β repertoires and identifying potential new invariant TCR α chains. BMC Immunol.; 2016; 17, 38. [DOI: https://dx.doi.org/10.1186/s12865-016-0177-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27729009]
25. Eisenhauer, E.A.; Therasse, P.; Bogaerts, J.; Schwartz, L.H.; Sargent, D.; Ford, R.; Dancey, J.; Arbuck, S.; Gwyther, S.; Mooney, M. et al. New response evaluation criteria in solid tumours: Revised RECIST guideline (version 1.1). Eur. J. Cancer; 2009; 45, pp. 228-247. [DOI: https://dx.doi.org/10.1016/j.ejca.2008.10.026] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19097774]
26. Kanda, Y. Investigation of the freely available easy-to-use software “EZR” for medical statistics. Bone Marrow Transpl.; 2013; 48, pp. 452-458. [DOI: https://dx.doi.org/10.1038/bmt.2012.244]
27. Kodama, Y.; Mashima, J.; Kosuge, T.; Katayama, T.; Fujisawa, T.; Kaminuma, E.; Ogasawara, O.; Okubo, K.; Takagi, T.; Nakamura, Y. The DDBJ Japanese genotype-phenotype archive for genetic and phenotypic human data. Nucleic Acids Res.; 2015; 43, pp. D18-D22. [DOI: https://dx.doi.org/10.1093/nar/gku1120]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
We investigated the tumor immune response in gastric cancer patients receiving third-line nivolumab monotherapy to identify immune-related biomarkers for better patient selection. Nineteen patients (10 males, median age 67 years) who received nivolumab as a third- or later-line therapy were enrolled. We analyzed the tumor immune response in durable clinical benefit (DCB) and non-DCB patients. Pre-treatment and early-on-treatment tumor transcriptomes were examined, and gene expression profiles, immunograms, and T cell receptor (TCR) repertoire were analyzed. DCB was observed in 15.8% of patients, with comparable secondary endpoints (ORR; objective response rate, OS; overall survival, PFS; progression-free survival) to previous trials. The immunograms of individual subjects displayed no significant changes before or early in the treatment, except for the regulatory T cell (Treg) score. Moreover, there were no consistent alterations observed among cases experiencing DCB. The intratumoral immune response was suppressed by previous treatments in most third- or later-line nivolumab recipients. TCR repertoire analysis revealed newly emerged clonotypes in early-on-treatment tumors, but clonal replacement did not impact efficacy. High T cell/Treg ratios and a low UV-radiation-response gene signature were linked to DCB and treatment response. This study emphasizes the tumor immune response’s importance in nivolumab efficacy for gastric cancer. High T cell/Treg ratios and specific gene expression signatures show promise as potential biomarkers for treatment response. The tumor-infiltrating immune response was compromised by prior treatments in third-line therapy, implying that, to enhance immunotherapeutic outcomes, commencing treatment at an earlier stage might be preferable. Larger cohort validation is crucial to optimize immune-checkpoint inhibitors in gastric cancer treatment.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details







1 Department of Gastrointestinal Surgery, Graduate School of Medicine, The University of Tokyo, Bunkyo-ku, Tokyo 113-8655, Japan;
2 Department of Gastrointestinal Surgery, Graduate School of Medicine, The University of Tokyo, Bunkyo-ku, Tokyo 113-8655, Japan;
3 Department of Immunotherapeutics, The University of Tokyo Hospital, Bunkyo-ku, Tokyo 113-8655, Japan;
4 cBioinformatics, Inc., Chiyoda-ku, Tokyo 101-0052, Japan;
5 Clinical Research Promotion Center, The University of Tokyo Hospital, Bunkyo-ku, Tokyo 113-8655, Japan;
6 Department of Gastrointestinal Surgery, Graduate School of Medicine, The University of Tokyo, Bunkyo-ku, Tokyo 113-8655, Japan;
7 Department of Gastrointestinal Surgery, Graduate School of Medicine, The University of Tokyo, Bunkyo-ku, Tokyo 113-8655, Japan;
8 Department of Immunotherapeutics, The University of Tokyo Hospital, Bunkyo-ku, Tokyo 113-8655, Japan;