A five-gene signature predicts overall survival of patients with papillary renal cell carcinoma

Identification of Cuproptosis-Associated Prognostic Gene Expression Signatures from 20 Tumor Types

We investigated the mRNA expression of 124 cuproptosis-associated genes in 7489 biopsies from 20 different tumor types of The Cancer Genome Atlas (TCGA). The KM plotter algorithm has been used to calculate Kaplan-Meier statistics and false discovery rate (FDR) corrections. Interaction networks have been generated using Ingenuity Pathway Analysis (IPA). High mRNA expression of 63 out of 124 genes significantly correlated with shorter survival times of cancer patients across all 20 tumor types. IPA analyses revealed that their gene products were interconnected in canonical pathways (e.g., cancer, cell death, cell cycle, cell signaling). Four tumor entities showed a higher accumulation of genes than the other cancer types, i.e., renal clear cell carcinoma (n = 21), renal papillary carcinoma (n = 13), kidney hepatocellular carcinoma (n = 13), and lung adenocarcinoma (n = 9). These gene clusters may serve as prognostic signatures for patient survival. These signatures were also of prognostic value for tumors with high mutational rates and neoantigen loads. Cuproptosis is of prognostic significance for the survival of cancer patients. The identification of specific gene signatures deserves further exploration for their clinical utility in routine diagnostics.

Discovery and Validation of Survival-Specific Genes in Papillary Renal Cell Carcinoma Using a Customized Next-Generation Sequencing Gene Panel

PURPOSE

Papillary renal cell carcinoma (PRCC), the second most common kidney cancer, is morphologically, genetically, and molecularly heterogeneous with diverse clinical manifestations. Genetic variations of PRCC and their association with survival are not yet well-understood. This study aimed to identify and validate survival-specific genes in PRCC and explore their clinical utility.

MATERIALS AND METHODS

Using machine learning, 293 patients from the Cancer Genome Atlas-Kidney Renal Papillary Cell Carcinoma (TCGA-KIRP) database were analyzed to derive genes associated with survival. To validate these genes, DNAs were extracted from the tissues of 60 Korean PRCC patients. Next generation sequencing was conducted using a customized PRCC gene panel of 202 genes, including 171 survival-specific genes. Kaplan-Meier and Log-rank tests were used for survival analysis. Fisher's exact test was performed to assess the clinical utility of variant genes.

RESULTS

A total of 40 survival-specific genes were identified in the TCGA-KIRP database through machine learning and statistical analysis. Of them, 10 (BAP1, BRAF, CFDP1, EGFR, ITM2B, JAK1, NODAL, PCSK2, SPATA13, and SYT5) were validated in the Korean-KIRP database. Among these survival gene signatures, three genes (BAP1, PCSK2, and SPATA13) showed survival specificity in both overall survival (OS) (p = 0.00004, p = 1.38 × 10-7, and p = 0.026, respectively) and disease-free survival (DFS) (p = 0.00002, p = 1.21 × 10-7, and p = 0.036, respectively). Notably, the PCSK2 mutation demonstrated survival specificity uniquely in both the TCGA-KIRP (OS: p = 0.010 and DFS: p = 0.301) and Korean-KIRP (OS: p = 1.38 × 10-7 and DFS: p = 1.21 × 10-7) databases.

CONCLUSIONS

We discovered and verified genes specific for the survival of PRCC patients in the TCGA-KIRP and Korean-KIRP databases. The survival gene signature, including PCSK2 commonly obtained from the 40 gene signature of TCGA and the 10 gene signature of the Korean database, is expected to provide insight into predicting the survival of PRCC patients and developing new treatment.

BUB1/KIF14 complex promotes anaplastic thyroid carcinoma progression by inducing chromosome instability

Chromosome instability (CIN) is a common contributor driving the formation and progression of anaplastic thyroid cancer (ATC), but its mechanism remains unclear. The BUB1 mitotic checkpoint serine/threonine kinase (BUB1) is responsible for the alignment of mitotic chromosomes, which has not been thoroughly studied in ATC. Our research demonstrated that BUB1 was remarkably upregulated and closely related to worse progression-free survival. Knockdown of BUB1 attenuated cell viability, invasion, migration and induced cell cycle arrests, whereas overexpression of BUB1 promoted the cell cycle progression of papillary thyroid cancer cells. BUB1 knockdown remarkably repressed tumour growth and tumour formation of nude mice with ATC xenografts and suppressed tumour metastasis in a zebrafish xenograft model. Inhibition of BUB1 by its inhibitor BAY-1816032 also exhibited considerable anti-tumour activity. Further studies showed that enforced expression of BUB1 evoked CIN in ATC cells. BUB1 induced CIN through phosphorylation of KIF14 at serine1292 (Ser1292 ). Overexpression of the KIF14ΔSer1292 mutant was unable to facilitate the aggressiveness of ATC cells when compared with that of the wild type. Collectively, these findings demonstrate that the BUB1/KIF14 complex drives the aggressiveness of ATC by inducing CIN.

Construction and validation of an autophagy-related long non-coding RNA signature to predict the prognosis of kidney renal papillary cell carcinoma

