Development Of A Three-Gene Prognostic Signature For Hepatitis B Virus Associated Hepatocellular Carcinoma Based On Integrated Transcriptomic Analysis

ECHDC2 inhibits the proliferation of gastric cancer cells by binding with NEDD4 to degrade MCCC2 and reduce aerobic glycolysis

BACKGROUND

The Enoyl-CoA hydratase/isomerase family plays a crucial role in the metabolism of tumors, being crucial for maintaining the energy balance and biosynthetic needs of cancer cells. However, the enzymes within this family that are pivotal in gastric cancer (GC) remain unclear.

METHODS

We employed bioinformatics techniques to identify key Enoyl-CoA hydratase/isomerase in GC. The expression of ECHDC2 and its clinical significance were validated through tissue microarray analysis. The role of ECHDC2 in GC was further assessed using colony formation assays, CCK8 assay, EDU assay, Glucose and lactic acid assay, and subcutaneous tumor experiments in nude mice. The mechanism of action of ECHDC2 was validated through Western blotting, Co-immunoprecipitation, and immunofluorescence experiments.

RESULTS

Our analysis of multiple datasets indicates that low expression of ECHDC2 in GC is significantly associated with poor prognosis. Overexpression of ECHDC2 notably inhibits aerobic glycolysis and proliferation of GC cells both in vivo and in vitro. Further experiments revealed that overexpression of ECHDC2 suppresses the P38 MAPK pathway by inhibiting the protein level of MCCC2, thereby restraining glycolysis and proliferation in GC cells. Ultimately, it was discovered that ECHDC2 promotes the ubiquitination and subsequent degradation of MCCC2 protein by binding with NEDD4.

CONCLUSIONS

These findings underscore the pivotal role of the ECHDC2 in regulating aerobic glycolysis and proliferation in GC cells, suggesting ECHDC2 as a potential therapeutic target in GC.

Identification of potential immunotherapeutic targets and prognostic biomarkers in Graves' disease using weighted gene co-expression network analysis

Graves' disease (GD) is an autoimmune disorder characterized by hyperthyroidism resulting from autoantibody-induced stimulation of the thyroid gland. Despite recent advancements in understanding GD's pathogenesis, the molecular processes driving disease progression and treatment response remain poorly understood. In this study, we aimed to identify crucial immunogenic factors associated with GD prognosis and immunotherapeutic response. To achieve this, we implemented a comprehensive screening strategy that combined computational immunogenicity-potential scoring with multi-parametric cluster analysis to assess the immunomodulatory genes in GD-related subtypes involving stromal and immune cells. Utilizing weighted gene co-expression network analysis (WGCNA), we identified co-expressed gene modules linked to cellular senescence and immune infiltration in CD4+ and CD8+ GD samples. Additionally, gene set enrichment analysis enabled the identification of hallmark pathways distinguishing high- and low-immune subtypes. Our WGCNA analysis revealed 21 gene co-expression modules comprising 1,541 genes associated with immune infiltration components in various stages of GD, including T cells, M1 and M2 macrophages, NK cells, and Tregs. These genes primarily participated in T cell proliferation through purinergic signaling pathways, particularly neuroactive ligand-receptor interactions, and DNA binding transcription factor activity. Three genes, namely PRSS1, HCRTR1, and P2RY4, exhibited robustness in GD patients across multiple stages and were involved in immune cell infiltration during the late stage of GD (p < 0.05). Importantly, HCRTR1 and P2RY4 emerged as potential prognostic signatures for predicting overall survival in high-immunocore GD patients (p < 0.05). Overall, our study provides novel insights into the molecular mechanisms driving GD progression and highlights potential key immunogens for further investigation. These findings underscore the significance of immune infiltration-related cellular senescence in GD therapy and present promising targets for the development of new immunotherapeutic strategies.

Machine Learning Methods for Survival Analysis with Clinical and Transcriptomics Data of Breast Cancer

Breast cancer is one of the most common cancers in women worldwide, which causes an enormous number of deaths annually. However, early diagnosis of breast cancer can improve survival outcomes enabling simpler and more cost-effective treatments. The recent increase in data availability provides unprecedented opportunities to apply data-driven and machine learning methods to identify early-detection prognostic factors capable of predicting the expected survival and potential sensitivity to treatment of patients, with the final aim of enhancing clinical outcomes. This tutorial presents a protocol for applying machine learning models in survival analysis for both clinical and transcriptomic data. We show that integrating clinical and mRNA expression data is essential to explain the multiple biological processes driving cancer progression. Our results reveal that machine-learning-based models such as random survival forests, gradient boosted survival model, and survival support vector machine can outperform the traditional statistical methods, i.e., Cox proportional hazard model. The highest C-index among the machine learning models was recorded when using survival support vector machine, with a value 0.688, whereas the C-index recorded using the Cox model was 0.677. Shapley Additive Explanation (SHAP) values were also applied to identify the feature importance of the models and their impact on the prediction outcomes.

Identification of diagnostic mRNA biomarkers in whole blood for ankylosing spondylitis using WGCNA and machine learning feature selection

Ankylosing spondylitis (AS) is a common inflammatory spondyloarthritis affecting the spine and sacroiliac joint that finally results in sclerosis of the axial skeleton. Aside from human leukocyte antigen B27, transcriptomic biomarkers in blood for AS diagnosis still remain unknown. Hence, this study aimed to investigate credible AS-specific mRNA biomarkers from the whole blood of AS patients by analyzing an mRNA expression profile (GSE73754) downloaded Gene Expression Omnibus, which includes AS and healthy control blood samples. Weighted gene co-expression network analysis was performed and revealed three mRNA modules associated with AS. By performing gene set enrichment analysis, the functional annotations of these modules revealed immune biological processes that occur in AS. Several feature mRNAs were identified by analyzing the hubs of the protein-protein interaction network, which was based on the intersection between differentially expressed mRNAs and mRNA modules. A machine learning-based feature selection method, SVM-RFE, was used to further screen out 13 key feature mRNAs. After verifying by qPCR, IL17RA, Sqstm1, Picalm, Eif4e, Srrt, Lrrfip1, Synj1 and Cxcr6 were found to be significant for AS diagnosis. Among them, Cxcr6, IL17RA and Lrrfip1 were correlated with severity of AS symptoms. In conclusion, our findings provide a framework for identifying the key mRNAs in whole blood of AS that is conducive for the development of novel diagnostic markers for AS.

miR-15a and miR-20b sensitize hepatocellular carcinoma cells to sorafenib through repressing CDC37L1 and consequent PPIA downregulation

