Expression pattern of COL5A3 in pan-cancer perspective
To assess the mRNA expression of COL5A3 in different cancer types, we analyzed datasets from 33 cancer types. As shown in Figure 1A, compared with normal tissues, COL5A3 was signifificantly upregulated in 28 of all 33 cancer types. This data indicated the mRNA expression of COL5A3 was differentially expressed in multiple tumor types.
COL5A3 mRNA and protein expression was upregulated in patients with pancreatic cancer
In order to determine the expression level of COL5A3mRNA and protein in pancreatic cancer, we analyzed the expression data of COL5A3 from TCGA data and HPA immunohistochemistry, and further analyzed and verified the expression of COL5A3 using GEO database. As shown in figure 1B, the expression level of COL5A3 in pancreatic cancer is significantly higher than that in normal tissues (P<0.001). Consistent with figure 1C, we can see that the expression of COL5A3 is up-regulated in the GSE16515 data set compared with normal tissues (P<0.05). Finally, as shown in figure 1D, the immunohistochemical staining of HPA also showed that the expression of COL5A3 protein was up-regulated in pancreatic cancer tissues. The above results showed that the mRNA and protein expression levels of COL5A3 in pancreatic cancer tissues were up-regulated.
The relationship between COL5A3 expression and clinical pathological characteristics in pancreatic cancer patients
We studied the clinicopathological features of pancreatic cancer patients with differential expression of COL5A3. As shown in Table 1, compared with the group with low expression of COL5A3, the group with high expression of COL5A3 had a significantly worse initial therapeutic effect (P=0.004), and the group with high expression of COL5A3 had a more significant history of alcohol (P=0.036). However, there was no statistical difference in other clinicopathological features, such as age, sex, TNM stage, pathological stage and so on. To sum up, the high expression of COL5A3 is closely related to some poor clinicopathological features, and also provides a new idea for the study of single gene molecules related to tumor drug resistance.
Table 1. Clinical characteristics of the pancreatic cancer patients (TCGA).
Characteristic
|
Low expression of COL5A3
|
High expression of COL5A3
|
p
|
n
|
89
|
89
|
|
T stage, n (%)
|
|
|
0.250
|
T1
|
6 (3.4%)
|
1 (0.6%)
|
|
T2
|
12 (6.8%)
|
12 (6.8%)
|
|
T3
|
68 (38.6%)
|
74 (42%)
|
|
T4
|
2 (1.1%)
|
1 (0.6%)
|
|
N stage, n (%)
|
|
|
0.754
|
N0
|
26 (15%)
|
24 (13.9%)
|
|
N1
|
59 (34.1%)
|
64 (37%)
|
|
M stage, n (%)
|
|
|
0.390
|
M0
|
31 (36.9%)
|
48 (57.1%)
|
|
M1
|
3 (3.6%)
|
2 (2.4%)
|
|
Pathologic stage, n (%)
|
|
|
0.758
|
Stage I
|
12 (6.9%)
|
9 (5.1%)
|
|
Stage II
|
71 (40.6%)
|
75 (42.9%)
|
|
Stage III
|
2 (1.1%)
|
1 (0.6%)
|
|
Stage IV
|
3 (1.7%)
|
2 (1.1%)
|
|
Radiation therapy, n (%)
|
|
|
0.058
|
No
|
66 (40.5%)
|
52 (31.9%)
|
|
Yes
|
17 (10.4%)
|
28 (17.2%)
|
|
Primary therapy outcome, n (%)
|
|
|
0.004**
|
PD
|
16 (11.5%)
|
33 (23.7%)
|
|
SD
|
8 (5.8%)
|
1 (0.7%)
|
|
PR
|
6 (4.3%)
|
4 (2.9%)
|
|
CR
|
40 (28.8%)
|
31 (22.3%)
|
|
Gender, n (%)
|
|
|
1.000
|
Female
|
40 (22.5%)
|
40 (22.5%)
|
|
Male
|
49 (27.5%)
|
49 (27.5%)
|
|
Race, n (%)
|
|
|
0.741
|
Asian
|
6 (3.4%)
|
5 (2.9%)
|
|
Black or African American
|
4 (2.3%)
|
2 (1.1%)
|
|
White
|
77 (44.3%)
|
80 (46%)
|
|
Age, n (%)
|
|
|
0.368
|
<=65
|
50 (28.1%)
|
43 (24.2%)
|
|
>65
|
39 (21.9%)
|
46 (25.8%)
|
|
Residual tumor, n (%)
|
|
|
0.720
|
R0
|
55 (33.5%)
|
52 (31.7%)
|
|
R1
|
24 (14.6%)
|
28 (17.1%)
|
|
R2
|
2 (1.2%)
|
3 (1.8%)
|
|
Histologic grade, n (%)
|
|
|
0.124
|
G1
|
21 (11.9%)
|
10 (5.7%)
|
|
G2
|
43 (24.4%)
|
52 (29.5%)
|
|
G3
|
22 (12.5%)
|
26 (14.8%)
|
|
G4
|
1 (0.6%)
|
1 (0.6%)
|
|
Anatomic neoplasm subdivision, n (%)
|
|
|
0.369
|
Head of Pancreas
|
66 (37.1%)
|
72 (40.4%)
|
|
Other
|
23 (12.9%)
|
17 (9.6%)
|
|
Smoker, n (%)
|
|
|
0.604
|
No
|
35 (24.3%)
|
30 (20.8%)
|
|
Yes
|
38 (26.4%)
|
41 (28.5%)
|
|
Alcohol history, n (%)
|
|
|
0.036*
|
No
|
25 (15.1%)
|
40 (24.1%)
|
|
Yes
|
57 (34.3%)
|
44 (26.5%)
|
|
History of diabetes, n (%)
|
|
|
0.774
|
No
|
52 (35.6%)
|
56 (38.4%)
|
|
Yes
|
20 (13.7%)
|
18 (12.3%)
|
|
History of chronic pancreatitis, n (%)
|
|
|
0.507
|
No
|
61 (43.3%)
|
67 (47.5%)
|
|
Yes
|
8 (5.7%)
|
5 (3.5%)
|
|
Family history of cancer, n (%)
|
|
|
0.271
|
No
|
26 (23.6%)
|
21 (19.1%)
|
|
Yes
|
27 (24.5%)
|
36 (32.7%)
|
|
Age, mean ± SD
|
63.6 ± 11.37
|
65.9 ± 10.13
|
0.155
|
*p<0.05, **p<0.01
Predictive value of COL5A3 for pancreatic cancer diagnosis and prognosis
In order to evaluate the effect of COL5A3 on the clinical prognosis of pancreatic cancer, we used K-M curve analysis to verify the prediction of clinical outcome of COL5A3. As shown in figure 2A-2B, the survival rate of OS in the high expression group was significantly lower than that in the low expression group of COL5A3 (p = 0.019), consistent with that in the DSS group, the survival rate in the high expression group was also significantly lower than that in the low expression group (p = 0.026). Next, we studied the clinical benefits of COL5A3, and we used ROC curves to demonstrate its value in the differential diagnosis of pancreatic cancer. As shown in figure 2C, the area under the curve (AUC) is 0.843, suggesting that COL5A3 has high sensitivity and specificity in the diagnosis of pancreatic cancer. The above results show that the up-regulation of COL5A3 indicates a worse prognosis and has high diagnostic value for clinical diagnosis.
COL5A3-associated PPI network and functional enrichment
In order to construct COL5A3-associated PPI networks and functional annotations, we used STRING database, GO and KEGG analysis. As shown in figure 3A, the PPI network shows that some genes are closely related to COL5A3, such as ADAMTS14, ADAMTS2, ADAMTS3, BMP1, COL11A1, COL1A2, LUM, P4HA3, PCOLCE and PCOLCE2。As shown in Figure 3B, COL5A3-related genes were involved in many biological processes (BP), cellular compositions (CC), and molecular functions (MF)。The results of GO and KEGG enrichment analysis of COL5A3-related genes were mainly involved in extracellular structure organization, extracellular matrix organization, collagen fibril organization, protein digestion and absorption. The GO and KEGG analysis demonstrated that COL5A3 might promote tumor proliferation, migration and invasion by affecting extracellular matrix. Analysis of the correlation between COL5A3 expression and co-expression genes in pancreatic cancer patients based on TCGA database, as shown in figure 4A-J.
Analysis of the relationship between the expression of COL5A3 and the level of immune cell infiltration in pancreatic carcinoma
We used TIMER database to analyze the expression of COL5A3 and tumor purity and the correlation between six kinds of infiltrating immune cells (B cells, CD8 T cells, CD4 T cells, macrophages, neutrophils and dendritic cells). As shown in figure 5A,The results displayed that the expression level of COL5A3 had obviously positive correlation with infiltrating levels of CD4+T cells(r=0.213、P=5.40e-03), macrophage cells(r=0.414、P=1.82e-08)、neutrophils(r=0.399、P=6.52e-08)、dendritic cells(r=0.375、P=4.45e-07) in pancreatic cancer, but no association with tumor purity and CD8+ T cells. P < 0.05 was defined as statistically significant. Figure 5B shows the relationship between COL5A3 and the expression of 28 kinds of tumor infiltrating lymphocytes in human cancer. As shown in figures 5C-H, the expression of COL5A3 is related to the abundance of Treg (r=0.37, P=4.43e-07), Tgd (r= 0.36, P=8.8e-07), Tcm CD4 (r= 0.35, P=1.83e-06), NK (r=0.348, P=2.15e-06), pDC (r=0.339, P=4.03e-06) and Tcm CD8 (r=0.257, P= 0.000526). These results suggest that the expression of COL5A3 may play an important role in the immune infiltration of pancreatic cancer.