Identification of 4-lncRNA prognostic signature in head and neck squamous cell carcinoma

Deregulated long noncoding RNAs (lncRNA) have been critically implicated in tumorigenesis and serve as novel diagnostic and prognostic biomarkers. Here we sought to develop a prognostic lncRNA signature in patients with head and neck squamous cell carcinoma (HNSCC). Original RNA-seq data of 499 HNSCC samples were retrieved from The Cancer Genome Atlas database, which was randomly divided into training and testing set. Univariate Cox regression survival analysis, robust likelihood-based survival model and random sampling iterations were applied to identify prognostic lncRNA candidates in the training cohort. A prognostic risk score was developed based on the Cox coefficient of four individual lncRNA imputed as follows: (0.14546 × expression level of RP11-366H4.1) + (0.27106 × expression level of LINC01123) + (0.54316 × expression level of RP11-110I1.14) + (−0.48794 × expression level of CTD-2506J14.1). Kaplan-Meier analysis revealed that patients with high-risk score had significantly reduced overall survival as compared with those with low-risk score when patients in training, testing, and validation cohorts were stratified into high- or low-risk subgroups. Multivariate survival analysis further revealed that this 4-lncRNA signature was a novel and important prognostic factor independent of multiple clinicopathological parameters. Importantly, ROC analyses indicated that predictive accuracy and sensitivity of this 4-lncRNA signature outperformed those previously well-established prognostic factors. Noticeably, prognostic score based on quantification of these 4-lncRNA via qRT-PCR in another independent HNSCC cohort robustly stratified patients into subgroups with high or low survival. Taken together, we developed a robust 4-lncRNA prognostic signature for HNSCC that might provide a novel powerful prognostic biomarker for precision oncology.     

A robust six-miRNA prognostic signature for head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains a major health problem worldwide. We aimed to identify a robust microRNA (miRNA)-based signature for predicting HNSCC prognosis. The miRNA expression profiles of HNSCC were obtained from The Cancer Genome Atlas (TCGA) database. The TCGA HNSCC cohort was randomly divided into the discovery and validation cohort. A miRNA-based prognostic signature was built up based on TGCA discovery cohort, and then further validated. The downstream targets of prognostic miRNAs were subjected to functional enrichment analyses. The role of miR-1229-3p, a prognosis-related miRNA, in tumorigenesis of HNSCC was further evaluated. A total of 305 significantly differentially expressed miRNAs were found between HNSCC samples and normal tissues. A six-miRNA prognostic signature was constructed, which exhibited a strong association with overall survival (OS) in the TCGA discovery cohort. In addition, these findings were successfully confirmed in TCGA validation cohort and our own independent cohort. The miRNA-based signature was demonstrated as an independent prognostic indicator for HNSCC. A risk signature-based nomogram model was constructed and showed good performance for predicting the OS for HNSCC. The functional analyses revealed that the downstream targets of these prognostic miRNAs were closely linked to cancer progression. Mechanistically, in vitro analysis revealed that miR-1229-3p played a tumor promoting role in HNSCC. In conclusion, our study has developed a robust miRNA-based signature for predicting the prognosis of HNSCC with high accuracy, which will contribute to improve the therapeutic outcome.     

Identification of a six-gene prognostic signature for oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of malignancies worldwide, and its morbidity and mortality have increased in the near term. Consequently, the purpose of the present study was to identify the notable differentially expressed genes (DEGs) involved in their pathogenesis to obtain new biomarkers or potential therapeutic targets for OSCC. The gene expression profiles of the microarray datasets GSE85195, GSE23558, and GSE10121 were obtained from the Gene Expression Omnibus (GEO) database. After screening the DEGs in each GEO dataset, 249 DEGs in OSCC tissues were obtained. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathway enrichment analysis was employed to explore the biological functions and pathways of the above DEGs. A protein–protein interaction network was constructed to obtain a central gene. The corresponding total survival information was analyzed in patients with oral cancer from The Cancer Genome Atlas (TCGA). A total of six candidate genes (CXCL10, OAS2, IFIT1, CCL5, LRRK2, and PLAUR) closely related to the survival rate of patients with oral cancer were identified, and expression verification and overall survival analysis of six genes were performed based on TCGA database. Time-dependent receiver operating characteristic curve analysis yields predictive accuracy of the patient's overall survival. At the same time, the six genes were further verified by quantitative real-time polymerase chain reaction using samples obtained from the patients recruited to the present study. In conclusion, the present study identified the prognostic signature of six genes in OSCC for the first time via comprehensive bioinformatics analysis, which could become potential prognostic markers for OCSS and may provide potential therapeutic targets for tumors.     