Sorafenib is a classical targeted drug for the treatment of advanced hepatocellular carcinoma (HCC), but intrinsic resistance severely limited its therapeutic effects. In the present study, we aimed to identify crucial genes in HCC cells that affect sorafenib resistance by a CRISPR/Cas9 genome-scale screening. The results indicated that the deficiency of miR-15a and miR-20b contributed to sorafenib resistance, whereas exogenous expression of miR-15a and miR-20b enhanced sorafenib sensitivity of HCC cells by cell viability, colony formation, and flow cytometry analyses. Further analyses revealed that cell division cycle 37 like 1 (CDC37L1) as a common target of miR-15a and 20b, was negatively regulated by the two miRNAs and could enhance sorafenib resistance of HCC cells in vitro and in vivo. Mechanistically, CDC37L1, as a cochaperone, effectively increased the expression of peptidylprolyl isomerase A (PPIA) through strengthening the binding between heat shock protein 90 (HSP90) and PPIA. The results from immunohistochemical staining of a HCC tissue microarray revealed a positive association between CDC37L1 and PPIA expression, and high expression of CDC37L1 and PPIA predicted worse prognosis of HCC patients after sorafenib therapy. Taken together, our findings reveal crucial roles of miR-15a, miR-20b, CDC37L1, and PPIA in sorafenib response of HCC cells. These factors may serve as therapeutic targets and predict prognosis for HCC treated with sorafenib.

Comprehensive Analysis Identified Mutation-Gene Signature Impacts the Prognosis Through Immune Function in Hepatocellular Carcinoma

Background

Hepatocellular carcinoma (HCC) is a life-threatening and refractory malignancy with poor outcome. Genetic mutations are the hallmark of cancer. Thus far, there is no comprehensive prognostic model constructed by mutation-gene transcriptome in HCC. The prognostic value of mutation-gene signature in HCC remains elusive.

Methods

RNA expression profiles and the corresponding clinical information were recruited from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. The least absolute shrinkage and selection operator (LASSO) Cox regression analysis was employed to establish gene signature. Kaplan–Meier curve and time-dependent receiver operating characteristic curve were implemented to evaluate the prognostic value. The Wilcoxon test was performed to analyze the expression of immune checkpoint genes, cell cycle genes, and tumor drug resistance genes in different risk groups. Finally, quantitative real-time PCR (qRT-RCR) and immunohistochemistry (IHC) were performed to validate the mRNA and protein expression between HCC and adjacent nontumorous tissues in an independent cohort.

Results

A prognostic model consisting of five mutated genes was established by LASSO Cox regression analysis. The prognostic model classified patients into high- and low-risk groups. Compared with the low‐risk group, patients in the high‐risk group had significantly worse survival results. The prognostic model can accurately predict the overall survival of HCC patients and predict overall survival more accurately when combined with stage. Furthermore, the immune checkpoint genes, cell cycle genes, and tumor drug resistance genes were higher expressed in the high-risk group compared in the low-risk group. In addition, the expression level of prognostic signature genes was validated in an independent sample cohort, which was consistent with RNA sequencing expression in the TCGA database.

Conclusion

The prediction model of HCC constructed using mutation-related genes is of great significance for clinical decision making and the personalized treatment of patients with HCC.

Comparative mRNA/micro-RNA co-expression network drives melanomagenesis by promoting epithelial-mesenchymal transition and vasculogenic mimicry signaling

This study aimed to identify a novel disease-associated differentially co-expressed mRNA-microRNA (miRNA) that is associated with vasculogenic mimicry (VM) and epithelial-to-mesenchymal transition (EMT) network at different stages of melanoma. By applying weighted gene co-expression network analysis, we constructed a VM+EMT biological network with the available microarray dataset downloaded from a public database. Quantitative real-time PCR, immunohistochemical staining, and CD31-periodic acid solution dual staining were performed to confirm the expression of genes associated with EMT and VM formation in subjects with malignant melanoma (n = 18) and primary melanoma (n = 13) and in healthy subjects (n = 10). Our findings suggested that phosphatidylserine-specific phospholipase A1-alpha (PLA1A) and dermokine (DMKN) genes function as oncogenes that trigger VM and EMT processes during melanomagenesis on interaction with miR-370, miR-563, and miR-770-5p. PLA1A and DMKN genes can be considered potential VM+EMT network-based diagnostic biomarkers for distinguishing between melanoma patients. We postulate that a network with altered PLA1A/miR-563 and DMNK/miR-770-5p/miR-370 may contribute to melanomagenesis by triggering the EMT signaling pathway and VM formation. This study provides a potentially valuable approach for the early diagnosis and prognosis of melanoma progression.

Accurate transcriptome assembly by Nanopore RNA sequencing reveals novel functional transcripts in hepatocellular carcinoma

The long reads of Nanopore sequencing permit accurate transcript assembly and ease in discovering novel transcripts with potentially important functions in cancers. The wide adoption of Nanopore sequencing for transcript quantification, however, is largely limited by high costs. To address this issue, we developed a bioinformatics software, NovelQuant, that can specifically quantify long-read-assembled novel transcripts with short-read sequencing data. Nanopore Direct RNA Sequencing was carried out on three hepatocellular carcinoma (HCC) patients' tumor, matched portal vein tumor thrombus, and peritumor to reconstruct the HCC transcriptome. Then, based on the reconstructed transcriptome, NovelQuant was applied on Illumina RNA sequencing data of 59 HCC patients' tumor and paired peritumor to quantify novel transcripts. Our further analysis revealed 361 novel transcripts dysregulated in HCC and that 101 of them were significantly associated with prognosis. There were 19 novel prognostic transcripts predicted to be long noncoding RNAs (lncRNAs), and some of them had regulatory targets that were reported to be associated with HCC. Additionally, 42 novel prognostic transcripts were predicted to be protein-coding mRNAs, and many of them could be involved in xenobiotic metabolism. Moreover, the tumor-suppressive roles of two representative novel prognostic transcripts, CDO1-novel (lncRNA) and CYP2A6-novel (protein-coding mRNA), were further functionally validated during HCC progression. Overall, the current study shows a possibility of combining long- and short-read sequencing to explore functionally important novel transcripts in HCC with accuracy and cost-efficiency, which expands the pool of molecular biomarkers that could enhance our understanding of the molecular mechanisms of HCC.

Robust rank aggregation and cibersort algorithm applied to the identification of key genes in head and neck squamous cell cancer

OBJECTIVE

