Identification of a Novel Prognostic Classifier for Platelet-related Genes Characteristics and Prognostic Significance of Aberrant Expression of PAFAH1B3 in Hepatocellular Carcinoma

DOI: https://doi.org/10.21203/rs.3.rs-910487/v1

Abstract

Background

Complex crosstalk between tumor cells and platelets is closely related to the development, relapse, and drug resistance of hepatocellular carcinoma (HCC). Therefore, an intensive analysis of the relationship between platelet-related genes and the effectiveness of immunotherapy is necessary for improving the poor prognosis of HCC patients.

Methods

Genes associated with platelets in the GeneCards database were collected and were used to identified molecular subtypes using a non-negative matrix decomposition algorithm (NMF) and constructed a platelet-related genes-based prognostic stratification model by the LASSO-Cox regression and stepwise Cox regression analysis. The effect of this feature on the immune microenvironment of HCC and the response to immune checkpoint inhibitors was also explored.

Results

After identifying two molecular subtypes, we constructed a platelet-related genes-based prognostic stratification model that can be effectively used for immune checkpoint inhibitor (PD1, PD-L1, PD-L2, and CTLA4) efficacy and prognosis prediction in HCC patients, which was subsequently validated using patient samples from ICGC, GSE14520 and a small sample size clinical cohort. We also found downregulation of PAFAH1B3 remarkably inhibited the proliferation and migration ability of Hep3B cells.

Conclusions

In conclusion, we constructed a prognostic classifier based on platelet-related genes that could effectively classify HCC patients for prognostic prediction provide new light on the selection of optimal individualized antiplatelet therapy for HCC patients in future clinical practice.

1. Introduction

Platelets are small, non-nucleated blood cells produced by megakaryocytes in the bone marrow, which play an important role in clotting and maintaining hemostasis after mechanical injury to blood vessels. In recent years, an increasing number of studies have demonstrated the strong link between platelets and tumors [1]. Tumors can cause abnormalities in platelet counts and functions of activation and aggregation. Platelet counts have been shown to be a prognostic predictor for patients with a variety of tumors, including lung, ovarian, prostate, colon cancers, and melanoma [2-5]. Abnormal activation and aggregation of platelets, especially tumor-activated platelets, may lead to thrombosis, especially cancer-related thrombosis, which may cause death in cancer patients [6, 7]. Interestingly, platelets contain a large number of bioactive molecules within them and express different receptors on their surface that can also, in turn, facilitate the growth and metastasis of tumors [8, 9]. Platelets can secrete transforming growth factor beta (TGF-β) to increase the proliferation of ovarian cancer cells, and the proliferation capacity is diminished with TGF-β1 blocking antibodies [10]. Platelet granules contain various proangiogenic factors, including vascular endothelial growth factor (VEGF), platelet-derived growth factor (PDGF), angiopoietin-1and IL-6, which can induce intratumor blood vessel formation and promote the proliferation of tumor-associated pericytes in the tumor microenvironment (TME) [11-13]. In addition, platelets can reduce E-cadherin levels and upregulate Snail, vimentin and fibronectin levels in tumor cells, inducing invasive EMT phenotype, and promoting tumor cell metastasis [13]. In short, there is a variety of crosstalk between platelets and cancer cells. In brief, there is a variety of crosstalk between platelets and cancer cells, and a deeper understanding of them can help clinicians develop strategies for the use of antiplatelet drugs in cancer prevention and treatment.

Considering that cancer and platelet crosstalk is a complex process involving multiple genes, we constructed and validated a prognostic stratification model in public datasets based on the platelet-related genes which could be efficiently applied for prognostic categorization of HCC patients. In addition, we found that platelet activating factor acetylhydrolase 1B catalytic subunit 3 (PAFAH1B3) was an independent prognostic factor in HCC patients and that reducing its expression inhibited tumor cell proliferation and migration.

2. Material And Methods

2.1 Public datasets and collection of genes related to platelet