Identification of autophagy-related long non-coding RNA prognostic signature for breast cancer

Autophagy-related long non-coding RNAs (lncRNAs) disorders are related to the occurrence and development of breast cancer. The purpose of this study is to explore whether autophagy-related lncRNA can predict the prognosis of breast cancer patients. The autophagy-related lncRNAs prognostic signature was constructed by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression. We identified five autophagy-related lncRNAs (MAPT-AS1, LINC01871, AL122010.1, AC090912.1, AC061992.1) associated with prognostic value, and they were used to construct an autophagy-related lncRNA prognostic signature (ALPS) model. ALPS model offered an independent prognostic value (HR = 1.664, 1.381-2.006), where this risk score of the model was significantly related to the TNM stage, ER, PR and HER2 status in breast cancer patients. Nomogram could be utilized to predict survival for patients with breast cancer. Principal component analysis and Sankey Diagram results indicated that the distribution of five lncRNAs from the ALPS model tends to be low-risk. Gene set enrichment analysis showed that the high-risk group was enriched in autophagy and cancer-related pathways, and the low-risk group was enriched in regulatory immune-related pathways. These results indicated that the ALPS model composed of five autophagy-related lncRNAs could predict the prognosis of breast cancer patients.     

Immune implication of an autophagy-related prognostic signature in uveal melanoma

Uveal Melanoma (UM) is a rare cancer deriving from melanocytes within the uvea. It has a high rate of metastasis, especially to the liver, and a poor prognosis thereafter. Autophagy, an intracellular programmed digestive process, has been associated with the development and progression of cancers, with controversial pro- and anti-tumour roles. Although previous studies have been conducted on autophagy-related genes (ARGs) in various cancer types, its role in UM requires a deeper understanding for improved diagnosis and development of novel therapeutics. In the present study, Zheng et al. used univariate Cox regression followed by least absolute shrinkage and selection operator (Lasso) regression to identify a robust 9-ARG signature prognostic of survival in a total of 230 patients with UM. The authors used the Cancer Genome Atlas (TCGA) UM cohort as a training cohort (n=80) to identify the signature and validated it in another four independent cohorts of 150 UM patients from the Gene Expression Omnibus (GEO) repository ({"type":"entrez-geo","attrs":{"text":"GSE22138","term_id":"22138"}}GSE22138, {"type":"entrez-geo","attrs":{"text":"GSE27831","term_id":"27831"}}GSE27831, {"type":"entrez-geo","attrs":{"text":"GSE44295","term_id":"44295"}}GSE44295 and {"type":"entrez-geo","attrs":{"text":"GSE84976","term_id":"84976"}}GSE84976). This 9-ARG signature was also significantly associated with the enrichment of cancer hallmarks, including angiogenesis, IL6-KJAK-STAT3 signalling, reactive oxygen species pathway and oxidative phosphorylation. More importantly, this signature is associated with immune-related functional pathways and immune cell infiltration. Thus, this 9-ARG signature predicts prognosis and provides deeper insights into the immune mechanisms in UM, with potential implications for future immunotherapy.     

Analysis of methylation-driven genes for predicting the prognosis of patients with head and neck squamous cell carcinoma