Although multiple hub genes have been identified in head and neck squamous cell cancer (HNSCC) in recent years, because of the limited sample size and inconsistent bioinformatics analysis methods, the results are not reliable. Therefore, it is urgent to use reliable algorithms to find new prognostic markers of HNSCC.

METHOD

The Robust Rank Aggregation (RRA) method was used to integrate 8 microarray datasets of HNSCC downloaded from the Gene Expression Omnibus (GEO) database to screen differentially expressed genes (DEGs). Later, Gene Ontology (GO) functional annotation together with Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was carried out to discover functions of those discovered DEGs. According to the KEGG results, those discovered DEGs showed tight association with the occurrence and development of HNSCC. Then cibersort algorithm was used to analyze the infiltration of immune cells of HNSCC and we found that the main infiltrated immune cells were B cells, dendritic cells and macrophages. A protein-protein interaction (PPI) network was established; moreover, key modules were also constructed to select 5 hub genes from the whole network using cytoHubba. 3 hub genes showed significant relationship with prognosis for TCGA-derived HNSCC patients.

RESULT

The potent DEGs along with hub genes were selected by the combined bioinformatic approach. AURKA, BIRC5 and UBE2C genes may be the potential prognostic biomarker and therapeutic targets of HNSCC.

CONCLUSIONS

The Robust Rank Aggregation method and cibersort algorithm method can accurately predict the potential prognostic biomarker and therapeutic targets of HNSCC through multiple GEO datasets.

CDC37L1 acts as a suppressor of migration and proliferation in gastric cancer by down-regulating CDK6

The co-chaperone protein CDC37 (Cell division cycle 37) is well known to regulate multiple protein kinases and involved in tumor progression. However to date, little is known about its analogue CDC37L1 (Cell division cycle 37 like 1) in tumorigenesis. This study aimed to explore the expression and function of CDC37L1 in gastric cancer (GC). The immunohistochemical staining in a tissue microarray showed a weak expression of CDC37L1 in high grade GC tissues compared with low grade tissues. Consistently, data from online database analysis demonstrated that CDC37L1 level was decreased in stage 4 patients and low expression of CDC37L1 indicated a poor prognosis. Functional studies revealed that CDC37L1 could inhibit GC cell proliferation and migration in CCK8, EdU incorporation, colony formation and transwell assays. In the meantime, CDC37L1 also inhibited the tumorigenicity of GC cells in nude mice. Mechanistically, we found that CDC37L1 had an impact on CDK6 protein expression by western blotting. Palbociclib, a specific CDK4/6 inhibitor, was discovered to block the rapid growth phenotype of GC cells induced by CDC37L1 silencing. Taken together, these findings unveiled a tumor-suppressive role of CDC37L1 in GC, and CDK6 may act as a downstream effector in this process.

System-Wide Pollution of Biomedical Data: Consequence of the Search for Hub Genes of Hepatocellular Carcinoma Without Spatiotemporal Consideration

Biomedical institutions rely on data evaluation and are turning into data factories. Big-data storage centers, supercomputing systems, and increased algorithmic efficiency allow us to analyze the ever-increasing amount of data generated every day in biomedical research centers. In network science, the principal intrinsic problem is how to integrate the data and information from different experiments on genes or proteins. Data curation is an essential process in annotating new functional data to known genes or proteins, undertaken by a biobank curator, which is then reflected in the calculated networks. We provide an example of how protein-protein networks today have space-time limits. The next step is the integration of data and information from different biobanks. Omics data and networks are essential parts of this step but also have flawed protocols and errors. Consider data from patients with cancer: from biopsy procedures to experimental tests, to archiving methods and computational algorithms, these are continuously handled so require critical and continuous "updates" to obtain reproducible, reliable, and correct results. We show, as a second example, how all this distorts studies in cellular hepatocellular carcinoma. It is not unlikely that these flawed data have been polluting biobanks for some time before stringent conditions for the veracity of data were implemented in Big data. Therefore, all this could contribute to errors in future medical decisions.

Identification of a Transcriptomic Prognostic Signature by Machine Learning Using a Combination of Small Cohorts of Prostate Cancer

Determining which treatment to provide to men with prostate cancer (PCa) is a major challenge for clinicians. Currently, the clinical risk-stratification for PCa is based on clinico-pathological variables such as Gleason grade, stage and prostate specific antigen (PSA) levels. But transcriptomic data have the potential to enable the development of more precise approaches to predict evolution of the disease. However, high quality RNA sequencing (RNA-seq) datasets along with clinical data with long follow-up allowing discovery of biochemical recurrence (BCR) biomarkers are small and rare. In this study, we propose a machine learning approach that is robust to batch effect and enables the discovery of highly predictive signatures despite using small datasets. Gene expression data were extracted from three RNA-Seq datasets cumulating a total of 171 PCa patients. Data were re-analyzed using a unique pipeline to ensure uniformity. Using a machine learning approach, a total of 14 classifiers were tested with various parameters to identify the best model and gene signature to predict BCR. Using a random forest model, we have identified a signature composed of only three genes (JUN, HES4, PPDPF) predicting BCR with better accuracy [74.2%, balanced error rate (BER) = 27%] than the clinico-pathological variables (69.2%, BER = 32%) currently in use to predict PCa evolution. This score is in the range of the studies that predicted BCR in single-cohort with a higher number of patients. We showed that it is possible to merge and analyze different small and heterogeneous datasets altogether to obtain a better signature than if they were analyzed individually, thus reducing the need for very large cohorts. This study demonstrates the feasibility to regroup different small datasets in one larger to identify a predictive genomic signature that would benefit PCa patients.

Integrated analysis identifies a long non-coding RNAs-messenger RNAs signature for prediction of prognosis in hepatitis B virus-hepatocellular carcinoma patients

