Next Article in Journal
Engagement in Lifestyle Activities is Associated with Increased Alzheimer’s Disease-Associated Cortical Thickness and Cognitive Performance in Older Adults
Next Article in Special Issue
Clinicopathological Significance of RUNX1 in Non-Small Cell Lung Cancer
Previous Article in Journal
MSI and EBV Positive Gastric Cancer’s Subgroups and Their Link with Novel Immunotherapy
Previous Article in Special Issue
Immune Checkpoint Inhibitors for Lung Cancer Treatment: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

HIP1R Expression and Its Association with PD-1 Pathway Blockade Response in Refractory Advanced NonSmall Cell Lung Cancer: A Gene Set Enrichment Analysis

1
Department of Pathology, Ajou University School of Medicine, Suwon 16499, Korea
2
Department of Thoracic and Cardiovascular Surgery, Ajou University School of Medicine, Suwon 16499, Korea
3
Department of Hematology-Oncology, Ajou University School of Medicine, Suwon 16499, Korea
*
Author to whom correspondence should be addressed.
J. Clin. Med. 2020, 9(5), 1425; https://doi.org/10.3390/jcm9051425
Submission received: 13 April 2020 / Revised: 8 May 2020 / Accepted: 8 May 2020 / Published: 11 May 2020

Abstract

:
Huntingtin-interacting protein 1-related protein (HIP1R) plays an important role in the regulation of programmed death-ligand 1 (PD-L1). The aim of this study was to investigate the expression of HIP1R and confirm its predictive or prognostic roles in anti-PD-1 therapy in nonsmall cell lung cancer (NSCLC) patients. HIP1R and PD-L1 immunohistochemical expression was examined in 52 refractory advanced NSCLC patients treated with anti-PD-1 inhibitors. We performed gene set enrichment analysis (GSEA) to detect HIP1R-specific gene sets. Patients in the PD-1 inhibitor responder group had lower HIP1R expression by univariate logistic regression analysis (odds ratio (OR) = 0.235, p = 0.015) and multivariate logistic regression analysis (OR = 0.209, p = 0.014). Patients with high HIP1R expression had poorer progression-free survival (PFS) than patients with low HIP1R expression in univariate analysis (p = 0.037) and multivariate Cox analysis (hazard ratio = 2.098, p = 0.019). The web-based mRNA dataset also showed that high HIP1R expression correlated with inferior overall survival in lung adenocarcinoma (p = 0.026). GSEA revealed that HIP1R levels correlate with a set of genes that reflect PD-L1-related immune pathways. HIP1R expression may be a promising predictor for determination of patient responses to anti-PD-1 treatment.

1. Introduction

Emergence of immune checkpoint inhibitors was a turning point in the treatment of advanced nonsmall cell lung cancer (NSCLC). Therapies targeting the programmed cell death protein 1 (PD-1) checkpoint, such as nivolumab and pembrolizumab, have yielded impressive responsive rates in advanced NSCLC patients otherwise refractory to multiples lines of therapy [1,2,3]. However, the overall response rate for PD-1 inhibitor therapy is approximately 15–20% in unselected patients with NSCLC, and between 15% and 45% in patients with PD-L1-expressing NSCLC [4]. We need a biomarker that can more accurately predict the response to PD-1 inhibitors.
Expression of programmed death ligand-1 (PD-L1), the PD-1 ligand, is currently the most widely used biomarker for PD-1 inhibition. To identify patients who preferentially respond to PD-1 blockade, we need to better understand how the PD-1 pathway is regulated. Recently, several mechanisms have been reported to underlie PD-1 pathway regulation. CKLF-like MARVEL transmembrane-domain-containing 6 (CMTM6) regulates the PD-1 pathway by maintaining the expression of PD-L1, and CMTM6 is a predictor of the response to PD-1 inhibitors [5,6]. F-box only protein 38 (FBXO38) mediates PD-1 ubiquitination of T cells, and knockout of FBXO38 in such cells induces tumor progression in a mouse model due to increased PD-1 expression by tumor-infiltrating T cells [7]. AXL expression displays a positive correlation with PD-L1 expression in lung adenocarcinoma with epidermal growth factor receptor (EGFR) mutation, and abolition of AXL kinase activity inhibits PD-L1 mRNA expression in a lung adenocarcinoma cell line with EGFR mutation [8].
Recent research has uncovered new strategies to remove specific unwanted proteins by using cellular protein degradation mechanisms, including lysosome-targeting molecules [9], proteolysis-targeting chimeras (PROTACs) [10], and tag-based degradation systems (dTAG) [11]. Wang et al. reported that Huntingtin-interacting protein 1-related protein (HIP1R) promotes lysosomal degradation of PD-L1, inhibits HIP1R-induced PD-L1 accumulation, and alters T cell–mediated cytotoxicity in a human colorectal cancer cell line [12]. A chimeric peptide including a lysosomal sorting signal and the HIP1R PD-L1-binding sequence significantly inhibits PD-L1 protein expression [12]. Although immune checkpoint inhibition is the most popular treatment for lung cancer, relationships involving HIP1R and immune checkpoint inhibitors in lung cancer have not been studied.
The present study was conducted to determine whether HIP1R protein expression affects the response of NSCLC patients to anti-PD-1 inhibitors and their prognosis. The relationship between HIP1R and PD-L1 was also evaluated, employing immunohistochemical and web-based mRNA expression data. In addition, we performed gene set enrichment analysis (GSEA) on RNA-sequencing data from The Cancer Genome Atlas (TCGA) to confirm the molecular pathways associated with HIP1R expression.

