CFHR3 is a potential novel biomarker for hepatocellular carcinoma

Complement factor H-related 3 (CFHR3) is a protein-coding gene acting in various diseases. However, its prognostic values of CFHR3 in hepatocellular carcinoma (HCC) are not understandable. Therefore, we present a further study on CFHR3 in HCC. CFHR3 expression data were acquired from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC). We compared the differential expression of CFHR3 between the low-stage (stage I and II) and high-stage (stage III and IV) patients with HCC in the TCGA and ICGC cohorts. Furthermore, we assessed the CFHR3 expression as a prognostic marker using the Kaplan-Meier survival analysis, univariate, and multivariate analysis. The Kaplan-Meier analysis declared that CFHR3 overexpression was correlated with a good prognosis for HCC patients. Multivariate analysis proved the prognostic significance of CFHR3 expression levels (P < .001 and .003 for TCGA and ICGC, respectively). Immune-related scores in low-risk cohorts were higher than high-risk cohorts. Gene set enrichment analysis implied that the low CFHR3 expression phenotype was significantly enriched in critical biological functions and pathways and was associated with tumorigenesis, such as regulation of cell activation cycle, and the WNT and NOTCH signal pathway. Above all, CFHR3 could be a novel prognostic biomarker and therapeutic target for HCC.     

Development and validation of a CTNNB1-associated metabolic prognostic model for hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a heterogeneous malignancy closely related to metabolic reprogramming. We investigated how CTNNB1 mutation regulates the HCC metabolic phenotype and thus affects the prognosis of HCC. We obtained the mRNA expression profiles and clinicopathological data from The Cancer Genome Atlas (TCGA), the International Cancer Genomics Consortium (ICGC) and the Gene Expression Omnibus database ( GSE14520 and GSE116174 ). We conducted gene set enrichment analysis on HCC patients with and without mutant CTNNB1 through TCGA dataset. The Kaplan-Meier analysis and univariate Cox regression analysis assisted in screening metabolic genes related to prognosis, and the prognosis model was constructed using the Lasso and multivariate Cox regression analysis. The prognostic model showed good prediction performance in both the training cohort (TCGA) and the validation cohorts (ICGC, GSE14520 , GSE116174 ), and the high-risk group presented obviously poorer overall survival compared with low-risk group. Cox regression analysis indicated that the risk score can be used as an independent predictor for the overall survival of HCC. The immune infiltration in different risk groups was also evaluated in this study to explore underlying mechanisms. This study is also the first to describe an metabolic prognostic model associated with CTNNB1 mutations and could be implemented for determining the prognoses of individual patients in clinical practice.     

Identification and Validation of a Prognostic Model Based on Three MVI-Related Genes in Hepatocellular Carcinoma

MVI has significant clinical value for treatment selection and prognosis evaluation in hepatocellular carcinoma (HCC). We aimed to construct a model based on MVI-Related Genes (MVIRGs) for risk assessment and prognosis prediction in patients with HCC. This study utilized various statistical analysis methods for prognostic model construction and validation in the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) cohorts, respectively. In addition, immunohistochemistry and qRT-PCR were used to analyze and identify the value of the model in our cohort. After the analyses, 153 differentially expressed MVIRGs were identified, and three key genes were selected to construct a prognostic model. The high-risk group showed significantly lower overall survival (OS), and this trend was observed in all subgroups: different age groups, genders, stages, and grades. Risk score was a risk factor independent of age, gender, stage, and grade. Moreover, the ICGC cohort validated the prognostic value of the model corresponding to the TCGA. In our cohort, qRT-PCR and immunohistochemistry showed that all three genes had higher expression levels in HCC samples than in normal controls. High expression levels of genes and high-risk scores showed significantly lower recurrence-free survival (RFS) and OS, especially in MVI-positive HCC samples. Therefore, the prognostic model constructed by three MVIRGs can reliably predict the RFS and OS of patients with HCC and is valuable for guiding clinical treatment selection and prognostic assessment of HCC.     

VNN3 is a potential novel biomarker for predicting prognosis in clear cell renal cell carcinoma