Hepatitis B virus (HBV) infection is a leading cause of hepatocellular carcinoma (HCC), but HBV-HCC related prognosis signature remains rarely investigated. This study was to identify an integrated long non-coding RNAs-messenger RNAs (lncRNA-mRNA) signature for prediction of overall survival (OS) and explore their underlying functions.One RNA-sequencing dataset (training set, n = 95) and one microarray dataset E-TABM-36 (validation set, n = 44) were collected. Least absolute shrinkage and selection operator analysis was performed to identify an lncRNA-mRNA prognosis signature. The OS difference of patients in the high-risk and low-risk risk groups was evaluated by Kaplan-Meier curve. Area under the receiver operating characteristic curve (AUC), Harrell concordance index (C-index) calculation, and multivariate analyses with clinical characteristics were used to determine the prognostic ability. Furthermore, a coexpression network was constructed to interpret the functions.Nine signature genes (3 lncRNAs and 6 mRNAs) were selected to generate the risk score model. Patients belonging to the high-risk group showed a significantly shorter survival than those of the low-risk group. The prediction accuracy of the risk score for 5-year OS was 0.936 and 0.905 for the training set and validation set, respectively. Also, this risk score was independent of various clinical variables for the prognosis prediction. Incorporation of the risk score remarkably increased the predictive power of the routine clinical prognostic factors (vascular invasion status, tumor recurrence status) (AUC = 0.942 vs 0.628; C-index = 0.7997 vs 0.6908). Furthermore, LncRNA insulin-like growth factor 2 antisense RNA (IGF2-AS) and long intergenic non-protein coding RNA 342 (LINC00342) were predicted to exert tumor suppression effects by regulating homeobox D1 (HOXD1) and secreted frizzled related protein 5 (SFRP5), respectively; while lncRNA rhophilin Rho GTPase binding protein 1 antisense RNA 1 (RHPN1-AS1) may possess carcinogenic potential by promoting the transcription of chromobox 2 (CBX2), cell division cycle 20 (CDC20), matrix metallopeptidase 12 (MMP12), stratifin (SFN), tripartite motif containing 16 (TRIM16), and uroplakin 3A (UPK3A). These mRNAs may be associated with cell proliferation or apoptosis related pathways.This study may provide a novel, effective prognostic biomarker, and some therapeutic targets for HBV-HCC patients.

A hypoxia-related signature for clinically predicting diagnosis, prognosis and immune microenvironment of hepatocellular carcinoma patients

BACKGROUND

Hypoxia plays an indispensable role in the development of hepatocellular carcinoma (HCC). However, there are few studies on the application of hypoxia molecules in the prognosis predicting of HCC. We aim to identify the hypoxia-related genes in HCC and construct reliable models for diagnosis, prognosis and recurrence of HCC patients as well as exploring the potential mechanism.

METHODS

Differentially expressed genes (DEGs) analysis was performed using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database and four clusters were determined by a consistent clustering analysis. Three DEGs closely related to overall survival (OS) were identified using Cox regression and LASSO analysis. Then the hypoxia-related signature was developed and validated in TCGA and International Cancer Genome Consortium (ICGC) database. The Gene Set Enrichment Analysis (GSEA) was performed to explore signaling pathways regulated by the signature. CIBERSORT was used for estimating the fractions of immune cell types.

RESULTS

A total of 397 hypoxia-related DEGs in HCC were detected and three genes (PDSS1, CDCA8 and SLC7A11) among them were selected to construct a prognosis, recurrence and diagnosis model. Then patients were divided into high- and low-risk groups. Our hypoxia-related signature was significantly associated with worse prognosis and higher recurrence rate. The diagnostic model also accurately distinguished HCC from normal samples and nodules. Furthermore, the hypoxia-related signature could positively regulate immune response. Meanwhile, the high-risk group had higher fractions of macrophages, B memory cells and follicle-helper T cells, and exhibited higher expression of immunocheckpoints such as PD1and PDL1.

CONCLUSIONS

Altogether, our study showed that hypoxia-related signature is a potential biomarker for diagnosis, prognosis and recurrence of HCC, and it provided an immunological perspective for developing personalized therapies.

Identification of Potential Therapeutic Targets and Immune Cell Infiltration Characteristics in Osteosarcoma Using Bioinformatics Strategy

Osteosarcoma is one of the most aggressive malignant bone tumors worldwide. Although great advancements have been made in its treatment owing to the advent of neoadjuvant chemotherapy, the problem of lung metastasis is a major obstacle in the improvement of survival outcomes. Thus, the aim of the present study is to screen novel and key biomarkers, which may act as potential prognostic markers and therapeutic targets in osteosarcoma. We utilized the robust rank aggregation (RRA) method to integrate three osteosarcoma microarray datasets downloaded from the Gene Expression Omnibus (GEO) database, and we identified the robust differentially expressed genes (DEGs) between primary and metastatic osteosarcoma tissues. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed to explore the functions of robust DEGs. The results of enrichment analysis showed that the robust DEGs were closely associated with osteosarcoma development and progression. Immune cell infiltration analysis was also conducted by CIBERSORT algorithm, and we found that macrophages are the most principal infiltrating immune cells in osteosarcoma, especially macrophages M0 and M2. Then, the protein–protein interaction network and key modules were constructed by Cytoscape, and 10 hub genes were selected by plugin cytoHubba from the whole network. The survival analysis of hub genes was also carried out based on the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. The integrated bioinformatics analysis was utilized to provide new insight into osteosarcoma development and metastasis and identified EGR1, CXCL10, MYC, and CXCR4 as potential biomarkers for prognosis of osteosarcoma.

Construction and Validation of Novel Diagnostic and Prognostic DNA Methylation Signatures for Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most prevalent life-threatening human cancers and the leading cause of cancer-related mortality, with increased global incidence within the last decade. Identification of effective diagnostic and prognostic biomarkers would enable reliable risk stratification and efficient screening of high-risk patients, thereby facilitating clinical decision-making. Herein, we performed a comprehensive, robust DNA methylation analysis based on genome-wide DNA methylation profiling. We constructed a diagnostic signature with five DNA methylation markers, which precisely distinguished HCC patients from normal controls. Cox regression and LASSO analysis were applied to construct a prognostic signature with four DNA methylation markers. A one-to-one correlation analysis was carried out between genes of the whole genome and our prognostic signature. Exploration of the biological function and the role of the underlying significantly correlated genes was conducted. A mixed dataset of 463 HCC patients and 253 normal controls, derived from six independent datasets, was used to valid the diagnostic signature. Results showed a specificity of 96.84% and sensitivity of 96.77%. Class scores for the diagnostic signature were significantly different between normal controls, individuals with liver diseases, and HCC patients. The present signature has the potential to serve as a biomarker to monitor health in normal controls. Additionally, HCC patients were successfully separated into low-risk and high-risk groups by the prognostic signature, with a better prognosis for patients in the low-risk group. Kaplan-Meier and ROC analysis confirmed that the prognostic signature performed well. We found eight of the top ten genes to positively correlate with risk scores of the prognostic signature, and to be involved in cell cycle regulation. This eight-gene panel also served as a prognostic signature. The robust evidence presented in this study therefore demonstrates the effectiveness of the prognostic signature. In summary, we constructed diagnostic and prognostic signatures, which have potential for use in diagnosis, surveillance, and prognostic prediction for HCC patients. Eight genes that were significantly and positively correlated with the prognostic signature were strongly associated with cell cycle processes. Therefore, the prognostic signature can be used as a guide by which to measure responsiveness to cell-cycle-targeting agents.