To help provide evidence for prognosis prediction and personalized targeted therapy for patients with head and neck squamous cell carcinoma (HNSCC), we investigated prognosis-specific methylation-driven genes in HNSCC. Survival time data, RNA sequencing data, and methylation data for HNSCC patients were downloaded from The Cancer Genome Atlas. The MethylMix R package based on the β mixture model was utilized to screen genes with different methylation statuses in tumor tissues and adjacent normal tissues, and a total of 182 HNSCC-related methylation-driven genes were then identified. A survival prediction scoring model based on multivariate Cox analysis was developed to screen the genes related to the prognosis of HNSCC, and a linear risk model of the methylation status of six genes (INA, LINC01354, TSPYL4, MAGEB2, EPHX3, and ZNF134) was constructed. The prognostic values of the six genes were further independently explored by survival analysis combined with methylation and gene expression analyses. The 5-year survival rate in the high-risk group of patients in the test set was 30.4% (95% CI: 22.7%-40.8%) and that in the low-risk group of patients was 65.5% (95% CI: 56.1%-76.5%). The area under the receiver operating characteristic curve for the model was 0.723, which further verified the specificity and sensitivity of the model. In addition, subsequent combined survival analysis revealed that all six genes could be used as independent prognostic markers and thus might be potential drug targets. The innovative method provides new insight into the molecular mechanism and prognosis of HNSCC.     

多发性骨髓瘤预后相关自噬基因筛选及建模

为探讨自噬相关基因(ARGs)在MM发生发展中的作用机制并建立相关的预后模型。基于MMRF与HADb数据库,通过R语言确定多发性骨髓瘤中自噬相关基因的差异表达,GO和KEGG分析自噬相关基因与多发性骨髓瘤发生发展的关系,使用COX回归算法建立多基因预后模型,Kaplan-Meier方法绘制生存曲线,ROC曲线评价预后模型的可靠性。最终从764例多发性骨髓瘤患者骨髓样本及4例正常骨髓样本中共发现104个基因的表达在多发性骨髓瘤样本中具有显著差异,其中上调基因46个,下调基因58个。GO富集主要集中在巨自噬、自噬调节、细胞对外部刺激的反应等本体学注释。KEGG富集主要集中在自噬、细胞凋亡、NOD样受体信号通路、PI3K-Akt信号通路。单因素COX分析发现33个自噬相关基因与多发性骨髓瘤患者整体生存明显相关。多因素COX回归筛选出13个预后相关自噬相关基因(NKX2-3、NCKAP1、BIRC5、PEX3、HGS、RUBCN、PARP1、ARSA、DNAJB9、HSP90AB1、EEF2、FKBP1B和CD46)建立多发性骨髓瘤自噬相关基因预后模型。Kaplan-Meier生存曲线分析显示...     

Identification of an autophagy-related gene signature predicting overall survival for hepatocellular carcinoma

The poor prognosis of hepatocellular carcinoma (HCC) calls for the development of accurate prognostic models. The growing number of studies indicating a correlation between autophagy activity and HCC indicates there is a commitment to finding solutions for the prognosis of HCC from the perspective of autophagy. We used a cohort in The Cancer Genome Atlas (TCGA) to evaluate the expression of autophagy-related genes in 371 HCC samples using univariate Cox and lasso Cox regression analysis, and the prognostic features were identified. A prognostic model was established by combining the expression of selected genes with the multivariate Cox regression coefficient of each gene. Eight autophagy-related genes were selected as prognostic features of HCC. We established the HCC prognostic risk model in TCGA dataset using these identified prognostic genes. The model’s stability was confirmed in two independent verification sets ({"type":"entrez-geo","attrs":{"text":"GSE14520","term_id":"14520"}}GSE14520 and {"type":"entrez-geo","attrs":{"text":"GSE36376","term_id":"36376"}}GSE36376). The model had a good predictive power for the overall survival (OS) of HCC (hazard ratio = 2.32, 95% confidence interval = 1.76–3.05, P<0.001). Moreover, the risk score computed by the model did not depend on other clinical parameters. Finally, the applicability of the model was demonstrated through a nomogram (C-index = 0.701). In the present study, we established an autophagy-related risk model having a high prediction accuracy for OS in HCC. Our findings will contribute to the definition of prognosis and establishment of personalized therapy for HCC patients.     