2. Materials and Methods

2.1. Patients

We retrospectively selected 52 advanced NSCLC patients who were administered PD-1 inhibitor from 2016 to 2019, and they previously received one or two lines of chemotherapy. This study was approved by the Institutional Review Board of Ajou University School of Medicine. Informed consent was waived due to the retrospective nature of the study (AJIRB-BMR-KSP-19-050 and 2019-03-26).
Patients treated with PD-1 inhibitor were assigned to either a responder group (complete response, partial response, or stable disease) or a nonresponder group (disease progression), according to the response evaluation criteria for solid tumors (RECIST) version 1.1 [13].

2.2. Immunohistochemical Staining and HIP1R Expression Scoring

One board-certified pathologist (YWK) reviewed hematoxylin and eosin (H&E)-stained tissue samples to determine a definitive pathologic diagnosis according to the 2015 World Health Organization Classification of Lung Tumors [14]. All patients were pathologically staged according to the eighth edition of the TNM classification.
HIP1R immunohistochemical (IHC) staining was performed with a Benchmark XT automatic IHC staining device (Ventana Medical Systems, Tucson, AZ, USA). The samples were incubated with an anti-HIP1R antibody (dilution 1:1000, 16814-1-AP, polyclonal, Proteintech, Rosemont, IL, USA). We used a human placenta tissue as positive control according to the manufacturer’s recommendations (Figure S1). We also evaluated the intensity of HIP1R staining on a four-point intensity scale: 0 (no staining), 1 (light yellow = faint staining), 2 (yellow-brown = moderate staining), and 3 (brown = strong staining) (Figure 1). We also evaluated the percentages (0–100%) of cytoplasmic versus membranous localization of HIP1R. We used H-scores to interpret HIP1R staining [15], where H-score = [1 × (% cells 1+) + 2 × (% cells 2+) + 3 × (% cells 3+)]. H-scores (0–300) were obtained by multiplying the percentage of cells by the intensity score.

2.3. Immunohistochemical Staining and PD-L1 Expression Scoring

Two PD-L1 antibodies (clone name SP263 or 22C3) were used to detect PD-L1 expression. Sp263 was a companion diagnostic assay for OPDIVO® (nivolumab), and 22c3 was a companion diagnostic assay for KEYTRUDA® (pembrolizumab). We performed SP263 and/or 22C3 assays prior to PD-1 inhibitor treatment for all NSCLC patients. Thirteen (25%) of the 52 specimens were tested for both SP263 and 22C3, 27 (51.9%) for only SP263, and 12 (23.1%) for only 22C3. Two PD-L1 tests used prediluted antibody (ready to use) according to the protocol. The SP263 assay was performed using a VENTANA BenchMark ULTRA instrument (Roche, Basel, Switzerland), and the 22C3 assay was conducted using the Dako Link-48 platform (Dako, Carpinteria, California, US), as recommended by the manufacturers [16]. PD-L1 intensity was also evaluated on a four-point intensity scale (0, none; 1, faint; 2, moderate; and 3, strong), and the percentage of membranous expression of PD-L1 was determined (Figures S2 and S3). When both the 22C3 and SP263 tests were conducted, mean values were used. High PD-L1 expression was defined as ≥ 50% of definitive tumor cells exhibiting PD-L1 staining, because 50% was the cut-off used for NSCLC [17].

2.4. Web-Based mRNA Profiling, GSEA, and Kaplan Meier Analysis

