Next Article in Journal
Impact of Center Volume on Cardiopulmonary and Mortality Outcomes after Immune-Checkpoint Inhibitors for Cancer: A Systematic Review and Meta-Analysis
Next Article in Special Issue
Predictive Signatures for Responses to Checkpoint Blockade in Small-Cell Lung Cancer in Second-Line Therapy Do Not Predict Responses in First-Line Patients
Previous Article in Journal
Evaluation of the Cost-Effectiveness of Evidence-Based Interventions to Increase Female Breast and Cervical Cancer Screens: A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contrasting Roles of Programmed Death-Ligand 1 Expression in Tumor and Stroma in Prognosis of Esophageal Squamous Cell Carcinoma

1
Department of Surgery, Hamamatsu University School of Medicine, 1-20-1 Handayama, Chuo-ku, Hamamatsu City 431-3192, Shizuoka, Japan
2
Department of Diagnostic Pathology; Hamamatsu University School of Medicine, 1-20-1 Handayama, Chuo-ku, Hamamatsu City 431-3192, Shizuoka, Japan
3
Department of Perioperative Functioning Care and Support, Hamamatsu University School of Medicine, 1-20-1 Handayama, Chuo-ku, Hamamatsu City 431-3192, Shizuoka, Japan
*
Author to whom correspondence should be addressed.
Cancers 2024, 16(6), 1135; https://doi.org/10.3390/cancers16061135
Submission received: 5 February 2024 / Revised: 6 March 2024 / Accepted: 11 March 2024 / Published: 13 March 2024
(This article belongs to the Special Issue Personalized Immunotherapy for Cancers)

Abstract

:

Simple Summary

Esophageal squamous cell carcinoma (ESCC) often recurs in advanced stages. With immune checkpoint inhibitors becoming a treatment option for recurrent ESCCs, assessing programmed death-ligand 1 (PD-L1) expression in tumor tissue has gained importance. However, challenges remain due to subjective interpretation methods, the unclear significance of staining intensity, and the contrasting roles of stromal and tumoral PD-L1 expression. Our study addresses these issues through a robust, machine learning-aided calculation of H-scores for PD-L1 expression, conducting separate analyses of tumoral and stromal areas in 194 cases of surgically resected ESCCs. Survival analysis revealed that high tumoral PD-L1 expression is an independent factor associated with prolonged survival. In contrast, high stromal PD-L1 expression correlated with less advanced pathological stages and a prolonged response to cytotoxic chemotherapy. This study highlights distinct PD-L1 roles in ESCC’s tumoral and stromal areas, underlining the importance of understanding these mechanisms in the tumor immune microenvironment to develop better treatments and improve patient outcomes.

Abstract

The assessment of programmed death-ligand 1 (PD-L1) expression in esophageal squamous cell carcinoma (ESCC) has become increasingly important with the rise of immune checkpoint inhibitors (ICIs). However, challenges persist, including subjective interpretation and the unclear significance of staining intensity, as well as contrasting roles in tumoral and stromal regions. Our study enhances the understanding of PD-L1 in ESCCs by analyzing its expression in tumors and stroma with H-scores, highlighting its distinct clinicopathological impacts. In a retrospective cohort of 194 ESCC specimens from surgical resection, we quantified PD-L1 expression in tumoral and stromal compartments using H-scores, analyzing whole slide images with digital pathology analysis software. Kaplan–Meier analysis demonstrated that higher PD-L1 expression is significantly associated with improved postoperative overall survival (OS) and recurrence-free survival (RFS) in both tumoral and stromal areas. Multivariable analysis identified high tumoral PD-L1 expression as an independent prognostic factor for prolonged OS and RFS (HR = 0.47, p = 0.007; HR = 0.54, p = 0.022, respectively). In a separate analysis, high stromal PD-L1 expression was found to correlate with less advanced pathological stages and a prolonged response to cytotoxic chemotherapy, with no similar correlation found for ICI treatment response. This study reveals PD-L1’s contrasting role in the ESCC tumor immune microenvironment, impacting prognosis, tumor stage, and treatment response.

Graphical Abstract

1. Introduction

Esophageal squamous cell carcinoma (ESCC) is a major global health concern, ranking as the sixth leading cause of cancer-related deaths worldwide [1]. While adenocarcinoma is more prevalent in Western countries, ESCC is particularly common in Asian countries, including Japan and China, accounting for up to 90% of esophageal cancers in these regions [2]. The standard treatment for stage II and III ESCC varies by region, with neoadjuvant cytotoxic chemotherapy (CTx) followed by surgery being common in Japan, and chemoradiotherapy (CRT) followed by surgery being common in Western countries [3,4,5]. A comprehensive nationwide database in Japan, the National Clinical Database (NCD), highlighted a significant reduction in 90-day mortality rates for esophagectomies, dropping from 3.2% in 2011 to 1.5% in 2020. Concurrently, the use of minimally invasive esophagectomy techniques rose markedly from 31.0% to 71.3%, reflecting a trend toward minimally invasive surgery [6]. Efforts to enhance postoperative outcomes and quality of life continue, focusing on tailoring treatment strategies and reconstructive procedures to individual patient and oncological needs for improved mid- to long-term outcomes [7]. However, recurrence rates remain high for stage II and III ESCCs, posing a significant challenge [8]. Current standard treatments for recurrent ESCC are diverse and evolving. They include CTx, CRT, and local resection. Recently, immune checkpoint inhibitors (ICIs) have been widely accepted as a first-line treatment option [3,4,5]. Studies such as the ATTRACTION-3 trial have demonstrated that ICIs have shown promising results in previously treated advanced ESCC patients, irrespective of tumor programmed death-ligand 1 (PD-L1) expression [9]. On the other hand, the Keynote-181 study suggested that pembrolizumab may offer a potential benefit in patients with high PD-L1 expression [10]. This indicates that the relationship between PD-L1 expression and ICI treatment response could be more complex than has been previously understood. The assessment of PD-L1 expression using combined positive score (CPS), as detailed in ESCC diagnostic guidelines, presents certain challenges [11]. These include the potential for inconsistency in the interpretation of 1+ staining intensity among pathologists, the unclear clinical significance of intensities exceeding 1+, and the lack of clinical implications of differential expression in tumor cells versus the stromal component. Among the various techniques for quantification, the H-score is particularly effective due to its unique scoring system that emphasizes the presence of strongly stained cells. This method enriches the evaluation of biomarker expression by combining a semi-quantitative assessment of staining intensity with the proportion of positively stained cells. Specifically, it assigns a score from 0 (no staining) to 3 (strong staining) based on the intensity, which is then integrated with the percentage of positive cells to compute the H-score, ranging from 0 to 300 [12]. Recent technological advancements have deepened our understanding of the complexity and diversity of the immune context within the tumor microenvironment. It is now recognized that the immune landscape of tumors is not uniform; rather, it comprises different subclasses of immune environments that significantly influence tumor initiation, progression, and response to therapy. The tumor immune microenvironment (TIME) encompasses a dynamic interplay of various cellular and molecular components, including immune cells, cytokines, and other factors that can either promote or inhibit tumor growth and metastasis [13]. The differential expression of PD-L1 in the tumor and stroma is not well understood at the basic research level. It is generally understood that PD-L1 expression can be induced in the tumor microenvironment under certain conditions, such as in the presence of inflammatory factors like interferons (IFNs), tumor necrosis factor-alpha (TNF-α), cell growth factors, hypoxic conditions, and exosomes [14,15]. Diverse findings exist regarding the role of PD-L1 expression in tumoral and stromal regions, with no consensus and limited research specifically in ESCCs [16,17,18,19,20]. This study delves into the clinicopathological significance of PD-L1 expression in ESCC, aiming to comprehensively understand its role in predicting responses to treatment modalities, including ICIs and CTx. Motivated by the ambiguous nature of PD-L1 as a predictive marker for ICI efficacy, our research seeks to broaden the investigation into its impact on disease progression and prognosis. This involves a detailed evaluation of PD-L1’s staining intensity and its variable implications in tumoral versus stromal areas, aiming to add another layer to current evaluation methods, which may overlook critical aspects of stromal expression. By adopting a holistic approach, our goal is to uncover the multifaceted role of PD-L1 in reflecting the status of the TIME, and its viability as a biomarker for customizing treatment strategies in ESCC. This investigation aims to clarify PD-L1’s prognostic value and refine treatment paradigms, contributing to enhanced patient outcomes in the evolving landscape of ESCC treatment.