Integrated bioinformatics analysis for the identification of potential key genes affecting the pathogenesis of clear cell renal cell carcinoma

Clear cell renal cell carcinoma (CCRCC) is a typical type of RCC with the worst prognosis among the common epithelial neoplasms of the kidney. However, its molecular pathogenesis remains unknown. Therefore, the aim of the present study was to screen for effective and potential pathogenic biomarkers of CCRCC. The gene expression profile of the GSE16441, GSE36895, GSE40435, GSE46699, GSE66270 and GSE71963 datasets were downloaded from the Gene Expression Omnibus database. First, the limma package in R language was used to identify differentially expressed genes (DEGs) in each dataset. The robust and strong DEGs were explored using the robust rank aggregation method. A total of 980 markedly robust DEGs were identified (429 upregulated and 551 downregulated). According to Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis, these DEGs exhibited an obvious enrichment in various cancer-related biological pathways and functions. The Search Tool for the Retrieval of Interacting Genes/Proteins database was used for the construction of a protein-protein interaction (PPI) network, the Cytoscape MCODE plug-in for module analysis and the cytoHubba plug-in to identify hub genes from the aforementioned DEGs. A total of four key modules were identified in the PPI network. A total of six hub genes, including C-X-C motif chemokine ligand 12, bradykinin receptor B2, adenylate cyclase 7, calcium sensing receptor (CASR), kininogen 1 and lysophosphatidic acid receptor 5, were identified. The DEG results of the hub genes were verified using The Cancer Genome Atlas database, and CASR was found to be significantly associated with the prognosis of patients with CCRCC. In conclusion, the present study provided new insight and potential biomarkers for the diagnosis and prognosis of CCRCC.

LINC02499, a novel liver-specific long non-coding RNA with potential diagnostic and prognostic value, inhibits hepatocellular carcinoma cell proliferation, migration, and invasion

AIM

Liver-specific non-coding RNAs have been reported to play crucial roles in hepatocellular carcinoma (HCC). We investigated the possible biological performance of a novel liver-specific long non-coding RNA, LINC02499, in HCC.

METHODS

The association between LINC02499 expression and HCC was evaluated based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and then confirmed in a HCC cohort by quantitative real-time polymerase chain reaction. The effects of LINC02499 on HCC cells were verified by gain- and loss-of-function assays. Pathway enrichment analyses were used to explore the potential mechanism of LINC02499 in HCC.

RESULTS

LINC02499 expression was remarkably decreased in HCC tissues compared to adjacent non-tumor tissues based on TCGA (P < 0.001) and GEO databases (P < 0.001) and our HCC cohort (P < 0.001). Decreased LINC02499 was also significantly associated with poorer overall survival in both the TCGA database (P = 0.009) and our HCC cohort (P = 0.002). Furthermore, the receiver operating characteristic analysis indicated that LINC02499 showed a good performance in HCC diagnosis (area under the curve = 0.879, P < 0.001), and both sensitivity and specificity were 83.8%. In addition, up- and downregulated LINC02499 significantly impacted proliferation, migration, and invasion abilities of HCC cells in vitro. Pathway enrichment analyses revealed that the potential target genes of LINC02499 were involved in "Complement and coagulation cascades" and "Butanoate metabolism" pathways.

CONCLUSION

LINC02499 could be a potential novel diagnostic and prognostic biomarker for HCC patients, and it could exert a tumor suppressor role in the progression of HCC.

Identification of hub genes in papillary thyroid carcinoma: robust rank aggregation and weighted gene co-expression network analysis

BACKGROUND

Papillary thyroid carcinoma (PTC), which is the most common endocrine malignancy, has been steadily increasing worldwide in incidence over the years, while mechanisms underlying the pathogenesis and diagnostic for PTC are incomplete. The purpose of this study is to identify potential biomarkers for diagnosis of PTC, and provide new insights into pathogenesis of PTC.

METHODS

Based on weighted gene co-expression network analysis, Robust Rank Aggregation, functional annotation, GSEA and DNA methylation, were employed for investigating potential biomarkers for diagnosis of PTC.

RESULTS

Black and turquoise modules were identified in the gene co-expression network constructed by 1807 DEGs that from 6 eligible gene expression profiles of Gene Expression Omnibus database based on Robust Rank Aggregation and weighted gene co-expression network analysis. Hub genes were significantly down-regulated and the expression levels of the hub genes were different in different stages in hub gene verification. ROC curves indicated all hub genes had good diagnostic value for PTC (except for ABCA6 AUC = 89.5%, the 15 genes with AUC > 90%). Methylation analysis showed that hub gene verification ABCA6, ACACB, RMDN1 and TFPI were identified as differentially methylated genes, and the decreased expression level of these genes may relate to abnormal DNA methylation. Moreover, the expression levels of 8 top hub genes were correlated with tumor purity and tumor-infiltrating immune cells. These findings, including functional annotations and GSEA provide new insights into pathogenesis of PTC.

CONCLUSIONS

The hub genes and methylation of hub genes may as potential biomarkers provide new insights for diagnosis of PTC, and all these findings may be the direction to study the mechanisms underlying of PTC in the future.

Identification of hub genes in hepatocellular carcinoma using integrated bioinformatic analysis

The molecular mechanisms underlying hepatocellular carcinoma (HCC) progression remain largely undefined. Here, we identified 176 commonly upregulated genes in HCC tissues based on three Gene Expression Omnibus datasets and The Cancer Genome Atlas (TCGA) cohort. We integrated survival and methylation analyses to further obtain 12 upregulated genes for validation. These genes were overexpressed in HCC tissues at the transcription and protein levels, and increased mRNA levels were related to higher tumor grades and cancer stages. The expression of all markers was negatively associated with overall and disease-free survival in HCC patients. Most of these hub genes can promote HCC proliferation and/or metastasis. These 12 hub genes were also overexpressed and had strong prognostic value in many other cancer types. Methylation and gene copy number analyses indicated that the upregulation of these hub genes was probably due to hypomethylation or increased gene copy numbers. Further, the methylation levels of three genes, KPNA2, MCM3, and LRRC1, were associated with HCC clinical features. Moreover, the levels of most hub genes were related to immune cell infiltration in HCC microenvironments. Finally, we identified three upregulated genes (KPNA2, TARBP1, and RNASEH2A) that could comprehensively and accurately provide diagnostic and prognostic value for HCC patients.