To identify the autophagy-related long non-coding RNAs (ARlncRNAs) associated with the prognosis of kidney renal papillary cell carcinoma (KIRP), thereby establishing a clinical prognostic model. The gene expression matrix and clinical survival information of patients with KIRP were downloaded from The Cancer Genome Atlas database, and were divided into the training and testing groups. ARlncRNAs associated with the KIRP prognosis were analyzed by univariate, Least Absolute Shrinkage and Selection Operator (LASSO(, and multivariate Cox regression to construct a signature. We combined clinical factors associated with the prognosis with ARlncRNAs to establish a prognostic model of patients with KIRP. A nomogram was established to predict 1-year, 3-year, and 5-year survival of patients with KIRP. Besides, we built the lncRNA-messenger RNA co-expression network and used Kyoto Encyclopedia of Genes and Genomes and Gene Set Enrichment Analysis to detect the biological functions of ARlncRNAs. LEF1-AS1, CU634019.6, C2orf48, AC027228.2, and AC107464.3 were identified. A prognosis-related ARlncRNAs signature was constructed in the training group and validated in the testing group. Patients with KIRP with a low risk score had significantly longer survival time than those with a high risk score. The risk score significantly affected the prognosis of patients, thereby being used for modeling. The area under the receiver operating characteristic curve values of 1-year, 3-year, and 5-year overall survival were 0.80, 0.78, and 0.84 in the training group, respectively. The signature had high concordance index and good accuracy in predicting the prognosis, which were confirmed by the nomogram. The prognosis-related ARlncRNAs signature we identified had a more accurate prediction for the prognosis of patients with KIRP.

Comprehensive Analyses of Mutation-Derived Long-Chain Noncoding RNA Signatures of Genome Instability in Kidney Renal Papillary Cell Carcinoma

Background: The role of long-chain noncoding RNA (lncRNA) in genomic instability has been demonstrated to be increasingly importance. Therefore, in this study, lncRNAs associated with genomic instability were identified and kidney renal papillary cell carcinoma (KIRP)-associated predictive features were analysed to classify high-risk patients and improve individualised treatment.

Methods: The training (n = 142) and test (n = 144) sets were created using raw RNA-seq and patient’s clinical data of KIRP obtained from The Cancer Genome Atlas (TCGA).There are 27 long-chain noncoding RNAs (lncRNAs) that are connected with genomic instability, these lncRNAs were identified using the ‘limma’ R package based on the numbers of somatic mutations and lncRNA expression profiles acquired from KIRP TCGA cohort. Furthermore, Cox regression analysis was carried out to develop a genome instability-derived lncRNA-based gene signature (GILncSig), whose prognostic value was confirmed in the test cohort as well as across the entire KIRP TCGA dataset.

Results: A GILncSig derived from three lncRNAs (BOLA3-AS1, AC004870, and LINC00839), which were related with poor KIRP survival, was identified, which was split up into high- and low-risk groups. Additionally, the GILncSig was found to be an independent prognostic predictive index in KIRP using univariate and multivariate Cox analysis. Furthermore, the prognostic significance and characteristics of GilncSig were confirmed in the training test and TCGA sets. GilncSig also showed better predictive performance than other prognostic lncRNA features.

Conclusion: The function of lncRNAs in genomic instability and the genetic diversity of KIRP were elucidated in this work. Moreover, three lncRNAs were screened for prediction of the outcome of KIRP survival and novel insights into identifying cancer biomarkers related to genomic instability were discussed.

Identification of the Ferroptosis-Related Long Non-Coding RNAs Signature to Improve the Prognosis Prediction in Papillary Renal Cell Carcinoma

Papillary renal cell carcinoma (pRCC) is one of the epithelial renal cell carcinoma (RCC) histological subtypes. Ferroptosis is a new iron-dependent form of cell death that has been seen in a variety of clinical situations. Using differentially expressed ferroptosis-related long non-coding RNAs (lncRNAs) from patients with pRCC in The Cancer Genome Atlas; we built a prognostic lncRNA-based signature. We discovered seven different lncRNAs that were strongly linked to the prognosis of patients with pRCC. High-risk scores were linked to a poor prognosis for pRCC, which was confirmed by the findings of Kaplan–Meier studies. In addition, the constructed lncRNA signature has a 1-year area under the curve (AUC) of 0.908, suggesting that it has a high predictive value in pRCC. In the high-risk group, Gene set enrichment analyses (GSEA) analysis identified immunological and tumor-related pathways. Furthermore, single-sample GSEA (ssGSEA) revealed significant differences in T cell functions checkpoint, antigen presenting cell (APC) co-stimulation, inflammation promoting, and para inflammation between the two groups with different risk scores. In addition, immune checkpoints like PDCD1LG2 (PD-L2), LAG3, and IDO1 were expressed differently in the two risk groups. In summary, a novel signature based on ferroptosis-related lncRNAs could be applied in predicting the prognosis of patients with pRCC.

A five-gene methylation signature predicts overall survival of patients with clear cell renal cell carcinoma

BACKGROUND

In this study, we aimed to screen methylation signatures associated with the prognosis of patients with clear cell renal cell carcinoma (ccRCC).

METHODS

Gene expression and methylation profiles of ccRCC patients were downloaded from publicly available databases, and differentially expressed genes (DEGs)-differentially methylated genes (DMGs) were obtained. Subsequently, gene set enrichment and transcription factor (TF) regulatory network analyses were performed. In addition, a prognostic model was constructed and the relationship between disease progression and immunity was analyzed.

RESULTS

A total of 23 common DEGs-DMGs were analyzed, among which 14 DEGs-DMGs were obtained with a cutoff value of PCC < 0 and p < 0.05. The enrichment analysis showed that the 14 DEGs-DMGs were enriched in three GO terms and three KEGG pathways. In addition, a total of six TFs were shown to be associated with the 14 DEGs-DMGs, including RP58, SOX9, NF-κB65, ATF6, OCT, and IK2. A prognostic model using five optimized DEGs-DMGs which efficiently predicted survival was constructed and validated using the GSE105288 dataset. Additionally, four types of immune cells (NK cells, macrophages, neutrophils, and cancer-associated fibroblasts), as well as ESTIMATE, immune, and stromal scores were found to be significantly correlated with ccRCC progression (normal, primary, and metastasis) in addition to the five optimized DEGs-DMGs.

CONCLUSION

A five-gene methylation signature with the predictive ability for ccRCC prognosis was investigated in this study, consisting of CCNB2, CDKN1C, CTSH, E2F2, and ERMP1. In addition, potential targets for methylation-mediated immunotherapy were highlighted.

A novel prognostic cancer-related lncRNA signature in papillary renal cell carcinoma

Background

Papillary renal cell carcinoma (pRCC) ranks second in renal cell carcinoma and the prognosis of pRCC remains poor. Here, we aimed to screen and identify a novel prognostic cancer-related lncRNA signature in pRCC.

Methods

The RNA-seq profile and clinical feature of pRCC cases were downloaded from TCGA database. Significant cancer-related lncRNAs were obtained from the Immlnc database. Differentially expressed cancer-related lncRNAs (DECRLs) in pRCC were screened for further analysis. Cox regression report was implemented to identify prognostic cancer-related lncRNAs and establish a prognostic risk model, and ROC curve analysis was used to evaluate its precision. The correlation between RP11-63A11.1 and clinical characteristics was further analyzed. Finally, the expression level and role of RP11-63A11.1 were studied in vitro.

Results

A total of 367 DECRLs were finally screened and 26 prognostic cancer-related lncRNAs were identified. Among them, ten lncRNAs (RP11-573D15.8, LINC01317, RNF144A-AS1, TFAP2A-AS1, LINC00702, GAS6-AS1, RP11-400K9.4, LUCAT1, RP11-63A11.1, and RP11-156L14.1) were independently associated with prognosis of pRCC. These ten lncRNAs were incorporated into a prognostic risk model. In accordance with the median value of the riskscore, pRCC cases were separated into high and low risk groups. Survival analysis indicated that there was a significant difference on overall survival (OS) rate between the two groups. The area under curve (AUC) in different years indicated that the model was of high efficiency in prognosis prediction. RP11-63A11.1 was mainly expressed in renal tissues and it correlated with the tumor stage, T, M, N classifications, OS, PFS, and DSS of pRCC patients. Consistent with the expression in pRCC tissue samples, RP11-63A11.1 was also down-regulated in pRCC cells. More importantly, up-regulation of RP11-63A11.1 attenuated cell survival and induced apoptosis.

Conclusions

Ten cancer-related lncRNAs were incorporated into a powerful model for prognosis evaluation. RP11-63A11.1 functioned as a cancer suppressor in pRCC and it might be a potential therapeutic target for treating pRCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02247-6.

NS398 as a potential drug for autosomal-dominant polycystic kidney disease: Analysis using bioinformatics, and zebrafish and mouse models

Autosomal‐dominant polycystic kidney disease (ADPKD) is characterized by uncontrolled renal cyst formation, and few treatment options are available. There are many parallels between ADPKD and clear‐cell renal cell carcinoma (ccRCC); however, few studies have addressed the mechanisms linking them. In this study, we aimed to investigate their convergences and divergences based on bioinformatics and explore the potential of compounds commonly used in cancer research to be repurposed for ADPKD. We analysed gene expression datasets of ADPKD and ccRCC to identify the common and disease‐specific differentially expressed genes (DEGs). We then mapped them to the Connectivity Map database to identify small molecular compounds with therapeutic potential. A total of 117 significant DEGs were identified, and enrichment analyses results revealed that they are mainly enriched in arachidonic acid metabolism, p53 signalling pathway and metabolic pathways. In addition, 127 ccRCC‐specific up‐regulated genes were identified as related to the survival of patients with cancer. We focused on the compound NS398 as it targeted DEGs and found that it inhibited the proliferation of Pkd1^−/− and 786‐0 cells. Furthermore, its administration curbed cystogenesis in Pkd2 zebrafish and early‐onset Pkd1‐deficient mouse models. In conclusion, NS398 is a potential therapeutic agent for ADPKD.

Inhibition of BUB1 suppresses tumorigenesis of osteosarcoma via blocking of PI3K/Akt and ERK pathways

Osteosarcoma (OS) is a primary malignant bone tumour that mainly affects teenagers, with patients displaying poor prognosis. Budding uninhibited by benzimidazoles 1 (BUB1), a type of serine/threonine kinase that is linked to pro-tumorigenic phenomena, has not been well studied in OS. Hence, this study aimed to explore the role of BUB1 in OS. The expression of BUB1 in OS specimens and cell lines was assessed using immunohistochemistry and Western blot analysis. Univariate and multivariate analyses were applied to evaluate the impact of BUB1 on patient survival. Cell counting kit-8, wound-healing and Transwell assays, as well as flow cytometry, were used to investigate the influence of BUB1 inhibition on OS in vitro. Moreover, a tumour xenograft model was established to investigate the in vivo effect of BUB1 inhibition on OS tumour growth. Results showed that BUB1 was overexpressed in OS specimens and cell lines. Furthermore, BUB1 overexpression was closely associated with the poor clinical outcomes of patients with OS. Inhibition of BUB1 markedly suppressed cell proliferation and tumour growth, cell migration, invasion and induced cell apoptosis of OS by blocking the PI3K/Akt and ERK signalling pathways. Thus, our study suggested that overexpression of BUB1 protein contributed to poor survival of OS patients and that inhibition of BUB1 resulted in considerable anti-tumour activity associated with proliferation, migration, invasion and apoptosis of OS.

An Immune-Related Signature Predicted Survival in Patients With Kidney Papillary Cell Carcinoma

Immune-related genes are important factors in tumor progression. The main aim of this study was to identify the immune-related genes in kidney papillary cell carcinoma (pRCC) patients. We downloaded RNAseq data and clinical information of pRCC patients from the TCGA database and retrieved the immune-related genes list from Immport. From the data, we mined out 2,468 differential expression genes (DEGs) and 183 immune-related DEGs. Four hub DEGs (NTS, BIRC5, ELN, and CHGA) were identified after conducting Cox analysis and LASSO analysis. Moreover, the prognostic value of the signature based on four hub DEGs was verified using Kaplan–Meier analysis (P = 0.0041 in the training set and p = 0.021 in the test set) and ROC analysis (AUC: 0.957 in 1 year, 0.965 in 2 years, and 0.901 in 3 years in the training set, and 0.963 in 1 year, 0.898 in 2 years, and 0.742 in 3 years in the test set). Furthermore, we found that the high-risk score group had a higher percentage of B cells in the immune component, a higher expression of immune-related genes (CTLA4, LAG3, PDCD1LG2, and TIGIT), and a better immunotherapy response.

In-Silico Evaluation of Genetic Alterations in Ovarian Carcinoma and Therapeutic Efficacy of NSC777201, as a Novel Multi-Target Agent for TTK, NEK2, and CDK1

Ovarian cancer is often detected at the advanced stages at the time of initial diagnosis. Early-stage diagnosis is difficult due to its asymptomatic nature, where less than 30% of 5-year survival has been noticed. The underlying molecular events associated with the disease’s pathogenesis have yet to be fully elucidated. Thus, the identification of prognostic biomarkers as well as developing novel therapeutic agents for targeting these markers become relevant. Herein, we identified 264 differentially expressed genes (DEGs) common in four ovarian cancer datasets (GSE14407, GSE18520, GSE26712, GSE54388), respectively. We constructed a protein-protein interaction (PPI) interaction network with the overexpressed genes (72 genes) and performed gene enrichment analysis. In the PPI networks, three proteins; TTK Protein Kinase (TTK), NIMA Related Kinase 2 (NEK2), and cyclin-dependent kinase (CDK1) with higher node degrees were further evaluated as therapeutic targets for our novel multi-target small molecule NSC777201. We found that the upregulated DEGs were enriched in KEGG and gene ontologies associated with ovarian cancer progression, female gamete association, otic vesicle development, regulation of chromosome segregation, and therapeutic failure. In addition to the PPI network, ingenuity pathway analysis also implicate TTK, NEK2, and CDK1 in the elevated salvage pyrimidine and pyridoxal pathways in ovarian cancer. The TTK, NEK2, and CDK1 are over-expressed, demonstrating a high frequency of genetic alterations, and are associated with poor prognosis of ovarian cancer cohorts. Interestingly, NSC777201 demonstrated anti-proliferative and cytotoxic activities (GI₅₀ = 1.6 µM~1.82 µM and TGI₅₀ = 3.5 µM~3.63 µM) against the NCI panels of ovarian cancer cell lines and exhibited a robust interaction with stronger affinities for TTK, NEK2, and CDK1, than do the standard drug, paclitaxel. NSC777201 displayed desirable properties of a drug-like candidate and thus could be considered as a novel small molecule for treating ovarian carcinoma.

Probabilistic ratiocination of hepatocellular carcinoma after resection: evaluation of expected to be promising approaches

Background

Precise prediction of survival after treatment is of great importance for patients with diseases with high mortality. RNA sequencing data and deep learning (DL) methods are expected to become promising approaches in the development of prediction models in the future. We aimed to evaluate the optimal covariates and methodology for patients with hepatocellular carcinoma (HCC) undergoing surgical resection.

Methods

The Cox proportional hazards regression model and the DL approach were used to develop prediction models incorporating clinical, genetic, and combined clinical and genetic variables for survival prediction in patients with HCC after resection. A total of 1,114 patients and 184 patients were enrolled in the present study from 2,163 and 601 patients from Eastern Hepatobiliary Surgery Hospital and Renji Hospital, respectively. The models were internally validated through random sampling and externally validated in clinical cohorts. Between-model comparisons were carried out in terms of the integrated discrimination improvement and net reclassification index.

Results

The Cox and DL clinical models were developed by adopting 7 independent prognostic factors (total bilirubin, prothrombin time, tumor size, tumor number, lymph node metastasis, and vascular invasion) and 22 clinical factors, respectively. Both the Cox clinical model and the DL clinical model showed excellent performances in the derivation [area under the curve (AUC): 0.75 vs. 0.77] and validation (AUC: 0.83 vs. 0.80) sets. The derived Cox genetic model with 6 significant prognostic genes was not as effective as the DL approach involving 686 genes. A combined clinical and genetic approach modified the performances of both the Cox and DL models. The integrated discrimination improvement and net reclassification index of the DL clinical model were generally better than those of the Cox clinical model.

Conclusions

Our Cox clinical model sufficiently provided precise survival prediction in patients with HCC after resection. It may serve as an accurate and cost-effective tool for predicting survival in such patients.

Prognostic and therapeutic implications of extracellular matrix associated gene signature in renal clear cell carcinoma

Complex interactions in tumor microenvironment between ECM (extra-cellular matrix) and cancer cell plays a central role in the generation of tumor supportive microenvironment. In this study, the expression of ECM-related genes was explored for prognostic and immunological implication in clear cell renal clear cell carcinoma (ccRCC). Out of 964 ECM genes, higher expression (z-score > 2) of 35 genes showed significant association with overall survival (OS), progression-free survival (PFS) and disease-specific survival (DSS). On comparison to normal tissue, 12 genes (NUDT1, SIGLEC1, LRP1, LOXL2, SERPINE1, PLOD3, ZP3, RARRES2, TGM2, COL3A1, ANXA4, and POSTN) showed elevated expression in kidney tumor (n = 523) compared to normal (n = 100). Further, Cox proportional hazard model was utilized to develop 12 genes ECM signature that showed significant association with overall survival in TCGA dataset (HR = 2.45; 95% CI [1.78-3.38]; p < 0.01). This gene signature was further validated in 3 independent datasets from GEO database. Kaplan-Meier log-rank test significantly associated patients with elevated expression of this gene signature with a higher risk of mortality. Further, differential gene expression analysis using DESeq2 and principal component analysis (PCA) identified genes with the highest fold change forming distinct clusters between ECM-rich high-risk and ECM-poor low-risk patients. Geneset enrichment analysis (GSEA) identified significant perturbations in homeostatic kidney functions in the high-risk group. Further, higher infiltration of immunosuppressive T-reg and M2 macrophages was observed in high-risk group patients. The present study has identified a prognostic signature with associated tumor-promoting immune niche with clinical utility in ccRCC. Further exploration of ECM dynamics and validation of this gene signature can assist in design and application of novel therapeutic approaches.

Development and Clinical Validation of a Seven-Gene Prognostic Signature Based on Multiple Machine Learning Algorithms in Kidney Cancer

About a third of patients with kidney cancer experience recurrence or cancer-related progression. Clinically, kidney cancer prognoses may be quite different, even in patients with kidney cancer at the same clinical stage. Therefore, there is an urgent need to screen for kidney cancer prognosis biomarkers. Differentially expressed genes (DEGs) were identified using kidney cancer RNA sequencing data from the Gene Expression Omnibus (GEO) database. Biomarkers were screened using random forest (RF) and support vector machine (SVM) models, and a multigene signature was constructed using the least absolute shrinkage and selection operator (LASSO) regression analysis. Univariate and multivariate Cox regression analyses were performed to explore the relationships between clinical features and prognosis. Finally, the reliability and clinical applicability of the model were validated, and relationships with biological pathways were identified. Western blots were also performed to evaluate gene expression. A total of 50 DEGs were obtained by intersecting the RF and SVM models. A seven-gene signature (RNASET2, EZH2, FXYD5, KIF18A, NAT8, CDCA7, and WNT7B) was constructed by LASSO regression. Univariate and multivariate Cox regression analyses showed that the seven-gene signature was an independent prognostic factor for kidney cancer. Finally, a predictive nomogram was established in The Cancer Genome Atlas (TCGA) cohort and validated internally. In tumor tissue, RNASET2 and FXYD5 were highly expressed and NAT8 was lowly expressed at the protein and transcription levels. This model could complement the clinicopathological characteristics of kidney cancer and promote the personalized management of patients with kidney cancer.

A gene-based risk score model for predicting recurrence-free survival in patients with hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) remains the most frequent liver cancer, accounting for approximately 90% of primary liver cancers worldwide. The recurrence-free survival (RFS) of HCC patients is a critical factor in devising a personal treatment plan. Thus, it is necessary to accurately forecast the prognosis of HCC patients in clinical practice.

METHODS

Using The Cancer Genome Atlas (TCGA) dataset, we identified genes associated with RFS. A robust likelihood-based survival modeling approach was used to select the best genes for the prognostic model. Then, the GSE76427 dataset was used to evaluate the prognostic model's effectiveness.

RESULTS

We identified 1331 differentially expressed genes associated with RFS. Seven of these genes were selected to generate the prognostic model. The validation in both the TCGA cohort and GEO cohort demonstrated that the 7-gene prognostic model can predict the RFS of HCC patients. Meanwhile, the results of the multivariate Cox regression analysis showed that the 7-gene risk score model could function as an independent prognostic factor. In addition, according to the time-dependent ROC curve, the 7-gene risk score model performed better in predicting the RFS of the training set and the external validation dataset than the classical TNM staging and BCLC. Furthermore, these seven genes were found to be related to the occurrence and development of liver cancer by exploring three other databases.

CONCLUSION

Our study identified a seven-gene signature for HCC RFS prediction that can be used as a novel and convenient prognostic tool. These seven genes might be potential target genes for metabolic therapy and the treatment of HCC.

Identification of the 11-lncRNA signatures associated with the prognosis of endometrial carcinoma

Endometrial carcinoma (EC) is the fourth most common cancer in women. Some long non-coding RNAs (lncRNAs) are regarded as potential prognostic biomarkers or targets for treatment of many types of cancers. We aim to screen prognostic-related lncRNAs and build a possible lncRNA signature which can effectively predict the survival of patients with EC. We obtained lncRNA expression profiling from the TCGA database. The patients were classified into training set and verification set. By performing Univariate Cox regression model, Robust likelihood-based survival analysis, and Cox proportional hazards model, we developed a risk score with the Cox co-efficient of individual lncRNAs in the training set. The optimum cut-off point was selected by ROC analysis. Patients were effectively divided into high-risk group and low-risk group according to the risk score. The OS of the low-risk patients was significantly prolonged compared with that of the high-risk group. At last, we validated this 11-lncRNA signature in the verification set and the complete set. We identified an 11-lncRNA expression signature with high stability and feasibility, which can predict the survival of patients with EC. These findings provide new potential biomarkers to improve the accuracy of prognosis prediction of EC.

cg04448376, cg24387542, cg08548498, and cg14621323 as a Novel Signature to Predict Prognosis in Kidney Renal Papillary Cell Carcinoma

INTRODUCTION

DNA methylation plays a vital role in prognosis prediction of cancers. In this study, we aimed to identify novel DNA methylation site biomarkers and create an efficient methylated site model for predicting survival in kidney renal papillary cell carcinoma (KIRP).

METHODS

DNA methylation and gene expression profile data were downloaded from The Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) database. Differential methylated genes (DMGs) and differential expression genes (DEGs) were identified and then searched for the hub genes. Cox proportional hazards regression was applied to identify DNA methylated site biomarkers from the hub genes. Kaplan-Meier survival and ROC analyses were used to validate the effective prognostic ability of the methylation gene site biomarker. The biomarker sites were validated in the GEO cohorts. The GO and KEGG annotation was done to explore the biological function of DNA methylated site signature.

