3.1 PLCXD2 mRNA expression in HNSC tissues
Here, we assessed the feasibility of PLCXD2 as a biomarker for HNSC. As anticipated, we found that PLCXD2 mRNA expressions in HNSC tissue were significantly elevated compared to non-cancerous tissues, with a notable statistical difference (log2FC < 1, P < 0.01) (Figure 1A). Additionally, high PLCXD2 mRNA expressions were strongly associated with poorer prognoses (HR = 1.3, P = 0.036) in patients with HNSC (Figure 1B).
3.2 Association between PLCXD2 protein levels and immune markers in HNSC tissues TME
We conducted mIHC staining on HNSC TMA sections to analyze the relationship between PLCXD2 protein expression and various TIICs, including CD4+ T cells, CD8+ T cells, CD56+ NK cells, CD66b+ neutrophils, CD68+CD86+ M1-like macrophages, CD68+CD163+ M2-like macrophages and Lamp3+ dendritic cells, and also examined their correlation with commonly used immune checkpoints PD-1, PD-L1, and CTLA-4.
PLCXD2 protein was primarily localized in the nuclear and cytoplasmic compartments of both HNSC tissues and benign squamous epithelial tissues. Notably, PLCXD2 protein expression varied across different cell types, being observed not only in cancer cells and non-tumor squamous epithelial cells but also in stromal cells. Wilcoxon rank sum test revealed that PLCXD2 protein expression was significantly higher in cancer cells (10.58 ± 14.07) compared to squamous epithelial cells (4.16 ± 5.24) (Z = -3.890, P < 0.001). Additionally, PLCXD2 protein expression in cancer stromal cells (21.51 ± 21.39) was significantly higher than in benign stromal cells (12.14 ± 13.24) (Z = -3.890, P = 0.007) (Figure 2A -B). These results aligned with the PLCXD2 mRNA expression findings from the GEPIA data. Furthermore, the Wilcoxon Signed Ranks Test demonstrated that PLCXD2 protein expression in cancer cells was significantly lower than in cancer stromal cells (Z = -6.498, P < 0.001).
We analyzed PLCXD2 protein expression in the TME of HNSC using Pearson's test and observed a correlation between PLCXD2 protein levels in cancer nests and that in cancer stromal cells (r = 0.216, P < 0.001). Additionally, PLCXD2 protein levels in cancer cells were positively correlated with CD4+ T cells (CD3+CD4+) (r = 0.236, P < 0.001), CD8+ T cells (CD3+CD8+) in the cancer stroma (r = 0.313, P < 0.001), and dendritic cells (Lamp3+) in both tumor intratumoral (r = 0.204, P = 0.008) and stroma (r = 0.212, P = 0.005) (Figure 2C, Figure 3A, Table S2-1). Furthermore, we observed CTLA-4 protein expression not only in TIICs but also in some cancer cells. PLCXD2 protein levels in cancer cells were correlated with CTLA4+ cells in cancer nests (r = 0.262, P < 0.001). However, no significant correlation was found between PLCXD2 protein expression in cancer stromal cells and TIICs or immune checkpoints (PD-1, PD-L1, and CTLA-4) in the TME of HNSC (Figure 2D, Figure 3B, Table S2-2).
Higher expression and lower expression of PLCXD2 protein in HNSC (A) compared to non-tumor squamous epithelial tissue (B). PLCXD2 protein expression is shown in purple. C: Immunostaining for CD3+ cells, CD4+ cells, and CD8+ cells in HNSC tissue. D: Immunostaining for Lamp3+ cells and CTLA-4+ cells in HNSC tissue. Cytokeratin staining is shown in green, and nuclei are stained with DAPI (blue).
3.3 Association between PLCXD2 protein levels and clinicopathological characteristics of HNSC patients
Based on the OS of HNSC patients, we determined a cutoff point for PLCXD2 protein expression in cancer cells (score 0 – 91.97). Scores equal to or less than 8.57 were classified as low-expression (n = 167), while scores ranging from 8.57 to 100 were categorized as high-expression (n = 108) groups. Chi-square test results confirmed that PLCXD2 protein expression in cancer cells was associated with age (χ2 = 4.331, P = 0.037), tumor location (χ2 = 12.261, P = 0.002), depth of tumor invasion (χ2 = 13.369, P = 0.010), lymph node metastasis (LNM) (χ2 = 5.041, P = 0.025), and TNM classification (χ2 = 15.543, P = 0.004) (Table 1).
However, no suitable cutoff value was identified to establish a relationship between the score of PLCXD2 protein expression in stromal cells (ranging from 0 to 92.59) and clinical outcomes. Even when considering the median or half of the population, no evidence was found to suggest a correlation between PLCXD2 protein expression in stromal cells and clinical features of HNSC (data not shown).
Table 1. The association between PLCXD2 protein expression in cancer cells and clinical characteristics of HNSC
Groups
|
n
|
PLCXD2 protein in cancer cells
|
χ2
|
p value
|
Low expression (%)
|
High expression (%)
|
Total
|
275
|
167(60.73)
|
108(39.27)
|
|
|
Age(years)
|
|
|
|
4.331
|
0.037*
|
≤ 60
|
118
|
80 (67.80)
|
38 (32.20)
|
|
|
> 60
|
157
|
87 (55.41)
|
70 (44.59)
|
|
|
Gender
|
|
|
|
0.677
|
0.411
|
Male
|
105
|
67 (63.81)
|
38 (36.19)
|
|
|
Female
|
170
|
100 (58.82)
|
70 (41.18)
|
|
|
Location
|
|
|
|
12.261
|
0.002*
|
Tongue
|
113
|
81 (71.68)
|
32 (28.32)
|
|
|
Buccal mucosa
|
86
|
51 (59.30)
|
35 (40.70)
|
|
|
Larynx
|
76
|
35 (46.05)
|
41 (53.95)
|
|
|
Differentiation
|
|
|
|
1.134
|
0.567
|
Well
|
150
|
94 (62.67)
|
56 (37.33)
|
|
|
Middle
|
115
|
66 (57.39)
|
49 (42.61)
|
|
|
Poor
|
10
|
7 (70.00)
|
3 (30.00)
|
|
|
T stage
|
|
|
|
13.369
|
0.010*
|
Tis
|
14
|
12 (85.71)
|
2 (14.29)
|
|
|
I
|
74
|
51 (68.92)
|
23 (31.08)
|
|
|
II
|
123
|
72 (58.54)
|
51 (41.46)
|
|
|
III
|
45
|
24 (53.33)
|
21 (46.67)
|
|
|
IV
|
7
|
1 (14.29)
|
6 (85.71)
|
|
|
Unknown
|
12
|
|
|
|
|
Lymph node metastasis
|
|
|
|
5.041
|
0.025*
|
Yes
|
51
|
24 (47.06)
|
27 (52.94)
|
|
|
No
|
212
|
136 (64.15)
|
76 (35.85)
|
|
|
Unknown
|
12
|
|
|
|
|
TNM stage
|
|
|
|
15.543
|
0.004*
|
Stage 0
|
18
|
14 (77.78)
|
4 (22.22)
|
|
|
Stage I
|
60
|
42 (70.00)
|
18 (30.00)
|
|
|
Stage II
|
94
|
59 (62.77)
|
35 (37.23)
|
|
|
Stage III
|
82
|
44 (53.66)
|
38 (46.34)
|
|
|
Stage IV
|
9
|
1 (11.11)
|
8 (88.89)
|
|
|
Unknown
|
12
|
|
|
|
|
* p < 0.05
3.4 The relation between the PLCXD2 protein expression and HNSC patient prognosis
To determine the relationship between PLCXD2 protein status and 5-year OS in patients, Cox regression analysis and K-M curves with log-rank tests were conducted. Univariate analysis revealed that PLCXD2 protein expression (HR: 2.232, P = 0.001), tumor differentiation (HR: 1.532, P = 0.032), depth of tumor invasion (HR: 1.482, P = 0.005), LNM (HR: 3.314, P < 0.001) and TNM classification (HR: 1.754, P < 0.001) significantly influenced patient survival. In multivariate analyses, PLCXD2 protein expression (HR: 1.955, P = 0.010) and TNM classification (HR: 1.617, P = 0.001) were identified as independently associated with the 5-year OS of HNSC patients (Table 2, Figure 4).
Table 2.Univariate and multivariate analysis of factors associated with HNSC patient survival.
|
Univariate analysis
|
Multivariate analysis
|
HR
|
P >|z|
|
95% CI
|
HR
|
P >|z|
|
95% CI
|
PLCXD2 protein expression in cancer cells
High vs low
|
2.332
|
0.001*
|
1.426
|
3.814
|
1.955
|
0.010*
|
1.175
|
3.253
|
Age (years)
≤60 vs >60
|
1.587
|
0.078
|
0.949
|
2.653
|
|
|
|
|
Gender
Male vs Female
|
1.441
|
0.180
|
0.844
|
2.460
|
|
|
|
|
Location
Tongue vs Buccal mucosa vs
|
1.125
|
0.432
|
0.839
|
1.508
|
|
|
|
|
Differentiation
Well and Middle vs Poor
|
1.532
|
0.032*
|
1.038
|
2.261
|
1.475
|
0.054
|
0.994
|
2.188
|
T
Tis vs T2 vs T3 vs T4
|
1.482
|
0.005*
|
1.125
|
1.953
|
|
|
|
|
Lymph node metastasis
Yes vs No
|
3.314
|
<0.001*
|
1.993
|
5.510
|
|
|
|
|
TNM stage
0 vs I vs II vs III and IV
|
1.754
|
<0.001*
|
1.333
|
2.309
|
1.617
|
0.001*
|
1.219
|
2.146
|
*p<0.05;