IRF1 association with tumor immune microenvironment and use as a diagnostic biomarker for colorectal cancer recurrence

Colorectal cancer (CRC) is considered to be one of the most lethal cancer types globally, and its recurrence is a major treatment challenge. Identifying the factors involved when determining the risk of CRC recurrence is required to improve personalized therapy for patients with CRC. Based on the GSE39582 dataset, the present study demonstrated that a higher ratio of M1 macrophages and activated memory CD4+ T cells indicated a better recurrence-free survival (RFS) time for CRC, using CIBERSORT and Pearson's correlation analysis. Through weighted correlation network analysis (WGCNA), an immune-associated module was identified that was significantly positively correlated with the ratio of M1 macrophages and activated memory CD4+ T cells. In this module, using WGCNA and a protein-protein interaction network, interferon regulatory factor 1 (IRF1), chemokine ligand 5, ubiquitin/ISG15-conjugating enzyme E2 L6, guanylate binding protein 1 and interleukin 2 receptor subunit beta were identified as hub genes. Among these genes, univariate Cox and multivariate Cox analysis revealed that IRF1 may be a potential diagnostic biomarker for RFS in patients with CRC. This was further validated using The Cancer Genome Atlas data. Gene set enrichment analysis demonstrated that IRF1 influenced the genes and pathways that are associated with immune cell recruitment and activation. Additionally, the DNA methylation of cg27587780 and cg15375424 CpG sites in the IRF1 gene region was indicated to be negatively correlated with IRF1 mRNA expression and positively correlated with the recurrence of CRC. Collectively, the results of the present study demonstrated that IRF1 may be a potential diagnostic biomarker for RFS in patients with CRC.

Identification of a Five-Autophagy-Related-lncRNA Signature as a Novel Prognostic Biomarker for Hepatocellular Carcinoma

Although great progresses have been made in the diagnosis and treatment of hepatocellular carcinoma (HCC), its prognostic marker remains controversial. In this current study, weighted correlation network analysis and Cox regression analysis showed significant prognostic value of five autophagy-related long non-coding RNAs (AR-lncRNAs) (including TMCC1-AS1, PLBD1-AS1, MKLN1-AS, LINC01063, and CYTOR) for HCC patients from data in The Cancer Genome Atlas. By using them, we constructed a five-AR-lncRNA prognostic signature, which accurately distinguished the high- and low-risk groups of HCC patients. All of the five AR lncRNAs were highly expressed in the high-risk group of HCC patients. This five-AR-lncRNA prognostic signature showed good area under the curve (AUC) value (AUC = 0.751) for the overall survival (OS) prediction in either all HCC patients or HCC patients stratified according to several clinical traits. A prognostic nomogram with this five-AR-lncRNA signature predicted the 3- and 5-year OS outcomes of HCC patients intuitively and accurately (concordance index = 0.745). By parallel comparison, this five-AR-lncRNA signature has better prognosis accuracy than the other three recently published signatures. Furthermore, we discovered the prediction ability of the signature on therapeutic outcomes of HCC patients, including chemotherapy and immunotherapeutic responses. Gene set enrichment analysis and gene mutation analysis revealed that dysregulated cell cycle pathway, purine metabolism, and TP53 mutation may play an important role in determining the OS outcomes of HCC patients in the high-risk group. Collectively, our study suggests a new five-AR-lncRNA prognostic signature for HCC patients.

A genomic-clinical nomogram predicting recurrence-free survival for patients diagnosed with hepatocellular carcinoma

Liver resection surgery is the most commonly used treatment strategy for patients diagnosed with hepatocellular carcinoma (HCC). However, there is still a chance for recurrence in these patients despite the survival benefits of this procedure. This study aimed to explore recurrence-related genes (RRGs) and establish a genomic-clinical nomogram for predicting postoperative recurrence in HCC patients. A total of 123 differently expressed genes and three RRGs (PZP, SPP2, and PRC1) were identified from online databases via Cox regression and LASSO logistic regression analyses and a gene-based risk model containing RRGs was then established. The Harrell’s concordance index (C-index), receiver operating characteristic (ROC) curves and calibration curves showed that the model performed well. Finally, a genomic-clinical nomogram incorporating the gene-based risk model, AJCC staging system, and Eastern Cooperative Oncology Group performance status was constructed to predict the 1-, 2-, and 3-year recurrence-free survival rates (RFS) for HCC patients. The C-index, ROC analysis, and decision curve analysis were good indicators of the nomogram’s performance. In conclusion, we identified three reliable RRGs associated with the recurrence of cancer and constructed a nomogram that performed well in predicting RFS for HCC patients. These findings could enrich our understanding of the mechanisms for HCC recurrence, help surgeons predict patients’ prognosis, and promote HCC treatment.

Identification of RAD54 homolog B as a promising therapeutic target for breast cancer

Breast cancer is a recognized threat to the health of women globally. Due to the lack of the knowledge about the molecular pathogenesis of breast cancer, therapeutic strategies remain inadequate, especially for aggressive breast cancer. In the present study, sequential bioinformatics analysis was performed using data from the GSE20711 dataset, and the results demonstrated that three genes may impact the survival of patients with breast cancer. One of these genes, RAD54 homolog B (RAD54B), may be a potential prognostic factor for breast cancer. A signature was established that could evaluate the overall survival for patients with breast cancer based on the risk score calculated from RAD54B expression and the Tumor-Node-Metastasis (TNM) stage [risk score=expRAD54B × 0.236 + TNM stage (I/II=0 or III/IV=1) ×1.025]. In addition, based on the GSE85871 dataset and inhibitory assay, the study identified a natural compound, Japonicone A, which may reduce the proliferation of breast cancer cells by inhibiting the expression of RAD54B. Overall, the present study identified a novel candidate gene and a candidate compound as promising therapeutic targets for the treatment of breast cancer.

NCF2, MYO1F, S1PR4, and FCN1 as potential noninvasive diagnostic biomarkers in patients with obstructive coronary artery: A weighted gene co-expression network analysis