RESULTS

Nine DMGs with opposite expression patterns containing 47 methylated sites were identified. Finally, four methylated sites were identified using the hazard regression model (cg04448376, cg24387542, cg08548498, and cg14621323) located in UTY, LGALS9B, SLPI, and PFN3, respectively. These sites classified patients into high- and low-risk groups in the training cohort. The 5-year survival rates for patients with low-risk and high-risk scores were 97.5% and 75.9% (P < 0.001). The prognostic accuracy and signature methylation sites were validated in the test (TCGA, n = 87) and GEO cohorts (n = 14). Multivariate regression analysis showed that the signature was an independent prediction prognostic factor for KIRP. Based on this analysis, we developed methylated site signature nomogram that predicts an individual's risk of survival. Functional analysis suggested that these signature genes are involved in the biological processes of protein binding.

CONCLUSIONS

Our study demonstrated that the methylated gene site signature might be a powerful prognostic tool for evaluating survival rate and guiding tailored therapy for KIRP patients.

Development and Validation of an Individualized Immune-Related Gene Pairs Prognostic Signature in Papillary Renal Cell Carcinoma

Papillary renal carcinoma (PRCC) is one of the important subtypes of kidney cancer, with a high degree of heterogeneity. At present, there is still a lack of robust and accurate biomarkers for the diagnosis, prognosis and treatment selection of PRCC. Considering the important role of tumor immunity in PRCC, we aim to construct a signature based on immune-related gene pairs (IRGPs) to estimate the prognostic of patients with PRCC. We obtained gene expression profiling and clinical information of patients with PRCC from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), which were divided into discovery (n = 287) and validation (n = 28) cohorts, respectively. By univariate analysis, multivariate Cox analysis, and least absolute shrinkage and selection operator (Lasso) analysis, we selected 14 IRGPs with a panel of 22 unique genes to construct the prognostic signature. According to the signature, we stratified patients into high-risk group and low-risk group. In both discovery and validation cohorts, the results of Kaplan-Meier analysis showed that there were significant differences in OS between the two groups (p < 0.001). Combined with multiple clinical and pathological factors, the results of multivariate analyses confirmed that this signature was an independent predictor of OS (HR, 3.548; 95%CI, 2.096-6.006; p < 0.001). The results of immune infiltration analysis demonstrated that the abundance of multiple tumor-infiltration lymphocytes such as CD8 + T cells, Tregs, and T follicular cell helper were significantly higher in the high-risk group. Functional analysis showed that multiple immune-related signaling pathways were enriched in the high-risk group. In conclusion, we successfully established an individualized prognostic IRGPs signature, which can accurately assess and predict the OS of patients with PRCC.