Although pathological observations provide approximate prognoses, it is difficult to achieve prognosis in patients with existing prognostic factors. Therefore, it is very important to find appropriate biomarkers to achieve accurate cancer prognosis. Renal cell carcinoma (RCC) has several subtypes, the discrimination of which is crucial for proper treatment. Here, we present a novel biomarker, VNN3, which is used to prognose clear cell renal cell carcinoma (ccRCC), the most common and aggressive subtype of kidney cancer. Patient information analyzed in our study was extracted from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) cohorts. VNN3 expression was considerably higher in stages III and IV than in stages I and II. Moreover, Kaplan–Meier curves associated high VNN3 expression with poor prognoses (TCGA, p?p?=?.00076), confirming that ccRCC prognosis can be predicted via VNN3 expression patterns. Consistent with all patient results, the prognosis of patients with higher VNN3 expression was worse in both low stage (I and II) and high stage (III and IV) (TCGA, p < 0.0001 in stage I and II; ICGC, p = 0.028 in stage I and II; TCGA, p = 0.005 in stage III and IV). Area under the curve and receiver operating characteristic curves supported our results that highlighted VNN3 expression as a suitable ccRCC biomarker. Multivariate analysis also verified the prognostic performance of VNN3 expression (TCGA, p?p?=?.017). Altogether, we suggest that VNN3 is applicable as a new biomarker to establish prognosis in patients with ccRCC.     

Integrated Machine Learning and Bioinformatic Analyses Constructed a Novel Stemness-Related Classifier to Predict Prognosis and Immunotherapy Responses for Hepatocellular Carcinoma Patients

Immunotherapy has made great progress in hepatocellular carcinoma (HCC), yet there is still a lack of biomarkers for predicting response to it. Cancer stem cells (CSCs) are the primary cause of the tumorigenesis, metastasis, and multi-drug resistance of HCC. This study aimed to propose a novel CSCs-related cluster of HCC to predict patients'' response to immunotherapy. Based on RNA-seq datasets from The Cancer Genome Atlas (TCGA) and Progenitor Cell Biology Consortium (PCBC), one-class logistic regression (OCLR) algorithm was applied to compute the stemness index (mRNAsi) of HCC patients. Unsupervised consensus clustering was performed to categorize HCC patients into two stemness subtypes which further proved to be a predictor of tumor immune microenvironment (TIME) status, immunogenomic expressions and sensitivity to neoadjuvant therapies. Finally, four machine learning algorithms (LASSO, RF, SVM-RFE and XGboost) were applied to distinguish different stemness subtypes. Thus, a five-hub-gene based classifier was constructed in TCGA and ICGC HCC datasets to predict patients'' stemness subtype in a more convenient and applicable way, and this novel stemness-based classification system could facilitate the prognostic prediction and guide clinical strategies of immunotherapy and targeted therapy in HCC.     

Identification of core genes and outcomes in hepatocellular carcinoma by bioinformatics analysis

Hepatocellular carcinoma (HCC) is the most common malignant liver disease in the world. However, the mechanistic relationships among various genes and signaling pathways are still largely unclear. In this study, we aimed to elucidate potential core candidate genes and pathways in HCC. The expression profiles GSE14520, GSE25097, GSE29721, and GSE62232, which cover 606 tumor and 550 nontumour samples, were downloaded from the Gene Expression Omnibus (GEO) database. Furthermore, HCC RNA-seq datasets were also downloaded from the Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were filtered using R software, and we performed gene ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis using the online databases DAVID 6.8 and KOBAS 3.0. Furthermore, the protein-protein interaction (PPI) network complex of these DEGs was constructed by Cytoscape software, the molecular complex detection (MCODE) plug-in and the online database STRING. First, a total of 173 DEGs (41 upregulated and 132 downregulated) were identified that were aberrantly expressed in both the GEO and TCGA datasets. Second, GO analysis revealed that most of the DEGs were significantly enriched in extracellular exosomes, cytosol, extracellular region, and extracellular space. Signaling pathway analysis indicated that the DEGs had common pathways in metabolism-related pathways, cell cycle, and biological oxidations. Third, 146 nodes were identified from the DEG PPI network complex, and two important modules with a high degree were detected using the MCODE plug-in. In addition, 10 core genes were identified, TOP2A, NDC80, FOXM1, HMMR, KNTC1, PTTG1, FEN1, RFC4, SMC4, and PRC1. Finally, Kaplan-Meier analysis of overall survival and correlation analysis were applied to these genes. The abovementioned findings indicate that the identified core genes and pathways in this bioinformatics analysis could significantly enrich our understanding of the development and recurrence of HCC; furthermore, these candidate genes and pathways could be therapeutic targets for HCC treatment.     