The mRNA sequencing data of 517 lung adenocarcinoma patients and 501 lung squamous cell carcinoma patients were downloaded from The Cancer Genome Atlas (TCGA) cBioportal (http://cbioportal.org) [18]. We conducted correlation analysis involving PD-L1 and HIP1R mRNA sequencing data.
GSEA is a method of analyzing associations between gene expression and biological information. We conducted GSEA using GSEA version 4.0.3 from the Broad Institute at MIT and Harvard (http://www.broadinstitute.org/gsea/index.jsp) [19]. TCGA mRNA sequencing data derived from lung adenocarcinoma and lung squamous cell carcinoma patients was used. Depending on the median value, it is divided into low and high HIP1R. Hallmark gene sets representing well-defined biological states or processes were used for GSEA. 1000 permutations were used for estimating nominal p values. If the p value was less than 0.05 and the False Discovery Rate (FDR) was less than 0.25, the findings were considered statistically significant.
We conducted survival analyses using an online Kaplan Meier plotter tool [20]. The online Kaplan Meier plotter tool provides mRNA expression data of cancer patients and allows for survival analysis. Survival analyses were performed in 719 patients with lung adenocarcinoma and 524 lung squamous cell carcinoma cases according to their HIP1R mRNA expression.

2.5. Statistical Analyses

Spearman’s rank order correlation coefficient analysis was used to measure monotonic relationships between continuous variables. Mann Whitney U-tests were used to compare differences between two independent groups. Univariate and multivariate logistic regression analyses were performed to determine factors that predicted a response to PD-1 inhibitors. Receiver operating curve (ROC) analysis was used to determine the cut-off values for HIP1R expression. The progression-free survival (PFS) or overall survival (OS) difference between the cohorts was determined using the log-rank test. Univariate and multivariate prognostic analyses were performed for PFS using a Cox proportional hazards regression model. SPSS Statistics for Windows (Version 25.0, IBM, Armonk, NY, USA) was used for all analyses, and p values <0.05 were considered to be statistically significant.

3. Results

3.1. Patient Demographics

Detailed patient and tumor characteristics are summarized in Table 1. Forty-six tissues were collected from lung lesions, and six were obtained from metastatic sites. Twenty-seven (51.9%) patients had been treated with nivolumab, and 25 (48.1%) patients received pembrolizumab. All patients were refractory to conventional treatments such as chemotherapy, radiotherapy, or target therapy. Therefore, they received PD-1 inhibitor as a second line or later setting. Twenty-seven (51.9%) patients were classified as responders, and 25 (48.1%) were classified as nonresponders. Four patients were treated with epidermal growth factor receptor (EGFR) inhibitors before PD-1 inhibitor administration. Patients with anaplastic lymphoma kinase (ALK) fusion were not identified in the present study.

3.2. Relationships Between HIP1R and PD-L1 Analyzed by IHC and mRNA Expression

We performed correlation analysis of HIP1R and PD-L1 expression using IHC techniques. There was no statistically significant correlation between HIP1R and PD-L1 expression (p = 0.905, Figure 2A).
Correlation analyses of HIP1R and PD-L1 expression were performed using mRNA data. From the TCGA dataset, HIP1R mRNA expression levels were negatively correlated with PD-L1 mRNA levels in adenocarcinoma and squamous cell carcinoma. (Spearman’s rho = −0.233, p < 0.001, Figure 2B; Spearman’s rho = −0.224, p < 0.001, Figure 2C).

3.3. Associations Involving HIP1R, PD-L1, Clinicopathologic Parameters, and Response to PD-1 Inhibitors

ROC analysis was performed to determine the cut-off value for HIP1R expression. Cut-off was determined as the value corresponding to the maximum joint sensitivity and specificity of the ROC curve. The area under the curve (AUC) was 0.659 for the expression of HIP1R, and the cut-off value was 180 (66% sensitivity and 68% specificity, Figure 3).
We explored the predictive capacity of HIP1R, PD-L1, and clinicopathologic factors in terms of responses to PD-1 inhibition. By univariate analysis, the expression of HIP1R was found to be a predictor of the response to anti-PD-1 therapy (OR) = 0.235, p = 0.015; Table 2). PD-L1 expression was also found to be a predictor of the response to anti-PD-1 therapy (OR = 4.062, p = 0.028; Table 2). By multivariate analysis, the expression of HIP1R was an independent predictor of anti-PD-1 therapy response (OR = 0.209, p = 0.014, Table 2).

3.4. GSEA According to HIP1R mRNA Expression

We performed GSEA to identify gene sets associated with HIP1R mRNA expression in the TCGA mRNA data of lung adenocarcinoma and lung squamous cell carcinoma cases. In lung adenocarcinoma, we identified the top 20 most prominent pathways that were upregulated in the low HIP1R mRNA expression group (Table S1). Four of the 20 were immune-related gene sets and were statistically significant (HALLMARK_ALLOGRAFT_REJECTION, HALLMARK_INFLAMMATORY_RESPONSE, HALLMARK_IL6_JAK_STAT3_SIGNALING, and HALLMARK_IL2_STAT5_SIGNALING) (Figure 4). HALLMARK_INTERFERON_GAMMA_RESPONSE is also upregulated in the low HIP1R mRNA expression group, although the statistical significance was marginal (p = 0.063). Core enrichment gene lists for HALLMARK_ALLOGRAFT_REJECTION, HALLMARK_INFLAMMATORY_RESPONSE, HALLMARK_IL6_JAK_STAT3_SIGNALING and HALLMARK_IL2_STAT5_SIGNALING are summarized in Tables S2–S5. In lung squamous cell carcinoma, there were no statistically significant immune-related gene sets associated with HIP1R mRNA expression (Table S6).

3.5. Prognostic Significance of HIP1R and PD-L1

Patients with high HIP1R expression had inferior PFS to patients with low HIP1R expression (p = 0.037, Figure 5A). Patients with high HIP1R expression also showed an inferior OS than patients with low HIP1R expression, however the statistical significance was not reached (p = 0.11, Figure 5B). Patients with high PD-L1 expression had superior PFS or OS than patients with low PD-L1 expression (p = 0.028, Figure 5C and p = 0.031, Figure 5D, respectively). Furthermore, patients with high HIP1R expression and low PD-L1 expression had lower PFS or OS than patients with other expression patterns (p < 0.001, Figure 5E and p = 0.001, Figure 5F, respectively). In multivariate analysis, high HIP1R expression was an independent prognostic factor for PFS (HR = 2.098, p = 0.019, Table 3).
We used a Kaplan Meier plotter tool and performed survival analysis according to HIP1R mRNA expression in lung adenocarcinoma and lung squamous cell carcinoma patients. The group with high HIP1R mRNA expression exhibited poorer OS in patients with adenocarcinoma (p = 0.026, Figure S4A). However, in lung squamous cell carcinoma, HIP1R mRNA expression was not correlated with OS (p = 0.63, Figure S4B)

4. Discussion

This study had several novel discoveries. First, we found that the expression of HIP1R was an independent predictive factor for anti-PD-1 treatment response by NSCLC patients. Second, the expression of HIP1R was an independent prognostic factor of PFS in patients treated with anti-PD-1 inhibitors. Third, GSEA revealed that HIP1R mRNA expression was tightly correlated with immune-related gene sets in lung adenocarcinoma. These GSEA results suggested that HIP1R mRNA expression plays an important role in regulating the expression of PD-L1.
GSEA revealed that low HIP1R mRNA expression was closely associated with allograft rejection, inflammatory responses, IL6-JAK-STAT3, IL2-STAT5, and interferon gamma response pathways in lung adenocarcinoma. PD-L1 expression is correlated with marked expression of adaptive immune responses (CD8+ T-cells) [21]. In our study, CD8 was also included in the core enrichment gene list of HALLMARK_ALLOGRAFT_REJECTION. Previous studies have also reported that the IL6-JAK-STAT3 pathway induces PD-L1 upregulation. IL-6 is positively correlated with PD-L1 expression in human hepatocellular carcinoma (HCC) cells, and IL-6 induces PD-L1 stability through glycosylation in a HCC cell line [22]. Glioblastoma-derived IL6 is required for up-regulation of myeloid PD-L1 in glioblastoma through a STAT3-dependent mechanism [23]. Combined blockade of IL6 and PD-L1 signaling achieves synergistic antitumor immune responses in colon carcinoma and murine melanoma models [24]. PD-L1 expression is also regulated by interferon gamma signaling in a melanoma cell line [25]. GSEA suggested that HIP1R expression plays an important role in adaptive immune responses associated with PD-L1.
In the present study, no correlation was identified between HIP1R and PD-L1 protein expression, However, HIP1R mRNA expression was negatively correlated with the mRNA expression level of PD-L1 in adenocarcinoma and squamous cell carcinoma. There are several possible explanations for this discrepancy. Post-transcriptional and post-translational modifications affect the level of protein expression [26]. Proteins can have significantly different half-lives in vivo [27]. There were no cases of surgery in our study, therefore the only sample we received was a biopsy. We cannot conduct additional experiments for mRNA testing of HIP1R and PD-L1, because very little tumor tissue remains in paraffin tissue.
Patients with high HIP1R mRNA expression exhibited poor clinical outcomes in web-based mRNA data of adenocarcinoma cases; however, HIP1R mRNA expression was not correlated with OS in squamous cell carcinoma. From our IHC data, HIP1R expression was correlated with poor clinical outcomes. However, we did perform subgroup analysis according to histologic type because of our small sample size. HIP1R levels also correlate with a set of genes that reflect PD-L1-related immune pathways in GSEA analysis of adenocarcinoma cases; however, there were no statistically significant immune-related gene sets associated with HIP1R mRNA expression in squamous cell carcinoma. Currently, lung adenocarcinoma and squamous cell carcinoma are known to involve different biologic mechanisms and prognoses. Therefore, the role of HIP1R in adaptive immune responses and its effect on clinical outcomes may vary depending on the histological type.
Despite some surprising discoveries, our study has certain limitations. First, our cohort was small. We performed multivariate logistic regression and prognostic analyses on only 52 samples. However, the web-based mRNA dataset also revealed results similar to ours. These results encourage further investigations involving larger populations. Second, we used an IHC method to detect HIP1R protein expression. There is no information regarding standardization, reliability, and reproducibility of IHC staining. We used the same antibody that Wang et al. used [12]. However, Wang et al. used HIP1R antibody (16814-1-AP) in Western Blott (WB) and immunofluorescence alone. Only recently has HIP1R attracted attention in cancer research, so few studies have been done on HIP1R. Therefore, there are no antibodies that are commonly used in immunohistochemistry. In the catalog of HIP1R antibody (16814-1-AP), it can be used in IHC, immunoprecipitation (IP), WB, and ELISA. According to the manufacturer’s guidelines, this antibody was validated by western blot in HeLa cells and human liver tissue. We used a human placenta tissue as positive control as recommended. An automatic IHC staining device (Benchmark XT) may improve the reproducibility of IHC staining. The H-scoring method is widely used for immunochemical staining, and is known to have relatively high reproducibility among pathologists [28,29]. Third, we examined the protein expression of HIP1R and PD-L1 in refractory advanced NSCLC; however, the mRNA profiles of HIP1R and PD-L1 were not evaluated. Because protein expression of HIP1R did not correlate with PD-L1 expression, the relationship between HIP1R and PD-L1 mRNA expression profile is very important. To verify the results of GSEA, we should evaluate the mRNA expression profiles of HIP1R. However, the sample we have is a small biopsy, and we have already performed several immunohistochemical stainings for diagnosis and ALK and EGFR mutation tests. Therefore, currently, very little tumor tissue remains in paraffin tissue and we cannot conduct additional experiments for mRNA testing. To confirm our experiments, future research should measure the mRNA expression level of HIP1R on many samples and investigate the relationship with the PD-1 blocker and PD-L1 expression.
In conclusion, we examined the expression of HIP1R in 52 refractory NSCLC samples from patients treated with PD-1 inhibitors. HIP1R expression was an independent biomarker predicting patient response to PD-1 inhibitors. High HIP1R expression was an independent predictor of poor PFS. In addition, HIP1R mRNA expression was significantly correlated with immune-related gene sets in lung adenocarcinoma. These immune-related gene sets are known to play important roles in PD-L1 regulation. Based on our findings, HIP1R expression may be a promising predictor for the therapeutic determination of responses to anti-PD-1 treatment.

Supplementary Materials

The following are available online at https://www.mdpi.com/2077-0383/9/5/1425/s1, Table S1. Gene sets within the top 20-ranked list related to low HIP1R in adenocarcinoma patients; Table S2. Core enrichment gene list of HALLMARK_ALLOGRAFT_REJECTION; Table S3. Core enrichment gene list of HALLMARK_INFLAMMATORY_RESPONSE; Table S4. Core enrichment gene list of HALLMARK_IL6_JAK_STAT3_SIGNALING; Table S5. Core enrichment gene list of HALLMARK_IL2_STAT5_SIGNALING; Table S6. Gene sets within the top 20-ranked list related to low HIP1R in squamous cell carcinoma patients; Figure S1. HIP1R expression in positive control. Positive HIP1R expression in placental tissue, X400; Figure S2. PD-L1 sp263 expression in nonsmall cell carcinoma. (A) No staining of PD-L1 sp263, x400. (B) Faint PD-L1 sp263 staining, x400. (C) Moderate PD-L1 sp263 staining, x400. (D) Strong PD-L1 sp263 staining, x400; Figure S3. PD-L1 22c3 expression in nonsmall cell carcinoma. (A) No staining of PD-L1 22c3, x400. (B) Faint PD-L1 22c3 staining, x400. (C) Moderate PD-L1 22c3 staining, x400. (D) Strong PD-L1 22c3 staining, x400; Figure S4. Comparison of survival rates, according to HIP1R mRNA expression in patients with nonsmall cell carcinoma. (A) Overall survival (OS) and HIP1R in lung adenocarcinoma. (B) OS and HIP1R in lung squamous cell carcinoma.

Author Contributions

Conceptualization, Y.W.K.; formal analysis, Y.W.K.; funding acquisition, Y.W.K.; investigation, Y.W.K., J.-H.H., S.H., and H.W.L.; methodology, Y.W.K., J.-H.H., S.H., and H.W.L.; resources, Y.W.K.; supervision, Y.W.K.; visualization, Y.W.K.; writing—original draft, Y.W.K., J.-H.H., S.H., and H.W.L.; writing—review and editing, Y.W.K., J.-H.H, S.H., and H.W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by Basic Science Research Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Science, ICT (NRF-2017R1C1B5076342 for Young Wha Koh) and the faculty research fund (Ajou translational research fund 2018) of Ajou University School of Medicine to Young Wha Koh and Seokjin Haam (M-2018-C0460-00035).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Garon, E.B.; Rizvi, N.A.; Hui, R.; Leighl, N.; Balmanoukian, A.S.; Eder, J.P.; Patnaik, A.; Aggarwal, C.; Gubens, M.; Horn, L.; et al. Pembrolizumab for the treatment of nonsmall-cell lung cancer. New Engl. J. Med. 2015, 372, 2018–2028. [Google Scholar] [CrossRef] [PubMed]
  2. Gettinger, S.; Rizvi, N.A.; Chow, L.Q.; Borghaei, H.; Brahmer, J.; Ready, N.; Gerber, D.E.; Shepherd, F.A.; Antonia, S.; Goldman, J.W.; et al. Nivolumab Monotherapy for First-Line Treatment of Advanced NonSmall-Cell Lung Cancer. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2016, 34, 2980–2987. [Google Scholar] [CrossRef] [PubMed]
  3. Herbst, R.S.; Baas, P.; Kim, D.W.; Felip, E.; Perez-Gracia, J.L.; Han, J.Y.; Molina, J.; Kim, J.H.; Arvis, C.D.; Ahn, M.J.; et al. Pembrolizumab versus docetaxel for previously treated, PD-L1-positive, advanced nonsmall-cell lung cancer (KEYNOTE-010): A randomised controlled trial. Lancet 2016, 387, 1540–1550. [Google Scholar] [CrossRef]
  4. Shukuya, T.; Carbone, D.P. Predictive Markers for the Efficacy of Anti-PD-1/PD-L1 Antibodies in Lung Cancer. J. Thorac. Oncol. Off. Publ. Int. Assoc. Study Lung Cancer 2016, 11, 976–988. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Mezzadra, R.; Sun, C.; Jae, L.T.; Gomez-Eerland, R.; de Vries, E.; Wu, W.; Logtenberg, M.E.W.; Slagter, M.; Rozeman, E.A.; Hofland, I.; et al. Identification of CMTM6 and CMTM4 as PD-L1 protein regulators. Nature 2017, 549, 106–110. [Google Scholar] [CrossRef]
  6. Koh, Y.W.; Han, J.H.; Haam, S.; Jung, J.; Lee, H.W. Increased CMTM6 can predict the clinical response to PD-1 inhibitors in nonsmall cell lung cancer patients. Oncoimmunology 2019, 8, e1629261. [Google Scholar] [CrossRef] [Green Version]
  7. Meng, X.; Liu, X.; Guo, X.; Jiang, S.; Chen, T.; Hu, Z.; Liu, H.; Bai, Y.; Xue, M.; Hu, R.; et al. FBXO38 mediates PD-1 ubiquitination and regulates anti-tumour immunity of T cells. Nature 2018, 564, 130–135. [Google Scholar] [CrossRef]
  8. Tsukita, Y.; Fujino, N.; Miyauchi, E.; Saito, R.; Fujishima, F.; Itakura, K.; Kyogoku, Y.; Okutomo, K.; Yamada, M.; Okazaki, T.; et al. Axl kinase drives immune checkpoint and chemokine signalling pathways in lung adenocarcinomas. Mol. Cancer 2019, 18, 24. [Google Scholar] [CrossRef] [Green Version]
  9. Bauer, P.O.; Goswami, A.; Wong, H.K.; Okuno, M.; Kurosawa, M.; Yamada, M.; Miyazaki, H.; Matsumoto, G.; Kino, Y.; Nagai, Y.; et al. Harnessing chaperone-mediated autophagy for the selective degradation of mutant huntingtin protein. Nat. Biotechnol. 2010, 28, 256–263. [Google Scholar] [CrossRef]
  10. Sakamoto, K.M.; Kim, K.B.; Kumagai, A.; Mercurio, F.; Crews, C.M.; Deshaies, R.J. Protacs: Chimeric molecules that target proteins to the Skp1-Cullin-F box complex for ubiquitination and degradation. Proc. Natl. Acad. Sci. USA 2001, 98, 8554–8559. [Google Scholar] [CrossRef] [Green Version]
  11. Nabet, B.; Roberts, J.M.; Buckley, D.L.; Paulk, J.; Dastjerdi, S.; Yang, A.; Leggett, A.L.; Erb, M.A.; Lawlor, M.A.; Souza, A.; et al. The dTAG system for immediate and target-specific protein degradation. Nat. Chem. Biol. 2018, 14, 431–441. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, H.; Yao, H.; Li, C.; Shi, H.; Lan, J.; Li, Z.; Zhang, Y.; Liang, L.; Fang, J.Y.; Xu, J. HIP1R targets PD-L1 to lysosomal degradation to alter T cell-mediated cytotoxicity. Nat. Chem. Biol. 2019, 15, 42–50. [Google Scholar] [CrossRef] [PubMed]
  13. 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, 228–247. [Google Scholar] [CrossRef]
  14. Travis, W.D.; Brambilla, E.; Nicholson, A.G.; Yatabe, Y.; Austin, J.H.; Beasley, M.B.; Chirieac, L.R.; Dacic, S.; Duhig, E.; Flieder, D.B.; et al. The 2015 World Health Organization Classification of Lung Tumors: Impact of Genetic, Clinical and Radiologic Advances Since the 2004 Classification. J. Thorac. Oncol. Official Publ. Int. Assoc. Study Lung Cancer 2015, 10, 1243–1260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. 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]
  16. Munari, E.; Rossi, G.; Zamboni, G.; Lunardi, G.; Marconi, M.; Sommaggio, M.; Netto, G.J.; Hoque, M.O.; Brunelli, M.; Martignoni, G.; et al. PD-L1 Assays 22C3 and SP263 are Not Interchangeable in NonSmall Cell Lung Cancer When Considering Clinically Relevant Cutoffs: An Interclone Evaluation by Differently Trained Pathologists. Am. J. Surg. Pathol. 2018, 42, 1384–1389. [Google Scholar] [CrossRef]
  17. Büttner, R.; Gosney, J.R.; Skov, B.G.; Adam, J.; Motoi, N.; Bloom, K.J.; Dietel, M.; Longshore, J.W.; López-Ríos, F.; Penault-Llorca, F.; et al. Programmed Death-Ligand 1 Immunohistochemistry Testing: A Review of Analytical Assays and Clinical Implementation in NonSmall-Cell Lung Cancer. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2017, 35, 3867–3876. [Google Scholar] [CrossRef]
  18. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio cancer genomics portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [Green Version]
  19. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  20. Gyorffy, B.; Surowiak, P.; Budczies, J.; Lanczky, A. Online survival analysis software to assess the prognostic value of biomarkers using transcriptomic data in nonsmall-cell lung cancer. PLoS ONE 2013, 8, e82241. [Google Scholar] [CrossRef] [Green Version]
  21. Huang, C.Y.; Wang, Y.; Luo, G.Y.; Han, F.; Li, Y.Q.; Zhou, Z.G.; Xu, G.L. Relationship between PD-L1 Expression and CD8+ T-cell Immune Responses in Hepatocellular Carcinoma. J. Immunother. 2017, 40, 323–333. [Google Scholar] [CrossRef] [PubMed]
  22. Chan, L.C.; Li, C.W.; Xia, W.; Hsu, J.M.; Lee, H.H.; Cha, J.H.; Wang, H.L.; Yang, W.H.; Yen, E.Y.; Chang, W.C.; et al. IL-6/JAK1 pathway drives PD-L1 Y112 phosphorylation to promote cancer immune evasion. J. Clin. Investig. 2019, 129, 3324–3338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Lamano, J.B.; Lamano, J.B.; Li, Y.D.; DiDomenico, J.D.; Choy, W.; Veliceasa, D.; Oyon, D.E.; Fakurnejad, S.; Ampie, L.; Kesavabhotla, K.; et al. Glioblastoma-Derived IL6 Induces Immunosuppressive Peripheral Myeloid Cell PD-L1 and Promotes Tumor Growth. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res. 2019, 25, 3643–3657. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Tsukamoto, H.; Fujieda, K.; Miyashita, A.; Fukushima, S.; Ikeda, T.; Kubo, Y.; Senju, S.; Ihn, H.; Nishimura, Y.; Oshiumi, H. Combined Blockade of IL6 and PD-1/PD-L1 Signaling Abrogates Mutual Regulation of Their Immunosuppressive Effects in the Tumor Microenvironment. Cancer Res. 2018, 78, 5011–5022. [Google Scholar] [CrossRef] [Green Version]
  25. Garcia-Diaz, A.; Shin, D.S.; Moreno, B.H.; Saco, J.; Escuin-Ordinas, H.; Rodriguez, G.A.; Zaretsky, J.M.; Sun, L.; Hugo, W.; Wang, X.; et al. Interferon Receptor Signaling Pathways Regulating PD-L1 and PD-L2 Expression. Cell Rep. 2017, 19, 1189–1201. [Google Scholar] [CrossRef] [Green Version]
  26. Greenbaum, D.; Colangelo, C.; Williams, K.; Gerstein, M. Comparing protein abundance and mRNA expression levels on a genomic scale. Genome Biol. 2003, 4, 117. [Google Scholar] [CrossRef] [Green Version]
  27. Nie, L.; Wu, G.; Zhang, W. Correlation of mRNA expression and protein abundance affected by multiple sequence features related to translational efficiency in Desulfovibrio vulgaris: A quantitative analysis. Genetics 2006, 174, 2229–2243. [Google Scholar] [CrossRef] [Green Version]
  28. Aviles-Salas, A.; Muniz-Hernandez, S.; Maldonado-Martinez, H.A.; Chanona-Vilchis, J.G.; Ramirez-Tirado, L.A.; HernaNdez-Pedro, N.; Dorantes-Heredia, R.; Rui, Z.M.J.M.; Motola-Kuba, D.; Arrieta, O. Reproducibility of the EGFR immunohistochemistry scores for tumor samples from patients with advanced nonsmall cell lung cancer. Oncol. Lett. 2017, 13, 912–920. [Google Scholar] [CrossRef] [Green Version]
  29. Detre, S.; Jotti, G.S.; Dowsett, M. A “quickscore” method for immunohistochemical semiquantitation: Validation for oestrogen receptor in breast carcinomas. J. Clin. Pathol. 1995, 48, 876–878. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Huntingtin-interacting protein 1-related protein (HIP1R) expression in nonsmall cell carcinoma. (A) No staining of HIP1R, x400. (B) Faint HIP1R staining, X400. (C) Moderate HIP1R staining, X400. (D) Strong HIP1R staining, X400.