This study aims to explore the predictive noninvasive biomarker for obstructive coronary artery disease (CAD). By using the data set GSE90074, weighted gene co-expression network analysis (WGCNA), and protein-protein interactive network, construction of differentially expressed genes in peripheral blood mononuclear cells was conducted to identify the most significant gene clusters associated with obstructive CAD. Univariate and multivariate stepwise logistic regression analyses and receiver operating characteristic analysis were used to predicate the diagnostic accuracy of biomarker candidates in the detection of obstructive CAD. Furthermore, functional prediction of candidate gene biomarkers was further confirmed in ST-segment elevation myocardial infarction (STEMI) patients or stable CAD patients by using the datasets of GSE62646 and GSE59867. We found that the blue module discriminated by WGCNA contained 13 hub-genes that could be independent risk factors for obstructive CAD (P < .05). Among these 13 hub-genes, a four-gene signature including neutrophil cytosol factor 2 (NCF2, P = .025), myosin-If (MYO1F, P = .001), sphingosine-1-phosphate receptor 4 (S1PR4, P = .015), and ficolin-1 (FCN1, P = .012) alone or combined with two risk factors (male sex and hyperlipidemia) may represent potential diagnostic biomarkers in obstructive CAD. Furthermore, the messenger RNA levels of NCF2, MYO1F, S1PR4, and FCN1 were higher in STEMI patients than that in stable CAD patients, although S1PR4 showed no statistical difference (P > .05). This four-gene signature could also act as a prognostic biomarker to discriminate STEMI patients from stable CAD patients. These findings suggest a four-gene signature (NCF2, MYO1F, S1PR4, and FCN1) alone or combined with two risk factors (male sex and hyperlipidemia) as a promising prognostic biomarker in the diagnosis of STEMI. Well-designed cohort studies should be implemented to warrant the diagnostic value of these genes in clinical purpose.

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

Background

CpG island methylator phenotype (CIMP), a common biological phenomenon characterized by a subset of concurrently methylated genes, can have an influence on the progression of multiple cancers. However, the potential mechanism of CIMP in hepatocarcinogenesis and its clinical relevance remains only partially understood.

Methods: We used a methylation array from the cancer genome atlas (TCGA) to stratify HCC patients into different CIMP subtypes, and evaluated their correlation with clinical characteristics. In addition, mutation, CNV, and transcriptome profiles were also utilized to evaluate the distinctive genomic patterns correlated with CIMP. Finally, a CIMP-associated prognostic model (CPM) was trained and validated using four independent datasets.

Findings

A subgroup of patients was identified as having CIMP-H, which was associated with worse OS and DFS. Gene enrichment analysis indicated that the terms “liver cancer with EPCAM up”, “tumor invasiveness up”, “methyltransferase complex”, and “translational initiation” were enriched in CIMP-H subgroup. Notably, somatic mutation analysis indicated that CIMP-H patients presented with a higher mutation burden of BRD4, DDIAS and NOX1. Moreover, four CPM associated genes could significantly categorize patients into low- and high-risk groups in the training dataset and another 3 independent validation datasets. Finally, a nomogram incorporating a classifier based on four mRNAs, pathological M stage and CIMP status was established, which showed a favorable discriminating ability and might contribute to clinical decision-making for HCC.

Interpretation

Our work highlights the potential clinical application value of CPM in predicting the overall survival of HCC patients and the mechanisms underlying the role of CIMP in hepatocarcinogenesis.

Fund

This work was supported by the State Key Project on Infectious Diseases of China (2018ZX10723204-003), the National Nature Science Foundation of China (Nos. 81874065, 81500565, 81874149, 81572427, and 81401997), the Hepato-Biliary-Pancreatic Malignant Tumor Investigation Fund of Chen Xiao-ping Foundation for the Development of Science and Technology of Hubei Province (CXPJJH11800001-2018356).

Identification of hub genes in prostate cancer using robust rank aggregation and weighted gene co-expression network analysis

The pathogenic mechanisms of prostate cancer (PCa) remain to be defined. In this study, we utilized the Robust Rank Aggregation (RRA) method to integrate 10 eligible PCa microarray datasets from the GEO and identified a set of significant differentially expressed genes (DEGs) between tumor samples and normal, matched specimens. To explore potential associations between gene sets and PCa clinical features and to identify hub genes, we utilized WGCNA to construct gene co-expression networks incorporating the DEGs screened with the use of RRA. From the key module, we selected LMNB1, TK1, ZWINT, and RACGAP1 for validation. We found that these genes were up-regulated in PCa samples, and higher expression levels were associated with higher Gleason scores and tumor grades. Moreover, ROC and K-M plots indicated these genes had good diagnostic and prognostic value for PCa. On the other hand, methylation analyses suggested that the abnormal up-regulation of these four genes likely resulted from hypomethylation, while GSEA and GSVA for single hub gene revealed they all had a close association with proliferation of PCa cells. These findings provide new insight into PCa pathogenesis, and identify LMNB1, TK1, RACGAP1 and ZWINT as candidate biomarkers for diagnosis and prognosis of PCa.

Identification of an eight-gene signature for survival prediction for patients with hepatocellular carcinoma based on integrated bioinformatics analysis

Background

Hepatocellular carcinoma (HCC) remains one of the leading causes of cancer-related death worldwide. Despite recent advances in imaging techniques and therapeutic intervention for HCC, the low overall 5-year survival rate of HCC patients remains unsatisfactory. This study aims to find a gene signature to predict clinical outcomes in HCC.

Methods

Bioinformatics analysis including Cox’s regression analysis, Kaplan-Meier (KM) and receiver operating characteristic curve (ROC) analysis and the random survival forest algorithm were performed to mine the expression profiles of 553 hepatocellular carcinoma (HCC) patients from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) public database.

Results

We selected a signature comprising eight protein-coding genes (DCAF13, FAM163A, GPR18, LRP10, PVRIG, S100A9, SGCB, and TNNI3K) in the training dataset (AUC = 0.77 at five years, n = 332). The signature stratified patients into high- and low-risk groups with significantly different survival in the training dataset (median 2.20 vs. 8.93 years, log-rank test P < 0.001) and in the test dataset (median 2.68 vs. 4.24 years, log-rank test P = 0.004, n = 221, GSE14520). Further multivariate Cox regression analysis showed that the signature was an independent prognostic factor for patients with HCC. Compared with TNM stage and another reported three-gene model, the signature displayed improved survival prediction power in entire dataset (AUC signature = 0.66 vs. AUC TNM = 0.64 vs. AUC gene model = 0.60, n = 553). Stratification analysis shows that it can be used as an auxiliary marker for many traditional staging models.

Conclusions

We constructed an eight-gene signature that can be a novel prognostic marker to predict the survival of HCC patients.

Development and validation of a TP53-associated immune prognostic model for hepatocellular carcinoma

BACKGROUND

TP53 mutation is the most common mutation in hepatocellular carcinoma (HCC), and it affects the progression and prognosis of HCC. We investigated how TP53 mutation regulates the HCC immunophenotype and thus affects the prognosis of HCC.