Genome-wide screening and cohorts validation identifying novel lncRNAs as prognostic biomarkers for clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the main subtype of renal cell carcinoma with varied prognosis. We aimed to identify and assess the possible prognostic long noncoding RNA (lncRNA) biomarkers. LncRNAs expression data and corresponding clinical information of 619 ccRCC patients were downloaded from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. Differentially expressed genes analysis, univariate Cox regression, the least absolute shrinkage and selection operator Cox regression model were utilized to identify hub lncRNAs. Multivariate Cox regression was used to establish the risk model. Statistical analysis was performed using R 3.5.3. The expression value of five lncRNAs and the risk-score levels were significantly associated with a survival prognosis of ccRCC patients (all P < .001). In the TCGA validation cohort, the area under the curve (AUC) for the integrated nomogram was 0.905 and 0.91 for 3-, 5-year prediction separately. The AUC reached up to 0.757 in an independent ICGC cohort. Besides, the calibration plots also illustrated well curve-fitting between observation values and predictive values. Weighted gene co-expression network analysis and subsequent pathway analysis revealed that the PI3K-Akt-mTOR and hypoxia-inducible factor signaling crosstalk might function as the most essential mechanisms related to the five-lncRNAs signature. Our study suggested that lncRNA AC009654.1, AC092490.2, LINC00524, LINC01234, and LINC01885 were significantly associated with ccRCC prognosis. The prognostic model based on this five lncRNA may predict the overall survival of ccRCC.     

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.     

A P53-related microRNA model for predicting the prognosis of hepatocellular carcinoma patients

Studies have shown that microRNAs (miRNAs) play a vital role in tumor progression and patients’ prognosis. Therefore, we aimed to construct a miRNA model for forecasting the survival of hepatocellular carcinoma (HCC) patients. The gene expression data of 433 patients with HCC from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus public databases were remined by survival analysis and receptor manipulation characteristic curve (ROC). A prognostic model including six miRNAs (hsa-mir-26a-1-3p, hsa-mir-188-5p, hsa-mir-212-5p, hsa-mir-149-5p, hsa-mir-105-5p, and hsa-mir-132-5p) were constructed in the training dataset (TCGA, n = 333). HCC patients were stratified into a high-risk group and a low-risk group with significantly different survival (median: 2.75 vs. 8.93 years, log-rank test p < .001). Then we proved its performance of stratification in another independent dataset (GSE116182, median: 2.55 vs 6.96 years, log-rank test p = .008). Cox regression analysis showed that the prognostic model was an independent prognostic indicator for HCC patients. Then time-dependent ROC analyses were performed to test the prognostic ability of the model with that of TNM staging, we found the model had a better performance, especially at 5 years (AUC = 0.76). Functional prediction showed that the genes targeted by the six prognostic miRNAs in the prognostic model were highly expressed in the P53-related pathway. In conclusion, we constructed a prognostic miRNA model that could indicate the survival of HCC patients.     

Hypermethylated and downregulated MEIS2 are involved in stemness properties and oxaliplatin-based chemotherapy resistance of colorectal cancer

The resistance against oxaliplatin (L-OHP) based regimens remains a major obstacle for its efficient usage in treating metastatic colorectal cancer (mCRC). In this study, we performed weighted gene coexpression network analysis (WGCNA) to systematically screen the relevant hub genes for L-OHP resistance using the raw microarray data of 30 consecutive mCRC samples from our earlier study (GSE69657). The results were further confirmed through datasets from Gene Expression Omnibus (GEO). From L-OHP resistance module, nine genes in both the coexpression and protein–protein interaction networks were chosen as hub genes. Among these genes, Meis Homeobox 2 (MEIS2) had the highest correlation with L-OHP resistance (r = −0.443) and was deregulated in L-OHP resistant tissues compared with L-OHP sensitive tissues in both our own dataset and GSE104645 testing dataset. The receiver operating characteristic curve validated that MEIS2 had a good ability in predicting L-OHP response in both our own dataset (area under the curve [AUC] = 0.802) and GSE104645 dataset (AUC = 0.746). Then, the down expression of MEIS2 was observed in CRC tissue compared with normal tissue in 12 GEO-sourced datasets and The Cancer Genome Atlas (TCGA) and was correlated with poor event-free survival. Furthermore, analyzing methylation data from TCGA showed that MEIS2 had increased promoter hypermethylation. In addition, MEIS2 expression was significantly decreased in CRC stem cells compared with nonstem cells in two GEO datasets (GSE14773 and GSE24747). Further methylation analysis from GSE104271 demonstrated that CRC stem cells had higher MEIS2 promoter methylation levels in cg00366722 and cg00610348 sites. Gene set enrichment analysis showed that MEIS2 might be involved in the Wnt/β-catenin pathway. In the overall view, MEIS2 had increased promoter hypermethylation and was downregulated in poor L-OHP response mCRC tissues. MEIS2 might be involved in the Wnt/β-catenin pathway to maintain CRC stemness, which leads to L-OHP resistance.     

