Research Paper Volume 13, Issue 6 pp 8706—8719

Predicting lung adenocarcinoma prognosis with a novel risk scoring based on platelet-related gene expression

Chengmao Zhou1,5, *, , Yongsheng Wang3, *, , Ying Wang1, , Lei Lei1, , Mu-Huo Ji1, , Guoren Zhou4, , Hongping Xia2,5,6, , Jian-Jun Yang1,5, ,

  • 1 Department of Anesthesiology, Pain and Perioperative Medicine, The First Affiliated Hospital of Zhengzhou University, Zhengzhou 450000, China
  • 2 Department of Pathology, School of Basic Medical Sciences & Key Laboratory of Antibody Technique of National Health Commission & Jiangsu Antibody Drug Engineering Research Center, Nanjing Medical University, Nanjing 211166, China
  • 3 Department of Respiratory Medicine, Affiliated Drum Tower Hospital of Nanjing University Medical School, Nanjing 210008, China
  • 4 Jiangsu Cancer Hospital, The Affiliated Cancer Hospital of Nanjing Medical University, Jiangsu Institute of Cancer Research, Nanjing 210009, China
  • 5 School of Medicine, Southeast University, Nanjing 210009, China
  • 6 Sir Run Run Hospital, Nanjing Medical University, Nanjing 211166, China
* Equal contribution

Received: May 6, 2020       Accepted: February 1, 2021       Published: February 22, 2021      

https://doi.org/10.18632/aging.202682
How to Cite

