A novel platelets-related gene signature for predicting prognosis, immune features and drug sensitivity in gastric cancer.

Platelets can dynamically regulate tumor development and progression. Nevertheless, research on the predictive value and specific roles of platelets in gastric cancer (GC) is limited. This research aims to establish a predictive platelets-related gene signature in GC with prognostic and therapeutic implications. We downloaded the transcriptome data and clinical materials of GC patients (n=378) from The Cancer Genome Atlas (TCGA) database. Prognostic platelets-related genes screened by univariate Cox regression were included in Least Absolute Shrinkage and Selection Operator (LASSO) analysis to construct a risk model. Kaplan-Meier curves and receiver operating characteristic curves (ROCs) were performed in the TCGA cohort and three independent validation cohorts. A nomogram integrating the risk score and clinicopathological features was constructed. Functional enrichment and tumor microenvironment (TME) analyses were performed. Drug sensitivity prediction was conducted through The Cancer Therapeutics Response Portal (CTRP) database. Finally, the expression of ten signature genes was validated by quantitative real-time PCR (qRT-PCR). A ten-gene (SERPINE1, ANXA5, DGKQ, PTPN6, F5, DGKB, PCDH7, GNG11, APOA1, and TF) predictive risk model was finally constructed. Patients were categorized as high- or low-risk using median risk score as the threshold. The area under the ROC curve (AUC) values for the 1-, 2-, and 3-year overall survival (OS) in the training cohort were 0.670, 0.695, and 0.707, respectively. Survival analysis showed a better OS in low-risk patients in the training and validation cohorts. The AUCs of the nomogram for predicting 1-, 2-, and 3-year OS were 0.708, 0.763, and 0.742, respectively. TME analyses revealed a higher M2 macrophage infiltration and an immunosuppressive TME in the high-risk group. Furthermore, High-risk patients tended to be more sensitive to thalidomide, MK-0752, and BRD-K17060750. The novel platelets-related genes signature we identified could be used for prognosis and treatment prediction in GC.

Li Q ，Zhang C ，Ren Y ，Qiao L ，Xu S ，Li K ，Liu Y ... -  《Frontiers in Immunology》

Identification of cuproptosis-related subtypes, construction of a prognosis model, and tumor microenvironment landscape in gastric cancer.

Cuproptosis is a novel identified regulated cell death (RCD), which is correlated with the development, treatment response and prognosis of cancer. However, the potential role of cuproptosis-related genes (CRGs) in the tumor microenvironment (TME) of gastric cancer (GC) remains unknown. Transcriptome profiling, somatic mutation, somatic copy number alteration and clinical data of GC samples were downloaded from the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database to describe the alterations of CRGs from genetic and transcriptional fields. Differential, survival and univariate cox regression analyses of CRGs were carried out to investigate the role of CRGs in GC. Cuproptosis molecular subtypes were identified by using consensus unsupervised clustering analysis based on the expression profiles of CRGs, and further analyzed by GO and KEGG gene set variation analyses (GSVA). Genes in distinct molecular subtypes were also analyzed by GO and KEGG gene enrichment analyses (GSEA). Differentially expressed genes (DEGs) were screened out from distinct molecular subtypes and further analyzed by GO enrichment analysis and univariate cox regression analysis. Consensus clustering analysis of prognostic DEGs was performed to identify genomic subtypes. Next, patients were randomly categorized into the training and testing group at a ratio of 1:1. CRG Risk scoring system was constructed through logistic least absolute shrinkage and selection operator (LASSO) cox regression analysis, univariate and multivariate cox analyses in the training group and validated in the testing and combined groups. Real-time quantitative polymerase chain reaction (RT-qPCR) was used to evaluate the expression of key Risk scoring genes. Sensitivity and specificity of Risk scoring system were examined by using receiver operating characteristic (ROC) curves. pRRophetic package in R was used to investigate the therapeutic effects of drugs in high- and low- risk score group. Finally, the nomogram scoring system was developed to predict patients' survival through incorporating the clinicopathological features and CRG Risk score. Most CRGs were up-regulated in tumor tissues and showed a relatively high mutation frequency. Survival and univariate cox regression analysis revealed that LIAS and FDX1 were significantly associated with GC patients' survival. After consensus unsupervised clustering analysis, GC patients were classified into two cuproptosis molecular subtypes, which were significantly associated with clinical features (gender, age, grade and TNM stage), prognosis, metabolic related pathways and immune cell infiltration in TME of GC. GO enrichment analyses of 84 DEGs, obtained from distinct molecular subtypes, revealed that DEGs primarily enriched in the regulation of metabolism and intracellular/extracellular structure in GC. Univariate cox regression analysis of 84 DEGs further screened out 32 prognostic DEGs. According to the expression profiles of 32 prognostic DEGs, patients were re-classified into two gene subtypes, which were significantly associated with patients' age, grade, T and N stage, and survival of patients. Nest, the Risk score system was constructed with moderate sensitivity and specificity. A high CRG Risk score, characterized by decreased microsatellite instability-high (MSI-H), tumor mutation burden (TMB) and cancer stem cell (CSC) index, and high stromal and immune score in TME, indicated poor survival. Four of five key Risk scoring genes expression were dysregulated in tumor compared with normal samples. Moreover, CRG Risk score was greatly related with sensitivity of multiple drugs. Finally, we established a highly accurate nomogram for promoting the clinical applicability of the CRG Risk scoring system. Our comprehensive analysis of CRGs in GC demonstrated their potential roles in TME, clinicopathological features, and prognosis. These findings may improve our understanding of CRGs in GC and provide new perceptions for doctors to predict prognosis and develop more effective and personalized therapy strategies.