Three-gene risk model in papillary renal cell carcinoma: a robust likelihood-based survival analysis

Background: Papillary renal cell carcinoma (PRCC) accounts for 15% of all renal cell carcinomas. The molecular mechanisms of renal papillary cell carcinoma remain unclear, and treatments for advanced disease are limited.

Result: We built the computing model as follows: Risk score = 1.806 * TPX2 - 0.355 * TXNRD2 - 0.805 * SLC6A20. The 3-year AUC of overall survival was 0.917 in the training set (147 PRCC samples) and 0.760 in the test set (142 PRCC samples). Based on the robust model, M2 macrophages showed positive correlation with risk score, while M1 macrophages were the opposite. PRCC patients with low risk score showed higher tumor mutation burden. TPX2 is a risk factor, and co-expression factors were enriched in cell proliferation and cancer-related pathways. Finally, the proliferation and invasion of PRCC cell line were decreased in the TPX2 reduced group, and the differential expression was identified. TPX2 is a potential risk biomarker which involved in cell proliferation in PRCC.

Conclusion: We conducted a study to develop a three gene model for predicting prognosis in patients with papillary renal cell carcinoma. Our findings may provide candidate biomarkers for prognosis that have important implications for understanding the therapeutic targets of papillary renal cell carcinoma.

