The CLEC10A mRNA expression in different cancers
We analyzed the mRNA expression of CLEC10A using the Xiantao database and TIMER website. The results showed that CLEC10A was significantly lower in most cancer tissues, than in corresponding normal tissues (Fig. 1A). In addition, we used the UALCAN database to validate the findings in the TIMER website and reported lower expression of CLEC10A in paired and nonpaired BC tissues than in normal tissues (Fig. 1B, C). In addition, CLEC10A expression was associated with pathologic stage, histological type, T stage and HER2 status (Fig. 1D-F, K), but not with N stage, M stage, ER status or PR status (Fig. 1G-J).
Prognostic Significance Of Clec10a Expression In Bc
We investigated the Kaplan-Meier plotter database for the prognostic significance of
CLEC10A expression in BC. Low levels of CLEC10A predicted poor prognosis in BC (Fig. 2A, C). As Kaplan-Meier plotter analyses only OS, disease-specific survival (DSS) and progression-free intervals (PFI) values, we assessed the multiple clinical prognostic value of CLEC10A in a variety of cancers by R project using the “survival” package. Forest plot showed CLEC10A as a risk factor for different prognoses in BC (Fig. 3A, C). These findings indicated that CLEC10A is a preventative factor in BC.
Clec10a Expression Is Correlated With Immune Infiltration In Bc
Tumor-infiltrating lymphocytes can independently predict sentinel lymphnode status and prognosis in various cancers [8]. Therefore, we used TIMER to analyze the correlation of CLEC10A levels with immune infiltration levels in BC. The results showed that CLEC10A expression was significantly positively correlated with the infiltration of B cells (r = 0.257, p = 3.65e-16), CD4 + T cells (r = 0.577, p = 1.21e-86), CD8 + T cells (r = 0.439, p = 3.42e-47), macrophages (r = 0.117, p = 2.30e-04), dendritic cells (r = 0.531, p = 2.04e-70) and neutrophils (r = 0.402, p = 2.94e-38) in BC (Fig. 4A). Then, we used the TISIDB database to further explore the relationship between CLEC10A levels and 28 tumor immune infiltrating cell subtypes. These results showed that CLEC10A was associated with all immune cell subtypes in BC (Fig. 4B, Table 1). In particular, activated CD8 T cells (Act_CD8) (r = 0.592, p < 2.2e-16), central memory CD8 + T cells (Tcm_CD8) (r = 0.134, p < 8.571e-6), activated CD4 T cells (Act_CD4) (r = 0.385, p < 2.2e-16), central memory CD4 + T cells (Tcm_CD4) (r = 0.296, p < 7.3e-24), type 2 T helper cells (r = 0.315, p < 2.2e-16), type 1 T helper cells (r = 0.69, p < 2.2e-16) and CLEC10A were moderately correlated with CLEC10A expression (Fig. 4C-H). Notably, macrophages (r = 0.585, p < 2.2e-16) also displayed a moderate correlation with CLEC10A expression (Fig. 4I). These results strongly implicate CLEC10A as a major tumor immune infiltration regulator in BC.
Clec10a Expression Is Associated With Immune Cell Type Markers
We assessed the correlation between CLEC10A expression and tumor-infiltrating immune cell gene marker levels in BC tissues by exploring the TIMER database. Our results showed that the CLEC10A level in BC tissues was strongly associated with immune markers of B cells, CD8 T cells, M2 macrophages, monocytes, natural killer (NK) cells, T cells (general), dendritic cells (DCs), T helper cells, Tregs, and T exhaustion cells (Table 1). Notably, we found that the CLEC10A level was significantly correlated with the levels of various subtypes of T cell markers, including CD8 + T markers, CD8A and CD8B; T cell (general) markers, CD3D, CD3E, and CD2; exhausted T cell markers CTLA4, HAVCR2, GZMB, LAG-3, and PD-1; Th2 markers; Th17 markers STAT6, STAT5A, and IL17A; Treg markers FOXP3, CCR8, STAT5B, and TGFB1; Tfh marker BCL6; and neutrophil markers ITGAM and CCR7,etc(Table 2). Significant correlations of CLEC10A levels with different macrophage markers (M2 macrophage markers MS4A4A, VSIG4 and CD163; monocyte markers CSF1R and CD86; tumor-associated macrophage (TAM) markers CD68, IL21, IL10, and CCL2; and B cell markers CD19 and CD79A) (Table 2). Furthermore, the expression of CLEC10A was not markedly related to marker genes for DC markers, CD1C, NOS2, CEACAM8, Th2, GATA3, Th17 cells and STATA3 in BC. These findings reveal that CLEC10A is involved in the regulation of tumor immune infiltration in BC.
Table 1 The correlation between CLEC10A expression and tumor lymphocyte infiltration in human cancer (TISIDB).
|
r
|
p
|
Activated CD8 T cell (Act _CD8)
|
0.59
|
<2.2e-16
|
Central memory CD8 T cell (Tcm _CD8)
|
0.134
|
8.57E-06
|
Effector memory CD8 T cell (Tem _CD8)
|
0.701
|
<2.2e-16
|
Activated CD4 T cell (Act _CD4)
|
0.382
|
<2.2e-16
|
Central memory CD4T cell (Tcm _CD4)
|
0.296
|
7.30E-24
|
Effector memory CD4 T cell (Tem _CD4)
|
0.392
|
<2.2e-16
|
T follicular helper cell(Tfh)
|
0.607
|
<2.2e-16
|
Gamma delta T cell (Tgd)
|
0.434
|
<2.2e-16
|
Type 1 T helper cell (Th1)
|
0.69
|
<2.2e-16
|
Type 17 T helper cell (Th17)
|
0.526
|
<2.2e-16
|
Type 2 T helper cell (Th2)
|
0.315
|
<2.2e-16
|
Regulatory T cell (Treg)
|
0.504
|
<2.2e-16
|
Activated B cell (Act _B)
|
0.782
|
<2.2e-16
|
Immature B cell (Imm _B)
|
0.688
|
<2.2e-16
|
Memory B cell (Mem _B)
|
0.381
|
<2.2e-16
|
natural killer cell (NK)
|
0.464
|
<2.2e-16
|
CD56bright natural killer cell (CD56bright)
|
0.388
|
<2.2e-16
|
CD56dim natural killer cell (CD56dim)
|
0.209
|
3.02E-10
|
Myeloid derived suppressor cell (MDSC)
|
0.566
|
<2.2e-16
|
Natural killer T cell (NKT)
|
0.576
|
<2.2e-16
|
Activated dendtritic cell (Act _DC)
|
0.326
|
<2.2e-16
|
Plasmacytoid dendtritic cell (pDC)
|
0.353
|
<2.2e-16
|
Immature dendtritic cell (iDC)
|
0.098
|
0.00111
|
Macrophage (Macrophage)
|
0.585
|
<2.2e-16
|
Eosinophi (Eosinophil)
|
0.601
|
<2.2e-16
|
Mast (Mast)
|
0.685
|
<2.2e-16
|
Monocyte (Monocyte)
|
0.483
|
<2.2e-16
|
Neutrophil (Neutrophil)
|
0.359
|
<2.2e-16
|
Table2 Correlation analysis between CLEC10A and relate genes and markers of immune cells in TIMER
Description
|
Gene markers
|
None
|
|
Purity
|
|
|
|
COR
|
P
|
COR
|
P
|
B CELL
|
CD19
|
0.655
|
***
|
0.532
|
***
|
|
CD79A
|
0.65
|
***
|
0.511
|
***
|
CD8 T CELL
|
CD8A
|
0.726
|
***
|
0.507
|
***
|
|
CD8B
|
0.677
|
***
|
0.565
|
***
|
DENDRITIC CELL
|
ITGAX
|
0.492
|
***
|
0.343
|
***
|
|
NRP1
|
0.273
|
***
|
0.121
|
***
|
|
CD1C
|
0.876
|
0.923
|
0.819
|
***
|
|
HLA-DPA1
|
0.642
|
***
|
0.511
|
***
|
|
HLA-DRA
|
0.648
|
***
|
0.514
|
***
|
|
HLA-DQB1
|
0.535
|
***
|
0.411
|
***
|
|
HLA-DPB1
|
0.715
|
***
|
0.586
|
***
|
M1 MACROPHAGE
|
PTGS2
|
0.415
|
***
|
0.284
|
***
|
|
IRF5
|
0.307
|
***
|
0.226
|
***
|
|
NOS2
|
0.011
|
0.72
|
0.022
|
0.493
|
M2 MACROPHAGE
|
MS4A4A
|
0.554
|
***
|
0.43
|
***
|
|
VSIG4
|
0.399
|
***
|
0.272
|
***
|
|
CD163
|
0.44
|
***
|
0.327
|
***
|
MONOCYTE
|
CSF1R
|
0.549
|
***
|
0.444
|
***
|
|
CD86
|
0.475
|
***
|
0.338
|
***
|
NATURAL KILLER CELL
|
KIR2DS4
|
0.332
|
***
|
0.234
|
***
|
|
KIR3DL3
|
0.235
|
***
|
0.181
|
***
|
|
KIR3DL2
|
0.469
|
***
|
0.356
|
***
|
|
KIR3DL1
|
0.432
|
***
|
0.34
|
***
|
|
KIR2DL4
|
0.388
|
***
|
0.279
|
***
|
|
KIR2DL3
|
0.356
|
***
|
0.266
|
***
|
|
KIR2DL1
|
0.34
|
***
|
0.255
|
***
|
NEUTROPHILS
|
CCR7
|
0.792
|
***
|
0.52
|
***
|
|
ITGAM
|
0.392
|
***
|
0.338
|
***
|
|
CEACAM8
|
0.041
|
0.178
|
0.005
|
0.886
|
T CELL GENERAL
|
CD3D
|
0.771
|
***
|
0.549
|
***
|
|
CD3E
|
0.782
|
***
|
0.56
|
***
|
|
CD2
|
0.735
|
***
|
0.528
|
***
|
T CELL EXHAUSTION
|
CTLA4
|
0.57
|
***
|
0.44
|
***
|
|
LAG3
|
0.387
|
***
|
0.286
|
***
|
|
HAVCR2
|
0.393
|
***
|
0.252
|
***
|
|
GZMB
|
0.625
|
***
|
0.508
|
***
|
|
PDCD1
|
0.675
|
***
|
0.552
|
***
|
TAM
|
CCL2
|
0.472
|
***
|
0.338
|
***
|
|
IL10
|
0.455
|
***
|
0.334
|
***
|
|
CD68
|
0.42
|
***
|
0.288
|
***
|
|
BCL6
|
0.126
|
***
|
0.083
|
**
|
|
IL21
|
0.367
|
***
|
0.279
|
***
|
TH1
|
TBX21
|
0.745
|
***
|
0.643
|
***
|
|
STAT4
|
0.686
|
***
|
0.554
|
***
|
|
STAT1
|
0.255
|
***
|
0.184
|
***
|
|
IFNG
|
0.52
|
***
|
0.404
|
***
|
|
IL13
|
0.273
|
***
|
0.216
|
***
|
TH2
|
GATA3
|
-0.198
|
***
|
-0.061
|
0.056
|
|
STAT6
|
0.249
|
***
|
0.223
|
***
|
|
STAT5A
|
0.424
|
***
|
0.31
|
***
|
TH17
|
STAT3
|
0.083
|
***
|
0.038
|
0.229
|
|
IL17A
|
0.203
|
***
|
0.102
|
***
|
TREG
|
FOXP3
|
0.525
|
***
|
0.384
|
***
|
|
CCR8
|
0.42
|
***
|
0.333
|
***
|
|
STAT5B
|
0.234
|
***
|
0.207
|
***
|
|
TGFB1
|
0.411
|
***
|
0.263
|
***
|
*p < 0.05, **p < 0.01, ***p < 0.001. |
Prognostic potential of CLEC10A in different tumors based on immune cells
This study showed that the CLEC10A level was associated with the immune infiltration of BC. Additionally, downregulated CLEC10A has a poor prognosis in BC patients. Thus, we propose a hypothesis that CLEC10A may affect the prognosis of BC patients partly through immune infiltration. We performed Kaplan-Meier plotter analyses of CLEC10A expression in BC and considered the infiltration level os B cells, CD4 + memory T cells, CD8 + T cells, macrophages, neutrophils, NK T cells, mesenchymal stem cells, regulatory T cells, Th1 cells, and Th2 cells. We found that BC samples with high CLEC10A levels were enriched in CD4 + memory T cells (p = 0.04), neutrophils (p = 0.032), macrophages (p = 0.002), natural killer T cells (p = 0.017), Treg cells (p = 0.04), and Th2 cells (p = 0.0018) (Fig. 5B, D, E, G, H, J). Unfortunately, there was no significant difference between the high and low CLEC10A expression groups in terms of OS for patients with high infiltration of B cells (p = 0.15), CD8 + T cells (p = 0.34), mesenchymal stem cells (p = 0.06), and Th1 cells (p = 0.1) (Fig. 5A, C, F, I). The above analysis suggested that immune infiltration may affect the prognosis of BC patients with low CLEC10A expression.