Wang J ，Qin D ，Tao Z ，Wang B ，Xie Y ，Wang Y ，Li B ，Cao J ，Qiao X ，Zhong S ，Hu X ... -  《Frontiers in Immunology》

A Comprehensive Prognostic and Immunological Analysis of a Six-Gene Signature Associated With Glycolysis and Immune Response in Uveal Melanoma.

Uveal melanoma (UM) is a subtype of melanoma with poor prognosis. This study aimed to construct a new prognostic gene signature that can be used for survival prediction and risk stratification of UM patients. In this work, transcriptome data from the Molecular Signatures Database were used to identify the cancer hallmarks most relevant to the prognosis of UM patients. Weighted gene co-expression network, univariate least absolute contraction and selection operator (LASSO), and multivariate Cox regression analyses were used to construct the prognostic gene characteristics. Kaplan-Meier and receiver operating characteristic (ROC) curves were used to evaluate the survival predictive ability of the gene signature. The results showed that glycolysis and immune response were the main risk factors for overall survival (OS) in UM patients. Using univariate Cox regression analysis, 238 candidates related to the prognosis of UM patients were identified (p < 0.05). Using LASSO and multivariate Cox regression analyses, a six-gene signature including ARPC1B, BTBD6, GUSB, KRTCAP2, RHBDD3, and SLC39A4 was constructed. Kaplan-Meier analysis of the UM cohort in the training set showed that patients with higher risk scores had worse OS (HR = 2.61, p < 0.001). The time-dependent ROC (t-ROC) curve showed that the risk score had good predictive efficiency for UM patients in the training set (AUC > 0.9). Besides, t-ROC analysis showed that the predictive ability of risk scores was significantly higher than that of other clinicopathological characteristics. Univariate and multivariate Cox regression analyses showed that risk score was an independent risk factor for OS in UM patients. The prognostic value of risk scores was further verified in two external UM cohorts (GSE22138 and GSE84976). Two-factor survival analysis showed that UM patients with high hypoxia or immune response scores and high risk scores had the worst prognosis. Moreover, a nomogram based on the six-gene signature was established for clinical practice. In addition, risk scores were related to the immune infiltration profiles. Taken together, this study identified a new prognostic six-gene signature related to glycolysis and immune response. This six-gene signature can not only be used for survival prediction and risk stratification but also may be a potential therapeutic target for UM patients.

Liu J ，Lu J ，Li W 《Frontiers in Immunology》

A nomogram model based on the number of examined lymph nodes-related signature to predict prognosis and guide clinical therapy in gastric cancer.