A panel of eight autophagy-related long non-coding RNAs is a good predictive parameter for clear cell renal cell carcinoma

Clear-cell renal cell carcinoma (ccRCC) carries a variable prognosis. Prognostic biomarkers can stratify patients according to risk, and can provide crucial information for clinical decision-making. We screened for an autophagy-related long non-coding lncRNA (lncRNA) signature to improve postoperative risk stratification in The Cancer Genome Atlas (TCGA) database. We confirmed this model in ICGC and SYSU cohorts as a significant and independent prognostic signature. Western blotting, autophagic-flux assay and transmission electron microscopy were used to verify that regulation of expression of 8 lncRNAs related to autophagy affected changes in autophagic flow in vitro. Our data suggest that 8-lncRNA signature related to autophagy is a promising prognostic tool in predicting the survival of patients with ccRCC. Combination of this signature with clinical and pathologic parameters could aid accurate risk assessment to guide clinical management, and this 8-lncRNAs signature related to autophagy may serve as a therapeutic target.     

The relationship among Girdin DNA methylation,its high expression,and immune infiltration in hepatocellular carcinoma: Clues from in silico analysis

Objective: The aim of the present study was to explore the relationship among Girdin DNA methylation, its high expression, and immune infiltration in human hepatocellular carcinoma (HCC).Materials and methods: The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and International Cancer Genome Consortium (ICGC) databases were used to compare Girdin mRNA expression between HCC tissues and normal tissues, and determine the relationship between Girdin expression and HCC prognosis. TCGA database was also used to analyze the expression of Girdin and its methylation status, as well as the relationship between Girdin DNA methylation and HCC prognosis. The Tumor IMmune Estimation Resource (TIMER) database was used to explore the correlation between Girdin expression and HCC immune infiltration.Results: Girdin expression was elevated in HCC tissues compared with that in normal tissues. The degree of methylation at cg03188526, a CpG site in the Girdin gene body, was positively correlated with Girdin mRNA expression, while high Girdin expression and cg03188526 hypermethylation were both correlated with poor HCC prognosis. Additionally, HCC tissue with high Girdin expression exhibited abundant immune infiltration, and the high Girdin expression was associated with a worse prognosis in macrophage-enriched HCC specimens.Conclusion: Our findings indicated that Girdin likely functions as an oncogene in HCC and that hypermethylation at cg03188526 in the Girdin gene body may explain the high Girdin expression levels in HCC tissue. Furthermore, we report for the first time that the adverse effects of high Girdin expression in HCC patients may be partially mediated by tumor macrophage infiltration.     

DNA methylation biomarkers for head and neck squamous cell carcinoma

DNA methylation plays an important role in the etiology and pathogenesis of head and neck squamous cell carcinoma (HNSCC). The current study aimed to identify aberrantly methylated-differentially expressed genes (DEGs) by a comprehensive bioinformatics analysis. In addition, we screened for DEGs affected by DNA methylation modification and further investigated their prognostic values for HNSCC. We included microarray data of DNA methylation (GSE25093 and GSE33202) and gene expression (GSE23036 and GSE58911) from Gene Expression Omnibus. Aberrantly methylated-DEGs were analyzed with R software. The Cancer Genome Atlas (TCGA) RNA sequencing and DNA methylation (Illumina HumanMethylation450) databases were utilized for validation. In total, 27 aberrantly methylated genes accompanied by altered expression were identified. After confirmation by The Cancer Genome Atlas (TCGA) database, 2 hypermethylated-low-expression genes (FAM135B and ZNF610) and 2 hypomethylated-high-expression genes (HOXA9 and DCC) were identified. A receiver operating characteristic (ROC) curve confirmed the diagnostic value of these four methylated genes for HNSCC. Multivariate Cox proportional hazards analysis showed that FAM135B methylation was a favorable independent prognostic biomarker for overall survival of HNSCC patients.     

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.     

Identification of an extracellular vesicle-related gene signature in the prediction of pancreatic cancer clinical prognosis