2. Materials and Methods

2.1. Patients and Samples

Between 23 January 2012, and 2 December 2020, our study enrolled 194 patients diagnosed with ESCC of stage pT1 or higher, and these cases were retrospectively reviewed. All of the patients in this study underwent surgical resection at Hamamatsu University School of Medicine. Patient characteristics are detailed in Table 1. The observation period for this study extended until June 2023, with a median follow-up duration of 60 months (IQR 42), excluding deceased cases.

2.2. Immunohistochemistry

Formalin-fixed paraffin-embedded blocks from the tumor center were selected based on the standard procedures of the Department of Diagnostic Pathology, focusing on the deepest tumoral regions, as advised by two pathologists. Sections from these blocks were prepared for immunohistochemical analysis. The sections underwent deparaffinization, the inhibition of endogenous peroxidase using 3% H2O2, and antigen retrieval using Dako Target Retrieval Solution, pH 9 (S236784-2). For the primary antibody, we used the research-grade rabbit monoclonal anti-PD-L1 antibody [28-8] (ab205921), which has been documented for use in various cancers, including esophageal carcinoma, in the literature [21,22]. The primary antibody was diluted at a 1:100 ratio and incubated overnight. The application of the subsequent secondary antibody and 3,3′-Diaminobenzidine (DAB) staining were performed using the EnVision Detection Systems (K500711-2, DAKO, Santa Clara, CA, USA). Counterstaining was achieved using Mayer’s hematoxylin solution. The slides were then dehydrated, cleared, and mounted with a cover slip. The pathologists involved in selection contributed to the confirmation of the stained samples.

2.3. Image Acquisition and Analysis

Whole slide image files were created by scanning the slides using NanoZoomer®-HT/NanoZoomer®2.0-HT (Hamamatsu Photonics, Hamamatsu, Japan) at 20× magnification in brightfield mode and saved as NDPi files. These files were analyzed using Qupath v0.4.3 [23]. Regions of interest (ROI) in the tumoral areas were established referencing the hematoxylin- and eosin-stained specimens that had been previously diagnosed by pathologists in the Department of Diagnostic Pathology. The positive cell selection feature was utilized to detect the sum of optical density (OD) for hematoxylin and DAB staining within the slides. The cell size was set within the range of 10–400 μm, with a cell expansion of 5μm. Stromal and tumoral regions were selected as annotations, and machine learning was applied to classify objects in the ROI as tumor or stroma based on their appearance to the researcher. Thresholds for 1+, 2+, and 3+ intensities were established by selecting the mean ± SD values of the compiled OD for stroma + tumor ROI (Figure S1A). Using these thresholds (0.100, 0.134, 0.168), H-scores for stroma, tumor, and stroma + tumor were calculated for the 194 slides (Figure S1B).

2.4. Statistical Analysis