Increasing evidence suggests that the number of examined lymph nodes (ELNs) is strongly linked to the survivorship of gastric cancer (GC). The goal of this study was to assess the prognostic implications of the ELNs number and to construct an ELNs-based risk signature and nomogram model to predict overall survival (OS) characteristics in GC patients. This inception cohort study included 19,317 GC patients from the U.S. Surveillance, Epidemiology, and End Results (SEER) database, who were separated into a training group and an internal validation group. The nomogram was built with the training set, then internally verified with SEER data, and externally validated with two different data sets. Based on the RNA-seq data, ELNs-related DERNAs (DElncRNAs, DEmiRNAs, andDEmRNAs) and immune cells were identified. The LASSO-Cox regression analysis was utilized to construct ELNs-related DERNAs and immune cell prognostic signature in The Cancer Genome Atlas (TCGA) cohort. The OS of subgroups with high- and low-ELN signature was compared using the Kaplan-Meier (K-M) analysis. A nomogram was successfully constructed based on the ELNs signature and other clinical characteristics. The concordance index (C-index), calibration plot, receiver operating characteristic curve, and decision curve analysis (DCA) were all used to evaluate the nomogram model. The meta-analysis, the Gene Expression Profiling Interactive Analysis database, and reverse transcription-quantitative PCR (RT-qPCR) were utilized to validate the RNA expression or abundance of prognostic genes and immune cells between GC tissues and normal gastric tissues, respectively. Finally, we analyzed the correlations between immune checkpoints, chemotherapy drug sensitivity, and risk score. The multivariate analysis revealed that the high ELNs improved OS compared with low ELNs (hazard ratio [HR] = 0.659, 95% confidence interval [CI]: 0.626-0.694, p < 0.0001). Using the training set, a nomogram incorporating ELNs was built and proven to have good calibration and discrimination (C-index [95% CI], 0.714 [0.710-0.718]), which was validated in the internal validation set (C-index [95% CI], 0.720 [0.714-0.726]), the TCGA set (C-index [95% CI], 0.693 [0.662-0.724]), and the Chinese set (C-index [95% CI], 0.750 [0.720-0.782]). An ELNs-related signature model based on ELNs group, regulatory T cells (Tregs), neutrophils, CDKN2B-AS1, H19, HOTTIP, LINC00643, MIR663AHG, TMEM236, ZNF705A, and hsa-miR-135a-5p was constructed by the LASSO-Cox regression analysis. The result showed that OS was remarkably lower in patients with high-ELNs signature compared with those with low-ELN signature (HR = 2.418, 95% CI: 1.804-3.241, p < 0.001). This signature performed well in predicting 1-, 3-, and 5-year survival (AUC [95% CI] = 0.688 [0.612-0.763], 0.744 [0.659-0.830], and 0.778 [0.647-0.909], respectively). The multivariate Cox analysis illustrated that the risk score was an independent predictor of survival for patients with GC. Moreover, the expression of prognostic genes (LINC00643, TMEM236, and hsa-miR-135a-5p) displayed differences between GC tissues and adjacent non-tumor tissues. The C-index of the nomogram that can be used to predict the OS of GC patients was 0.710 (95% CI: 0.663-0.753). Both the calibration plots and DCA showed that the nomogram has good predictive performance. Moreover, the signature was significantly correlated with the N stage and T stage. According to our analysis, GC patients in the low-ELN signature group may have a better immunotherapy response and OS outcome. We explored the prognostic role of ELNs in GC and successfully constructed an ELNs signature linked to the GC prognosis in TCGA. The findings manifested that the signature is a powerful predictive indicator for patients with GC. The signature might contain potential biomarkers for treatment response prediction for GC patients. Additionally, we identified a novel and robust nomogram combining the characteristics of ELNs and clinical factors for predicting 1-, 3-, and 5-year OS in GC patients, which will facilitate personalized survival prediction and aid clinical decision-making in GC patients.

Li H ，Lin D ，Yu Z ，Li H ，Zhao S ，Hainisayimu T ，Liu L ，Wang K ... -  《Frontiers in Immunology》

Construction of a novel mRNA-signature prediction model for prognosis of bladder cancer based on a statistical analysis.

Bladder cancer (BC) is a common malignancy neoplasm diagnosed in advanced stages in most cases. It is crucial to screen ideal biomarkers and construct a more accurate prognostic model than conventional clinical parameters. The aim of this research was to develop and validate an mRNA-based signature for predicting the prognosis of patients with bladder cancer. The RNA-seq data was downloaded from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Differentially expressed genes (DEGs) were screened in three datasets, and prognostic genes were identified from the training set of TCGA dataset. The common genes between DEGs and prognostic genes were narrowed down to six genes via Least Absolute Shrinkage and Selection Operator (LASSO) regression, and stepwise multivariate Cox regression. Then the gene-based risk score was calculated via Cox coefficient. Time-dependent receiver operating characteristic (ROC) and Kaplan-Meier (KM) survival analysis were used to assess the prognostic power of risk score. Multivariate Cox regression analysis was applied to construct a nomogram. Decision curve analysis (DCA), calibration curves, and time-dependent ROC were performed to assess the nomogram. Finally, functional enrichment of candidate genes was conducted to explore the potential biological pathways of candidate genes. SORBS2, GPC2, SETBP1, FGF11, APOL1, and H1-2 were screened to be correlated with the prognosis of BC patients. A nomogram was constructed based on the risk score, pathological stage, and age. Then, the calibration plots for the 1-, 3-, 5-year OS were predicted well in entire TCGA-BLCA patients. Decision curve analysis (DCA) indicated that the clinical value of the nomogram was higher than the stage model and TNM model in predicting overall survival analysis. The time-dependent ROC curves indicated that the nomogram had higher predictive accuracy than the stage model and risk score model. The AUC of nomogram time-dependent ROC was 0.763, 0.805, and 0.806 for 1-year, 3-year, and 5-year, respectively. Functional enrichment analysis of candidate genes suggested several pathways and mechanisms related to cancer. In this research, we developed an mRNA-based signature that incorporated clinical prognostic parameters to predict BC patient prognosis well, which may provide a novel prognosis assessment tool for clinical practice and explore several potential novel biomarkers related to the prognosis of patients with BC.

Li J ，Cao J ，Li P ，Yao Z ，Deng R ，Ying L ，Tian J ... -  《BMC CANCER》