MDK and PTN exhibited abnormal expression levels in CRC samples
We first analyzed data from the TCGA cohort and found that MDK was significantly upregulated (P<0.0001), while PTN was significantly downregulated (P<0.0001) in CRC samples versus normal tissues (Fig. 1a and 1b). Pearson correlation analysis showed no statistically significant correlation between the mRNA expression levels of PTN and MDK (r=0.050, P=0.279), suggesting that the expression of PTN and MDK may be subjected to unrelated regulatory mechanisms during the occurrence and development of CRC (Fig. 2).
MDK and PTN were associated with multiple biological pathways in CRC
To further elucidate the biological functions of MDK and PTN in CRC, we performed GO and KEGG signaling pathway enrichment analyses based on the LinkedOmics database. The results showed that MDK was mainly enriched with immune response and inflammatory response pathways in CRC and was also involved in RNA metabolic process and cell cycle regulation (Fig. 1c and 1d). Meanwhile, PTN was mainly enriched with metabolic process and DNA replication pathways (GO pathways) in CRC (Fig. 1e and 1f).
MDK and PTN were correlated with immune functions in CRC
Based on the LinkedOmics database, we performed GSEA on MDK and PTN in CRC. The results showed that the expression of MDK was significantly positively correlated with that of genes involved in humoral immune response, adaptive immune response, T cell activation, lymphocyte-mediated immunity, cell killing, regulation of inflammatory response, natural killer (NK) cell-mediated cytotoxicity, differentiation of Th1 and Th2 cells, intestinal immune network for IgA production, Th17 cell differentiation, inflammatory bowel disease, and autoimmune thyroid disease, and so on, in CRC (p<0.05, false discovery rate [FDR]<0.05; Fig. 3). Furthermore, the expression of PTN was significantly negatively correlated with that of genes implicated in NADH dehydrogenase complex assembly, tricarboxylic acid metabolic process, mitochondrial RNA metabolic process, protein localization to chromosome, cytoplasmic translation, DNA damage response, detection of DNA damage, pentose phosphate pathway, glyoxylate and dicarboxylate metabolisms, pyruvate metabolism, DNA replication, fructose and mannose metabolisms, citrate cycle (TCA cycle), and so on, in CRC (p<0.05, FDR<0.05; Fig. 4).
Associations of MDK and PTN with immune infiltration
Since our GSEA found that MDK and PTN played important immunological roles in CRC samples, we next sought to determine the relationships between MDK and PTN expression and the levels of infiltrating immune cells. Based on the TIMER database, we discovered that the expression of MDK was negatively correlated with tumor purity, and positively correlated with B cell, neutrophil and dendritic cell (DC) infiltration in CRC tissues (Fig. 5a). Meanwhile, PTN expression was also negatively correlated with tumor purity, but positively correlated with B cell, CD8+ T cell, CD4+ T cell, macrophage, neutrophil and dendritic cell infiltration in CRC tissues (Fig. 5a). Based on the “GSVA” package (1.34.0 version) of R (3.6.3 version), the correlations of MDK and PTN expression with the infiltration of various subtypes of aDCs (activated DCs), B cells, CD8+ T cells, cytotoxic cells, DCs, eosinophils, iDCs (immature DCs), macrophages, mast cells, neutrophils, NK cells, CD56bright cells [26], NK CD56dim cells [26], NK cells, pDCs (plasmacytoid DCs), T cells, T helper cells, Tcm (T central memory) cells, Tem (T effector memory) cells, Tfh (T follicular helper) cells, Tgd (T gamma delta) cells, Th1 cells, Th17 cells, Th2 cells, and Tregs (regulatory T cells) were comprehensively analyzed. The results (Fig. 5b, Table 1) were consistent with those obtained based on the TIMER database (Fig. 5a). We also calculated the immune score, stromal score, and ESTIMATE score for each CRC patient by using the “estimate” package (version 1.0.13) in R (version 3.6.3). Mann-Whitney U analysis revealed that the expression levels of MDK and PTN were both positively correlated with immune score, stromal score, and ESTIMATE score of the CRC patients (Fig. 5c).
Table 1
Correlation analysis between MDK and PTN expression and the infiltration of immune cells in CRC
| PTN | MDK |
Cell type | Correlation coefficient (Pearson) | p value (Pearson) | Correlation coefficient (Pearson) | p value (Pearson) |
aDCs | 0.192 | <0.001 | 0.275 | <0.001 |
B cells | 0.298 | <0.001 | 0.157 | <0.001 |
CD8 T cells | 0.238 | <0.001 | 0.119 | 0.011 |
Cytotoxic cells | 0.170 | <0.001 | 0.335 | <0.001 |
DCs | 0.332 | <0.001 | 0.248 | <0.001 |
Eosinophils | 0.355 | <0.001 | 0.190 | <0.001 |
iDCs | 0.459 | <0.001 | 0.328 | <0.001 |
Macrophages | 0.455 | <0.001 | 0.178 | <0.001 |
Mast cells | 0.562 | <0.001 | 0.199 | <0.001 |
Neutrophils | 0.262 | <0.001 | 0.168 | <0.001 |
NK CD56bright cells | -0.109 | 0.021 | 0.301 | <0.001 |
NK CD56dim cells | 0.096 | 0.041 | 0.242 | <0.001 |
NK cells | 0.454 | <0.001 | 0.238 | <0.001 |
pDCs | 0.197 | <0.001 | 0.084 | 0.074 |
T cells | 0.205 | <0.001 | 0.280 | <0.001 |
T helper cells | 0.222 | <0.001 | -0.139 | 0.003 |
Tcm cells | 0.153 | 0.001 | -0.127 | 0.007 |
Tem cells | 0.268 | <0.001 | 0.244 | <0.001 |
Tfh cells | 0.405 | <0.001 | 0.214 | <0.001 |
Tgd cells | 0.378 | <0.001 | 0.037 | 0.432 |
Th1 cells | 0.292 | <0.001 | 0.252 | <0.001 |
Th17 cells | -0.180 | <0.001 | -0.101 | 0.031 |
Th2 cells | 0.084 | 0.073 | -0.097 | 0.038 |
Tregs | 0.290 | <0.001 | 0.300 | <0.001 |
Prognostic analysis of MDK and PTN in CRC samples
We analyzed the overall survival (OS), disease-free survival (DSS) and progression-free interval (PFI) of CRC patients and found that low MDK expression was associated with a better prognosis in terms of OS (hazard ratio [HR]=1.58, 95% confidence interval (CI): 1.00-2.49; p=0.052) (Fig. 6a) and DSS (HR=2.26, 95% CI: 1.18-4.33; p=0.016) (Fig. 6b), but was not significantly associated with the PFI of CRC patients (HR=1.34, 95% CI: 0.93-1.94; p=0.116) (Fig. 6c). We also analyzed the associations between MDK expression and OS of CRC patients in the T1/T2 and T3/T4 stages. The results indicated that MDK expression correlated negatively with the OS of T1/T2-stage CRC patients (HR=2.99, 95% CI: 0.57-15.55; p = 0.194) (Fig. 6d) but positively with the OS of T3/T4-stage CRC patients (HR=1.62, 95% CI: 1.01-2.61; p=0.048) (Fig. 6e). Meanwhile, high PTN expression was associated with a better prognosis in terms of OS (HR=0.65, 95% CI: 0.44-0.96; p=0.031) (Fig.e 6f), but was not significantly associated with the DSS (HR=0.66, 95% CI: 0.40-1.09; p=0.103) (Fig. 6g) and PFI (HR=0.85, 95% CI: 0.58-1.23; p=0.378) (Fig. 6h) of CRC patients. We also analyzed the correlations between PTN expression and OS of CRC patients in the T1/T2 and T3/T4 stages. The findings revealed that PTN expression had no significant correlation with the OS of patients with T1/T2-stage CRC (HR=0.74, 95% CI: 0.14-3.80; p=0.715) (Fig. 6i) but had a negative correlation with the OS of patients with T3/T4-stage CRC (HR=0.67, 95% CI: 0.44-1.00; p=0.052) (Fig. 6j). In addition, the area under the ROC curve (AUC) values for the expression of MDK and PTN were 0.682 and 0.970, respectively, in CRC patients (Fig. 6k and 6l), indicating that the expression levels of MDK and PTN were closed related with the prognosis of CRC. Therefore, these two genes are potential biomarkers for CRC.