FOXD1 is a prognostic biomarker and correlated with macrophages infiltration in head and neck squamous cell carcinoma

Background: Forkhead Box D1 (FOXD1) is differentially expressed in various tumors. However, its role and correlation with immune cell infiltration remains uncertain in head and neck squamous cell carcinoma (HNSC).Methods: FOXD1 expression was analyzed in The Cancer Genome Atlas (TCGA) pan-cancer data. The clinical prognosis influence of FOXD1 was evaluated by clinical survival data of TCGA. Enrichment analysis of FOXD1 was performed using R packages ‘clusterProfiler’. We downloaded the immune cell infiltration score of TCGA samples from published articles, and analyzed the correlation between immune cell infiltration level and FOXD1 expression.Results: FOXD1 was highly expressed and associated with poorer overall survival (OS, P<0.0001), disease-specific survival (DSS, P=0.00011), and progression-free interval (PFI, P<0.0001) in HNSC and some other tumors. In addition, FOXD1 expression was significantly correlated with infiltration of immune cells. Tumor-associated macrophages (TAMs) infiltration increased in tissues with high FOXD1 expression in HNSC. Immunosuppressive genes such as PD-L1, IL-10, TGFB1, and TGFBR1 were significantly positively correlated with FOXD1.Conclusions: Our study suggests FOXD1 to be an oncogene and act as an indicator of poor prognosis in HNSC. FOXD1 might contribute to the TAM infiltration in HNSC. High FOXD1 may be associated with tumor immunosuppression status.     

Identification of CCT3 as a prognostic factor and correlates with cell survival and invasion of head and neck squamous cell carcinoma

Background: Recurrent locally advanced or metastatic head and neck squamous cell carcinoma (HNSCC) is associated with dismal prognosis because of its highly invasive behavior and resistance to conventional intensive chemotherapy. The identification of effective markers for early diagnosis and prognosis is important for reducing mortality and ensuring that therapy for HNSCC is effective. Chaperonin-containing TCP-1 3 (CCT3) folds cancer-related proteins to control carcinogenesis. The prognostic value and growth association of CCT3 and HNSCC remain unknown.Methods: The GEO, Oncomine and UALCAN databases were used to examine CCT3 expression in HNSCC. A few clinical HNSCC samples with normal tissues were used to detect CCT3 expression by using immunohistochemistry method. The TCGA-HNSC dataset was used to evaluate the association between expression of CCT3 and prognosis. The molecular mechanism was investigated with gene set enrichment analysis (GSEA). CCK-8 and wound healing assays were used to detect cell growth and invasion of HNSCC, respectively.Results: CCT3 expression was significantly up-regulated in HNSCC in both mRNA and protein levels. In addition, up-regulated CCT3 expression was associated with various clinicopathological parameters. High expression of CCT3 was significantly correlated with inferior survival of HNSCC patients. Knockdown of CCT3 significantly inhibited cell growth and invasion of HNSCC cell lines. GSEA analysis indicated that CCT3 was closely correlated with tumor-related signaling pathways and HNSCC cell survival.Conclusion: Our findings suggest that CCT3 is a biomarker of poor prognosis and related to the process of HNSCC.     

Systematic profiling of invasion-related gene signature predicts prognostic features of lung adenocarcinoma

