Bei Wang1, Limiao Zhao1, Dandan Chen2. 1. Institute of Integration of Traditional Chinese and Western Medicine, the Hospital Affiliated to Jiangnan University, Wuxi, People's Republic of China. 2. Department of Education, School of Humanities, Jiangnan University, Wuxi, People's Republic of China.
Abstract
Background: Coactosin-like protein (COTL1) was first identified as protein that binds 5-lipoxygenase and F-actin; its functions in tumors remain unknown. COTL1 could inhibit the proliferation of breast cancer (BRCA) in vivo and in vitro; however, online public databases including UALCAN and Kaplan-Meier plotter showed high COTL1 expression in breast cancer tissue, which was correlated with poor prognosis. Therefore, we studied the role of COTL1 expression in human breast cancer and its use in determining clinical prognosis. Methods: We first used the UALCAN database and immunohistochemical analysis to elucidate COTL1 expression in BRCA. We then performed Kaplan-Meier plotting and immunohistochemical analysis to assess prognosis in BRCA in relation to COTL1 expression. Finally, we used the CancerSEA and LinkedOmics databases to evaluate the function of COTL1 in BRCA. The TIMER and TISIDB databases were used to evaluate the association between COTL1 expression and immune cell infiltration in BRCA. Results: UALCAN and immunohistochemical analysis showed that COTL1 was highly expressed in breast cancer. Furthermore, high COTL1 expression was correlated with poor prognosis in BRCA. We also found that COTL1 is involved in immune response via the CancerSEA and LinkedOmics databases. The TIMER and TISIDB databases showed that high COTL1 expression was correlated with immune cell infiltration. Conclusion: COTL1 expression was higher in breast cancer tissues than in normal tissues, and high COTL1 expression was correlated with poor prognosis and immune cell infiltration. These results provide a basis for the development of applications of COTL1 in determining the prognosis of breast cancer and its treatment.
Background: Coactosin-like protein (COTL1) was first identified as protein that binds 5-lipoxygenase and F-actin; its functions in tumors remain unknown. COTL1 could inhibit the proliferation of breast cancer (BRCA) in vivo and in vitro; however, online public databases including UALCAN and Kaplan-Meier plotter showed high COTL1 expression in breast cancer tissue, which was correlated with poor prognosis. Therefore, we studied the role of COTL1 expression in human breast cancer and its use in determining clinical prognosis. Methods: We first used the UALCAN database and immunohistochemical analysis to elucidate COTL1 expression in BRCA. We then performed Kaplan-Meier plotting and immunohistochemical analysis to assess prognosis in BRCA in relation to COTL1 expression. Finally, we used the CancerSEA and LinkedOmics databases to evaluate the function of COTL1 in BRCA. The TIMER and TISIDB databases were used to evaluate the association between COTL1 expression and immune cell infiltration in BRCA. Results: UALCAN and immunohistochemical analysis showed that COTL1 was highly expressed in breast cancer. Furthermore, high COTL1 expression was correlated with poor prognosis in BRCA. We also found that COTL1 is involved in immune response via the CancerSEA and LinkedOmics databases. The TIMER and TISIDB databases showed that high COTL1 expression was correlated with immune cell infiltration. Conclusion: COTL1 expression was higher in breast cancer tissues than in normal tissues, and high COTL1 expression was correlated with poor prognosis and immune cell infiltration. These results provide a basis for the development of applications of COTL1 in determining the prognosis of breast cancer and its treatment.
Coactosin-like protein (COTL1) was first identified as a binding protein of 5-lipoxygenase and F-actin;1–3 its functions in tumors remain unclear.4–7 Xia et al recently reported that COTL1 could inhibit the proliferation of breast cancer in vivo and in vitro.6 Previous studies have focused on the function of COTL1 in cell and animal models, but not in humans. Online public databases including UALCAN and Kaplan-Meier plotter show that COTL1 is highly expressed in breast cancer tissue and high COTL1 expression correlates with poor prognosis. Therefore, we aimed to determine the role of COTL1 expression in human breast cancer, and the related clinical prognosis.Breast cancer is one of the most common malignancies among women, and is a major cause of cancer-related deaths worldwide.8 The immune system plays two roles in breast cancer. It can promote tumor formation through inflammatory pathways, and can suppress tumor formation through active immune surveillance.9 In this study, we used immunohistochemistry and public online databases including UALCAN and Kaplan-Meier plotter to evaluate COTL1 expression and related prognosis in breast cancer. In addition, we explored the correlation between COTL1 expression and immune cell infiltration in breast cancer by using the CancerSEA, LinkedOmics, and TISIDB public online databases. Our results could provide a basis for the clinical application of COTL1 in breast cancer in the future.
Materials and Methods
Immunohistochemistry
We used a tissue microarray (obtained from XinChao Biological Technology Co., Ltd), which included 170 samples of breast cancer tissues and sixty-two samples from para-carcinoma tissues. Before dewaxing, the tissue microarray was incubated for 30 minutes at 60 °C. The slides were repaired with citrate antigen retrieval solution (Beyotime, #P0083). After cooling down, non-specific binding was blocked with endogenous peroxidase blocking buffer (Beyotime, #P0100A) and immunol staining blocking buffer (Beyotime, #P0102). Primary antibody (COTL1 (1: 400, Proteintech, #17119-1-AP)) was added dropwise and incubated overnight at 4°C. The increasing agent was added dropwise and the samples were placed in a humidified incubator at room temperature for 20 min. The enzyme-labeled secondary antibody (MXB, #KIT 9902) was added dropwise and the samples were placed in a humidified incubator at room temperature for 30 min. DAB horseradish peroxidase color development kit (Beyotime, #P0203) and hematoxylin contrast (Beyotime, #C0107) were used for staining. Routine alcohol gradient dehydration, transparent xylene, and neutral resin mounting were performed. The negative control was exposed to PBS instead of the primary antibody. Staining was assessed by scanning the entire tissue specimens. An immunoreactivity score system was adopted based on the scale and the intensity of positively-stained target cells. The following standards were used: 1) the number of cells with positive staining (≤ 5% = 0; 6% - 25% = 1; 26% - 50% = 2; 51% - 75% = 3; and > 75% = 4); 2) the staining intensity (colorless = 0; pallide-flavens = 1; yellow = 2; brown = 3). We multiplied the scores of 1 and 2 in the standard, and the dyeing grade was divided into none (0 points), weak (1–4 points), medium (5–8 points), or strong (9–12 points). Immunohistochemistry staining and evaluation were performed as described previously.10 This work was approved by the affiliated hospital of Jiangnan University (#LS2021076).
UALCAN
UALCAN ()11 is a comprehensive, user-friendly, and interactive web resource for analyzing COTL1 expression in breast cancer.
Kaplan-Meier Plotter Analysis
Kaplan-Meier plotter (),12 a database that integrates gene expression data and clinical data, was used to analyze the prognostic value of COTL1 in breast cancer.
LinkedOmics
LinkedOmics (),13 was used to study the associated RNA Sequence and for enrichment analysis.
CancerSEA
CancerSEA (),14 was used to comprehensively analyze COTL1 function in BRCA at single-cell resolution.
Timer
TIMER (),15 was used to comprehensively analyze tumor immune cells and COTL1 expression in BRCA.
TISIDB
TISIDB (),16 was used to study the correlation between COTL1 expression and immune system in BRCA.
Statistical Analysis
Chi-square tests were used to analyze immunohistochemical score levels in different clinical pathological feature groups. Univariate and multivariate analyses were used to evaluate the effect of clinical variables on survival. A two-tailed p-value <0.05 was considered statistically significant.
Results
High COTL1 Expression in BRCA
To elucidate the differences in COTL1 expression in tumor and normal tissues, the COTL1 mRNA levels in breast tumors and normal tissues were analyzed using the UALCAN database. COTL1 mRNA expression was significantly higher in breast cancer tissues than in normal tissues (Figure 1A). In addition, subgroup analysis based on stage, race, gender, age, subclasses, menopause status, histologic subtype, metastasis, and TP53 mutation status indicated that COTL1 mRNA expression was significantly higher in breast cancer tissues than in normal tissues (Figure 1B–J). COTL1 expression was much higher in triple negative breast cancer compared to other subtypes. Triple negative breast cancer has a worse clinical prognosis than does non-triple negative breast cancer.17 Further, the TISIDB database showed that COTL1 expression was associated with molecular and immune subtypes (Figure 1K and L). We then studied COTL1 protein expression by using a tissue microarray, which included 170 breast cancer tissue samples and sixty-two para-carcinoma tissue samples. The result of immunohistochemical analysis was consistent with that of UALCAN analysis (Figure 1M and N). To further evaluate the correlation between COTL1 expression and clinical characteristics, COTL1 expression and clinical data were downloaded from the cancer genome atlas (TCGA). The results indicated that COTL1 expression was associated with clinical parameters such as laterality, ER, PR, and radiation therapy (Table 1). The results from our immunohistochemical analysis demonstrated that COTL1 expression was correlated with clinical parameters at T, N, and HER2. (Table 2). These contradictory findings could be due to our small sample size, or because COTL1 expression is unrelated to clinical parameters.
Figure 1
Expression of COTL1 in breast cancer. (A) COTL1 mRNA expression (from The Cancer Genome Atlas) in breast cancer was significantly higher than that in normal tissue. (B-L) Differences in COTL1 mRNA expression depending on stage, race, gender, age, subclasses, menopause status, histologic subtype, metastasis, TP53 mutation status, molecular and immune subtypes. (M-N) COTL1 protein expression in breast cancer tissues was higher than that in para-carcinoma tissues. Data are shown as mean ± SD, *<0.05, **<0.01, ***p < 0.001.
Table 1
The Correlation of COTL1 Expression and Clinical Characteristics from TCGA (n=1104)
Characteristics
N
COTL1 Low
COTL1 High
P value
Gender
Male
12
7
5
0.568
Female
1089
545
544
Age
<60
589
291
298
0.603
≥60
512
261
251
Laterality
Right
445
243
202
0.004
Left
556
253
303
Stage
I
182
86
96
0.781
II
625
313
312
III
252
131
121
IV
28
15
13
T
T1
281
134
147
0.588
T2
639
321
318
T3
138
71
67
T4
41
24
17
N
N0
515
240
275
0.12
N1
367
201
166
N2
120
59
61
N3
79
38
41
M
M0
916
475
441
0.803
M1
22
12
10
Grade
1
210
105
105
0.145
2
469
254
215
3
79
47
32
4
30
21
9
Radiation therapy
No
445
243
202
0.004
Yes
556
253
303
ER
Negative
255
54
201
2.15E-26
Positive
796
473
323
PR
Negative
259
95
164
1.47E-11
Positive
522
325
197
HER2
Negative
843
423
420
0.675
Positive
146
76
70
Notes: P value from Chi-square test; Bold values indicate p<0.05.
Abbreviations: PR, progesterone receptor; ER, Estrogen receptor; HER2, Human epidermal growth factor receptor 2.
Table 2
The Correlation of COTL1 Expression and Clinical Characteristics in Breast Cancer Patients (n=138)
Clinicopathological Characteristics
N
COTL1 High Expression
COTL1 Low Expression
P value
Age (years)
≦50
64
40
24
0.849
>50
74
44
30
Tumor size (cm)
≦5
114
69
45
0.857
>5
24
15
9
Stage
I, IIa
35
17
18
0.127
IIb, III
103
67
36
Grade
I, II
91
50
41
0.072
III
47
34
13
T
T1
30
13
17
0.044
T2-T3
108
71
37
N
N0-N1
94
51
43
0.032
N2-N3
44
33
11
M
M0
138
84
54
/
M1
0
0
0
ER
Negative
45
28
17
0.821
Positive
93
56
37
PR
Negative
56
39
17
0.081
Positive
82
45
37
HER2
Negative
90
49
41
0.034
Positive
48
35
13
Note: Bold values indicate p<0.05.
Abbreviations: PR, progesterone receptor; ER, Estrogen receptor; HER2, Human epidermal growth factor receptor 2.
The Correlation of COTL1 Expression and Clinical Characteristics from TCGA (n=1104)Notes: P value from Chi-square test; Bold values indicate p<0.05.Abbreviations: PR, progesterone receptor; ER, Estrogen receptor; HER2, Human epidermal growth factor receptor 2.The Correlation of COTL1 Expression and Clinical Characteristics in Breast Cancer Patients (n=138)Note: Bold values indicate p<0.05.Abbreviations: PR, progesterone receptor; ER, Estrogen receptor; HER2, Human epidermal growth factor receptor 2.Expression of COTL1 in breast cancer. (A) COTL1 mRNA expression (from The Cancer Genome Atlas) in breast cancer was significantly higher than that in normal tissue. (B-L) Differences in COTL1 mRNA expression depending on stage, race, gender, age, subclasses, menopause status, histologic subtype, metastasis, TP53 mutation status, molecular and immune subtypes. (M-N) COTL1 protein expression in breast cancer tissues was higher than that in para-carcinoma tissues. Data are shown as mean ± SD, *<0.05, **<0.01, ***p < 0.001.
Clinical Prognosis of COTL1 in Patients with BRCA
We determined whether COTL1 expression is correlated with prognosis in patients with breast cancer. The influence of COTL1 expression on survival rates was first evaluated using the Kaplan-Meier plotter database. The results showed that COTL1 mRNA level was positively correlated with poor prognosis (Figure 2A). We then analyzed GSE20486 and GSE12276 from the Gene Expression Omnibus (GEO) database, and found that higher COTL1 mRNA levels corresponded with poorer overall survival (OS; Figure 2B and C). To confirm the relationship between COTL1 and clinical prognosis, the overall survival analysis from our tissue microarray results demonstrated that high COTL1 expression in breast cancer tissues was associated with poor prognosis (Figure 2D). To determine whether COTL1 is an independent risk factor for overall survival in patients with breast cancer, univariate and multivariate Cox analyses were performed using the SPSS software. In these analyses, estrogen receptor (ER), and COTL1 expression were independent risk factors for OS (Table 3). Thus, COTL1 could be an independent prognostic biomarker for BRCA.
Figure 2
Clinical prognosis of COTL1 in patients with breast cancer. (A-C) Survival curve based on COTL1 expression was analyzed using the Kaplan-Meier Plotter and PrognoScan database. (D) Survival curve based on COTL1 expression was analyzed by GraphPad Prism. COTL1 expression observed by IHC staining analysis was quantified using the H-score.
Table 3
Univariate and Multivariate Analysis of the Correlation of COTL1 Expression with Overall Survival Among Breast Cancer Patients (n=139)
Parameter
Univariate Analysis
Multivariate Analysis
HR(95% CI)
P value
HR(95% CI)
P value
Age
1.45(0.78–2.69)
0.246
/
Grade
2.39(1.29–4.41)
0.005
1.87(0.99–3.54)
0.055
Stage
1.09(0.53–2.23)
0.811
/
PR
0.58(0.30–1.09)
0.093
/
ER
0.45(0.24–0.83)
0.011
0.45(0.24–0.85)
0.014
HER2
0.98(0.51–1.90)
0.960
/
COTL1
3.33(1.54–7.20)
0.002
3.02(1.38–6.61)
0.006
Note: Bold values indicate p<0.05.
Abbreviations: PR, progesterone receptor; ER, Estrogen receptor; HER2, Human epidermal growth factor receptor 2.
Univariate and Multivariate Analysis of the Correlation of COTL1 Expression with Overall Survival Among Breast Cancer Patients (n=139)Note: Bold values indicate p<0.05.Abbreviations: PR, progesterone receptor; ER, Estrogen receptor; HER2, Human epidermal growth factor receptor 2.Clinical prognosis of COTL1 in patients with breast cancer. (A-C) Survival curve based on COTL1 expression was analyzed using the Kaplan-Meier Plotter and PrognoScan database. (D) Survival curve based on COTL1 expression was analyzed by GraphPad Prism. COTL1 expression observed by IHC staining analysis was quantified using the H-score.
COTL1 Function in BRCA
COTL1 can interact with F-actin in a calcium-independent manner and influence ALOX5 stability and activity in leukotriene synthesis as a chaperone for ALOX52. We used the CancerSEA database to investigate the new function of COTL1 in BRCA at the single-cell resolution. The results indicated that COTL1 was positively associated with differentiation, inflammation, metastasis, angiogenesis, apoptosis, and quiescence and negatively associated with DNA repair, invasion, stemness, and DNA damage (Figure 3A). According to data from Braune EB (n=369) and Jordan NV(n=70), COTL1 expression was positively correlated with angiogenesis, differentiation, inflammation and quiescence18,19 (Figure 3B and C). To elucidate the function of COTL1 in BRCA, the LinkedOmics database was used, and the results showed that COTL1 was positively and negatively correlated with the top fifty significant genes (Figure 4A). Results of gene set enrichment analysis (GSEA), including details of biological process, cellular components, and molecular function of COTL1 are shown in Figure 4B. Overrepresentation enrichment analysis (ORA) showed that COTL1 is involved in and regulates multiple immune responses (Figure 4C). Specifically, COTL1 increased interferon-gamma and interleukin-1 production and activated natural killer cells (Figure 4D). In addition, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway showed that COTL1 is involved in multiple signal pathways such as the IL-17, TNF, and B cell receptor signaling pathways (Figure 4E). These results suggest that COTL1 plays a role in immune processes in breast cancer.
Figure 3
Function of COTL1 in breast cancer determined using the CancerSEA database. (A) Analysis from the CancerSEA database at single-cell resolution indicated that COTL1 was primarily involved in differentiation, inflammation, metastasis, and angiogenesis. (B and C) According to data from Braune EB and Jordan NV, COTL1 expression was significantly positively correlated with angiogenesis, differentiation, inflammation, and quiescence. ***p < 0.001, **p < 0.01.
Figure 4
Function of COTL1 in breast cancer from LinkedOmics. (A) Top 50 positively and negatively regulated genes highly co-expressed with COTL1 in TCGA from the LinkedOmics database. (B) Gene set enrichment analysis (GSEA) showed the biological processes, cellular components, and molecular functions of COTL1 in BRCA. (C) Overrepresentation enrichment analysis (ORA) demonstrated that COTL1 is involved in and regulates multiple immune response. (D) COTL1 increased the immune associated factors. E. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway showed that COTL1 is involved in multiple signal pathways such as the IL-17, TNF, and B cell receptor signaling pathways.
Function of COTL1 in breast cancer determined using the CancerSEA database. (A) Analysis from the CancerSEA database at single-cell resolution indicated that COTL1 was primarily involved in differentiation, inflammation, metastasis, and angiogenesis. (B and C) According to data from Braune EB and Jordan NV, COTL1 expression was significantly positively correlated with angiogenesis, differentiation, inflammation, and quiescence. ***p < 0.001, **p < 0.01.Function of COTL1 in breast cancer from LinkedOmics. (A) Top 50 positively and negatively regulated genes highly co-expressed with COTL1 in TCGA from the LinkedOmics database. (B) Gene set enrichment analysis (GSEA) showed the biological processes, cellular components, and molecular functions of COTL1 in BRCA. (C) Overrepresentation enrichment analysis (ORA) demonstrated that COTL1 is involved in and regulates multiple immune response. (D) COTL1 increased the immune associated factors. E. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway showed that COTL1 is involved in multiple signal pathways such as the IL-17, TNF, and B cell receptor signaling pathways.
Correlations Between COTL1 and Immune Cell Infiltration in BRCA
Since COTL1 expression was associated with immunity, we investigated the correlation between COTL1 expression and immune cell infiltration in BRCA. The TIMER database was used to assess the correlation between COTL1 expression and the infiltration levels of six immune cells. COTL1 expression was strongly positively correlated with immune cell infiltration (Purity: cor = 0.415, p = 1.2e-42; B Cell: cor = 0.435, p = 2.23e-46; CD8+ T cell: cor = 0.283, p = 1.83e-71; CD4+T Cell: cor = 0.532, p = 2.08e-71; Macrophage: cor = 0.156, p = 8.84e-07; Neutrophil: cor = 0.523, p = 9.34e-68, Dendritic cell: cor = 0.622, p = 4.66–103; Figure 5). To confirm these results, we used the TISIDB database to determine the correlation between COTL1 expression and twenty-eight immune cells in BRCA. COTL1 expression in BRCA tissues was significantly correlated with immune cells including helper T cells, NK cells, NKT cells, and regulatory T cells in addition to the aforementioned six immune cell types (Figure 6A). Furthermore, TIMER was utilized to determine the correlation between COTL1 expression and immune marker genes of different immune cells including CD8+T cells, T cells, B cells, monocytes, TAMs, M1 and M2 macrophages, neutrophils, NK cells, DCs, Th1, Th2, Th17, Treg, and T cell exhaustion in BRCA. COTL1 expression was strongly correlated with 94.8% immune cell gene markers (Table 4). The TISIDB database was then utilized to determine the correlation between COTL1 expression and immune inhibitors. COTL1 expression was significantly correlated with immune inhibitors. The T cell exhaustion markers: CTLA4, LAG3, PD1, and PDCD1 were positively correlated with COTL1 (Figure 6B), which was consistent with the TIMER data (Table 4). These results strongly suggest that COTL1 is closely correlated with immune cell infiltration in BRCA.
Figure 5
Correlation of COTL1 expression with tumor immune cell infiltration (purity, B cell, CD8+ T cell, CD4+ T cell, macrophage, neutrophil, and DCs) in the TIMER database.
Figure 6
Correlation of COTL1 expression and lymphocyte and immune inhibitors in the TISIDB database. (A) Correlation between COTL1 expression and lymphocyte in BRCA. (B) Correlation between COTL1 expression and immune inhibitors in BRCA.
Table 4
Correlation Analysis Between COTL1 and Relate Genes and Markers of Immune Cells in BRCA by TIMER
Gene Markers
BRCA
None
Purity
Cor
P
Cor
P
CD8+T cell
CD8A
0.52
***
0.418
***
CD8B
0.582
***
0.499
***
T cell(general)
CD3D
0.61
***
0.516
***
CD3E
0.597
***
0.5
***
CD2
0.591
***
0.498
***
B cell
CD19
0.502
***
0.379
***
CD79A
0.489
***
0.353
***
Monocyte
CD86
0.58
***
0.514
***
CD115 (CSF1R)
0.557
***
0.474
***
TAM
CCL2
0.482
***
0.382
***
CD68
0.507
***
0.438
***
IL10
0.442
***
0.352
***
M1 Macrophage
INOS(NOS2)
0.173
***
0.158
***
IRF5
0.358
***
0.318
***
COX2(PTGS2)
0.346
***
0.237
***
M2 Macrophage
CD163
0.476
***
0.408
***
VSIG4
0.401
***
0.318
***
MS4A4A
0.471
***
0.375
***
Neutrophils
CD66b(CEACAM8)
0.102
***
0.118
***
CD11B(ITGAM)
0.429
***
0.35
***
CCR7
0.527
***
0.411
***
Natural killer cell
KIR2DL1
0.312
***
0.247
***
KIR2DL3
0.315
***
0.234
***
KIR2DL4
0.403
***
0.334
***
KIR3DL1
0.363
***
0.283
***
KIR3DL2
0.405
***
0.317
***
KIR3DL3
0.243
***
0.202
***
KIR2DS4
0.302
***
0.233
***
Dendritic cell
HLA-DPB1
0.599
***
0.507
***
HLA-DQB1
0.519
***
0.425
***
HLA-DRA
0.578
***
0.489
***
HLA-DPA1
0.533
***
0.431
***
BDCA-1(CD1C)
0.46
***
0.322
***
BDCA-4(NRP1)
0.217
***
0.109
***
CD11c(ITGAX)
0.572
***
0.487
***
Th1
T-bet(TBX21)
0.576
***
0.477
***
STAT4
0.523
***
0.408
***
STAT1
0.291
***
0.25
***
IFN-γ(IFNG)
0.513
***
0.435
***
TNF-α(TNF)
0.417
***
0.371
***
Th2
GATA3
−0.516
***
−0.464
***
STAT6
−0.027
0.362
−0.073
*
STAT5A
0.309
***
0.223
***
IL13
0.263
***
0.209
***
Tfh
BCL6
0.009
0.765
−0.039
0.215
IL21
0.347
***
0.281
***
Th17
STAT3
0.002
0.945
−0.021
0.499
IL17A
0.215
***
0.159
***
Treg
FOXP3
0.563
***
0.476
***
CCR8
0.444
***
0.387
***
STAT5B
−0.014
0.641
−0.059
0.063
TGFβ(TGFB1)
0.32
***
0.201
***
T cell exhaustion
PD-1(PDCD1)
0.582
***
0.492
***
CTLA4
0.594
***
0.515
***
TIM-3(HAVCR2)
0.511
***
0.439
***
GZMB
0.582
***
0.5
***
LAG3
0.512
***
0.457
***
PDL1(CD274)
0.35
***
0.262
***
Notes: *P < 0.01; *** P < 0.0001.
Abbreviations: TAM, tumor-associated macrophage; Th, T helper cell; Tfh, Follicular helper T cell; Treg, regulatory T cell; Cor, R value of Spearman correlation; None, correlation without adjustment. Purity, correlation adjusted by purity.
Correlation Analysis Between COTL1 and Relate Genes and Markers of Immune Cells in BRCA by TIMERNotes: *P < 0.01; *** P < 0.0001.Abbreviations: TAM, tumor-associated macrophage; Th, T helper cell; Tfh, Follicular helper T cell; Treg, regulatory T cell; Cor, R value of Spearman correlation; None, correlation without adjustment. Purity, correlation adjusted by purity.Correlation of COTL1 expression with tumor immune cell infiltration (purity, B cell, CD8+ T cell, CD4+ T cell, macrophage, neutrophil, and DCs) in the TIMER database.Correlation of COTL1 expression and lymphocyte and immune inhibitors in the TISIDB database. (A) Correlation between COTL1 expression and lymphocyte in BRCA. (B) Correlation between COTL1 expression and immune inhibitors in BRCA.
Discussion
COTL1 was first identified in pancreatic cancer tissues, where it was highly expressed compared to normal tissues.20 Subsequently, Jeong et al found that COTL1 was highly expressed in small cell lung cancer (SCLC).21 COTL1 was correlated with low metastatic potential in human non-small cell lung cancer (NSCLC), which suggested that COTL1 inhibits lung cancer metastasis.4 A recent study showed that COTL1 promotes glioblastoma (GBM) cell growth in vitro and in vivo, and immunohistochemical assays showed that COTL1 is highly expressed in human GBM tissues.5 Xia et al found that COTL1 inhibited breast cancer cell growth in vitro and in vivo,6 However, multiple online databases show that COTL1 is highly expressed in breast cancer tissue compared to normal tissues. Therefore, we studied COTL1 expression and associated prognosis in human breast cancer tissues by using a tissue microarray with 170 breast cancer tissue samples and sixty-two para-carcinoma tissue samples.Analysis using the TCGA database and immunohistochemistry showed that COTL1 expression was higher in breast cancer tissues than in normal tissues. We then determined whether COTL1 expression was correlated with prognosis in patients with breast cancer. The results from the Kaplan-Meier Plotter database were inconsistent: COTL1 mRNA level was positively correlated with poor prognosis (Figure 2A) while COTL1 protein level was negatively correlated with poor prognosis (). This could be because of the different data sources and processing methods of the Kaplan-Meier plotter database. Our tissue microarray data showed that high expression of COTL1 in breast cancer tissues was associated with poor prognosis (Figure 2D). Although only 139 cases were included in our study, our data strongly demonstrated that ER, and COTL1 were independent risk factors for OS (Table 3).Breast cancer is a heterogeneous disease with unique histological, molecular, and clinical phenotypes. The expression of estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) are considered the major determinants of breast cancer and can be used to refine molecular and prognostic taxonomy of breast cancer.22 Triple negative breast cancer (TNBC) accounts for 10–15% breast cancer cases and lacks ER and PR, and expresses low levels of HER29. TNBC involves immune activation and is associated with tumor-infiltrating lymphocytes (TILs).23 Immune checkpoint antagonists have been applied in the treatment of various tumors, specifically for CTLA-4, PD-1, and PD-L1.24 COTL1 expression is the highest in triple negative breast cancer and correlates with immune checkpoint antagonists, which suggests that COTL1 is a potential therapeutic target for breast cancer, especially triple negative breast cancer.Xia et al found that COTL1 could inhibit the proliferation of breast cancer. In current study, it was found that the expression of COTL1 was related to the high infiltration levels of lymphocytes immune cells. Therefore, we hypothesized that COTL1 may inhibit the proliferation of breast cancer by recruiting CD8+ T lymphocytes. However, high CD8+T cell infiltration was found in breast cancer tissues with high COTL1 expression, which was associated with poor OS in BRCA. This observation appears to be contradictory assuming that the presence of CD8+ T cells improves the patient survival. However, T cell exhaustion is widely described as a mechanism that inhibits CD8+ T cells proliferation and kills tumor cells. In fact, we found that high expression of T cell exhaustion markers PD-1, CTLA4, TIM-3, GZMB, LAG-3, and PDL1 were all associated with high expression of COTL1, supporting that T cell exhaustion could suppress T cell functions in BRCA. Therefore, these findings may explain why the high level of CD8+ T cells failed to induce a survival benefit in BRCA. Xia et al proposed that COTL1 inhibited breast cancer growth in the early stage and lost its growth inhibition function in the late stage, in current study, we found COTL1 expression was higher in breast cancer tissues than in normal tissues, and high COTL1 expression was correlated with poor prognosis. Xia et al speculated that COTL1 may play a bidirectional role in tumor like TGFβ1. The mechanism by which COTL1 functions as a tumor suppressor in the early stage of breast cancer progression and as an oncogene in the late stage merits further study.In conclusion, we used multiple online public databases and conducted immunohistochemical analysis of breast cancer tissues, and found that COTL1 is highly expressed in breast cancer and is correlated with poor prognosis. However, our study has limitations, such as small sample size and data from online public database. Future studies with larger samples can be used to verify the clinical expression of COTL1 in breast cancer and associated prognosis. Our results provide a basis for the development of clinical applications of COTL1 in breast cancer treatment.
Conclusions
We found that COTL1 expression was higher in breast cancer tissues than in normal tissues, and that high COTL1 expression was correlated with poor prognosis and immune cell infiltration. The current study provided novel insights for the future application of COTL1 in the treatment of breast cancer.
Authors: Taiwen Li; Jingyu Fan; Binbin Wang; Nicole Traugh; Qianming Chen; Jun S Liu; Bo Li; X Shirley Liu Journal: Cancer Res Date: 2017-11-01 Impact factor: 12.701
Authors: Darshan S Chandrashekar; Bhuwan Bashel; Sai Akshaya Hodigere Balasubramanya; Chad J Creighton; Israel Ponce-Rodriguez; Balabhadrapatruni V S K Chakravarthi; Sooryanarayana Varambally Journal: Neoplasia Date: 2017-07-18 Impact factor: 5.715
Authors: Nicole Vincent Jordan; Aditya Bardia; Ben S Wittner; Cyril Benes; Matteo Ligorio; Yu Zheng; Min Yu; Tilak K Sundaresan; Joseph A Licausi; Rushil Desai; Ryan M O'Keefe; Richard Y Ebright; Myriam Boukhali; Srinjoy Sil; Maristela L Onozato; Anthony J Iafrate; Ravi Kapur; Dennis Sgroi; David T Ting; Mehmet Toner; Sridhar Ramaswamy; Wilhelm Haas; Shyamala Maheswaran; Daniel A Haber Journal: Nature Date: 2016-08-24 Impact factor: 49.962