Method: Gene expression matrix and clinical data were obtained from TCGA (The Cancer Genome Atlas), GSE26574, GSE2048, and GSE7023. Prognostic factors were identified using “survival” and “rbsurv” packages, and a risk score was constructed using Multivariate Cox regression analysis. The co-expression networks of the factors in model were constructed using the “WGCNA” package. The co-expression genes of factors were enriched and displayed the biological process. Based on this robust risk model, immune cells infiltration proportions and tumor mutation burdens were compared between risk groups. Subsequently, using the PRCC cell line, the role of TPX2 was determined by Cell proliferation assay, 5-Ethynyl-20-deoxyuridine assay and Transwell assay.

Mitophagy-associated genes PINK1 and PARK2 are independent prognostic markers of survival in papillary renal cell carcinoma and associated with aggressive tumor behavior

The aim of this study was to investigate the mitophagy-related genes PINK1 and PARK2 in papillary renal cell carcinoma and their association with prognosis. In silico data of PINK1 and PARK2 were analyzed in TCGA cohorts of papillary renal cell carcinoma comprising 290 tumors and 33 corresponding non-neoplastic renal tissues. Protein expression data from a cohort of 95 papillary renal cell carcinoma patients were analyzed and associated with clinical-pathological parameters including survival. PINK1 and PARK2 were significantly downregulated in papillary renal cell carcinoma at transcript and protein levels. Reduced transcript levels of PINK1 and PARK2 were negatively associated with overall survival (p < 0.05). At the protein level, PARK2 and PINK1 expression were positively correlated (correlation coefficient 0.286, p = 0.04) and reduced PINK1 protein expression was prognostic for shorter survival. Lower PINK1 protein levels were found in tumors with metastases at presentation and in tumors of higher pT-stages. The multivariate analysis revealed mRNA expression of PINK1 and PARK2 as well as PINK1 protein expression as independent prognostic factors for shorter overall survival. The downregulation of PINK1 is a strong predictor of poor survival in papillary renal cell carcinoma. Immunohistochemical PINK1 expression in resected pRCC should be considered as an additional prognostic marker for routine practice.