Figure 1. Huntingtin-interacting protein 1-related protein (HIP1R) expression in nonsmall cell carcinoma. (A) No staining of HIP1R, x400. (B) Faint HIP1R staining, X400. (C) Moderate HIP1R staining, X400. (D) Strong HIP1R staining, X400.
Jcm 09 01425 g001
Figure 2. Correlation analyses involving Huntingtin-interacting protein 1-related protein (HIP1R) and programmed death-ligand 1 (PD-L1) expression. (A) Correlation between HIP1R and PD-L1 detected immunohistochemically. (B) Correlation between HIP1R and PD-L1 mRNA expression in lung adenocarcinoma. (C) Correlation between HIP1R and PD-L1 mRNA expression in lung squamous cell carcinoma.
Figure 2. Correlation analyses involving Huntingtin-interacting protein 1-related protein (HIP1R) and programmed death-ligand 1 (PD-L1) expression. (A) Correlation between HIP1R and PD-L1 detected immunohistochemically. (B) Correlation between HIP1R and PD-L1 mRNA expression in lung adenocarcinoma. (C) Correlation between HIP1R and PD-L1 mRNA expression in lung squamous cell carcinoma.
Jcm 09 01425 g002
Figure 3. Receiver operating characteristic (ROC) and area under the curve (AUC) analysis of HIP1R expression.
Figure 3. Receiver operating characteristic (ROC) and area under the curve (AUC) analysis of HIP1R expression.
Jcm 09 01425 g003
Figure 4. Gene set enrichment analysis (GSEA) according to HIP1R mRNA expression. (A) HALLMARK_ALLOGRAFT_REJECTION pathway; (B) HALLMARK_INFLAMMATORY_RESPONSE pathway; (C) HALLMARK_IL6_JAK_STAT3_SIGNALING pathway; (D) HALLMARK_IL2_STAT5_SIGNALING pathway.
Figure 4. Gene set enrichment analysis (GSEA) according to HIP1R mRNA expression. (A) HALLMARK_ALLOGRAFT_REJECTION pathway; (B) HALLMARK_INFLAMMATORY_RESPONSE pathway; (C) HALLMARK_IL6_JAK_STAT3_SIGNALING pathway; (D) HALLMARK_IL2_STAT5_SIGNALING pathway.
Jcm 09 01425 g004
Figure 5. Comparison of survival rates according to HIP1R and PD-L1 expression. (A) Progression-free survival (PFS) and expression of HIP1R. (B) Overall survival (OS) and expression of HIP1R (C) PFS and PD-L1. (D) OS and PD-L1. (E) PFS, HIP1R, and PD-L1. (F) OS, HIP1R, and PD-L1.
Figure 5. Comparison of survival rates according to HIP1R and PD-L1 expression. (A) Progression-free survival (PFS) and expression of HIP1R. (B) Overall survival (OS) and expression of HIP1R (C) PFS and PD-L1. (D) OS and PD-L1. (E) PFS, HIP1R, and PD-L1. (F) OS, HIP1R, and PD-L1.
Jcm 09 01425 g005
Table 1. Demographic and clinical characteristics of patients.
Table 1. Demographic and clinical characteristics of patients.
VariableNumber (%)
Age, Median (Range) (Years)64 (38–85)
Male Sex43 (82.7%)
Smoking Sistory31 (73.8%)
Histologic Subtype
Adenocarcinoma22 (42.3%)
Squamous Cell Carcinoma19 (36.5%)
Pleomorphic Carcinoma4 (7.7%)
NSCLC, NOS7 (13.5%)
Clinical Stage at Diagnosis
III13 (25%)
IV39 (75%)
Genetic Alteration Status
EGFR-Mutated4 (9.1%)
ALK-Rearranged0 (0%)
Wild Type44 (92.3%)
Type of PD-1 Blockade
Nivolumab27 (51.9%)
Pembrolizumab25 (48.1%)
PD-L1 Expression
Low (<50%)17 (32.7%)
High (≥50%)35 (67.3%)
Response to PD-1 Blockade
Responder27 (51.9%)
Nonresponder25 (48.1%)
Smoking history was collected for 42 patients. EGFR test was performed in 44 patients. ALK test was performed in 47 patients. Abbreviations: epidermal growth factor receptor; EGFR, nonsmall-cell lung cancer—not otherwise specified; NSCLC, NOS, programmed cell death protein 1, PD-1; programmed death-ligand 1, PD-L1.
Table 2. Univariate and multivariate logistic regression analysis for predicting clinical response to PD-1 blockade.
Table 2. Univariate and multivariate logistic regression analysis for predicting clinical response to PD-1 blockade.
Univariate AnalysisMultivariate Analysis
CovariateOR95% CIp-Value OR95% CIp-Value
Age (≥65 Years vs. <65 Years)1.5910.532–4.7570.406
Sex (Male vs. Female)2.5260.558–11.440.229
Smoking History (+ vs. −)3.2380.720–14.560.126
Presence of EGFR Mutation (+ vs. −)0.2220.021–2.3300.210
Type of PD-1 Blockade (Nivolumab vs. Pembrolizumab)1.8750.622–5.6490.264
PD-L1 (>50% vs. ≤50%)4.0621.166–14.150.0284.6641.198–18.150.026
HIP1R (>180 vs. ≤180)0.2350.074–0.7510.0150.2090.060–0.7310.014
logistic regression analysis. Abbreviations: confidence interval; CI, epidermal growth factor receptor; EGFR, Huntingtin Interacting Protein 1 Related; HIP1R, odd ratio; OR, programmed cell death protein 1, PD-1; programmed death-ligand 1, PD-L1.
Table 3. Univariate and multivariate analyses of progression-free survival.
Table 3. Univariate and multivariate analyses of progression-free survival.
Univariate AnalysisMultivariate Analysis
CovariateHR95% CIp-ValueHR95% CIp-Value
Age (≥65 Years vs. <65 Years)1.1200.597–2.1010.724
Sex (Male vs. Female)0.5530.251–1.2180.141
Smoking History (+ vs. −)1.1210.517–2.4290.773
Presence of EGFR Mutation (+ vs. −)1.6030.482–5.3290.441
Type of PD-1 Blockade (Nivolumab vs. Pembrolizumab)1.4820.785–2.8000.225
PD-L1 (>50% vs. ≤50%)0.4890.254–0.9420.0320.4320.222–0.8440.014
HIP1R (>180 vs. ≤180)1.9351.027–3.6480.0412.0981.136–4.1330.019
Cox proportional hazards regression model. Abbreviations: CI, confidence interval; EGFR, epidermal growth factor receptor; HIP1R, Huntingtin Interacting Protein 1 Related; HR, hazard ratio; PD-1, programmed cell death protein 1; PD-L1, programmed death-ligand 1.