Due to the high heterogeneity of lung adenocarcinoma (LUAD), molecular subtype based on gene expression profiles is of great significance for diagnosis and prognosis prediction in patients with LUAD. Invasion-related genes were obtained from the CancerSEA database, and LUAD expression profiles were downloaded from The Cancer Genome Atlas. The ConsensusClusterPlus was used to obtain molecular subtypes based on invasion-related genes. The limma software package was used to identify differentially expressed genes (DEGs). A multi-gene risk model was constructed by Lasso-Cox analysis. A nomogram was also constructed based on risk scores and meaningful clinical features. 3 subtypes (C1, C2 and C3) based on the expression of 97 invasion-related genes were obtained. C3 had the worst prognosis. A total of 669 DEGs were identified among the subtypes. Pathway enrichment analysis results showed that the DEGs were mainly enriched in the cell cycle, DNA replication, the p53 signalling pathway and other tumour-related pathways. A 5-gene signature (KRT6A, MELTF, IRX5, MS4A1 and CRTAC1) was identified by using Lasso-Cox analysis. The training, validation and external independent cohorts proved that the model was robust and had better prediction ability than other lung cancer models. The gene expression results showed that the expression levels of MS4A1 and KRT6A in tumour tissues were higher than in normal tissues, while CRTAC1 expression in tumour tissues was lower than in normal tissues. The 5-gene signature prognostic stratification system based on invasion-related genes could be used to assess prognostic risk in patients with LUAD.     

A two-CpG-based prognostic signature for oral squamous cell carcinoma overall survival

Oral squamous cell carcinoma (OSCC) represents one of the most common head and neck cancer that with dire prognosis due partly to the lack of reliable prognostic biomarker. Here, we aimed to develop a CpG site–based prognostic signature through which we could accurately predict overall survival (OS) of patients with OSCC. We obtained OSCC-related DNA methylation and gene expression data sets from the public accessible Gene Expression Omnibus. Correlations between methylation level of CpG sites and OS of patients with OSCC were assessed by univariate Cox regression analysis followed by robust likelihood-based survival analysis on those CpG sites with permutation P < 0.05 for further screening the optimal CpG sites for OSCC OS prediction based on the risk score formula that composed of the methylation level of optimal CpG sites weighted by their regression coefficients. Besides, differential expression genes (DEGs) and differential methylation genes (DMGs) in OSCC samples compared with normal samples were obtained and shared genes were considered as vital genes in OSCC tumorgenesis and progression. As a result, two CpG sites including cg17892178 and cg17378966 that located in NID2 and IDO1, respectively, were identified as the optimal prognostic signatures for OSCC OS. In addition, 12 overlapping genes between DEGs and DMGs that closely associated with inflammation or blood and tissue development–related biological processes were obtained. In conclusions, this study should provide valuable signatures for OSCC diagnosis and treatment.     

Plectin promotes migration and invasion of cancer cells and is a novel prognostic marker for head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is usually found at a late stage and distant metastasis occurs at high frequency; therefore, novel prognostic markers are needed. This study was aimed to identify novel tumor markers in HNSCC. We identified 65 proteins which were significantly increased or decreased in the tumors by 2D-DIGE using 12 HNSCC and adjacent non-cancer tissues. Western blotting and immunohistochemical analysis confirmed that the expression of plectin was significantly increased in most cancer tissues as compared with non-cancer tissues. Strikingly, the suppression of endogenous plectin using siRNA inhibited the proliferation, migration and invasion of HNSCC cells and down-regulated Erk 1/2 kinase. Furthermore, immunohistochemistry using paraffin-embedded tissues from 62 patients showed not only that the frequency of recurrence was correlated with the plectin expression but that the prognosis of patients with a high plectin was extremely poor. Moreover, the survival rate of patients with a high plectin was significantly lower than that of patients with low E-cadherin levels, which is known to correlate with the poor prognosis of HNSCC. Our findings suggest that plectin promotes the migration and invasion of HNSCC cells through activation of Erk 1/2 kinase and is a potential prognostic biomarker of HNSCC.     