A risk stratification model based on four novel biomarkers predicts prognosis for patients with renal cell carcinoma

BACKGROUND

Accurate prediction of the prognosis of RCC using a single biomarker is challenging due to the genetic heterogeneity of the disease. However, it is essential to develop an accurate system to allow better patient selection for optimal treatment strategies. ARL4C, ECT2, SOD2, and STEAP3 are novel molecular biomarkers identified in earlier studies as survival-related genes by comprehensive analyses of 43 primary RCC tissues and RCC cell lines.

METHODS

To develop a prognostic model based on these multiple biomarkers, the expression of four biomarkers ARL4C, ECT2, SOD2, and STEAP3 in primary RCC tissue were semi-quantitatively investigated by immunohistochemical analysis in an independent cohort of 97 patients who underwent nephrectomy, and the clinical significance of these biomarkers were analyzed by survival analysis using Kaplan-Meier curves. The prognostic model was constructed by calculation of the contribution score to prognosis of each biomarker on Cox regression analysis, and its prognostic performance was validated.

RESULTS

Patients whose tumors had high expression of the individual biomarkers had shorter cancer-specific survival (CSS) from the time of primary nephrectomy. The prognostic model based on four biomarkers segregated the patients into a high- and low-risk scored group according to defined cut-off value. This approach was more robust in predicting CSS compared to each single biomarker alone in the total of 97 patients with RCC. Especially in the 36 metastatic RCC patients, our prognostic model could more accurately predict early events within 2 years of diagnosis of metastasis. In addition, high risk-scored patients with particular strong SOD2 expression had a much worse prognosis in 25 patients with metastatic RCC who were treated with molecular targeting agents.

CONCLUSIONS

Our findings indicate that a prognostic model based on four novel biomarkers provides valuable data for prediction of clinical prognosis and useful information for considering the follow-up conditions and therapeutic strategies for patients with primary and metastatic RCC.

Bioinformatics profiling integrating a four immune-related long non-coding RNAs signature as a prognostic model for papillary renal cell carcinoma

Background: Papillary renal cell carcinoma (pRCC) was the 2^nd most common subtype, accounting for approximately 15% incidence of renal cell carcinoma (RCC). Immune related long non-coding RNAs (IR-lncRs) plentiful in immune cells and immune microenvironment (IME) are potential in evaluating prognosis and assessing the effects of immunotherapy. A completed and meaningful IR-lncRs analysis based on abundant pRCC gene samples from The Cancer Genome Atlas (TCGA) will provide insight in this field.

Results: 17 IR-lncRs were selected by Pearson correlation analysis of immune score and the lncRNA expression level, and 5 sIRlncRs were significantly correlated with the OS of pRCC patients. 4 sIRlncRs (AP001267.3, AC026471.3, SNHG16 and ADAMTS9-AS1) with the most remarkable prognostic values were identified to establish the IRRS model and the OS of the low-risk group was longer than that in the high-risk group. The IRRS was certified as an independent prognosis factor and correlated with the OS. The high-risk group and low-risk group showed significantly different distributions and immune status through PCA and GSEA. In addition, we further found the expression levels of SNHG16 was remarkably enhanced in female patients with more advanced T-stages, but ADAMTS9-AS1 showed the opposite results.

