2.1 Overexpression of PRKDC mRNA in NSCLC
We analyzed the mRNA and protein expression of PRKDC from the ONCOMINE(www.oncomine.org),the UALCAN(http://ualcan.path.uab.edu)and the human protein atlas(https://www.proteinatlas.org).We measured PRKDC mRNA expression patterns in 20 cancers in the ONCOMINE database and compared them with normal tissues. The results showed that PRKDC was overexpressed in 13 of the 19 Lung Cancer Studies (see figure 1).At the same time, PRKDC mRNA expression was significantly increased in Luad and LUSC samples from multiple datasets.In the Bhattacharjee [22] lung dataset, the mRNA expression of PRKDC in the LUAD sample increased 3.022 fold (p = 0.002).Likewise,in the Garber [23] , Stearman [24] , Landi [25] and Hou [26]lung dataset ,the mRNA expression of PRKDC increased 1.905-fold (p = 9.00E-04) , 1.836-fold (p = 1.13E-07) , 1.787-fold (p = 3.11E-12) and 1.73-fold (p = 3.09E-09) , respectively.In addition, significant increases in PRKDC mRNA expression were also observed in the LUSC sample.In Bhattacharjee's [22]lung dataset, PRKDC mRNA expression increased 8.885-fold (P = 3.87E-06). Similarly, in Talbot [27] , Garber [23] and Hou [26] lung datasets, the mRNA expression of PRKDC increased 2.484-fold (p = 1.19E-12) , 2.365-fold (p = 1.63E-04) and 2.224-fold (6.45E-13) , respectively (see table 1).PRKDC mRNA expression was again measured using the UALCAN database, which is based on clinical data from the level 3 RNA-seq and TCGA database for 31 cancer types, unlike the ORCOMINE database.The results showed that the mRNA expression of PRKDC was significantly overexpressed in both primary Luad and LUSC tissues (p = 1.62E-12)(see figures 2A and 2B).
Table 1
Relationship between mRNA expression of PRKDC in different lung cancer data sets
|
Types of NSCLC VS Lung
|
Fold Change
|
P value
|
t-test
|
Ref
|
PRKDC
|
Lung Adenocarcinoma
|
3.022
|
0.002
|
3.316
|
Bhattacharjee Lung
|
|
Lung Adenocarcinoma
|
1.905
|
0.0009
|
4.187
|
Garber Lung
|
|
Lung Adenocarcinoma
|
1.836
|
1.13E-07
|
7.124
|
Stearman Lung
|
|
Lung Adenocarcinoma
|
1.787
|
3.11E-12
|
7.754
|
Landi Lung
|
|
Lung Adenocarcinoma
|
1.730
|
3.9E-09
|
6.75
|
Hou Lung
|
PRKDC
|
Lung Squamous cell Carcinoma
|
8.885
|
3.87-06E
|
5.338
|
Bhattacharjee Lung
|
|
Lung Squamous cell Carcinoma
|
2.482
|
1.19E-012
|
8.775
|
Talbot Lung
|
|
Lung Squamous cell Carcinoma
|
2.365
|
0.000163
|
4.677
|
Garber Lung
|
|
Lung Squamous cell Carcinoma
|
2.224
|
6.45E-13
|
10.594
|
Hou Lung
|
2.2 Analysis of relationship between clinical pathological parameters and PRKD
UALCAN(http://ualcan.path.uab.edu)was used to analyze the relationship between PRKDC mRNA expression and clinicopathologic parameters of NSCLC patients, including tumor stage, smoking habit, sex, race and lymph node metastasis.In general, patients with advanced cancer tend to have higher mRNA expression at the tumor stage. In patients with Luad and LUSC, the highest mRNA expression of PRKDC was at stage 3(see figures 3A and 3B).In terms of smoking habits, PRKDC mRNA expression tends to be higher as the duration of smoking increases.In LUAD patients, the mRNA expression of PRKDC was significantly higher in smokers and smokers over 15 years than in non-smokers and smokers under 15 years,and there are statistical differences (figure 4A).In patients with LUSC, there was no significant difference in PRKDC mRNA expression between smoking and non-smoking patients, and between patients over 15 years of age and under 15 years of age.For gender differences, the highest mRNA expression of PRKDC was found in male in both Luad and LUSC patients (figure 4B). In terms of racial differences,among patients with Luad and LUSC, the highest mRNA expression levels of PRKDC were found in Asians (see figures 4C and 4D).For lymph node metastasis, PRKDC mRNA expression was highest in patients with Luad at stage N2 and highest in patients with LUSC at stage N1. The highest expression level of PRKDC mRNA was not at stage N3, possibly due to the small sample size (N = 2 in Luad and N = 5 in LUSC)(see figures 4E and 4F) .
In short, these results suggest that PRKDC mRNA expression is significantly correlated with clinicopathologic parameters of NSCLC.
2.3 Prognostic value of PRKDC mRNA expression in NSCLC patients
Kaplan Myer plotter(http://kmplot.com/analysis)was used to analyze the prognostic value of PRKDC mRNA expression in NSCLC patients.The expression of PRKDC mRNA was significantly correlated with the prognosis of NSCLC.The results showed that the high mRNA expression level of PRKDC (HR = 1.76,95% CI: 1.48-2.10, p = 1.4E-10) was significantly associated with poor prognosis (see figure 5).Therefore, PRKDC may be used as a biomarker for predicting survival in NSCLC patients.
2.4 Analysis of PRKDC gene mutation and prognosis
We analyzed the mutations in the PRKDC gene using the CBIOPORTAL database.The results showed that 119 out of 510 Luad patients were found to have the gene change, with a mutation rate of 23% (see figure 6A).Of 484 patients with LUSC sequencing, 90 LUAD patients were found to have gene changes, with a mutation rate of 19% (see figure 6B).In addition, the Kaplan-Meier and log-rank tests showed that the gene alterations in PRKDC were significantly different from those in Luad patients with shorter OS (figure 6C, p = 5.194E-3) and DFS (figure 6D, P = 5.256E-4) ,were not statistically significant compared with patients with LUSC with Short Os (figure 6E, p = 0.0835) and DFS (figure 6F, p = 0.364) .These results suggest that PRKDC has a high mutation rate in NSCLC patients, and PRKDC gene mutation may influence the prognosis of LUAD patients.
2.5 Function prediction and enrichment analysis of PRKDC and similar genes
We searched GEPIA2.0(http://gepia2.cancer-pku.cn)for 100 adjacent genes significantly associated with the PRKDC mutation and constructed a complete network by STRING11.0(https://string-db.org)(figure 7A).In addition, GO and KEGG analyzed the function of PRKDC and its 100 mutated neighboring genes in Metascape(http://www.metascape.org).See figure 7B-7D,biological Processes such as DNA repair (GO: 0006281) and regulation of DNA metabolic processes (GO: 0051052) are significantly regulated by the PRKDC mutation in LUAD. Cellular components such as the Chromosomal region (GO: 0098687) and the DNA damage site (GO: 0090734) are significantly regulated by the PRKDC mutation in LUAD.In addition, PRKDCA mutations also significantly affected molecular functions, such as impaired DNA binding (GO: 0003684) , DNA replicator binding (GO: 0003688) and so on.In KEGG analysis, 99 pathways including KO04110 cell cycle and KO03030DNA replication were associated with the function of PRKDC in LUAD (figure 7E).