To determine the most significant H-score cutoff affecting postoperative overall survival (OS) and postoperative recurrence-free survival (RFS), the minimum p-value method was employed with the logrank function from the “cutoff” package in R 4.3.0 (R Core Team, 2023. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL https://www.R-project.org/ (accessed on 9 July 2023) via RStudio Version 2023.06.1+524 (RStudio Team, 2023) [24]. Survival analysis was performed for OS and RFS using the Kaplan–Meier method and analyzed using the logrank test. Cox regression analysis was performed using various clinicopathological factors, including age, gender, differentiation status, pT, pN, lymphatic invasion, venous invasion, intramural metastasis status, and the H-score status of PD-L1 in the stroma, tumor, and stroma + tumor. In the context of this analysis for OS, the PD-L1 H-score cutoff values for stroma, tumor, and stroma + tumor were set based on the value yielding the minimum p-value specifically for OS, thereby categorizing patients into two distinct groups. Separately, a univariate analysis was performed, setting the significance level at p < 0.05. Variables of clinical significance identified from the univariate analysis were then subjected to a multivariate analysis. For the RFS analysis, an identical approach was adopted, with the only difference using H-score cutoff values that were set based on the minimum p-value for RFS. Survival analyses were performed using SPSS Version 29 (IBM Corp., Armonk, NY, USA). Correlation analysis was conducted to assess the relationship between clinicopathological factors and PD-L1 expression status in the stroma, tumor, and stroma + tumor regions. The cutoff values for categorizing PD-L1 expression as high or low were determined based on the H-scores that yielded the minimum p-value in OS analysis. Chi-square test analysis was performed using SPSS Version 29 (IBM Corp., Armonk, NY, USA). Statistical significance was set at p < 0.05. The H-scores for PD-L1 expression in stroma, tumor, and stroma + tumor were compared using the Kruskal–Wallis and Dunn’s multiple comparison test. For the analysis of treatment duration and PD-L1 H-score, the Mann–Whitney U test was utilized. These analyses were conducted using GraphPad Prism (version 9.0, GraphPad Software, LLC., San Diego, CA, USA).

3. Results

3.1. PD-L1 Expression in Stroma, Tumor, and Stroma + Tumor

To gain an overview of the distribution of PD-L1 expression in tumor and stroma, we analyzed the median H-scores in our cohort of 194 cases. The results revealed that median H-scores for PD-L1 expression in the tumor were significantly higher than those in the stroma (median 155.0 vs. 97.4; p < 0.001, Figure 1A). Representative images of PD-L1 staining using the quartile H-scores of the tumor and stroma as cutoff values are shown (Figure 1B).

3.2. Survival Analysis Based on PD-L1 H-Score

To evaluate the clinical significance of PD-L1 expression in tumor, stroma, and stroma + tumor on postoperative prognosis, we performed a survival analysis using PD-L1 H-Scores. Kaplan–Meier survival curves generated using quartile H-scores as cutoff values for stroma, tumor, and stroma + tumor showed a trend between higher H-score and improved prognosis (Figure S1C–H). Next, we attempt to determine the optimal H-Score cutoff values that will give the most clinical impact in the analyzed cases. The minimum p-value method for the logrank test was used to identify the optimal cutoff. The identified optimal cutoff values for OS for stroma, tumor, and stroma + tumor were 114.86 (p = 0.000020), 159.83 (p = 0.002498), and 165.22 (p = 0.003303), respectively (Figure S2A–C). The optimal cutoff values for RFS were 114.86 (p = 0.000080), 164.22 (p = 0.000556), and 185.45 (p = 0.004172), respectively (Figure S2D–F). Kaplan–Meier survival curves were generated using the optimal H-score cutoff values. For OS, groups with high H-scores in stroma, tumor, and stroma + tumor showed significantly longer survival (p < 0.001, p = 0.005, and p = 0.007 respectively; Figure 1C–E). Similarly, for RFS, longer survival durations were observed in groups with high H-scores in stroma, tumor, and stroma + tumor (p < 0.001, p = 0.008, and p = 0.009 respectively; Figure 1F–H).

3.3. Multivariate Analysis of Clinicopathological Factors

To determine the significance of PD-L1 expression alongside other clinicopathological factors in influencing OS in ESCC, we performed a Cox regression analysis. In the univariate analysis to identify clinicopathological factors affecting OS, preoperative treatment, pT2 or higher, pN positive, lymphatic invasion positive, venous invasion positive, intramural metastasis positive, and low PD-L1 expression in stroma, tumor, and stroma + tumor were identified as significant factors for shorter OS (Table 2). In the multivariate analysis, pN positive (HR = 2.70, 95% CI 1.12–6.52, p = 0.027) and low PD-L1 expression in the tumor (HR = 0.47, 95% CI 0.27–0.83, p = 0.010) remained significant to predict poor OS. Next, Cox regression analysis was performed to identify clinicopathological factors affecting RFS (Table 3). In the univariate analysis, the inclusion of a poorly differentiated tumor, preoperative treatment, pT2 or higher, pN positive, lymphatic invasion positive, venous invasion positive, intramural metastasis positive, low PD-L1 expression in stroma, tumor, and stroma + tumor were identified as significant factors for shorter RFS. In the multivariate analysis, pN positive (HR = 2.56, 95% C.I. 1.17–5.58, p = 0.018) and low PD-L1 expression in the tumor (HR = 0.54, 95%C.I. 0.32–0.91, p = 0.022) remained significant to predict poor RFS.

3.4. Correlation Analysis of PD-L1 Expression and Clinicopathological Factors

Given that only tumoral PD-L1 expression was identified as a significant factor in the multivariate analysis for survival, we investigated whether specific clinicopathological factors differentially influence PD-L1 expression in the stroma, tumor, and stroma + tumor compartments. Clinicopathological factors correlated with PD-L1 expression in stroma, tumor, and stroma + tumor were analyzed (Table 4). Higher pT Stage, higher pN Stage, positive lymphatic invasion, and positive venous invasion were associated with lower stroma PD-L1(p < 0.001, respectively). However, no such correlation was observed for the tumor. For stroma + tumor, significant correlations were observed with pN stage (p = 0.041) and venous invasion (p = 0.002).

3.5. PD-L1 Expression and Its Association with Recurrence Patterns and Metastatic Sites

To further investigate the role of PD-L1 expression in disease progression of ESCC, we conducted an analysis focused on recurrence patterns, total recurrence types, and total metastatic sites in relation to PD-L1 expression in the stroma, tumor, and stroma + tumor regions. This analysis was specifically performed in cases of recurrence (Table 5). Our findings indicate that there were no statistically significant differences across these variables in relation to PD-L1 expression.

3.6. Comparison of Survival Times Based on PD-L1 Expression Patterns in Stroma and Tumor

Given that our previous analyses suggested the presence of a PD-L1 expression pattern with distinct prognostic implications, we sought to assess the impact of various combinations of PD-L1 expression in stromal and tumoral regions. These were characterized as either high or low in both compartments, leading us to conduct a comparative analysis of OS and RFS. The minimum p-value cutoff for OS and RFS was utilized to categorize PD-L1 expression in stroma and tumor into four bins, resulting in distinct patterns of expression (Figure 2A). When we generated Kaplan-Meier survival curves for these four groups, we found a significant difference in both OS and RFS between the group with both low stromal and tumoral PD-L1 expression and the group with high expression in both regions (p < 0.001, respectively Figure 2B,C).

3.7. PD-L1 Expression and Post-Recurrence Treatment Efficacy: CTx and Immune Checkpoint Inhibitors

Our observations of varying recurrence rates linked to PD-L1 expression levels led us to examine their potential impact on post-recurrence treatment outcomes, focusing on the duration of CTx and responses to ICIs without progression. Among the 85 cases that recurred, 70 received CTx for recurrence treatment. Our analysis indicated a significantly higher stromal H-score in cases with a CTx treatment period of 8 months or more (median 53.8 vs. 98.3, p = 0.028, Figure 2D). However, no significant difference in PD-L1 expression was observed for the tumor (median 133,8 vs. 154.6, p = 0.200, Figure 2E) and stroma + tumor (median 113.7 vs. 122.3, p = 0.355, Figure 2F). Furthermore, among the recurrent cases, 24 received ICIs. Analysis showed no significant difference in stroma, tumor, or stroma + tumor PD-L1 expression between cases that were able to continue ICIs without progression for 90 days or more and those that did not (median 102.4 vs. 106.6, p = 0.450; median 148.5 vs. 158.7, p = 0.568; median 108.1 vs. 122.3, p = 0.287, respectively; Figure 2G–I).

4. Discussion

This study quantified PD-L1 expression in 194 surgical specimens of ESCC using digital pathology, classifying the whole tumor within the slide into tumoral and stromal compartments enhanced using AI. By employing Kaplan–Meier survival analysis and categorizing H-scores into quartiles, we found that higher scores were associated with better OS and RFS. Utilizing the minimum p value method for logrank test, we determined statistically optimal H-score cutoff values for OS and RFS, dividing the cohort into groups to assess the impact on survival. In the multivariate analysis regarding survival, tumoral PD-L1 expression emerged as an independent predictor of prognosis for both OS and RFS. Our correlation analysis with clinicopathological variables revealed that high stromal PD-L1 expression was strongly linked to less advanced pathological stages. Furthermore, when comparing H-scores with the duration of CTx and ICI treatment, we discovered that cases treated for more than 8 months exhibited higher stromal H-scores. However, no difference in H-scores between compartments was observed in cases treated with ICI for more than 90 days versus those treated for shorter periods. Our approach to evaluating PD-L1 expression in this study marks a significant shift from traditional methods by leveraging digital pathology for quantitative analysis of PD-L1 in ESCC, thus deepening our grasp of its clinical-pathological significance. Enhanced by AI-based tissue classification, this method efficiently differentiates between tumoral and stromal regions, spotlighting the underexplored area of stromal PD-L1 expression and its relative clinical-pathological importance in comparison to tumoral PD-L1 expression in ESCC. This pioneering approach not only improves the precision of PD-L1 measurement but also expands our insight into its role within the TIME, thereby offering a more nuanced understanding of the interplay between tumoral and stromal components in ESCC. By providing a detailed analysis of PD-L1 expression in both tumoral and stromal compartments, our study contributes to resolving the conflicting perspectives in the broader context of ESCC research, where previous studies have yielded mixed results regarding the prognostic significance of PD-L1 expression. Certain investigations indicated that higher PD-L1 expression in tumor cells was linked to shorter disease-free survival and OS [25]. In contrast, there are studies that found higher PD-L1 expression to be associated with better prognosis [26,27]. The distinguishing of PD-L1 expression in tumoral and stromal compartments is seldom addressed. Some studies have evaluated PD-L1 expression in both compartments but not in ESCC [16,17,18]. Previous research on stromal PD-L1 expression has primarily focused on immune cells within the tumor tissue. Some studies have associated high stromal PD-L1 expression with worse prognosis [17,28] while others have linked it to better prognosis [16,19,20]. In our study, we could not differentiate the stromal expression of PD-L1 in immune cells from other stromal cells upon the classification procedure due to the complexity and heterogeneity of stromal components. It is anticipated that future advancements in image analysis and machine learning will enable the differentiation of subtle features in the TIME, thereby facilitating a better understanding of PD-L1 dynamics in ESCC. Although, this brings us to the intriguing possibility that TIME components such as cancer associated fibroblasts (CAFs), could also have an influence on PD-L1 expression. With regard to PD-L1 expression on CAFs, there are only a handful of studies, and their findings remain inconclusive. One study in nonmetastatic non-small cell lung cancer suggested that PD-L1 expression on CAFs was associated with a better prognosis [20], whereas another study in colorectal cancer reported that CAFs upregulate PD-L1 expression in tumor cells, contributing to worse prognosis [29]. The molecular mechanisms underlying PD-L1 expression in CAFs remain largely unexplored to date. Thus, further investigations are necessary to elucidate the role of CAFs in PD-L1 expression and to determine its prognostic significance. Our study revealed that elevated stromal PD-L1 expression was linked to a longer treatment period (8 months or more) in patients undergoing CTx, implying that stromal PD-L1 expression might serve as a predictive biomarker for a favorable response to CTx. This echoes a study on triple-negative breast cancer, which demonstrated the potential of PD-L1 as a response biomarker for CTx [30]. Our research offers an innovative perspective by examining PD-L1’s influence on CTx effectiveness in ESCC, moving beyond the conventional focus on its correlation with ICI responsiveness. By understanding PD-L1 status from resected specimens, we can predict the responsiveness to CTx, which plays a significant role from the first line recurrence treatment, paving the way for more precise and personalized medicine. As discussed so far, the role of PD-L1 in either promoting progression or exerting a suppressive effect on tumors remains controversial in the literature. Our findings, which indicate that high tumoral PD-L1 expression is associated with better prognosis, align with some previous studies in ESCC. Additionally, our report that higher stromal PD-L1 expression is linked to less advanced disease progression and our discovery that PD-L1 expression influences responsiveness to CTx in ESCC both represent first time in the context of ESCC. However, contrary to the KEYNOTE181 trials, our study did not observe a correlation between PD-L1 expression and an improved response to immunotherapy. This discrepancy could potentially be attributed to the smaller sample size in our study compared to these large clinical trials. It is also possible that other factors, such as variations in treatment protocols or patient characteristics, could have influenced this outcome. Thus, while our findings contribute to the evolving understanding of PD-L1’s role in immunotherapy, more comprehensive studies are needed to definitively ascertain this relationship. PD-L1 expression in cells are known to be influenced by factors like inflammatory cytokines, hypoxic conditions, and exosomes [14,15]. TIME components and the altered cytokine milieu could further contribute to this complexity [31,32]. These factors might be involved in mechanisms that render tumors more responsive to CTx and potentially immunotherapy. Further research is necessary to unravel these complex dynamics and could lead to strategies for effectively targeting stromal PD-L1 to enhance the response to therapy. Building upon our foundational work in delineating PD-L1 expression dynamics within the stromal and tumor cells of ESCC, future research will delve further into the TIME. This exploration includes assessing T-cell and tumor-associated macrophage markers, along with inflammatory cytokines that potentially exacerbate ESCC prognosis through mechanisms related to cachexia. Each step forward in deciphering the intricate interactions within the TIME not only aids in predicting treatment responses but also strives to identify novel therapeutic targets, moving us closer to improving patient outcomes in ESCC. There are several methodological limitations in this study. Future directions will focus on expanding the study scale and diversity by undertaking a multi-institutional approach to compile a larger, more diverse dataset. This expansion should enhance the persuasiveness, generalizability, and credibility of the results by incorporating varied patient demographics and treatment histories. Insights into PD-L1 expression in patients treated with ICIs as a first-line treatment can offer new perspectives on mechanisms of resistance and significantly contribute to our understanding of PD-L1′s role in ESCCs. Integration of multi-omics data from public databases and existing resources for genomic, transcriptomic, and proteomic analyses can deepen our understanding of PD-L1 dynamics and interactions within TIME. Longitudinal analysis will allow us to monitor changes in PD-L1 levels over time, for example, by using endoscopic biopsy specimens before and after treatment could provide insights into PD-L1′s prognostic significance, inform treatment planning, and possibly reveal aspects of treatment responsiveness. Regarding digital pathology analysis, expanding the dataset to include slides evaluated with CPS and TPS will enable a comparison between traditional and digital pathology methods, potentially validating the clinical utility of digital pathology. To sum, this research advances the field by elucidating the complex role of PD-L1 in both tumoral and stromal compartments of ESCC, employing quantification methods to capture staining intensity nuances not considered by CPS and TPS. The digital pathology-based H-score quantification methods offer a more nuanced view of PD-L1’s impact on cancer progression and therapeutic avenues. This work underscores the necessity of a holistic view of PD-L1 expression within the TIME, contributing significantly to the understanding of its role in ESCC and guiding future therapeutic strategies.

5. Conclusions

Our study demonstrates that high tumoral PD-L1 expression is an independent prognostic factor for longer survival. Conversely, high stromal PD-L1 expression correlates with less advanced pT and pN stages, negative lymphatic and venous invasion, and prolonged duration of CTx treatment. Contrary to expectations, PD-L1 expression in the stromal, tumor, and stroma + tumor compartments did not affect the duration of ICI treatment.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers16061135/s1. Figure S1: PD-L1 staining intensity thresholds determination and H-score calculation workflow, followed by survival analysis based on quartile H-scores across stroma, tumor, and stroma + tumor compartments in ESCC. Figure S2: Identification of optimal cutoff values for OS and RFS using the minimum p-value method for the logrank test in ESCC.

Author Contributions

Conceptualization T.M. (Tomohiro Murakami); methodology, T.M. (Tomohiro Murakami) and S.F.; validation, T.M. (Tomohiro Murakami), E.B., M.F. and S.B.; formal analysis, T.M. (Tomohiro Murakami) and S.F.; investigation, T.M. (Tomohiro Murakami); resources, H.T. and S.B.; data curation, Y.S., K.S., R.H., T.M. (Tomohiro Murakami) and T.M. (Tomohiro Matsumoto); writing—original draft preparation, T.M. (Tomohiro Murakami) and E.B.; writing—review and editing, T.M. (Tomohiro Murakami), E.B., Y.M., H.K. and H.T.; visualization, T.M. (Tomohiro Murakami) and S.F.; supervision, E.B., H.K., Y.H. and H.T.; project administration, T.M. (Tomohiro Murakami) and H.K.; funding acquisition, T.M. (Tomohiro Murakami) All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI, Grant Number 22K16485.

Institutional Review Board Statement

All procedures were conducted in accordance with institutional and national standards on human experimentation, as confirmed by the ethics committee of Hamamatsu University School of Medicine (IRB No. 17-323; approved 1 March 2018), and with the Declaration of Helsinki of 1964 and its later versions.

Informed Consent Statement

The IRB waived the requirement for written patient consent for the use of clinicopathological data, and all patients agreed to participate through an opt-out approach.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, [T.M.], upon reasonable request.

Acknowledgments

The authors would like to express their sincere gratitude to all the patients and medical staff at the institutions who contributed to this study. We also extend our thanks to Yuki Kurita for his technical expertise and support in the implementation of immunostaining. Additionally, the graphical abstract was created using Canva (https://www.canva.com/, accessed on January 13 2024).

Conflicts of Interest

The authors declare no conflicts of interest.

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, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Arnold, M.; Abnet, C.C.; Neale, R.E.; Vignat, J.; Giovannucci, E.L.; McGlynn, K.A.; Bray, F. Global Burden of 5 Major Types of Gastrointestinal Cancer. Gastroenterology 2020, 159, 335–349.e15. [Google Scholar] [CrossRef] [PubMed]
  3. Kitagawa, Y.; Ishihara, R.; Ishikawa, H.; Ito, Y.; Oyama, T.; Oyama, T.; Kato, K.; Kato, H.; Kawakubo, H.; Kawachi, H.; et al. Esophageal cancer practice guidelines 2022 edited by the Japan esophageal society: Part 1. Esophagus 2023, 20, 343–372. [Google Scholar] [CrossRef] [PubMed]
  4. Kitagawa, Y.; Ishihara, R.; Ishikawa, H.; Ito, Y.; Oyama, T.; Oyama, T.; Kato, K.; Kato, H.; Kawakubo, H.; Kawachi, H.; et al. Esophageal cancer practice guidelines 2022 edited by the Japan Esophageal Society: Part 2. Esophagus 2023, 20, 373–389. [Google Scholar] [CrossRef]
  5. Obermannová, R.; Alsina, M.; Cervantes, A.; Leong, T.; Lordick, F.; Nilsson, M.; van Grieken, N.C.T.; Vogel, A.; Smyth, E.C. Oesophageal cancer: ESMO Clinical Practice Guideline for diagnosis, treatment and follow-up. Ann. Oncol. 2022, 33, 992–1004. [Google Scholar] [CrossRef]
  6. Kajiwara, Y.; Takahashi, A.; Ueno, H.; Kakeji, Y.; Hasegawa, H.; Eguchi, S.; Goi, T.; Saiura, A.; Sasaki, A.; Takiguchi, S.; et al. Annual report on National Clinical Database 2020 for gastroenterological surgery in Japan. Ann. Gastroenterol. Surg. 2023, 7, 367–406. [Google Scholar] [CrossRef]
  7. Booka, E.; Takeuchi, H.; Morita, Y.; Hiramatsu, Y.; Kikuchi, H. What is the best reconstruction procedure after esophagectomy? A meta-analysis comparing posterior mediastinal and retrosternal approaches. Ann. Gastroenterol. Surg. 2023, 7, 553–564. [Google Scholar] [CrossRef]
  8. Zhou, S.; Liu, S.; Zhang, L.; Guo, S.; Shen, J.; Li, Q.; Yang, H.; Feng, Y.; Liu, M.; Lin, S.H.; et al. Recurrence Risk Based on Pathologic Stage after Neoadjuvant Chemoradiotherapy in Esophageal Squamous Cell Carcinoma: Implications for Risk-Based Postoperative Surveillance Strategies. Ann. Surg. Oncol. 2018, 25, 3639–3646. [Google Scholar] [CrossRef]
  9. Kato, K.; Cho, B.C.; Takahashi, M.; Okada, M.; Lin, C.Y.; Chin, K.; Kadowaki, S.; Ahn, M.J.; Hamamoto, Y.; Doki, Y.; et al. Nivolumab versus chemotherapy in patients with advanced oesophageal squamous cell carcinoma refractory or intolerant to previous chemotherapy (ATTRACTION-3): A multicentre, randomised, open-label, phase 3 trial. Lancet Oncol. 2019, 20, 1506–1517. [Google Scholar] [CrossRef]
  10. Kojima, T.; Shah, M.A.; Muro, K.; Francois, E.; Adenis, A.; Hsu, C.H.; Doi, T.; Moriwaki, T.; Kim, S.B.; Lee, S.H.; et al. Randomized Phase III KEYNOTE-181 Study of Pembrolizumab Versus Chemotherapy in Advanced Esophageal Cancer. J. Clin. Oncol. 2020, 38, 4138–4148. [Google Scholar] [CrossRef]
  11. PD-L1 IHC 22C3 pharmDx Interpretation Manual—Esophageal Cancer. Available online: https://www.agilent.com/cs/library/usermanuals/public/29439-d67239-pd-l1-ihc22c3-ec-kn590-int-man-en.pdf (accessed on 25 July 2023).
  12. McCarty, K.S., Jr.; Szabo, E.; Flowers, J.L.; Cox, E.B.; Leight, G.S.; Miller, L.; Konrath, J.; Soper, J.T.; Budwit, D.A.; Creasman, W.T.; et al. Use of a monoclonal anti-estrogen receptor antibody in the immunohistochemical evaluation of human tumors. Cancer Res. 1986, 46, 4244s–4248s. [Google Scholar] [PubMed]
  13. Binnewies, M.; Roberts, E.W.; Kersten, K.; Chan, V.; Fearon, D.F.; Merad, M.; Coussens, L.M.; Gabrilovich, D.I.; Ostrand-Rosenberg, S.; Hedrick, C.C.; et al. Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat. Med. 2018, 24, 541–550. [Google Scholar] [CrossRef]
  14. Han, Y.; Liu, D.; Li, L. PD-1/PD-L1 pathway: Current researches in cancer. Am. J. Cancer Res. 2020, 10, 727–742. [Google Scholar] [PubMed]
  15. Jiang, X.; Wang, J.; Deng, X.; Xiong, F.; Ge, J.; Xiang, B.; Wu, X.; Ma, J.; Zhou, M.; Li, X.; et al. Role of the tumor microenvironment in PD-L1/PD-1-mediated tumor immune escape. Mol. Cancer 2019, 18, 10. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, C.Q.; Xu, J.; Zhou, Z.G.; Jin, L.L.; Yu, X.J.; Xiao, G.; Lin, J.; Zhuang, S.M.; Zhang, Y.J.; Zheng, L. Expression patterns of programmed death ligand 1 correlate with different microenvironments and patient prognosis in hepatocellular carcinoma. Br. J. Cancer 2018, 119, 80–88. [Google Scholar] [CrossRef]
  17. Thompson, E.D.; Zahurak, M.; Murphy, A.; Cornish, T.; Cuka, N.; Abdelfatah, E.; Yang, S.; Duncan, M.; Ahuja, N.; Taube, J.M.; et al. Patterns of PD-L1 expression and CD8 T cell infiltration in gastric adenocarcinomas and associated immune stroma. Gut 2017, 66, 794–801. [Google Scholar] [CrossRef]
  18. Zhai, Q.; Fan, J.; Lin, Q.; Liu, X.; Li, J.; Hong, R.; Wang, S. Tumor stromal type is associated with stromal PD-L1 expression and predicts outcomes in breast cancer. PLoS ONE 2019, 14, e0223325. [Google Scholar] [CrossRef]
  19. Sumitomo, R.; Hirai, T.; Fujita, M.; Murakami, H.; Otake, Y.; Huang, C.L. PD-L1 expression on tumor-infiltrating immune cells is highly associated with M2 TAM and aggressive malignant potential in patients with resected non-small cell lung cancer. Lung Cancer 2019, 136, 136–144. [Google Scholar] [CrossRef]
  20. Teramoto, K.; Igarashi, T.; Kataoka, Y.; Ishida, M.; Hanaoka, J.; Sumimoto, H.; Daigo, Y. Clinical significance of PD-L1-positive cancer-associated fibroblasts in pN0M0 non-small cell lung cancer. Lung Cancer 2019, 137, 56–63. [Google Scholar] [CrossRef]
  21. Wu, X.; Ke, X.; Ni, Y.; Kuang, L.; Zhang, F.; Lin, Y.; Lin, W.; Xiong, X.; Huang, H.; Lin, X.; et al. Tumor-Infiltrating Immune Cells and PD-L1 as Prognostic Biomarkers in Primary Esophageal Small Cell Carcinoma. J. Immunol. Res. 2020, 2020, 8884683. [Google Scholar] [CrossRef]
  22. Wang, P.; Chen, Y.; Long, Q.; Li, Q.; Tian, J.; Liu, T.; Wu, Y.; Ding, Z. Increased coexpression of PD-L1 and TIM3/TIGIT is associated with poor overall survival of patients with esophageal squamous cell carcinoma. J. Immunother. Cancer 2021, 9, PMC8504357. [Google Scholar] [CrossRef] [PubMed]
  23. Bankhead, P.; Loughrey, M.B.; Fernández, J.A.; Dombrowski, Y.; McArt, D.G.; Dunne, P.D.; McQuaid, S.; Gray, R.T.; Murray, L.J.; Coleman, H.G.; et al. QuPath: Open source software for digital pathology image analysis. Sci. Rep. 2017, 7, 16878. [Google Scholar] [CrossRef] [PubMed]
  24. Furuhashi, S.; Morita, Y.; Matsumoto, A.; Ida, S.; Muraki, R.; Kitajima, R.; Takeda, M.; Kikuchi, H.; Hiramatsu, Y.; Takeuchi, H. Tenascin C in pancreatic cancer-associated fibroblasts enhances epithelial mesenchymal transition and is associated with resistance to immune checkpoint inhibitor. Am. J. Cancer Res. 2023, 13, 5641–5655. [Google Scholar] [PubMed]
  25. Okadome, K.; Baba, Y.; Nomoto, D.; Yagi, T.; Kalikawe, R.; Harada, K.; Hiyoshi, Y.; Nagai, Y.; Ishimoto, T.; Iwatsuki, M.; et al. Prognostic and clinical impact of PD-L2 and PD-L1 expression in a cohort of 437 oesophageal cancers. Br. J. Cancer 2020, 122, 1535–1543. [Google Scholar] [CrossRef]
  26. Jiang, C.; Zhu, Y.; Tang, S.; Zhang, G.; Lin, Q.; Xu, Y.; Shang, J. High PD-L1 expression is associated with a favorable prognosis in patients with esophageal squamous cell carcinoma undergoing postoperative adjuvant radiotherapy. Oncol. Lett. 2019, 17, 1626–1634. [Google Scholar] [CrossRef] [PubMed]
  27. Guo, W.; Zhang, F.; Shao, F.; Wang, P.; Li, Z.; Yang, X.; He, Z.; Shi, S.; Gao, Y.; He, J. PD-L1 expression on tumor cells associated with favorable prognosis in surgically resected esophageal squamous cell carcinoma. Hum. Pathol. 2019, 84, 291–298. [Google Scholar] [CrossRef] [PubMed]
  28. Takahashi, H.; Watanabe, H.; Hashimura, M.; Matsumoto, T.; Yokoi, A.; Nakagawa, M.; Ishibashi, Y.; Ito, T.; Ohhigata, K.; Saegusa, M. A combination of stromal PD-L1 and tumoral nuclear β-catenin expression as an indicator of colorectal carcinoma progression and resistance to chemoradiotherapy in locally advanced rectal carcinoma. J. Pathol. Clin. Res. 2022, 8, 458–469. [Google Scholar] [CrossRef]
  29. Gao, Y.; Sun, Z.; Gu, J.; Li, Z.; Xu, X.; Xue, C.; Li, X.; Zhao, L.; Zhou, J.; Bai, C.; et al. Cancer-Associated Fibroblasts Promote the Upregulation of PD-L1 Expression Through Akt Phosphorylation in Colorectal Cancer. Front. Oncol. 2021, 11, 748465. [Google Scholar] [CrossRef]
  30. Yoshikawa, K.; Ishida, M.; Yanai, H.; Tsuta, K.; Sekimoto, M.; Sugie, T. Prognostic significance of PD-L1-positive cancer-associated fibroblasts in patients with triple-negative breast cancer. BMC Cancer 2021, 21, 239. [Google Scholar] [CrossRef]
  31. Farlow, J.L.; Brenner, J.C.; Lei, Y.L.; Chinn, S.B. Immune deserts in head and neck squamous cell carcinoma: A review of challenges and opportunities for modulating the tumor immune microenvironment. Oral. Oncol. 2021, 120, 105420. [Google Scholar] [CrossRef]
  32. Chen, D.S.; Mellman, I. Elements of cancer immunity and the cancer–immune set point. Nature 2017, 541, 321–330. [Google Scholar] [CrossRef] [PubMed]
Figure 1. PD-L1 expression and its impact on survival in ESCC patients. (A) Scatter bar graph showing the median H-scores for PD-L1 expression in the stroma (S), tumor (T), and stroma + tumor (S + T) in a cohort of 194 ESCC cases. Significant differences among the groups were observed, as determined by the Kruskal–Wallis and Dunn’s multiple comparison test (*** p < 0.001). (B) Representative PD-L1 staining images selected based on quartile H-scores for tumor and stroma. Scale bar: 50 μm. (CH) Kaplan–Meier survival curves illustrating OS and RFS stratified by optimal PD-L1 H-score cutoffs in stroma, tumor, and stroma + tumor. Survival differences between groups were analyzed using the logrank test.
Figure 1. PD-L1 expression and its impact on survival in ESCC patients. (A) Scatter bar graph showing the median H-scores for PD-L1 expression in the stroma (S), tumor (T), and stroma + tumor (S + T) in a cohort of 194 ESCC cases. Significant differences among the groups were observed, as determined by the Kruskal–Wallis and Dunn’s multiple comparison test (*** p < 0.001). (B) Representative PD-L1 staining images selected based on quartile H-scores for tumor and stroma. Scale bar: 50 μm. (CH) Kaplan–Meier survival curves illustrating OS and RFS stratified by optimal PD-L1 H-score cutoffs in stroma, tumor, and stroma + tumor. Survival differences between groups were analyzed using the logrank test.
Cancers 16 01135 g001
Figure 2. Impact of PD-L1 expression patterns in stroma and tumor on survival times and post-recurrence treatment efficacy. (A) Representative image showing four distinct PD-L1 expression patterns in stroma and tumor, categorized based on the minimum p-value cutoffs for OS. Scale bar: 50 μm. (B,C) Kaplan–Meier survival curves comparing OS and RFS among four groups with distinct PD-L1 expression patterns in stroma and tumor, highlighting the impact of PD-L1 expression on survival times. Groups were categorized using minimum p-value cutoffs for OS and RFS, respectively. The patterns are abbreviated as follows: stroma low/tumor low (S-L/T-L), stroma low/tumor high (S-L/T-H), stroma high/tumor low (S-H/T-L), and stroma high/tumor high (S-H/T-H). Survival differences between groups were analyzed using the logrank test. (DI) Scatter bar graphs illustrating H-scores for PD-L1 expression in the stroma, tumor, and stroma + tumor compartments, focusing on the correlation with post-recurrence treatment efficacy in CTx and ICIs (* p < 0.05, NS -not significant). (DF) Graphs compare cases that experienced recurrence and received CTx for 8 months or longer versus those treated for less than 8 months. (GI) Graphs compare cases that experienced recurrence and received ICIs for 90 days or longer versus those treated for less than 90 days. Median and interquartile ranges are indicated on the graphs. The Mann–Whitney U test was used to assess statistical significance for both comparisons.
Figure 2. Impact of PD-L1 expression patterns in stroma and tumor on survival times and post-recurrence treatment efficacy. (A) Representative image showing four distinct PD-L1 expression patterns in stroma and tumor, categorized based on the minimum p-value cutoffs for OS. Scale bar: 50 μm. (B,C) Kaplan–Meier survival curves comparing OS and RFS among four groups with distinct PD-L1 expression patterns in stroma and tumor, highlighting the impact of PD-L1 expression on survival times. Groups were categorized using minimum p-value cutoffs for OS and RFS, respectively. The patterns are abbreviated as follows: stroma low/tumor low (S-L/T-L), stroma low/tumor high (S-L/T-H), stroma high/tumor low (S-H/T-L), and stroma high/tumor high (S-H/T-H). Survival differences between groups were analyzed using the logrank test. (DI) Scatter bar graphs illustrating H-scores for PD-L1 expression in the stroma, tumor, and stroma + tumor compartments, focusing on the correlation with post-recurrence treatment efficacy in CTx and ICIs (* p < 0.05, NS -not significant). (DF) Graphs compare cases that experienced recurrence and received CTx for 8 months or longer versus those treated for less than 8 months. (GI) Graphs compare cases that experienced recurrence and received ICIs for 90 days or longer versus those treated for less than 90 days. Median and interquartile ranges are indicated on the graphs. The Mann–Whitney U test was used to assess statistical significance for both comparisons.
Cancers 16 01135 g002
Table 1. Detailed characteristics of the 194 patients diagnosed with ESCC of stage pT1 or higher analyzed in this study.
Table 1. Detailed characteristics of the 194 patients diagnosed with ESCC of stage pT1 or higher analyzed in this study.
VariablesN(%)
Age
  under 70 y.o114(59%)
  over and 70 y.o80(41%)
Sex
  male169(88%)
  female25(13%)
Tumor Location
  Ce4(2%)
  Ut22(11%)
  Mt95(49%)
  Lt72(37%)
  Ae1(1%)
Differentiation
  differentiated140(73%)
  include poorly differentiated41(21%)
  NA13(7%)
Preoperative treatment
  None97(50%)
  Chemotherapy87(45%)
  Chemotherapy+Radiation10(5%)
pT
  T194(49%)
  T226(13%)
  T374(38%)
pN
  N071(37%)
  N165(34%)
  N234(18%)
  N324(12%)
Lymphatic invasion
  negative70(36%)
  positive124(64%)
Venous invasion
  negative55(28%)
  positive139(72%)
Intramural metastasis
  negative176(91%)
  positive18(9%)
pStage *
  I55(28%)
  II43(22%)
  III59(31%)
  IVA37(19%)
Total recurrence cases85(44%)
  Lymphogenous metastasis
    negative16(19%)
    positive69(81%)
  Hematogenous metastasis
    negative33(39%)
    positive52(61%)
  Locoregional recurrence
    negative76(89%)
    positive9(11%)
  Pleural/peritoneal dissemination
    negative66(78%)
    positive19(22%)
* TNM-8—TNM Classification of Malignant Tumours—8th edition; Ce—Cervical part of the esophagus; Ut—Upper thoracic part of the esophagus; Mt—Middle thoracic part of the esophagus; Lt—Lower thoracic part of the esophagus; Ae—Abdominal part of the esophagus.
Table 2. Results of the Cox regression analysis identifying clinicopathological factors affecting OS.
Table 2. Results of the Cox regression analysis identifying clinicopathological factors affecting OS.
VariablesNUnivariableMultivariable
HR95%CIp-ValueHR95%CIp-Value
Age
  under 70 y.o1141.000.481.300.355
  over and 70 y.o800.79
Sex
  male 1691.000.311.470.321
  female250.67
Differentiation
  differentiated1401.000.952.790.077
  include poorly diff.411.63
Preoperative treatment
  none971.001.774.97<0.0011.000.953.080.075
  chemotherapy(+/-Radiation)972.971.71
pT *
  <T2941.002.547.82<0.0011.000.953.650.069
  T2 and over 1004.461.86
pN *
  negative 711.002.8012.28<0.0011.001.126.520.027
  positive1235.862.70
Lymphatic invasion
  negative 701.001.776.11<0.0011.000.472.270.943
  positive1243.281.03
Venous invasion
  negative 551.001.969.35<0.0011.000.613.890.357
  positive1394.231.55
Intramural metastasis
  negative 1761.001.585.28<0.0011.000.802.930.201
  positive182.881.53
PD-L1 Stroma
  low1171.000.170.55<0.0011.000.211.170.108
  high770.300.49
PD-L1 Tumor
  low1021.000.300.820.0061.000.270.830.010
  high920.490.47
PD-L1 Stroma+Tumor
  low1431.000.200.800.0101.000.453.840.616
  high510.401.32
* TNM-8—TNM Classification of Malignant Tumours—8th edition.
Table 3. Results of the Cox regression analysis identifying clinicopathological factors affecting RFS.
Table 3. Results of the Cox regression analysis identifying clinicopathological factors affecting RFS.
VariablesNUnivariableMultivariable
HR95%CIp-ValueHR95%CIp-Value
Age
  under 70 y.o1141.000.521.260.344
  over and 70 y.o800.81
Sex
  male 1691.000.251.180.121
  female250.54
Differentiation
  differentiated1401.001.052.740.0311.000.711.920.534
  include poorly diff.411.701.17
Preoperative treatment
  none971.001.914.79<0.0011.000.992.980.051
  chemotherapy(+/−Radiation)973.031.72
pT *
  <T2941.002.366.13<0.0011.000.892.910.114
  T2 and over 1003.801.61
pN *
  negative 711.002.689.15<0.0011.001.175.580.018
  positive1234.962.56
Lymphatic invasion
  negative 701.001.664.82<0.0011.000.532.150.850
  positive1242.831.07
Venous invasion
  negative 551.001.715.80<0.0011.000.623.030.436
  positive1393.141.37
Intramural metastasis
  negative 1761.001.444.560.0011.000.702.440.397
  positive182.561.31
PD-L1 Stroma
  low1171.000.240.65<0.0011.000.431.480.465
  high770.390.79
PD-L1 Tumor
  low1081.000.350.870.0101.000.320.910.022
  high870.550.54
PD-L1 Stroma+Tumor
  low1631.000.150.800.0131.000.332.770.930
  high310.350.95
* TNM-8—TNM Classification of Malignant Tumours.
Table 4. Correlation analysis between PD-L1 expression in stroma, tumor, and stroma + tumor, with various clinicopathological factors using the chi-square test.
Table 4. Correlation analysis between PD-L1 expression in stroma, tumor, and stroma + tumor, with various clinicopathological factors using the chi-square test.
VariablesPD-L1 StromaPD-L1 TumorPD-L1 Stroma + Tumor
Low
N = 117
High
N = 77
p-ValueLow
N = 102
High
N = 92
p-ValueLow
N = 143
High
N = 51
p-Value
Age
  under 70 y.o37%22%0.33331%27% 0.756 44%14% 0.514
  over and 70 y.o23%18%21%20%29%12%
Sex
  male55%32%0.07447%40% 0.357 65%22% 0.237
  female6%7%6%7%8%5%
Tumor Location
  Ce2%1%0.7362%1% 0.298 2%1% 0.498
  Ut7%4%7%5%8%4%
  Mt29%20%28%21%37%12%
  Lt22%15%16%21%27%10%
  Ae0%1%0%1%0%1%
Differentiation
  differentiated42%30%0.70136%37% 0.267 53%19% 0.962
  include poorly differentiated14%7%12%9%15%6%
  NA4%3%5%2%5%2%
Preoperative treatment
  None28%22%0.39428%22% 0.538 35%15% 0.510
  Chemotherapy29%16%22%23%35%10%
  Chemotherapy+Radiation4%2%3%2%4%1%
pT *
  T122%27%<0.00125%23% 0.988 32%16% 0.051
  T210%4%7%6%10%3%
  T329%9%20%18%31%7%
pN *
  N015%21%<0.00118%19% 0.362 23%14% 0.041
  N122%11%21%13%26%7%
  N213%5%8%9%14%3%
  N310%3%6%6%10%2%
Lymphatic invasion
  negative16%20%<0.00119%18% 0.810 24%12% 0.057
  positive44%20%34%30%50%14%
Venous invasion
  negative10%19%<0.00113%15% 0.352 16%12% 0.002
  positive51%21%39%32%57%14%
Intramural metastasis
  negative54%37%0.27847%43% 0.790 67%24% 0.880
  positive7%3%5%4%7%3%
pStage *
  I8%20%<0.00114%14% 0.658 16%12% 0.013
  II17%5%11%11%17%5%
  III20%10%18%12%24%7%
  IVA15%4%10%9%16%3%
* TNM-8—TNM Classification of Malignant Tumours—8th edition; Ce—Cervical part of the esophagus; Ut—Upper thoracic part of the esophagus; Mt—Middle thoracic part of the esophagus; Lt—Lower thoracic part of the esophagus; Ae—Abdominal part of the esophagus.
Table 5. Correlation analysis of recurrence patterns, total recurrence types, and total metastatic sites with PD-L1 expression in stroma, tumor, and stroma + tumor regions using the chi-square test.
Table 5. Correlation analysis of recurrence patterns, total recurrence types, and total metastatic sites with PD-L1 expression in stroma, tumor, and stroma + tumor regions using the chi-square test.
PD-L1 StromaPD-L1 TumorPD-L1 Stroma + Tumor
VariablesLow
N = 65
High
N = 20
p-ValueLow
N = 52
High
N = 33
p-ValueLow
N = 79
High
N = 6
p-Value
Lymphogenous metastasis
  negative22%10%0.24823%12% 0.339 20%0% 0.135
  positive78%90%77%88%80%100%
Hematogenous metastasis
  negative35%50%0.24140%36% 0.247 38%50% 0.127
  positive65%50%60%64%62%50%
Locoregional recurrence
  negative86%100%0.07888%91% 0.386 89%100% 0.289
  positive14%0%12%9%11%0%
Pleural/peritoneal dissemination
  negative78%75%0.74579%76%0.73978%67%0.503
  positive22%25%21%24%22%33%
Total recurrence types
145%55%0.51546%48% 0.139 43%50% 0.742
237%25%38%27%32%17%
314%20%15%15%13%33%
45%0%0%9%3%0%
Total metastatic sites
140%55%0.35244%42%0.31841%50%0.115
222%10%25%9%19%0%
323%20%19%27%19%17%
49%5%8%9%8%17%
56%5%4%9%4%17%
60%5%0%3%0%0%
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

Murakami, T.; Booka, E.; Furuhashi, S.; Sakai, Y.; Sekimori, K.; Haneda, R.; Fujihiro, M.; Matsumoto, T.; Morita, Y.; Kikuchi, H.; et al. Contrasting Roles of Programmed Death-Ligand 1 Expression in Tumor and Stroma in Prognosis of Esophageal Squamous Cell Carcinoma. Cancers 2024, 16, 1135. https://doi.org/10.3390/cancers16061135

AMA Style

Murakami T, Booka E, Furuhashi S, Sakai Y, Sekimori K, Haneda R, Fujihiro M, Matsumoto T, Morita Y, Kikuchi H, et al. Contrasting Roles of Programmed Death-Ligand 1 Expression in Tumor and Stroma in Prognosis of Esophageal Squamous Cell Carcinoma. Cancers. 2024; 16(6):1135. https://doi.org/10.3390/cancers16061135

Chicago/Turabian Style

Murakami, Tomohiro, Eisuke Booka, Satoru Furuhashi, Yuki Sakai, Kenichi Sekimori, Ryoma Haneda, Mayu Fujihiro, Tomohiro Matsumoto, Yoshifumi Morita, Hirotoshi Kikuchi, and et al. 2024. "Contrasting Roles of Programmed Death-Ligand 1 Expression in Tumor and Stroma in Prognosis of Esophageal Squamous Cell Carcinoma" Cancers 16, no. 6: 1135. https://doi.org/10.3390/cancers16061135

APA Style

Murakami, T., Booka, E., Furuhashi, S., Sakai, Y., Sekimori, K., Haneda, R., Fujihiro, M., Matsumoto, T., Morita, Y., Kikuchi, H., Hiramatsu, Y., Baba, S., & Takeuchi, H. (2024). Contrasting Roles of Programmed Death-Ligand 1 Expression in Tumor and Stroma in Prognosis of Esophageal Squamous Cell Carcinoma. Cancers, 16(6), 1135. https://doi.org/10.3390/cancers16061135

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