Conclusion: The IRRS model based on the identified 4 sIRlncRs showed the significant values on forecasting prognoses of pRCC patients, with the longer OS in the low-risk group.

Methods: We integrated the expression profiles of LncRNA and overall survival (OS) in the 322 pRCC patients based on the TCGA dataset. The immune scores calculated on account of the expression level of immune-related genes were used to verify the most relevant IR-lncRs. Survival-related IR-lncRs (sIRlncRs) were estimated by COX regression analysis in pRCC patients. The high-risk group and low-risk group were identified by the median immune-related risk score (IRRS) model established by the screened sIRlncRs. Functional annotation was displayed by gene set enrichment analysis (GSEA) and principal component analysis (PCA), and the immune composition and purity of the tumor were evaluated through microenvironment cell count records. The expression levels of sIRlncRs of pRCC samples were verified by real-time quantitative PCR.

Cell Cycle Changes, DNA Ploidy, and PTTG1 Gene Expression in HTLV-1 Patients

Human T-cell lymphotropic virus type-1 (HTLV-1) is a pathogenic retrovirus that is associated with adult T-cell leukemia/lymphoma (ATL). Genetic instability is the hallmark of ATL. Cell cycle progression is needed for virus particle reproduction. HTLV-1 encoded Tax protein ultimately disrupts the mitotic spindle checkpoint, leading to incorrect chromosome segregation, resulting in aneuploidy. Cell cycle abnormalities have been described in T cells transfected with HTLV-1 virus in vitro, but not in HTLV-1 asymptomatic carriers. PTTG1 and HTLV-1 viral protein Tax exhibit a cooperative transforming activity. Overexpressed PTTG1 results in chromosome instability and aneuploidy, which has been suggested as a mechanism underlying PTTG1 transforming activity. Here we aimed to investigate cell cycle, DNA ploidy and PTTG1 mRNA expression in CD4⁺ and CD8⁺ T cells in healthy subjects (HS), HTLV-1 asymptomatic carriers and ATL patients. We have identified that HTLV-1 asymptomatic carriers have shown DNA aneuploidy and cell cycle arrest at cell cycle phase G₀/G₁ in CD4⁺ T cells. CD8⁺ T cells of HTLV-1 asymptomatic carriers also demonstrated DNA aneuploidy but without alteration in cell cycle. In ATL, CD4⁺ and CD8⁺ T cells present a higher number of cells in cell cycle S-phase and PTTG1 overexpression. These studies provide insight into malignant transformation of HTLV-1 asymptomatic carriers to ATL patients.

Identification of Crucial Genes Associated With Immune Cell Infiltration in Hepatocellular Carcinoma by Weighted Gene Co-expression Network Analysis

The dreadful prognosis of hepatocellular carcinoma (HCC) is primarily due to the low early diagnosis rate, rapid progression, and high recurrence rate. Valuable prognostic biomarkers are urgently needed for HCC. In this study, microarray data were downloaded from GSE14520, GSE22058, International Cancer Genome Consortium (ICGC), and The Cancer Genome Atlas (TCGA). Differentially expressed genes (DEGs) were identified among GSE14520, GSE22058, and ICGC databases. Weighted gene co-expression network analysis (WGCNA) was used to establish gene co-expression modules of DEGs, and genes of key modules were examined to identify hub genes using univariate Cox regression in the ICGC cohort. Expression levels and time-dependent receiver operating characteristic (ROC) and area under the curve (AUC) were determined to estimate the prognostic competence of the hub genes. These hub genes were also validated in the Gene Expression Profiling Interactive Analysis (GEPIA) and TCGA databases. TIMER algorithm and GSCALite database were applied to analyze the association of the hub genes with immunocytotic infiltration and their pathway enrichment. Altogether, 276 DEGs were identified and WGCNA described a unique and significantly DEGs-associated co-expression module containing 148 genes, with 10 hub genes selected by univariate Cox regression in the ICGC cohort (BIRC5, FOXM1, CENPA, KIF4A, DTYMK, PRC1, IGF2BP3, KIF2C, TRIP13, and TPX2). Most of the genes were validated in the GEPIA databases, except IGF2BP3. The results of multivariate Cox regression analysis indicated that the abovementioned hub genes are all independent predictors of HCC. The 10 genes were also confirmed to be associated with immune cell infiltration using the TIMER algorithm. Moreover, four-gene signature was developed, including BIRC5, CENPA, FOXM1, DTYMK. These hub genes and the model demonstrated a strong prognostic capability and are likely to be a therapeutic target for HCC. Moreover, the association of these genes with immune cell infiltration improves our understanding of the occurrence and development of HCC.

Identification of key genes in non-small cell lung cancer by bioinformatics analysis

Background

Non-small cell lung cancer (NSCLC) is one of the most common malignant tumors in the world, and it has become the leading cause of death of malignant tumors. However, its mechanisms are not fully clear. The aim of this study is to investigate the key genes and explore their potential mechanisms involving in NSCLC.

Methods

We downloaded gene expression profiles GSE33532, GSE30219 and GSE19804 from the Gene Expression Omnibus (GEO) database and analyzed them by using GEO2R. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes were used for the functional and pathway enrichment analysis. We constructed the protein-protein interaction (PPI) network by STRING and visualized it by Cytoscape. Further, we performed module analysis and centrality analysis to find the potential key genes. Finally, we carried on survival analysis of key genes by GEPIA.

Results

In total, we obtained 685 DEGs. Moreover, GO analysis showed that they were mainly enriched in cell adhesion, proteinaceous extracellular region, heparin binding. KEGG pathway analysis revealed that transcriptional misregulation in cancer, ECM-receptor interaction, cell cycle and p53 signaling pathway were involved in. Furthermore, PPI network was constructed including 249 nodes and 1,027 edges. Additionally, a significant module was found, which included eight candidate genes with high centrality features. Further, among the eight candidate genes, the survival of NSCLC patients with the seven high expression genes were significantly worse, including CDK1, CCNB1, CCNA2, BIRC5, CCNB2, KIAA0101 and MELK. In summary, these identified genes should play an important role in NSCLC, which can provide new insight for NSCLC research.