Copyright: © 2021 Zhou et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Lung adenocarcinoma is the most common subtype of non-small cell lung cancer, and platelet receptor-related genes are related to its occurrence and progression. A new prognostic indicator based on platelet receptor-related genes was developed with multivariate COX analysis. Prognostic markers based on platelet-related risk score perform moderately in prognosis prediction. The functional annotation of this risk model in high-risk patients shows that the pathways related to cell cycle, glycolysis and platelet-derived related factors are rich. It is worth noting that somatic mutation analysis shows that TTN and MUC16 have higher mutation burdens in high-risk patients. Moreover, the differential genes of high- and low-risk groups are regulated by copy number variation and miRNA. And we provide a free online nomogram web tool based on clinical factors and the risk score (https://wsxzaq.shinyapps.io/wsxzaq_nomogram/). The score has been verified among three independent external cohorts (GSE13213, GSE68465 and GSE72094), and is still an independent risk factor for lung adenocarcinoma. In addition, among the other 6 cancers, the OS prognosis of high and low-risk groups of PRS is different (P < 0.05). Our research results have screened multiple platelet differential genes with clinical significance and constructed a meaningful prognostic risk score (PRS).

Introduction

Global data show that approximately 1.5 million patients die of lung cancer each year, with a mortality rate over 25% [1]. Lung cancer is divided into small cell lung cancer and non-small cell lung cancer (NSCLC), with the latter accounting for 85% of total cases [2]. Lung adenocarcinoma (LUAD) accounts for 30–35% of all lung cancers, and about 500,000 people die of LUAD each year [3]. Although targeted drugs such as ALK and EGFR-TKI have survival benefits for LUAD patients with sensitive gene mutations, the survival time for advanced LUAD patients without sensitive gene mutations remains short. Thus, more treatment targets need to be identified to provide advanced LUAD patients with more effective treatment options and to prolong survival time.

In recent years, the incidence of adenocarcinoma has proliferated, and adenocarcinoma has replaced squamous cell carcinoma as the most prevalent type of NSCLC [4]. LUAD progresses rapidly with a high degree of malignancy. Patients who show clinical symptoms and signs are often at an advanced stage, causing an onerous burden on both themselves and society. Studies [5, 6] have shown that platelets play a non-negligible role in tumor metastasis, and there is a correlation between platelet count and cancer prognosis. Platelets are known to be involved in pathophysiological processes such as inflammation, vascular hemostasis, thrombosis and atherosclerosis [7]. There is also considerable evidence that platelets aggregate and activate malignant tumors [8]. They are also associated with distant tumor metastasis. Antiplatelet drugs (such as heparin) inhibit tumor metastasis, suggesting that inhibiting platelet activation may become a target for tumor treatment.

Potential mechanisms regarding platelet-associated receptors and LUAD prognosis may involve epigenetic changes and genetic mutations. Therefore, it is of clinical significance to build molecular models of platelet-associated receptor genes that predict LUAD prognosis and have the potential to guide personalized clinical treatment strategies. In this study, we developed a specific, accurate, and reliable prognostic risk scoring system based on platelet-associated receptor genes for LUAD. We also performed bioinformatics analysis to explore potential regulatory mechanisms.

Results

Identifying differentially expressed genes (TCGA)

The EdgeR algorithm can identify 188 differentially expressed genes, 83 up-regulated genes and 105 down-regulated genes (Figure 1A). As expected, gene function enrichment analysis showed that the platelet pathway was the most common. Platelet activation is one of the most common biological terms in biological processes (Figure 1B). And platelet activation is most often enriched by differentially expressed genes. (Figure 1C).

Differentially expressed platelet receptor-related genes, Gene functional enrichment and identification of hub platelet receptor-related genes. (A) Volcano plot of differentially expressed platelet receptor-related genes (DEPG) in the TCGA. (B) GO enrichment analysis of DEPG. (C) KEGG enrichment analysis of DEPG. (D) Single-factor analysis for DEPG. (E) Diagnosing lung adenocarcinoma patients.

Figure 1. Differentially expressed platelet receptor-related genes, Gene functional enrichment and identification of hub platelet receptor-related genes. (A) Volcano plot of differentially expressed platelet receptor-related genes (DEPG) in the TCGA. (B) GO enrichment analysis of DEPG. (C) KEGG enrichment analysis of DEPG. (D) Single-factor analysis for DEPG. (E) Diagnosing lung adenocarcinoma patients.

Identifying differential genes related to survival

After a single factor analysis, we found that 18 differences are related to lung adenocarcinoma patients’ overall survival rates, and these results are statistically significant (Figure 1D). The forest map shows that most of the identified central genes in the lung adenocarcinoma samples are up-regulated. The forest map of hazard ratios shows that most of these genes are risk factors. Nine genes were selected by multifactor regression, but because the validation group did not have the expression value of TREML1, we finally included only the other eight genes into the study. And we constructed a prognostic marker and divided lung adenocarcinoma patients into two groups. The calculation formula is as follows: Platelet-related risk score (PRS) = RAB27B *0.017537331+CCNA2*0.020495528+PCDH7 *0.038831123+HAS2*0.062731768+P2RX1*0.173961877+PFKP*0.006157281+ENTPD2*0.065271557+TIMP1*0.000897402

Platelet-related receptor genes based on multivariate Cox may be an important tool for diagnosing lung adenocarcinoma patients according to potential, and the ROC curve of CCNA2 is 0.970. (Figure 1E).

Survival analysis showed that the high-risk group’s survival rate was lower than that of the low-risk group, and these results were statistically significant. Figure 2A represents a risk map, including signature-based, individual survival status between groups and the expression levels and distributions of the included modeling genes. It shows a clear difference in survival status between the risk groups, with red dots indicating death and green dots indicating survival. Many deaths occurred in the high-risk group, while most patients in the low-risk group survived through follow-up. In addition, PRS scores are verified in three other independent queues (GSE13213, GSE68465 and GSE72094), and the results are consistent with those of TCGA (Figure 2B2D). Multivariate Cox regression analysis showed that after adjusting for other parameters like stage and T, the prognostic index can become an independent predictor (Figure 2E). Moreover, multiple Cox regression analysis showed that in GSE13213, GSE68465 and GSE72094 data sets, the prognosis index could also become an independent prediction index (Figure 2F2H).

The prognostic value of PRS (platelet receptor-related risk score). (A–D) Kaplan-Meier curves of overall survival based on the PRS in the TCGA, GSE13213, GSE68465 and GSE72094 cohort (E–H) Multivariate Cox regression analysis in the TCGA, GSE13213, GSE68465 and GSE72094 cohort.

Figure 2. The prognostic value of PRS (platelet receptor-related risk score). (A–D) Kaplan-Meier curves of overall survival based on the PRS in the TCGA, GSE13213, GSE68465 and GSE72094 cohort (E–H) Multivariate Cox regression analysis in the TCGA, GSE13213, GSE68465 and GSE72094 cohort.

The ROC curve was 0.732, indicating that there may be a medium potential for PRS based prognosis in survival monitoring. (Figure 3A) There are still survival differences between the high-risk and low-risk groups in the other three datasets (GSE13213, GSE68465 and GSE72094), and the ROC curve values formed by PRS scores are 0.791, 0.737 and 0.606, respectively (Figure 3B3D). It was also found that for different lung cancer patient subtypes, the platelet-associated gene risk score was a moderately viable indicator for prognosis. (Supplementary Figure 1).

Development of the nomogram based on PRS and clinical features. (A–D) The receiver operating characteristic (ROC) curve of OS in TCGA, GSE13213, GSE68465 and GSE72094 cohort. (E) ROC curve of the nomogram based on PRS in TCGA. (F) The nomogram based on PRS in TCGA. (G–I) 1, 3 and 5-year calibration curves of nomogram in TCGA.

Figure 3. Development of the nomogram based on PRS and clinical features. (A–D) The receiver operating characteristic (ROC) curve of OS in TCGA, GSE13213, GSE68465 and GSE72094 cohort. (E) ROC curve of the nomogram based on PRS in TCGA. (F) The nomogram based on PRS in TCGA. (G–I) 1, 3 and 5-year calibration curves of nomogram in TCGA.

Using the above eight genes and clinical factors for lung cancer, we developed a nomogram to predict the 1-year, 3-year, and 5-year overall survival. The ROC analysis showed sufficient discrimination with an AUC as 0.746, 0.746 and 0.695, indicating that the nomogram has medium prediction performance (Figure 3E and Figure 3F). Besides, the calibration chart showed the best prediction accuracy, and the predicted survival rate was equal to the actual survival rate (Figure 3G3I). Moreover, we constructed the nomogram by combining PRS with clinical factors, and provided a free nomogram prediction web tool for clinicians (https://wsxzaq.shinyapps.io/wsxzaq_nomogram/) (Supplementary Figure 2).

Correlation analysis

Correlation analysis between PRS and immune cell infiltration showed that PRS and B_cell or CD4 T_cell (Cor = –0.229, P = 4.952e-07 and Cor = –0.102, P = 0.026) were related, and correlation analysis between PRS and immune checkpoints showed that PRS and PD-L1 (Cor = 0.183, P = 5.575e-05) were related. PRS was negatively correlated with ESTIMATE score (Cor = –0.112, P = 0.014), indicating that PRS was positively correlated with tumor purity of LUAD. PRS was positively correlated with the MKI67 score, which proves that the high-risk group may be closely related to the regulation of cancer cell cycle and cancer cell proliferation (Cor = 0.295, P = 4.865e-11). (Figure 4).

The relationships between PRS and other prognostic features (Immune infiltration, immune checkpoints and cell proliferation-related genes, etc.).

Figure 4. The relationships between PRS and other prognostic features (Immune infiltration, immune checkpoints and cell proliferation-related genes, etc.).

Functional annotation of established eight gene signatures

GSEA was performed in TCGA high-risk populations to study key biological and cellular processes related to prognosis. There were clear abundant GO pathways in high-risk patients, including cell cycle, glycolysis and platelet-derived related factors (Figure 5A). Most of these pathways were involved in cell proliferation and migration and may lead to lung cancer metastasis and recurrence. Previous studies have shown that platelet-derived growth factors regulate glycolysis and therefore affect the cell cycle of cancer cells. This may also be a potential mechanism for high-risk PRS, suggesting poor LUAD prognosis.

Genome-wide analysis of PRS risk score. (including mutation, copy number change and microRNA change, etc.). (A) GSEA was performed in TCGA high-risk populations. (B) Mutation of 8 Factors Constructed into Risk Score in CBioportal. (C) Differential copy number in PRS high and low group. (D) Distribution of the First 30 Mutant Genes in High and Low-Risk Groups. (E) Overlapping of Up-regulated and Down-regulated Genes for Copy Number Change Genes and Differences between High and Low-Risk Groups (F) Downregulated MicRNA target control gene and overlapping gene of up-regulated and down-regulated genes with a high and low-risk difference.

Figure 5. Genome-wide analysis of PRS risk score. (including mutation, copy number change and microRNA change, etc.). (A) GSEA was performed in TCGA high-risk populations. (B) Mutation of 8 Factors Constructed into Risk Score in CBioportal. (C) Differential copy number in PRS high and low group. (D) Distribution of the First 30 Mutant Genes in High and Low-Risk Groups. (E) Overlapping of Up-regulated and Down-regulated Genes for Copy Number Change Genes and Differences between High and Low-Risk Groups (F) Downregulated MicRNA target control gene and overlapping gene of up-regulated and down-regulated genes with a high and low-risk difference.

Analyzing mutations and CNV in patients with high and low PRS expression

CBioportal database (http://www.cbioportal.org/) shows that among the 8 genes, HAS2 has the highest mutation rate of 10%. (Figure 5B) TTN and MUC16 were detected with high mutation frequencies in tumors expressed in the high-risk group, and they are well-known tumor-associated genes (Figure 5D). This demonstrates that PRS distinguishes mutant genes with poor prognosis. Figure 5C shows the differential copy number in PRS high and low risk, which get copy number in chromosomes 10, 11, 12 and 20, and mitigated in 8, 10, 11 and 12. Among them, the copy numbers of 36 and 20 genes in the high and low-risk groups, respectively, also changed (Figure 5E). The differential expression genes of PRS between high and low expression levels are affected by copy number variation, which indicates that there is no independence between the abnormal copy number expression genes and the differential expression genes of PRS.

The correlation between miRNA and risk score differential expression

MiRNA is a small non-coding RNA, which can be involved in targeted mRNA cleavage and post-transcriptional regulation. Previous studies have shown that miRNA expression imbalance is a key step in lung adenocarcinoma carcinogenesis. Our work evaluated the relationship between miRNA target-controlled down-regulated genes and risk score differentially expressed genes. Therefore, we used TargetScan to predict the relationship between miRNA and its target genes. A total of 120 up-regulated miRNA and 14 down-regulated miRNA were detected in the high-risk group (the low-risk group served as a reference). We identified 233 pairs of down/miRNA-mRNA interactions by TargetScan (i.e. the miRNA and target control gene scores were greater than or equal to 0.6), of which 10 pairs overlapped with risk score differential expression (Figure 5F). Overlapping genes, including XCT and UCK2, are adverse prognostic factors of lung adenocarcinoma, and the mechanism may be related to regulating the cell cycle of cancer cells. These results indicate that miRNA in the down-regulated group has feedback loops in the carcinogenesis processes of the high-risk and low-risk groups.

External data support

By analyzing the HPA (http://www.proteinatlas.org) and TCPA (http://tcpaportal.org) protein databases, we found that the protein expression of related genes constructed by PRS was also consistent with the previous results. (Figure 6A and 6B) Finally, we verified the relationship between PRS and TP53 mutations and KRAS mutations on GSE13213 data sets and found that the risk score was not affected by these two gene mutations. (Figure 6C and 6D).

External data support. (A) Protein Expression of Risk Score Components in TCPA Database. (B) Immunohistochemical Results of Risk Score Components in HPA Database. (C) The relationship between PRS and TP53 mutations and KRAS mutations on GSE13213.

Figure 6. External data support. (A) Protein Expression of Risk Score Components in TCPA Database. (B) Immunohistochemical Results of Risk Score Components in HPA Database. (C) The relationship between PRS and TP53 mutations and KRAS mutations on GSE13213.

In addition, among the other 6 cancers, the high-risk scoring group of PRS also has a significantly worse prognosis than the low-risk group, and the difference is statistically significant (P < 0.05). (Supplementary Figure 3).

Discussion

This study has identified and verified platelet receptor-related gene risk scoring systems for predicting LUAD’s OS. After a comprehensive analysis, 188 genes were selected as the final specific genes. Finally, eight prognostic gene signatures (PRS: RAB27B, CCNA2, PCDH7, HAS2, P2RX1, PFKP, ENTPD2 and TIMP1) were constructed by COX multivariate regression model and verified in three other independent datasets (GSE13213, GSE68465 and GSE72094). Besides, GSEA enrichment analysis showed that high-risk group scores may affect cell cycle, glycolysis and platelet-derived related factors. Moreover, the potential mechanism of PRS risk score in lung adenocarcinoma prognosis was comprehensively analyzed.

At present, studies on these gene-related cancers have been reported. Studies [9] have shown that Rab27B regulates the aggressive growth and metastasis of ER-positive breast cancer cell lines, and its increased expression is associated with poor prognosis. Similarly, a high expression of Rab27b is associated with poor LUAD prognosis, which may be a potential indicator of LUAD metastasis and poor prognosis [10]. The SNP (rs769236) on the CCNA2 promoter may be associated with an increased risk of colon, liver, and lung cancer [11]. Additionally, HAS2 may be involved in the aggressive phenotypes of primary breast cancer [12]. The simultaneous silencing of HA synthase 2 (HAS2) and HAS3 or CD44 and RHAMM inhibits cell proliferation and survival, as well as the EGFR/AKT/ERK signaling pathway [13]. It has been shown that PFKP is one of the differential lncRNA-miRNA-mRNA networks in NSCLC, and it is related to glycolysis [14]. In hepatocellular carcinoma, hypoxia induces the expression of ENTPD2 on cancer cells; This increases extracellular 5'-AMP, and in turn, promotes MDSC maintenance by preventing its differentiation [15]; moreover, studies have shown that ENTPD2 is of great significance in the diagnosis, monitoring and prognosis of LUAD [16]. Inhibiting TIMP1 expression reduces proliferation and metastasis and increases apoptosis. TIMP1 mRNA expressed in the platelets of colorectal cancer patients can be transported into colorectal cancer cells, thereby promoting tumor growth [17]. The meta-analysis also shows that highly-expressed TIMP-1 is associated with poor prognosis in NSCLC patients [18]. In addition, studies have shown that P2RX1 is one of the prognostic genes of LUAD [19].

Gene enrichment analysis of high-risk patients also shows that the pathways related to cell cycle, platelet-derived growth factor and glycolysis were enriched. Studies have shown that CCNA2 [20], HAS2 [21] and TIMP1 [22] can regulate the cell cycle of cancer cells. Besides, studies [23] have also shown that PFKP can affect the proliferation of cancer cells through the glycolytic pathway. Thus, platelet-related receptor genes in lung adenocarcinoma may affect the proliferation and metastasis of lung adenocarcinoma by regulating cell cycle and glycolysis. And the potential mechanism of distinguishing lung adenocarcinoma by this risk score may be related to copy number variation and negative feedback by miRNA and may be related to regulating the cell cycle.

Many studies have shown that these eight genes play an important role in the biological function of platelets. Rab27b has been a key modulator of the secretion of dense particles in platelets [24]. Studies have shown that there is a close relationship between CCNA2 and plasma release of platelets [25]. Platelet-derived growth factors modulate the activity of HAS2 [26]. Studies also have shown that P2RX1 is closely related to the platelet activation pathway [27]. Phosphofructokinase platelet-type (PFKP) is a rate-limiting enzyme involved in glycolysis and is closely related to the progression of non-small cell lung cancer [28]. Excised nucleoside triphosphate diphosphate hydrolase (ENTPD) is involved in inflammation and platelet aggregation [29]. Studies have shown a high correlation between TIMP1 hemoglobin levels and in vitro platelet reactivity [30].

Correlation analysis shows that the high-risk group is positively correlated with tumor mutation load and PD-L1.

However, certain limitations of this study should also be noted. First, this study is only a retrospective analysis, and a prospective multicenter study should also be conducted. Second, we deleted the gene TREML1 because the validation group did not have it. However, in the forest plot, TREML1 has very strong statistical significance. And some studies have shown that TREML1 is a prognostic gene of prostate cancer. However, there is no research on TREML1-associated lung adenocarcinoma. Therefore, it is necessary to further experiment with the role of TREML1 in the tumor in future research.

Conclusions

In conclusion, we have constructed and confirmed a prognostic risk scoring system for lung adenocarcinoma platelet receptor-related genes consisting of eight genes. And we provided a free nomogram prediction web tool for clinicians (https://wsxzaq.shinyapps.io/wsxzaq_nomogram/). According to the results of functional annotation, the high-risk group participates in the cell cycle, and the pathways related to platelet-derived growth factor and glycolysis are abundant. This may help to explain the molecular mechanism of lung adenocarcinoma occurrence and development. In addition, our data provide new and promising evidence for predicting lung adenocarcinoma biomarkers and targeted therapy.

Materials and Methods

Clinical samples and data collection

RNA sequencing and clinical data were available from The Cancer Genome Atlas (TCGA: https://portal.gdc.cancer.gov/) and the Gene Expression Omnibus (GEO: https://www.ncbi.nlm.nih.gov/geo/) data portals. We adopted three datasets (GSE13213, GSE72094 and GSE68465) from the GEO database. We also populated a list of platelet-associated receptor genes from the GO Consortium database (http://geneontology.org/). The mutation spectrum and copy number variation (CNV) data for the entire TCGA dataset study population were also downloaded. Samples with less than 30 days of survival and samples with 0 expression values higher than 10% were removed.

Analysis of differentially expressed genes

We performed differential gene analysis on all transcription data, setting the false discovery rate (FDR) to < 0.05, and |fold change| > 1.5 as the cutoff value [31]. The differentially expressed genes were then extracted from all the differentially expressed genes. Functional enrichment analysis was performed through the GO and KEGG pathways.

Survival analysis

Survival related genes were selected by univariate and multivariate COX analysis, which was performed with R’s survival software package.

Building a prognostic index based on platelet-related genes

We performed a multivariate analysis of core platelet-related genes to construct a platelet risk score (PRS). This was based on the expression data multiplied by the Cox regression coefficient. Patients were divided into high-risk and low-risk groups based on the median. The prognostic value of PRS was evaluated in patients with different subtypes. The TIMER database was used to analyze the abundance of tumor-infiltrating immune cells. The association was then calculated between PRS and LUAD immune cell infiltration.

Gene set enrichment analysis

Gene set enrichment analysis (GSEA) is a computational method to identify over-represented biological processes and pathways.

Mutation and copy number variation analysis

The genes were analyzed with different mutations in different risk score subtypes (high and low risks). The mutation spectrum of the entire study population from the TCGA dataset was determined as previously described. Copy number variation (CNV) data were downloaded from GDAC Firehose and divided into datasets according to the risk score expression.

Gene expression and miRNA integration

To study microRNAs’ (miRNAs’) potential gene regulation, we focused on the abnormally expressed miRNAs (adjusted p-value < 0.05, log2 | multiple changes | > 1) and the RNAseq genes with statistically significant differences selected from high and low-risk patients. The correlation between miRNA and regulatory genes was obtained through TargetScan analysis (http://www.targetscan.org/vert_72/).

Statistical analysis

The AUC for the survival ROC curve was calculated by the survivalROC R software package to verify the prognostic signature’s performance. Differences between clinical parameters were tested using independent t-tests. A p-value of less than 0.05 was considered statistically significant.

Ethics approval and consent to participate

Our study did not require the approval of an ethics committee, as it was a secondary analysis of a public database.

Data availability statement

The data that support the findings of this study are openly available in TCGA at https://cancergenome.nih.gov/, and GEO at https://www.ncbi.nlm.nih.gov/geo/, reference number: GSE13213, GSE68465 and GSE72094.

Supplementary Materials

Supplementary Figures

Author Contributions

All of the authors worked on this study. All authors have read and approved the final manuscript.

Conflicts of Interest

None.

Funding

This study was supported by grants from the Postgraduate Research & Practice Innovation Program of Jiangsu Province (No. KYCX19_0113), and the grants from The Recruitment Program of Overseas High-Level Young Talents, “Innovative and Entrepreneurial Team” (No. (2018)2015), Science and Technology Grant (BE2019758) of Jiangsu Province.

References

  • 1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018; 68:7–30. https://doi.org/10.3322/caac.21442 [PubMed]
  • 2. Samarghandian S, Azimi-Nezhad M, Farkhondeh T. Thymoquinone-induced antitumor and apoptosis in human lung adenocarcinoma cells. J Cell Physiol. 2019; 234:10421–10431. https://doi.org/10.1002/jcp.27710 [PubMed]
  • 3. Yang QS, Li B, Xu G, Yang SQ, Wang P, Tang HH, Liu YY. Long noncoding RNA LINC00483/microRNA-144 regulates radiosensitivity and epithelial-mesenchymal transition in lung adenocarcinoma by interacting with HOXA10. J Cell Physiol. 2019; 234:11805–11821. https://doi.org/10.1002/jcp.27886 [PubMed]
  • 4. Chen W. Cancer statistics: updated cancer burden in China. Chin J Cancer Res. 2015; 27:1. https://doi.org/10.3978/j.issn.1000-9604.2015.02.07 [PubMed]
  • 5. Menter DG, Tucker SC, Kopetz S, Sood AK, Crissman John D, Honn Kenneth V. Platelets and cancer: a casual or causal relationship: revisited. Cancer Metastasis Rev. 2014; 33:231–269. https://doi.org/10.1007/s10555-014-9498-0 [PubMed]
  • 6. Zhou C, Lei L, Ji MH, Yang J, Xia H. The prognostic significance of CLEC-2-related gene signal in lung adenocarcinoma: A multicenter development and validation cohort. Clin Transl Med. 2020; 10:e100. https://doi.org/10.1002/ctm2.100 [PubMed]
  • 7. Thomas MR, Storey RF. The role of platelets in inflammation. Thromb Haemost. 2015; 114:449–458. https://doi.org/10.1160/th14-12-1067 [PubMed]
  • 8. Zhou C, Wang Y, Lei L, Ji MH, Yang JJ, Xia H. Identifying Common Genes Related to Platelet and Immunity for Lung Adenocarcinoma Prognosis Prediction. Front Mol Biosci. 2020; 7:563142. https://doi.org/10.3389/fmolb.2020.563142 [PubMed]
  • 9. Hendrix A, Maynard D, Pauwels P, Braems G, Denys H, Van den BR, Lambert J, Van BS, Cocquyt V, Gespach C, Bracke M, Seabra MC, Gahl WA, et al. Effect of the secretory small GTPase Rab27B on breast cancer growth, invasion, and metastasis. J Natl Cancer Inst. 2010; 102:866–880. https://doi.org/10.1093/jnci/djq153 [PubMed]
  • 10. Zhang L, Fan W, Xu L, Mao Q, Chen Y, Mao Y, Xu L, Wang J. Rab27b Is a Potential Indicator for Lymph Node Metastasis and Unfavorable Prognosis in Lung Adenocarcinoma. Dis Markers. 2018; 2018:7293962. https://doi.org/10.1155/2018/7293962 [PubMed]
  • 11. Kim DH, Park SE, Kim M, Ji YI, Kang MY, Jung EH, Ko E, Kim Y, Kim S, Shim YM, Park J. A functional single nucleotide polymorphism at the promoter region of cyclin A2 is associated with increased risk of colon, liver, and lung cancers. Cancer. 2011; 117:4080–4091. https://doi.org/10.1002/cncr.25930 [PubMed]
  • 12. Lien HC, Lee YH, Jeng YM, Lin CH, Lu YS, Yao YT. Differential expression of hyaluronan synthase 2 in breast carcinoma and its biological significance. Histopathology. 2014; 65:328–339. https://doi.org/10.1111/his.12390 [PubMed]
  • 13. Song JM, Im J, Nho RS, Han YH, Upadhyaya P, Kassie F. Hyaluronan-CD44/RHAMM interaction-dependent cell proliferation and survival in lung cancer cells. Mol Carcinog. 2019; 58:321–333. https://doi.org/10.1002/mc.22930 [PubMed]
  • 14. He Z, Dang J, Song A, Cui X, Ma Z, Zhang Z. Identification of LINC01234 and MIR210HG as novel prognostic signature for colorectal adenocarcinoma. J Cell Physiol. 2019; 234:6769–6777. https://doi.org/10.1002/jcp.27424 [PubMed]
  • 15. Chiu DK, Tse AP, Xu IM, Di CJ, Lai RK, Li LL, Koh HY, Tsang FH, Wei LL, Wong CM, Ng IO, Wong CC. Hypoxia inducible factor HIF-1 promotes myeloid-derived suppressor cells accumulation through ENTPD2/CD39L1 in hepatocellular carcinoma. Nat Commun. 2017; 8:517. https://doi.org/10.1038/s41467-017-00530-7 [PubMed]
  • 16. Zhang J, Zhang J, Yuan C, Luo Y, Li Y, Dai P, Sun W, Zhang N, Ren J, Zhang J, Gong Y, Xie C. Establishment of the Prognostic Index Reflecting Tumor Immune Microenvironment of Lung Adenocarcinoma Based on Metabolism-Related Genes. J Cancer. 2020; 11:101–7115. https://doi.org/10.7150/jca.49266 [PubMed]
  • 17. Yang L, Jiang Q, Li DZ, Zhou X, Yu DS, Zhong J. TIMP1 mRNA in tumor-educated platelets is diagnostic biomarker for colorectal cancer. Aging (Albany NY). 2019; 11:8998–9012. https://doi.org/10.18632/aging.102366 [PubMed]
  • 18. Kang X, Kong F, Huang K, Li L, Li Z, Wang X, Zhang W, Wu X. LncRNA MIR210HG promotes proliferation and invasion of non-small cell lung cancer by upregulating methylation of CACNA2D2 promoter via binding to DNMT1. Onco Targets Ther. 2019; 12:3779–3790. https://doi.org/10.2147/ott.s189468 [PubMed]
  • 19. Ma C, Luo H, Cao J, Zheng X, Zhang J, Zhang Y, Fu Z. Identification of a Novel Tumor Microenvironment-Associated Eight-Gene Signature for Prognosis Prediction in Lung Adenocarcinoma. Front Mol Biosci. 2020; 7:571641. https://doi.org/10.3389/fmolb.2020.571641 [PubMed]
  • 20. Li J, Ying Y, Xie H, Jin K, Yan H, Wang S, Xu M, Xu X, Wang X, Yang K, Zheng X, Xie L. Dual regulatory role of CCNA2 in modulating CDK6 and MET-mediated cell-cycle pathway and EMT progression is blocked by miR-381-3p in bladder cancer. FASEB J. 2019; 33:1374–1388. https://doi.org/10.1096/fj.201800667r [PubMed]
  • 21. Zhang HY, Liang F, Wang F, Zhang JW, Wang L, Kang XG, Wang J, Duan QL. In Vitro Effects of HAS-2 Gene Silencing on the Proliferation and Apoptosis of the MCF-7 Human Breast Cancer Cell Line. Cell Physiol Biochem. 2016; 40:807–817. https://doi.org/10.1159/000453140 [PubMed]
  • 22. Hekmat O, Munk S, Fogh L, Yadav R, Francavilla C, Horn H, Würtz SØ, Schrohl AS, Damsgaard B, Rømer MU, Belling KC, Jensen NF, Gromova I, et al. TIMP-1 increases expression and phosphorylation of proteins associated with drug resistance in breast cancer cells. J Proteome Res. 2013; 12:4136–4151. https://doi.org/10.1021/pr400457u [PubMed]
  • 23. Zhou K, Yao YL, He ZC, Chen C, Zhang XN, Yang KD, Liu YQ, Liu Q, Fu WJ, Chen YP, Niu Q, Ma QH, Zhou R, et al. VDAC2 interacts with PFKP to regulate glucose metabolism and phenotypic reprogramming of glioma stem cells. Cell Death Dis. 2018; 9:988. https://doi.org/10.1038/s41419-018-1015-x [PubMed]
  • 24. Tolmachova T, Abrink M, Futter C, Authi K, Seabra M. Rab27b regulates number and secretion of platelet dense granules. Proc Natl Acad Sci U S A. 2007; 104:5872–5877. https://doi.org/10.1073/pnas.0609879104 [PubMed]
  • 25. Tsai WC, Yu TY, Lin LP, Lin MS, Wu YC, Liao CH, Pang JH. Platelet rich plasma releasate promotes proliferation of skeletal muscle cells in association with upregulation of PCNA, cyclins and cyclin dependent kinases. Platelets. 2017; 28:491–497. https://doi.org/10.1080/09537104.2016.1227061 [PubMed]
  • 26. Kastner S, Thomas G, Jenkins R, Davies M, Steadman R. Hyaluronan induces the selective accumulation of matrix- and cell-associated proteoglycans by mesangial cells. Am J Pathol. 2007; 171:1811–1821. https://doi.org/10.2353/ajpath.2007.070085 [PubMed]
  • 27. Bongiovanni D, Santamaria G, Klug M, Santovito D, Felicetta A, Hristov M, von SM, Aslani M, Cibella J, Weber C, Moretti A, Laugwitz KL, Peano C, et al. Transcriptome Analysis of Reticulated Platelets Reveals a Prothrombotic Profile. Thromb Haemost. 2019; 119:1795–1806. https://doi.org/10.1055/s-0039-1695009 [PubMed]
  • 28. Wang F, Li L, Zhang Z. Platelet isoform of phosphofructokinase promotes aerobic glycolysis and the progression of non-small cell lung cancer. Mol Med Rep. 2021; 23:74. https://doi.org/10.3892/mmr.2020.11712 [PubMed]
  • 29. Pelikant-Małecka I, Sielicka A, Kaniewska E, Smoleński RT, Słomińska EM. Influence of 4-pyridone-3-carboxamide-1Β-D-ribonucleoside (4PYR) on activities of extracellular enzymes in endothelial human cells. Nucleosides Nucleotides Nucleic Acids. 2016; 35:732–736. https://doi.org/10.1080/15257770.2016.1174263 [PubMed]
  • 30. Shah P, Yang W, Sun S, Pasay J, Faraday N, Zhang H. Platelet glycoproteins associated with aspirin-treatment upon platelet activation. Proteomics. 2017; 17:1600199. https://doi.org/10.1002/pmic.201600199 [PubMed]
  • 31. Yan H, Qu J, Cao W, Liu Y, Zheng G, Zhang E, Cai Z. Identification of prognostic genes in the acute myeloid leukemia immune microenvironment based on TCGA data analysis. Cancer Immunol Immunother. 2019; 68:1971–1978. https://doi.org/10.1007/s00262-019-02408-7 [PubMed]