Identification of potential core genes and pathways predicting pathogenesis in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the most common subtype of head and neck cancer; however, its pathogenesis and potential therapeutic targets remain largely unknown. In the present study, we analyzed three gene expression profiles and screened differentially expressed genes (DEGs) between HNSCC and normal tissues. The DEGs were subjected to gene ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), protein–protein interaction (PPI), and survival analyses, while the connectivity map (CMap) database was used to predict candidate small molecules that may reverse the biological state of HNSCC. Finally, we measured the expression of the most relevant core gene in vitro and examined the effect of the top predicted potential drug against the proliferation of HNSCC cell lines. Among the 208 DEGs and ten hub genes identified, CDK1 and CDC45 were associated with unfavorable HNSCC prognosis, and three potential small molecule drugs for treating HNSCC were identified. Increased CDK1 expression was confirmed in HNSCC cells, and menadione, the top predicted potential drug, exerted significant inhibitory effects against HNSCC cell proliferation and markedly reversed CDK1 expression. Together, the findings of the present study suggest that the ten hub genes and pathways identified may be closely related to HNSCC pathogenesis. In particular, CDK1 and CDC45 overexpression could be reliable biomarkers for predicting unfavorable prognosis in patients with HNSCC, while the new candidate small molecules identified by CMap analysis provide new avenues for the development of potential drugs to treat HNSCC.     

Identification of a prognostic alternative splicing signature in oral squamous cell carcinoma

Alternative splicing (AS) is critically associated with tumorigenesis and patient's prognosis. Here, we systematically analyzed survival-associated AS signatures in oral squamous cell carcinoma (OSCC) and evaluated their prognostic predictive values. Survival-related AS events were identified by univariate and multivariate Cox regression analyses using OSCC data from the TCGA head neck squamous cell carcinoma data set. The Percent Spliced In calculated by SpliceSeq from 0 to 1 was used to quantify seven types of AS events. A predictive model based on AS events was constructed by least absolute shrinkage and selection operator Cox regression assay and further validated using a training-testing cohort design. Patient survival was estimated using the Kaplan–Meier method and compared with Log-rank test. The receiver operating characteristics curve area under the curves was used to evaluate the predictive abilities of these predictive models. Furthermore, gene–gene interaction networks and the splicing factors (SFs)-AS regulatory network was generated by Cytoscape. A total of 825 survival-related AS events within 719 genes were identified in OSCC samples. The integrative predictive model was better at predicting outcomes of patients as compared to those models built with the individual AS event. The predictive model based on three AS-related genes also effectively predicted patients’ survival. Moreover, seven survival-related SFs were detected in OSCC including RBM4, HNRNPD, and HNRNPC, which have been linked to tumorigenesis. The SF-AS network revealed a significant correlation between survival-related AS genes and these SFs. Our findings revealed a systemic portrait of survival-associated AS events and the splicing network in OSCC, suggesting that AS events might serve as novel prognostic biomarkers and therapeutic targets for OSCC.     

Construction and validation of a seven-gene signature for predicting overall survival in patients with kidney renal clear cell carcinoma via an integrated bioinformatics analysis

ABSTRACT

Kidney renal clear cell carcinoma (KIRC) remains a significant challenge worldwide because of its poor prognosis and high mortality rate, and accurate prognostic gene signatures are urgently required for individual therapy. This study aimed to construct and validate a seven-gene signature for predicting overall survival (OS) in patients with KIRC. The mRNA expression profile and clinical data of patients with KIRC were obtained from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). Prognosis-associated genes were identified, and a prognostic gene signature was constructed. Then, the prognostic efficiency of the gene signature was assessed. The results obtained using data from the TCGA were validated using those from the ICGC and other online databases. Gene set enrichment analyses (GSEA) were performed to explore potential molecular mechanisms. A seven-gene signature (PODXL, SLC16A12, ZIC2, ATP2B3, KRT75, C20orf141, and CHGA) was constructed, and it was found to be effective in classifying KIRC patients into high- and low-risk groups, with significantly different survival based on the TCGA and ICGC validation data set. Cox regression analysis revealed that the seven-gene signature had an independent prognostic value. Then, we established a nomogram, including the seven-gene signature, which had a significant clinical net benefit. Interestingly, the seven-gene signature had a good performance in distinguishing KIRC from normal tissues. GSEA revealed that several oncological signatures and GO terms were enriched. This study developed a novel seven-gene signature and nomogram for predicting the OS of patients with KIRC, which may be helpful for clinicians in establishing individualized treatments.