Although extracellular vesicles (EVs) in body fluid have been considered to be ideal biomarkers for cancer diagnosis and prognosis, it is still difficult to distinguish EVs derived from tumor tissue and normal tissue. Therefore, the prognostic value of tumor-specific EVs was evaluated through related molecules in pancreatic tumor tissue. NA sequencing data of pancreatic adenocarcinoma (PAAD) were acquired from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). EV-related genes in pancreatic cancer were obtained from exoRBase. Protein–protein interaction (PPI) network analysis was used to identify modules related to clinical stage. CIBERSORT was used to assess the abundance of immune and non-immune cells in the tumor microenvironment. A total of 12 PPI modules were identified, and the 3-PPI-MOD was identified based on the randomForest package. The genes of this model are involved in DNA damage and repair and cell membrane-related pathways. The independent external verification cohorts showed that the 3-PPI-MOD can significantly classify patient prognosis. Moreover, compared with the model constructed by pure gene expression, the 3-PPI-MOD showed better prognostic value. The expression of genes in the 3-PPI-MOD had a significant positive correlation with immune cells. Genes related to the hypoxia pathway were significantly enriched in the high-risk tumors predicted by the 3-PPI-MOD. External databases were used to verify the gene expression in the 3-PPI-MOD. The 3-PPI-MOD had satisfactory predictive performance and could be used as a prognostic predictive biomarker for pancreatic cancer.     

A novel autophagy-related lncRNA survival model for lung adenocarcinoma

Long non-coding RNA (lncRNA) is an important regulatory factor in the development of lung adenocarcinoma, which is related to the control of autophagy. LncRNA can also be used as a biomarker of prognosis in patients with lung adenocarcinoma. Therefore, it is important to determine the prognostic value of autophagy-related lncRNA in lung adenocarcinoma. In this study, autophagy-related mRNAs-lncRNAs were screened from lung adenocarcinoma and a co-expression network of autophagy-related mRNAs-lncRNAs was constructed by using The Cancer Genome Atlas (TCGA). The univariate and multivariate Cox proportional hazard analyses were used to evaluate the prognostic value of the autophagy-related lncRNAs and finally obtained a survival model composed of 11 autophagy-related lncRNAs. Through Kaplan-Meier analysis, univariate and multivariate Cox regression analysis and time-dependent receiver operating characteristic (ROC) curve analysis, it was further verified that the survival model was a new independent prognostic factor for patients with lung adenocarcinoma. In addition, based on the survival model, gene set enrichment analysis (GSEA) was used to illustrate the function of genes in low-risk and high-risk groups. These 11 lncRNAs were GAS6-AS1, AC106047.1, AC010980.2, AL034397.3, NKILA, AL606489.1, HLA-DQB1-AS1, LINC01116, LINC01806, FAM83A-AS1 and AC090559.1. The hazard ratio (HR) of the risk score was 1.256 (1.196-1.320) (P < .001) in univariate Cox regression analysis and 1.215 (1.149-1.286) (P < .001) in multivariate Cox regression analysis. And the AUC value of the risk score was 0.809. The 11 autophagy-related lncRNA survival models had important predictive value for the prognosis of lung adenocarcinoma and may become clinical autophagy-related therapeutic targets.     

环氧化物水解酶2在肝细胞癌中的表达水平与功能作用

本研究通过公共数据和实验数据,全面分析环氧化物水解酶2(epoxide hydrolase 2, EPHX2)在肝细胞癌中的表达情况、功能作用以及预后意义。利用GEO和MitoCarta数据集,筛选肝细胞癌中呈差异表达的线粒体相关基因;利用TCGA数据库分析EPHX2及其相关基因在肝细胞癌中的表达水平;运行R包绘制Kaplan-Meier生存曲线和功能富集分析;基于STRING和GSEA构建蛋白质互作网络和基因集富集分析;荧光定量PCR和GEO数据集验证EPHX2在肝细胞癌中的表达水平。本研究共筛选得到15个在肝细胞癌中呈差异表达的线粒体相关基因。EPHX2在肝细胞癌组织中的表达水平显著降低(P<0.01)。EPHX2表达水平与肝癌患者性别、分期和级别有关,而与年龄、T分期等因素无关。与EPHX2低表达组肝癌患者相比,EPHX2高表达组肝癌患者预后较好。功能富集结果显示,EPHX2与补体途径、脂肪酸降解等信号通路有关。蛋白质互作网络结果显示,EPHX2与HAO1、AGXT、ACOX1、GSTκ1、SCP-2、CAT、CYP2C8,CYP2C9,CYP2B6,和CYP2J2等密切相关。GSEA结果显示,EPHX2低表达组与肝癌细胞增殖、肝癌复发等基因集正相关。荧光定量PCR和GEO数据集验证结果显示,EPHX2在肝细胞癌组织和肝癌细胞株中呈显著低表达。EPHX2在肝细胞癌中呈显著低表达,提示其可能在肝细胞癌发生发展过程中发挥抑癌基因作用,但具体作用机制还需进一步验证。     

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 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.