METHODS

We investigated TP53 mutation status and RNA expression in different populations and platforms and developed an immune prognostic model (IPM) based on immune-related genes that were differentially expressed between TP53^WT and TP53^MUT HCC samples. Then, the influence of the IPM on the immune microenvironment in HCC was comprehensively analysed.

FINDINGS

TP53 mutation resulted in the downregulation of the immune response in HCC. Thirty-seven of the 312 immune response-related genes were differentially expressed based on TP53 mutation status. An IPM was established and validated based on 865 patients with HCC to differentiate patients with a low or high risk of poor survival. A nomogram was also established for clinical application. Functional enrichment analysis showed that the humoral immune response and immune system diseases pathway represented the major function and pathway, respectively, related to the IPM genes. Moreover, we found that the patients in the high-risk group had higher fractions of T cells follicular helper, T cells regulatory (Tregs) and macrophages M0 and presented higher expression of CTLA-4, PD-1 and TIM-3 than the low-risk group.

INTERPRETATION

TP53 mutation is strongly related to the immune microenvironment in HCC. Our IPM, which is sensitive to TP53 mutation status, may have important implications for identifying subgroups of HCC patients with low or high risk of unfavourable survival. FUND: This work was supported by the International Science and Technology Cooperation Projects (2016YFE0107100), the Capital Special Research Project for Health Development (2014-2-4012), the Beijing Natural Science Foundation (L172055 and 7192158), the National Ten Thousand Talent Program, the Fundamental Research Funds for the Central Universities (3332018032), and the CAMS Innovation Fund for Medical Science (CIFMS) (2017-I2M-4-003 and 2018-I2M-3-001).

Identification of special key genes for alcohol-related hepatocellular carcinoma through bioinformatic analysis

Background

Alcohol-related hepatocellular carcinoma (HCC) was reported to be diagnosed at a later stage, but the mechanism was unknown. This study aimed to identify special key genes (SKGs) during alcohol-related HCC development and progression.

Methods

The mRNA data of 369 HCC patients and the clinical information were downloaded from the Cancer Genome Atlas project (TCGA). The 310 patients with certain HCC-related risk factors were included for analysis and divided into seven groups according to the risk factors. Survival analyses were applied for the HCC patients of different groups. The patients with hepatitis B virus or hepatitis C virus infection only were combined into the HCC-V group for further analysis. The differentially expressed genes (DEGs) between the HCCs with alcohol consumption only (HCC-A) and HCC-V tumors were identified through limma package in R with cutoff criteria│log2 fold change (logFC)|>1.0 and p < 0.05. The DEGs between eight alcohol-related HCCs and their paired normal livers of GSE59259 from the Gene Expression Omnibus (GEO) were identified through GEO2R (a built-in tool in GEO database) with cutoff criteria |logFC|> 2.0 and adj.p < 0.05. The intersection of the two sets of DEGs was considered SKGs which were then investigated for their specificity through comparisons between HCC-A and other four HCC groups. The SKGs were analyzed for their correlations with HCC-A stage and grade and their prognostic power for HCC-A patients. The expressional differences of the SKGs in the HCCs in whole were also investigated through Gene Expression Profiling Interactive Analysis (GEPIA). The SKGs in HCC were validated through Oncomine database analysis.

Results

Pathological stage is an independent prognostic factor for HCC patients. HCC-A patients were diagnosed later than HCC patients with other risk factors. Ten SKGs were identified and nine of them were confirmed for their differences in paired samples of HCC-A patients. Three (SLC22A10, CD5L, and UROC1) and four (SLC22A10, UROC1, CSAG3, and CSMD1) confirmed genes were correlated with HCC-A stage and grade, respectively. SPP2 had a lower trend in HCC-A tumors and was negatively correlated with HCC-A stage and grade. The SKGs each was differentially expressed between HCC-A and at least one of other HCC groups. CD5L was identified to be favorable prognostic factor for overall survival while CSMD1 unfavorable prognostic factor for disease-free survival for HCC-A patients and HCC patients in whole. Through Oncomine database, the dysregulations of the SKGs in HCC and their clinical significance were confirmed.

Conclusion

The poor prognosis of HCC-A patients might be due to their later diagnosis. The SKGs, especially the four stage-correlated genes (CD5L, SLC22A10, UROC1, and SPP2) might play important roles in HCC development, especially alcohol-related HCC development and progression. CD5L might be useful for overall survival and CSMD1 for disease-free survival predication in HCC, especially alcohol-related HCC.

Comprehensive analysis of potential prognostic genes for the construction of a competing endogenous RNA regulatory network in hepatocellular carcinoma

BACKGROUND

Hepatocellular carcinoma (HCC) is an extremely common malignant tumor with worldwide prevalence. The aim of this study was to identify potential prognostic genes and construct a competing endogenous RNA (ceRNA) regulatory network to explore the mechanisms underlying the development of HCC.

METHODS

Integrated analysis was used to identify potential prognostic genes in HCC with R software based on the GSE14520, GSE17548, GSE19665, GSE29721, GSE60502, and the Cancer Genome Atlas databases. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway-enrichment analyses were performed to explore the molecular mechanisms of potential prognostic genes. Differentially expressed miRNAs (DEMs) and lncRNAs (DELs) were screened based on the Cancer Genome Atlas database. An lncRNA-miRNA-mRNA ceRNA regulatory network was constructed based on information about interactions derived from the miRcode, TargetScan, miRTarBase, and miRDB databases.

RESULTS

A total of 152 potential prognostic genes were screened that were differentially expressed in HCC tissue and significantly associated with overall survival of HCC patients. There were 13 key potential prognostic genes in the ceRNA regulatory network: eleven upregulated genes (CCNB1, CEP55, CHEK1, EZH2, KPNA2, LRRC1, PBK, RRM2, SLC7A11, SUCO, and ZWINT) and two downregulated genes (ACSL1 and CDC37L1) whose expression might be regulated by eight DEMs and 61 DELs. Kaplan-Meier curve analysis showed that nine DELs (AL163952.1, AL359878.1, AP002478.1, C2orf48, C10orf91, CLLU1, CLRN1-AS1, ERVMER61-1, and WARS2-IT1) in the ceRNA regulatory network were significantly associated with HCC-patient prognoses.

CONCLUSION

This study identified potential prognostic genes and constructed an lncRNA- miRNA-mRNA ceRNA regulatory network of HCC, which not only has important clinical significance for early diagnoses but also provides effective targets for HCC treatments and could provide new insights for HCC-interventional strategies.