A Gene Signature of Survival Prediction for Kidney Renal Cell Carcinoma by Multi-Omic Data Analysis

Kidney renal cell carcinoma (KIRC), which is the most common subtype of kidney cancer, has a poor prognosis and a high mortality rate. In this study, a multi-omics analysis is performed to build a multi-gene prognosis signature for KIRC. A combination of a DNA methylation analysis and a gene expression data analysis revealed 863 methylated differentially expressed genes (MDEGs). Seven MDEGs (BID, CCNF, DLX4, FAM72D, PYCR1, RUNX1, and TRIP13) were further screened using LASSO Cox regression and integrated into a prognostic risk score model. Then, KIRC patients were divided into high- and low-risk groups. A univariate cox regression analysis revealed a significant association between the high-risk group and a poor prognosis. The time-dependent receiver operating characteristic (ROC) curve shows that the risk group performs well in predicting overall survival. Furthermore, the risk group is contained in the best multivariate model that was obtained by a multivariate stepwise analysis, which further confirms that the risk group can be used as a potential prognostic biomarker. In addition, a nomogram was established for the best multivariate model and shown to perform well in predicting the survival of KIRC patients. In summary, a seven-MDEG signature is a powerful prognosis factor for KIRC patients and may provide useful suggestions for their personalized therapy.

Identification of novel prognostic alternative splicing signature in papillary renal cell carcinoma

Papillary renal cell carcinoma (pRCC) is a heterogeneous disease containing multifocal or solitary tumors with an aggressive phenotype. Increasing evidence has indicated the involvement of aberrant splicing variants in renal cell cancer, while systematic profiling of aberrant alternative splicing (AS) in pRCC was lacking and largely unknown. In the current study, comprehensive profiling of AS events were performed based on the integration of pRCC cohort from the Cancer Genome Atlas database and SpliceSeq software. With rigorous screening and univariate Cox analysis, a total of 2077 prognoses AS events from 1642 parent genes were identified. Then, stepwise least absolute shrinkage and selection operator method-penalized Cox regression analyses with 10-fold cross-validation followed by multivariate Cox regression were used to construct the prognostic AS signatures within each AS type. And a final 21 AS event-based signature was proposed which showed potent prognostic capability in stratifying patients into low- and high-risk subgroups (P < .0001). Furthermore, time-dependent receiver operating characteristics curves confirmed that the final AS signature was effective and robust in predicting overall survival for pRCC patients with the area under the curve above 0.9 from 1 to 5 years. In addition, splicing correlation network was built to uncover the potential regulatory pattern among prognostic splicing factors and candidate AS events. Besides, gene set enrichment analysis revealed the involvement of these candidates AS events in tumor-related pathways including extracellular matrix organization, oxidative phosphorylation, and P53 signaling pathways. Taken together, our results could contribute to elucidating the underlying mechanism of AS in the oncogenesis process and broaden the novel field of prognostic and clinical application of molecule-targeted approaches in pRCC.

Identification of seven-gene signature for prediction of lung squamous cell carcinoma

Background and aim: Lung squamous cell carcinoma (LUSC), is a pathological subtype of lung cancer, accounting for 30% of the lung cancers. A reliable model was constructed, based on the whole gene expression profiles, to predict the prognosis of patients with LUSC. Methods: The RNA-Seq data of LUSC was downloaded from the TCGA database, and differentially expressed genes (p<0.05, |log2fold change| >1) were screened out. By univariate and multivariate Cox regression analysis, we identified seven prognosis-related genes. Then, we established a risk score staging system to predict the prognosis of patients with LUSC. Compared with other clinical parameters, the risk score was an independent prognostic factor and had a better performance in predicting prognosis. Finally, GSEA analysis was carried out to determine the enrichment pathway significantly. The risk score models were established by Cox proportional hazard regression analysis; the ROC curve was applied to test the performance of risk score model. All the statistical analysis was accomplished by R packages. Results: In this study, a model was constructed to predict prognosis, which contains seven genes: CSRNP1, CLEC18B, MIR27A, AC130456.4, DEFA6, ARL14EPL, and ZFP42. Based on the model, the risk score of each patient was calculated with LUSC (hazard ratio [HR]=2.673, 95% CI=1.871-3.525). It was found that the risk score can distinguish high-risk and low-risk groups in prognosis of LUSC patients, independently. Furthermore, the model was validated by ROC curves in the testing dataset and the whole dataset. Lastly, by gene set enrichment analysis (GSEA), we showed the main enrichment pathways were DNA damage stimulus, DNA repair, and DNA replication. It was suggested that the risk score may provide a new and reliable method for prognosis prediction. Conclusion: The results of this study suggested that the risk score based on seven-genes could indicate a promising and independent prognostic biomarker for LUSC patients.

Serine/threonine kinase BUB1 promotes proliferation and radio-resistance in glioblastoma

BUB1 (Budding uninhibited by benzimidazoles 1), a mitotic checkpoint serine/threonine kinase, has been linked in numerous cancers to pro-tumorigenic phenomena including elevation of cellular proliferation, tumor growth, metastatic potential, and poorer patient prognosis. However, the role of BUB1 in glioblastoma remains poorly investigated. In this study, clinical analyses determined significant enrichment of BUB1 in glioblastoma with direct correlation of elevated expression to poorer prognosis in glioma patients. Genetic inhibition of BUB1 in glioblastoma tumor cells via shRNA silencing diminished both proliferative ability and tumorigenicity in vitro and in vivo. Silencing of BUB1 was additionally determined to promote the cytotoxic effect of irradiation on glioblastoma tumor cells, and investigation of the underlying pathways revealed the roles of DNA mismatch repair, spliceosome and c-Myc pathways. Mechanistically, FOXM1 was determined to positively regulate transcription of BUB1 via direct promoter region binding. For validation, pharmacologic inhibition through administration of a BUB1 inhibitor demonstrated attenuated glioblastoma cellular proliferation in vitro as well as delayed tumor growth with prolonged survival in vivo. Collectively, this study demonstrates a novel therapeutic target for glioblastoma in the form of BUB1, which plays a pivotal role in GBM proliferative and radio-resistance capacities in a FOXM1-dependant manner.