Share and Cite

MDPI and ACS Style

Koh, Y.W.; Han, J.-H.; Haam, S.; Lee, H.W. HIP1R Expression and Its Association with PD-1 Pathway Blockade Response in Refractory Advanced NonSmall Cell Lung Cancer: A Gene Set Enrichment Analysis. J. Clin. Med. 2020, 9, 1425. https://doi.org/10.3390/jcm9051425

AMA Style

Koh YW, Han J-H, Haam S, Lee HW. HIP1R Expression and Its Association with PD-1 Pathway Blockade Response in Refractory Advanced NonSmall Cell Lung Cancer: A Gene Set Enrichment Analysis. Journal of Clinical Medicine. 2020; 9(5):1425. https://doi.org/10.3390/jcm9051425

Chicago/Turabian Style

Koh, Young Wha, Jae-Ho Han, Seokjin Haam, and Hyun Woo Lee. 2020. "HIP1R Expression and Its Association with PD-1 Pathway Blockade Response in Refractory Advanced NonSmall Cell Lung Cancer: A Gene Set Enrichment Analysis" Journal of Clinical Medicine 9, no. 5: 1425. https://doi.org/10.3390/jcm9051425

APA Style

Koh, Y. W., Han, J. -H., Haam, S., & Lee, H. W. (2020). HIP1R Expression and Its Association with PD-1 Pathway Blockade Response in Refractory Advanced NonSmall Cell Lung Cancer: A Gene Set Enrichment Analysis. Journal of Clinical Medicine, 9(5), 1425. https://doi.org/10.3390/jcm9051425

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