This study includes transcriptomic expression data and clinical characteristics of HCC patients from three publicly available datasets, including TCGA-LIHC, ICGC (LIRI-JP), and GSE14520, and the clinical features of HCC patients in the three cohorts were showed in Table S1. Genes associated with platelets in the GeneCards database (https://www.genecards.org/) were collected [14]. Genes with a correlation score >12 were considered platelet-associated genes using "platelet" as the keyword for the search.

2.2 Molecular subtypes identification by non-negative matrix factorization (NMF) algorithm

Based on the platelet-associated genes, the HCC samples were clustered using the NMF algorithm with the criterion "brunet" and 50 iterations. The number of clusters (K) ranged from 2 to 6 while the optimal number of clusters was determined based on cophenetic, dispersion, and profile. Kaplan-Meier survival examination was further performed to assess survival differences between different subtypes identified by NMF.

2.3 DEG Identification and Functional Analysis

The expression levels of platelet-associated genes were explored in Gene Expression Profiling Interactive Analysis (GEPIA) [15] (http://gepia.cancer-pku.cn/) and genes with statistically significant differences were defined as differential expression genes (DEGs). Then functional annotation of DEGs for GO and KEGG pathway analysis were explored in the Metascape database [16] (http://metascape.org/). P<0.05 was considered statistically significant.

2.4 TME and tumor-infiltrating immune cells Analysis

Stromal and immune scores in tumor tissue were estimated by ESTIMATE based on gene expression profiles of HCC samples to assess the abundance of stromal and immune cells within the tumor [17]. In addition, abundance ratios of 22 tumor-infiltrating immune cells (TIICs) in HCC TME were measured by CIBERSORT [18] (http://cibersort.stanford.edu/).

2.5 Prognostic risk score model construction

Univariate Cox regression was performed to calculate the association between platelet-associated genes and survival outcomes of HCC patients. Overlapping genes between prognosis-related genes and DEGs were used for further analysis. Then LASSO-Cox regression strategy and stepwise Cox regression examination were performed to investigate the quality and establish the classifier associated with prognosis. A risk score was at last set up based on the premise of directly combining the equation underneath with the mRNA expression level duplicated the multivariate Cox relapse coefficient (β) demonstrate. Risk score = (βmRNA1 × mRNA1) + (βmRNA2 × mRNA2) +…+ (βmRNAn × mRNAn). We stratified patients in TCGA dataset into two subgroups due to the ideal hazard score edge. The prescient control and autonomy of the prognostic signature in TCGA were evaluated by ROC examination, Kaplan-Meier survival examination and cox relative risks relapse investigation. Gene set enrichment analysis (GSEA) between the two subgroups was performed to distinguish the altogether cautioned GO and KEGG items with FDR < 0.05.

2.6 Clinical specimens and Quantitative real-TME PCR (qRT-PCR) analysis

Fresh frozen tumor tissues from previously collected HCC patients were selected as a small sample size clinical validation cohort [19] and the clinical characteristics of HCC patients were showed in Table S1. qRT-PCR was used to detect the mRNA levels of genes in the model [20]. After the relative mRNAs expression levels were normalized to β-ACTIN and log2 transformed, patients were stratified into two subgroups according to the above formula. Primer sequences are showed in Table S2.

2.7 Genetic alterations and Immune checkpoint genes analysis

The mutation and CNA data of 350 HCC patients were downloaded from TCGA to analyze the difference of genetic alterations between different subgroups with R package “maftools”. Furthermore, to explore the impacts of the prognostic classifier on immunotherapies, we calculated the relationship between risk score and 4 potentially available targeted immune checkpoint genes including PD1, PD-L1, PD-L2, and CTLA4.

2.8 Drug susceptibility analysis

The association between anticancer drug sensitivity and mRNA molecules in our risk model was directly explored in the CellMiner database[21]. 574 in advanced clinical trials and 216 Food and Drug Administration (FDA)-approved drugs were used for follow-up analyses. Drugs with adjusted P value <0.001 and Pearson correlation coefficient >0.3 as cut-off criteria were considered tumor-sensitive drugs.

2.9 Cell culture and cell infection

Hep3B cells purchased from Cell Bank of the Shanghai Institute of Cells, Chinese Academy of Science (Shanghai, China) were cultured in recommended DMEM medium (Sangon Biotech, China) with 10% fetal bovine serum (FBS, Sangon Biotech, Shanghai, China) in 100% humidity at 37 °C with 5% CO2. Small interfering RNA (siRNA) (5′-CCUCUGCAUGCACUUAACUTT-3′) and negative control (5′-UUCUCCGAACGUGUCACGUTT-3′) were selected for reducing PAFAH1B3 expression in Hep3B cells using lipofectamine TM 3000 transfection reagent (Invitrogen, Carlsbad, USA) according to the manufacturer’s instructions.

2.10 Cell proliferation and migration assay

The viability of cells was assessed through the cell counting kit-8 (CCK-8, Sangon Biotech, Shanghai, China) under the manufacturer’s instructions. For the CCK-8 assay, in short, cells were seeded in 96-well cell culture clusters at a density of 1 x 105 cells per well and cultured for 1 day, 2 days, and 3 days, respectively. After culturing, 10 μL CCK-8 solution was added into each well, and then the absorbance was detected at a wavelength of 450 nm within 4 h with a microplate reader. For the colony-formation assay, 100 Hep3B cells were plated into six-well plates and cultured for 10 days. Then cell colonies were fixed and stained with crystal violet (Sangon Biotech, Shanghai, China) in 10% ethanol for 5 minutes. Finally, cell colonies were imaged and counted. For the wound-healing assay, cells were seeded in 6-well plates and cultured to approximately 80% confluence in serum-free medium and then cell monolayers were scratched with a sterile pipette tip. Cells were cultured in DMEM medium with 10% FBS for the next 24h after removal of cell debris by PBS washing. The area of the wound width was measured after photographing the wound width of the cell monolayers. Three independent duplicates needed to ensure the accuracy in this assay and the wound closure rate was calculated as [1 - (wound area / original wound area)] from photographs. For the transwell assay, after transwell filters were coated with Matrigel, cultured cells were resuspended in 200 μL serum-free DMEM at a density of 1 x 104 cells per mL, and plated into the transwell inserts while the wells were filled with 500 μL DMEM supplemented with 10% FBS. After incubated at 37℃ for 48h, cells attached to the downside of the transwell filters were stained with 0.1% crystal violet in PBS for 15min and counted under microscopy at 200 x magnification.

2.11 Statistical analysis

Quantitative variables were analyzed using the independent-samples t test. ROC curve analysis and Kaplan-Meier survival analysis were performed to assess the prediction performance of survival outcomes with R software (Version 4.0.3). Cox proportional model was performed to analyze the relationship between a prognostic classifier and survival outcomes, together with other clinical features. Results were considered statistically significant when P value <0.05.

3. Results

3.1 Molecular subtypes identification

A total of 160 platelet-associated genes were collected in the GeneCards database. The optimal number of clusters was identified as two based on cophenetic, dispersion, and profile (Figure 1A, Supplementary Figure S1). HCC samples were divided into two subcategories based on the platelet-associated genes and significant differences between patients in the two subgroups were observed (Figure 1B). Furthermore, patients in Cluster 2 had significantly higher PD1, PD-L1, PD-L2, and CTLA4 expression when compared with patients in Cluster 1 (Figure 1C-F). In addition, TME characteristics among the two subgroups were also compared. As shown in Figure 2A, samples from Cluster 1 had lower immune, stromal, and ESTIMATE scores compared with samples from Cluster 2. According to the CIBISORT, five types of TIICs, including T cell, CD4+ T cell, Dendritic cell, Mast cell, and NK cell, differed between Cluster 1 and Cluster 2 subgroups (Figure 2B). GO items with FDR < 0.05 were identified based on gene set enrichment analysis (GSEA) performed between Cluster 1 and Cluster 2 subgroups, which were mainly involved in homeostasis (Figure S2). Subsequently, genetic alteration analysis demonstrated that the mutation rates of the top 10 most significantly mutated genes were significantly different in the two subgroups. CTNNB1 (34%) was the most commonly mutated gene in the Cluster 1 subgroup, while TP53 (37%) was the most commonly mutated gene in the Cluster 2 subgroup (Figure S3). All above results demonstrated that prognosis, TME characteristics, TIICs infiltration, functional alterations, and somatic mutations were different between patients in the two subgroups, and patients in the Cluster 2 subgroup might be more effective against immune checkpoint blockade therapy.

3.2 Establishment of a signature for Cancer Stem Cell Characteristics in HCC

After 24 DEGs between normal and tumor tissues were identified from the GEPIA database, 52 prognosis-related genes were calculated by univariate Cox regression. Ten overlapping genes between prognosis-related genes and DEGs were identified (Figure 3A, B) and their enriched GO and KEGG items were shown in Figure S4.  Then the ten platelet-associated genes were selected by performing the LASSO-Cox regression model based on the minimum value of λ and 7 genes including DTNBP1, PAFAH1B3, SERPINE1, VEGFA, AP3B1, CD63 and PDGFA were screened as shown in Figure 3C. These 7 genes were then placed into a stepwise Cox proportional model and finally, a prognostic six-gene signature was identified. Risk score = (0. 0.2579817×PAFAH1B3) + (0.1354325×SERPINE1) + (0.3748835×AP3B1). Risk scores for HCC patients were calculated with the above formula, and patients were stratified into high- or low-risk subgroups with an optimal risk score threshold (Figure 3D). The association between risk score and clinical characteristics including age, gender, grade, TNM stage, vascular invasion, T stage, value of AFP, cirrhosis, HBV infection status, and tumor status were evaluated. The results revealed that higher risk scores were linked to vascular invasion, advanced TNM stage, AFP values, later grade, later T stage, and recurrence (Supplementary Figure S5). Kaplan-Meier survival analysis revealed that patients with higher risk score were significantly relevant to poorer survival outcomes (Figure 3E). Furthermore, ROC analysis revealed that this classifier had a good prognostic performance with AUCs at 1-, 2-, 3-, 5-year of 0.727, 0.654, 0.683, 0.644 (Figure 3F). The decision curve analysis (DCA) of the model was shown in Figure 3G. Finally, to explore whether this classifier could be acted as an independent prognostic model for HCC patients, univariable and multivariate Cox analyses were performed. Univariable Cox regression analysis revealed that this classifier was statistically associated with survival outcomes for HCC patients (HR=4.263, 95%CI 2.478-7.335, P < 0.001). Then statistically significant variables obtained above were entered into multivariate Cox regression analysis, which revealed that this classifier could be served as an independent prognostic factor for HCC patients (HR=2.677, 95%CI 1.538-4.667, P < 0.001) after adjusting for other clinical features. The expression levels of the three platelet-associated genes between normal and tumor tissues in the GEPIA database were shown in Figure S6A. The results of the survival variance analysis showed that all of the three genes were significantly related to overall survival (Figure S6B) while only PAFAH1B3 and AP3B1 were significantly related to disease-free survival (Figure S6C). Moreover, the protein expression in tissues (Figure S6D) and mRNA expression in HCC cell lines (Figure S6E) of the three genes were explored in Human Protein Atlas database [22] (HPA, www.proteinatlas.org) and Cancer Cell Line Encyclopedia (CCLE) database [23] (https://portals.broadinstitute.org/ccle), respectively. However, SERPINE1 protein was not detected in HCC.

3.3 Functional, Genetic alterations and Immune checkpoint genes analysis 

GO and KEGG items with FDR < 0.05 were identified based on gene set enrichment analysis (GSEA) performed between high- and low-risk score subgroups, which were mainly involved in the catabolic and metabolism process (Figure S7A, B). Subsequently, genetic alteration analysis demonstrated that the mutation rates of the top 10 most significantly mutated genes were significantly different in the two subgroups. TP53 (41%) was the most commonly mutated gene in the high-risk score subgroup, while CTNNB1 (31%) was the most commonly mutated gene in the low-risk score subgroup (Figure S8). Finally, we calculated the relationship between risk scores and expression levels of PD1, PD-L1, PD-L2, and CTLA4 genes to explore the impact of the prognostic model on immunotherapy and observed significant differences and positive correlations between them (Figure 4A, B). According to the CIBISORT, five types of TIICs, including T cell, CD4+ T cell, Dendritic cell, Mast cell, and NK cell, differed between high- and low- risk score subgroups (Figure S9). All above results demonstrated that functional alterations and somatic mutations were different between patients with higher and lower risk scores and patients in the high-risk score subgroup might be more effective against immune checkpoint blockade therapy.

3.4 Verification of the signature in ICGC and GSE14520 cohort

To validate the signature, ICGC and GSE14520 datasets were applied as validation cohorts. Risk scores of patients were calculated with the same formula, and patients were stratified into high- or low-risk subgroups in the ICGC cohort (Figure 5A) and GSE14520 cohort (Figure 5D). Kaplan-Meier survival analysis revealed that patients with higher risk scores were prominently relevant to poorer OS rates in the ICGC cohort (Figure 5B), as well as in the GSE14520 cohort (Figure 5E). ROC analysis revealed that this signature had a good prognostic performance with AUCs at 1-, 2-, 3-year of 0.675, 0.647, 0.652 in the ICGC cohort (Figure 5C) and 0.543, 0.608, 0.604 in the GSE14520 cohort (Figure 5F), respectively. Univariable and multivariate Cox analysis were performed to explore whether the six-gene signature could be acted as an independent prognostic factor for HCC patients. Univariable Cox regression analysis revealed that this six-gene signature was statistically associated with survival outcomes for HCC patients (IGGC: HR=5.228, 95%CI 1.985-13.770, P < 0.001; GSE14520: HR=2.529, 95%CI 1.333-4.798, P = 0.004). Then statistically significant variables obtained above were entered into multivariate Cox regression analyses, which revealed that this classifier could be served as an independent prognostic factor for HCC patients (IGGC: HR=3.372, 95%CI 1.258-90.38, P = 0.015; GSE14520: HR=1.982, 95%CI 1.055-3.711, P = 0.032) after adjusting for other clinical features. 

3.5 Establishment of a nomogram model in TCGA

To explore the coefficient prediction efficiency of this classifier, a nomogram model was built in TCGA dataset, and the result revealed that the nomogram with a C-index of 0.708 could help us provide a quantitative method for predicting the 1-, 2-, 3-, 5-year survival rate accurately (Figure 6A). The overlap between the forecasted and actual probabilities of 1-, 2-, 3-, 5-year survival rates in the calibration curves indicated good agreement (Figure 6B).

3.6 Drug susceptibility analysis

Among the 574 in advanced clinical trials and 216 Food and Drug Administration (FDA)-approved drugs, 60 were considered tumor-sensitive drugs (Table S3), and the top 16 most significant tumor-sensitive drugs were shown in Figure S10.

3.7  Expression levels of genes in this classifier

The expression levels of the three genes in our model were explored by qRT-PCR in an independent clinical cohort. All three genes showed differential expression in normal and tumor tissues (Figure 7A), while only PAFAH1B3 was significantly associated with overall survival in HCC patients (Figure 7B).  Furthermore, risk scores of patients were calculated with the same formula, and patients were stratified into high- or low-risk subgroups after normalization (Figure 7C). Kaplan-Meier survival analysis revealed that patients with higher risk scores were prominently relevant to poorer OS rates (Figure 7D). ROC analysis revealed that this signature had a good prognostic performance with AUCs at 1-, 2-, 3-year of 0.754, 0.657, 0.688 (Figure 7E).  All above demonstrated that our model could act as a good classification indicator for HCC prognosis.

3.8 Downregulation of PAFAH1B3 remarkably inhibited the proliferation and migration ability of Hep3B cells.

Considering that only PAFAH1B3 was associated with overall survival among the three genes in our collected clinical samples, we further explored the relationship between PAFAH1B3 expression and cell proliferation and migratory capacity. PAFAH1B3 siRNA transfection could significantly reduce PAFAH1B3 expression in Hep3B cells (Figure 8A). Cell proliferation capacity was significantly inhibited in Hep3B cells transfected with PAFAH1B3 siRNA according to the result of CCK-8 assay (Figure 8B) and colony-formation assay (Figure 8C). In addition, the proportion of wound closure (Figure 8D) and the number of cells crossing the membrane (Figure 8E) in the siPAFAH1B3 group were significantly lower than in the negative control group, indicating that cell migration capacity was significantly inhibited in Hep3B cells by downregulating the expression of PAFAH1B3.

4. Discussion

There is growing evidence of complex mutual crosstalk between tumor cells and platelets not only in the blood but also in the TME [1, 24]. Tumor cells can transfer mutated RNA into platelets by releasing tumor-derived exosomes, thus altering the platelet RNA profile [25]. Tumor cells may also promote increased platelet counts in patients by secreting thrombopoietic cytokines such as interleukin-1 (IL1), IL3, IL6, and IL11 [26]. In turn, platelets can promote tumor cells to maintain proliferative activity, assist tumor cells to resist cell death, and induce tumor angiogenesis [10, 12]. Platelets can also facilitate tumor invasion and migration by inducing an aggressive epithelial-mesenchymal transition (EMT) phenotype in tumor cells [27]. Tumor-activated platelets can release large amounts of TGF-β and downregulate NKG2D expression, thereby protecting tumor cells from NK cell-mediated lysis [28]. Therefore, tumor-associated platelet protein and RNA profiles may be useful biomarkers for cancer detection, prognosis prediction, and clinical treatment monitoring of HCC patients. In this study, after platelet-associated genes were collected in the GeneCards database, we first utilized a non-negative matrix decomposition algorithm to classify HCC patients into two molecular subtypes and significant differences in survival outcomes, expression of ICIs genes, and TME of HCC patients between the two subtypes were found. Subsequently, we constructed and validated a prognostic stratification model based on platelet-related genes that can be effectively applied to the prognostic classification of HCC patients. Finally, we found that downregulation of PAFAH1B3 remarkably inhibited the proliferation and migration ability of Hep3B cells. Our findings may potentially guide the selection of optimal individualized antiplatelet therapy for HCC patients in future clinical practice.

Serine family E member 1 (SERPINE1) is a serine protease inhibitor that regulates urokinase and plasminogen activator with the function of reducing excessive degradation of the extracellular matrix, regulating cell adhesion balance, and stimulating cell growth and vascular generation [29]. It has been demonstrated that overexpressed SERPINE1 could facilitate the aggressiveness and invasiveness of tumor cells in colorectal cancer [30], esophageal cancer [31], gastric cancer [32], and breast cancer [33]. SERPINE1 was also involved as a downstream of ELK3 in the regulation of prostate cancer proliferation, migration, and apoptosis [34]. Adaptor protein complex 3 subunit β3A (AP3B1) was a protein-coding gene whose mutations could cause a rare autosomal recessive primary immunodeficiency called Hermansky Pudlak type 2 syndrome (HPS2) [35] and Hemophagocytic lymphohistiocytosis (HLH) [36]. AP3B1 might also be involved in regulating the prognosis of HCC patients as a target of microRNA-9 [37]. In our study, we found SERPINE1 and AP3B1 expression significantly differed between normal and HCC tissues, while both of them were not associated with the overall survival of HCC patients, which may be related to our relatively small sample size, and follow-up experiments are still needed to verify it. Platelet-activating factor acetylhydrolase 1B catalytic subunit 3 (PAFAH1B3) is a member of the PAFAH1B complex, which could regulate the deacetylation activity of platelet-activating factor (PAF). In recent years, aberrant expression of PAFAH1B3 has been detected as a metabolic oncogene in different kinds of tumors, including gastric cancer [38], hypopharyngeal squamous cell carcinoma [39], and osteosarcoma [40]. When analyzed using pan-cancer TCGA data, upregulated PAFAH1B3 was strongly associated with poor prognosis in a variety of tumors, including HCC [41]. Consistent with previous studies, we found that PAFAH1B3 expression levels were significantly higher in HCC tissues compared to normal tissues and were independently associated with patient prognosis, and that reducing PAFAH1B3 with siRNA significantly inhibited the proliferative and migratory properties of Hep3B cells. All these results suggested that PAFAH1B3 could be utilized as a biomarker for HCC prognosis.

Undeniably, our study has some limitations. the diversity and individual variability of HCC patients may reduce the performance of this feature. In addition, the small sample size limits the validation of the model, and future multicenter randomized controlled studies are needed to evaluate this feature. Finally, the specific mechanism of the PAFAH1B3 gene in HCC needs to be further evaluated by more in vivo and in vitro experiments in the future.

In summary, we constructed and validated a prognostic stratification model based on platelet-related genes that can be effectively applied to the prognostic classification of HCC patients, and we found that downregulation of PAFAH1B3 remarkably inhibited the proliferation and migration ability of Hep3B cells. Our findings may potentially provide new light on the selection of optimal individualized antiplatelet therapy for HCC patients in future clinical practice.

Declarations

ACKNOWLEDGMENTS

Not applicable.

AUTHOR CONTRIBUTIONS

ZGH designed the study, downloaded, analyzed the data and wrote the manuscript and the manuscript. The final manuscript has been approved by all authors.

FUNDING

No funding. 

DATA AVAILABILITY STATEMENT

We can find the datasets analyzed in this study in the https://xena.ucsc.eduhttps://dcc.icgc.org/projects/LIRI-JP and https://www.ncbi.nlm.nih.gov/geo.

ETHICS APPROVAL AND CONSENT TO PARTICIPATE

This study was supported by the Ethics Committees of Zhengzhou University. 

CONSENT FOR PUBLICATION

Not applicable.

COMPETING INTERESTS

The authors declare no conficts of interest regarding this work. 

AUTHOR DETAILS

1 Department of Blood transfusion, the First Affiliated Hospital of Zhengzhou University, Zhengzhou, China.

References

1.     Xu XR, Yousef GM, Ni H: Cancer and platelet crosstalk: opportunities and challenges for aspirin and other antiplatelet agents. Blood 2018, 131(16):1777-1789.

2.     Sylman JL, Boyce HB, Mitrugno A, Tormoen GW, Thomas IC, Wagner TH, Lee JS, Leppert JT, McCarty OJT, Mallick P: A Temporal Examination of Platelet Counts as a Predictor of Prognosis in Lung, Prostate, and Colon Cancer Patients. Sci Rep 2018, 8(1):6564.

3.     Zu R, Yu S, Yang G, Ge Y, Wang D, Zhang L, Song X, Deng Y, He Q, Zhang K et al: Integration of platelet features in blood and platelet rich plasma for detection of lung cancer. Clin Chim Acta 2020, 509:43-51.

4.     Rachidi S Md P, Kaur MM, Lautenschlaeger TM, Li Z Md P: Platelet count correlates with stage and predicts survival in melanoma. Platelets 2019, 30(8):1042-1046.

5.     Hu Q, Hada A, Han L: Platelet count as a biomarker for monitoring treatment response and disease recurrence in recurrent epithelial ovarian cancer. J Ovarian Res 2020, 13(1):78.

6.     Meikle CK, Kelly CA, Garg P, Wuescher LM, Ali RA, Worth RG: Cancer and Thrombosis: The Platelet Perspective. Front Cell Dev Biol 2016, 4:147.

7.     Schulman S: How I treat recurrent venous thromboembolism in patients receiving anticoagulant therapy. Blood 2017, 129(25):3285-3293.

8.     Leblanc R, Peyruchaud O: Metastasis: new functional implications of platelets and megakaryocytes. Blood 2016, 128(1):24-31.

9.     Li N: Platelets in cancer metastasis: To help the "villain" to do evil. Int J Cancer 2016, 138(9):2078-2087.

10.    Cho MS, Bottsford-Miller J, Vasquez HG, Stone R, Zand B, Kroll MH, Sood AK, Afshar-Kharghan V: Platelets increase the proliferation of ovarian cancer cells. Blood 2012, 120(24):4869-4872.

11.    Battinelli EM, Markens BA, Kulenthirarajan RA, Machlus KR, Flaumenhaft R, Italiano JE, Jr.: Anticoagulation inhibits tumor cell-mediated release of platelet angiogenic proteins and diminishes platelet angiogenic response. Blood 2014, 123(1):101-112.

12.    Farooqi AA, Siddik ZH: Platelet-derived growth factor (PDGF) signalling in cancer: rapidly emerging signalling landscape. Cell Biochem Funct 2015, 33(5):257-265.

13.    Holmes CE, Levis JE, Schneider DJ, Bambace NM, Sharma D, Lal I, Wood ME, Muss HB: Platelet phenotype changes associated with breast cancer and its treatment. Platelets 2016, 27(7):703-711.

14.    Paolacci S, Precone V, Acquaviva F, Chiurazzi P, Fulcheri E, Pinelli M, Buffelli F, Michelini S, Herbst KL, Unfer V et al: Genetics of lipedema: new perspectives on genetic research and molecular diagnoses. Eur Rev Med Pharmacol Sci 2019, 23(13):5581-5594.

15.    Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z: GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res 2017, 45(W1):W98-w102.

16.    Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, Benner C, Chanda SK: Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun 2019, 10(1):1523.

17.    Yoshihara K, Shahmoradgoli M, Martínez E, Vegesna R, Kim H, Torres-Garcia W, Treviño V, Shen H, Laird PW, Levine DA et al: Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun 2013, 4:2612.

18.    Gentles AJ, Newman AM, Liu CL, Bratman SV, Feng W, Kim D, Nair VS, Xu Y, Khuong A, Hoang CD et al: The prognostic landscape of genes and infiltrating immune cells across human cancers. Nat Med 2015, 21(8):938-945.

19.    Su L, Zhang G, Kong X: Prognostic Significance of Pregnancy Zone Protein and Its Correlation with Immune Infiltrates in Hepatocellular Carcinoma. Cancer Manag Res 2020, 12:9883-9891.

20.    Zhang G: Expression and Prognostic Significance of BANF1 in Triple-Negative Breast Cancer. Cancer Manag Res 2020, 12:145-150.

21.    Reinhold WC, Sunshine M, Liu H, Varma S, Kohn KW, Morris J, Doroshow J, Pommier Y: CellMiner: a web-based suite of genomic and pharmacologic tools to explore transcript and drug patterns in the NCI-60 cell line set. Cancer Res 2012, 72(14):3499-3511.

22.    Colwill K, Gräslund S: A roadmap to generate renewable protein binders to the human proteome. Nat Methods 2011, 8(7):551-558.

23.    Nusinow DP, Szpyt J, Ghandi M, Rose CM, McDonald ER, 3rd, Kalocsay M, Jané-Valbuena J, Gelfand E, Schweppe DK, Jedrychowski M et al: Quantitative Proteomics of the Cancer Cell Line Encyclopedia. Cell 2020, 180(2):387-402.e316.

24.    Schlesinger M: Role of platelets and platelet receptors in cancer metastasis. J Hematol Oncol 2018, 11(1):125.

25.    Best MG, Sol N, Kooi I, Tannous J, Westerman BA, Rustenburg F, Schellen P, Verschueren H, Post E, Koster J et al: RNA-Seq of Tumor-Educated Platelets Enables Blood-Based Pan-Cancer, Multiclass, and Molecular Pathway Cancer Diagnostics. Cancer Cell 2015, 28(5):666-676.

26.    Lin RJ, Afshar-Kharghan V, Schafer AI: Paraneoplastic thrombocytosis: the secrets of tumor self-promotion. Blood 2014, 124(2):184-187.

27.    Labelle M, Begum S, Hynes RO: Direct signaling between platelets and cancer cells induces an epithelial-mesenchymal-like transition and promotes metastasis. Cancer Cell 2011, 20(5):576-590.

28.    Placke T, Salih HR, Kopp HG: GITR ligand provided by thrombopoietic cells inhibits NK cell antitumor activity. J Immunol 2012, 189(1):154-160.

29.    Azimi I, Petersen RM, Thompson EW, Roberts-Thomson SJ, Monteith GR: Hypoxia-induced reactive oxygen species mediate N-cadherin and SERPINE1 expression, EGFR signalling and motility in MDA-MB-468 breast cancer cells. Sci Rep 2017, 7(1):15140.

30.    Mazzoccoli G, Pazienza V, Panza A, Valvano MR, Benegiamo G, Vinciguerra M, Andriulli A, Piepoli A: ARNTL2 and SERPINE1: potential biomarkers for tumor aggressiveness in colorectal cancer. J Cancer Res Clin Oncol 2012, 138(3):501-511.

31.    Klimczak-Bitner AA, Kordek R, Bitner J, Musiał J, Szemraj J: Expression of MMP9, SERPINE1 and miR-134 as prognostic factors in esophageal cancer. Oncol Lett 2016, 12(5):4133-4138.

32.    Liao P, Li W, Liu R, Teer JK, Xu B, Zhang W, Li X, McLeod HL, He Y: Genome-scale analysis identifies SERPINE1 and SPARC as diagnostic and prognostic biomarkers in gastric cancer. Onco Targets Ther 2018, 11:6969-6980.

33.    Zhang Q, Lei L, Jing D: Knockdown of SERPINE1 reverses resistance of triple‑negative breast cancer to paclitaxel via suppression of VEGFA. Oncol Rep 2020, 44(5):1875-1884.

34.    Mao Y, Li W, Hua B, Gu X, Pan W, Chen Q, Xu B, Wang Z, Lu C: Silencing of ELK3 Induces S-M Phase Arrest and Apoptosis and Upregulates SERPINE1 Expression Reducing Migration in Prostate Cancer Cells. Biomed Res Int 2020, 2020:2406159.

35.    Lorenzi L, Tabellini G, Vermi W, Moratto D, Porta F, Notarangelo LD, Patrizi O, Sozzani S, de Saint Basile G, Latour S et al: Occurrence of nodular lymphocyte-predominant hodgkin lymphoma in hermansky-pudlak type 2 syndrome is associated to natural killer and natural killer T cell defects. PLoS One 2013, 8(11):e80131.

36.    Miao Y, Zhu HY, Qiao C, Xia Y, Kong Y, Zou YX, Miao YQ, Chen X, Cao L, Wu W et al: Pathogenic Gene Mutations or Variants Identified by Targeted Gene Sequencing in Adults With Hemophagocytic Lymphohistiocytosis. Front Immunol 2019, 10:395.

37.    Zhang J, Cheng J, Zeng Z, Wang Y, Li X, Xie Q, Jia J, Yan Y, Guo Z, Gao J et al: Comprehensive profiling of novel microRNA-9 targets and a tumor suppressor role of microRNA-9 via targeting IGF2BP1 in hepatocellular carcinoma. Oncotarget 2015, 6(39):42040-42052.

38.    Xie T, Guo X, Wu D, Li S, Lu Y, Wang X, Chen L: PAFAH1B3 Expression Is Correlated With Gastric Cancer Cell Proliferation and Immune Infiltration. Front Oncol 2021, 11:591545.

39.    Xu J, Zang Y, Cao S, Lei D, Pan X: Aberrant expression of PAFAH1B3 associates with poor prognosis and affects proliferation and aggressiveness in hypopharyngeal squamous cell carcinoma. Onco Targets Ther 2019, 12:2799-2808.

40.    Fan J, Yang Y, Qian JK, Zhang X, Ji JQ, Zhang L, Li SZ, Yuan F: Aberrant Expression of PAFAH1B3 Affects Proliferation and Apoptosis in Osteosarcoma. Front Oncol 2021, 11:664478.

41.    Kumari R, Majumder MM, Lievonen J, Silvennoinen R, Anttila P, Nupponen NN, Lehmann F, Heckman CA: Prognostic significance of esterase gene expression in multiple myeloma. Br J Cancer 2021, 124(8):1428-1436.