Five Core Genes Related to the Progression and Prognosis of Hepatocellular Carcinoma Identified by Analysis of a Coexpression Network

Identification of RFC4 as a potential biomarker for pan-cancer involving prognosis, tumour immune microenvironment and drugs

RFC4 is required for DNA polymerase δ and DNA polymerase ε to initiate DNA template expansion. Downregulated RFC4 inhibits tumour proliferation by causing S-phase arrest and inhibiting mitosis, resulting in the reduction of tumour cells. RFC4 has been implicated that it plays an important role in the initiation and progression of cancers, but a comprehensive analysis of the role of RFC4 in cancer has not been performed. We comprehensively analysed the expression, prognosis, methylation level, splicing level, relationship of RFC4 and immune infiltration, and pan-cancer immunotherapy response used various databases (including TCGA, GTEx, UALCAN, Oncosplicing, TIDE, TISCH, HPA and CAMOIP), and experimented its biological function in HCC. Through pan-cancer analysis, we found that RFC4 is significantly upregulated in most tumours. The tumour patients with high expression of RFC4 have poor prognosis. The methylation level and variable splicing level of RFC4 were abnormal in most tumours compared with the adjacent tissues. Furthermore, RFC4 was closely associated with immune cell infiltration in various cancers. RFC4 was significantly co-expressed with immune checkpoints and other immune-related genes. The expression of RFC4 could indicate the immunotherapy efficacy of some tumours. The RFC4 expression was associated with sensitivity to specific small molecule drugs. Cell experiments have shown that downregulated RFC4 can inhibit cell cycle and tumour cell proliferation. We conducted a systematic pan-cancer analysis of RFC4, and the results showed that RFC4 can serve as a biomarker for cancer diagnosis and prognosis. These findings open new perspectives for precision medicine.

Inhibition of CDK9 exhibits anticancer activity in hepatocellular carcinoma cells via targeting ribonucleotide reductase

Cyclin-dependent kinase 9 (CDK9) inhibitors are a novel category of anticancer treatment for cancers. However, their effects on hepatocellular carcinoma (HCC) are rarely investigated. Human ribonucleotide reductase (RR, which consists of RRM1 and RRM2 subunits) catalyzes the conversion of ribonucleoside diphosphate into 2'-deoxyribonucleoside diphosphate to maintain the homeostasis of nucleotide pools, which play essential roles in DNA synthesis and DNA repair. In this study, we identified that CDK9 protein expression in adjacent non-tumor tissues predicted HCC patients' overall and progression-free survivals. The anticancer activity of a CDK9-selective inhibitor, LDC000067, on HCC cells was positively associated with its ability to inhibit the expression of RRM1 and RRM2. LDC000067 downregulated RRM1 and RRM2 expression through post-transcriptional pathway. Specifically, LDC000067 triggered RRM2 protein degradation via multiple pathways, including proteasome-, lysosome-, and calcium-dependent pathways. Furthermore, CDK9 positively correlates with RRM1 or RRM2 expression in HCC patients, and the expressions of these three genes were associated with the higher infiltration of immune cells in HCC. Taken together, this study identified the prognostic relevance of CDK9 in HCC and the molecular mechanism for the anticancer effect of CDK9 inhibitors on HCC.

Screening of Prognostic Markers for Hepatocellular Carcinoma Patients Based on Multichip Combined Analysis

Methods

GSE (14520, 36376, 57957, 76427) datasets were accessed from GEO database. 55 differential mRNAs (DEGs) were obtained by differential analysis based on the datasets. GO and KEGG analysis results indicated that the DEGs were enriched in xenobiotic metabolic process and other pathways. Expression profiles and clinical data of TCGA-LIHC mRNAs were from TCGA database. We established a prognostic model of HCC through univariate and multivariate Cox risk regression analyses. ROC curve analysis was used to examine the prognostic model performance. GSEA analysis was performed between the high- and low-risk score sample groups.

Results

A 4-gene HCC prognostic model was constructed, in which the gene expressions correlated to HCC patients' survival. The AUC value presented 0.734 in the ROC analysis for the prognostic model.

Conclusion

The four-gene model could be introduced as an independent prognostic factors to assess HCC patients' survival status.

High Expression of RRM1 Mediated by ncRNAs Correlates with Poor Prognosis and Tumor Immune Infiltration of Hepatocellular Carcinoma

Introduction

Hepatocellular carcinoma (HCC) is one of several tumors with poor prognosis and causes a significant social burden. A growing number of studies have shown that RRM1 plays a crucial role in the development and progression of multiple human cancers. However, the specific role and mechanism of RRM1 have not been fully defined in HCC.

Methods

TCGA and GTEx data were used for the first time to conduct a pan-cancer analysis of RRM1 expression and prognosis, and identified RRM1 as a possible potential oncogene in HCC. At the same time, a combination of analyses (including expression analysis, correlation analysis or survival analysis) identified non-coding RNAs (ncRNAs) that contribute to RRM1 overexpression.

Results

MIR4435-2HG/miR-22-3p and SNHG6/miR-101-3p were identified as the most promising RRM1 upstream ncRNA-related pathways in HCC. In addition, RRM1 levels were significantly and positively correlated with tumor immune cell infiltration, immune cell biomarker or immune checkpoint expression.

Conclusion

These results suggest that high expression of RRM1 mediated by ncRNAs is associated with poor prognosis and tumor immune infiltration in HCC.

A Novel Four-Gene Signature as a Potential Prognostic Biomarker for Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a common malignant tumor with high incidence and mortality rates. However, a reliable prognostic signature has not yet been confirmed. Thus, the purpose of the present study was to develop a biomarker with high specificity and sensitivity for the diagnosis and prognosis of patients with HCC. The mRNA expression profiles of HCC were obtained from the GSE19665, GSE41804, and TCGA databases. Subsequently, 193 differentially expressed genes (DEGs) were identified from the intersection of the data from the three datasets. Bioinformatics analysis showed that the identified DEGs are related to the cell cycle, oocyte meiosis, and p53 signaling pathway, among other factors, in cancers. A protein-protein interaction (PPI) and a functional analysis were performed to investigate the biological function of the DEGs and obtain the candidate genes using the MCODE of Cytoscape. The candidate genes were introduced into the TCGA database for survival analysis, and the four candidate genes that were hub genes and meaningful for survival were retained for further verification. We validated the gene and protein expression and determined the prognosis of our patient cohort. In addition, we evaluated the biological functions regulating tumor cell proliferation and metastasis in vitro. According to the ROC curve analysis of gene expression in clinical samples, it was found that the four genes can be used to predict the diagnosis. A survival analysis based on data from the TCGA database and clinical samples showed that the four genes may be used as biomarkers for providing prognoses for patients. The cell functional experiments revealed that these four genes were related to tumor proliferation, migration, and invasion. In conclusion, the genes identified in the present study could be used as markers to diagnose and predict the prognosis of patients with HCC and guide targeted therapy.

Screening Hub Genes of Hepatocellular Carcinoma Based on Public Databases

Tumor recurrence and metastasis often occur in HCC patients after surgery, and the prognosis is not optimistic. Hence, searching effective biomarkers for prognosis of is of great importance. Firstly, HCC-related data was acquired from the TCGA and GEO databases. Based on GEO data, 256 differentially expressed genes (DEGs) were obtained firstly. Subsequently, to clarify function of DEGs, clusterProfiler package was used to conduct functional enrichment analyses on DEGs. Protein-protein interaction (PPI) network analysis screened 20 key genes. The key genes were filtered via GEPIA database, by which 11 hub genes (F9, CYP3A4, ASPM, AURKA, CDC20, CDCA5, NCAP, PRC1, PTTG1, TOP2A, and KIFC1) were screened out. Then, univariate Cox analysis was applied to construct a prognostic model, followed by a prediction performance validation. With the risk score calculated by the model and common clinical features, univariate and multivariate analyses were carried out to assess whether the prognostic model could be used independently for prognostic prediction. In conclusion, the current study screened HCC prognostic gene signature based on public databases.

Identification of Atrial Fibrillation-Associated Genes ERBB2 and MYPN Using Genome-Wide Association and Transcriptome Expression Profile Data on Left-Right Atrial Appendages

More reliable methods are needed to uncover novel biomarkers associated with atrial fibrillation (AF). Our objective is to identify significant network modules and newly AF-associated genes by integrative genetic analysis approaches. The single nucleotide polymorphisms with nominal relevance significance from the AF-associated genome-wide association study (GWAS) data were converted into the GWAS discovery set using ProxyGeneLD, followed by merging with significant network modules constructed by weighted gene coexpression network analysis (WGCNA) from one expression profile data set, composed of left and right atrial appendages (LAA and RAA). In LAA, two distinct network modules were identified (blue: p = 0.0076; yellow: p = 0.023). Five AF-associated biomarkers were identified (ERBB2, HERC4, MYH7, MYPN, and PBXIP1), combined with the GWAS test set. In RAA, three distinct network modules were identified and only one AF-associated gene LOXL1 was determined. Using human LAA tissues by real-time quantitative polymerase chain reaction, the differentially expressive results of ERBB2, MYH7, and MYPN were observed (p < 0.05). This study first demonstrated the feasibility of fusing GWAS with expression profile data by ProxyGeneLD and WGCNA to explore AF-associated genes. In particular, two newly identified genes ERBB2 and MYPN via this approach contribute to further understanding the occurrence and development of AF, thereby offering preliminary data for subsequent studies.

Replication Factor C4 in human hepatocellular carcinoma: A potent prognostic factor associated with cell proliferation

Replication Factor c4 (RFC4) has been found to play important roles in many carcinomas and is correlated with poor prognosis. The present study was performed to investigate the specific role of RFC4 in hepatocellular carcinoma (HCC) and the underlying molecular mechanism. Public datasets including TCGA and GTEx were applied to explore the expression of RFC4 in HCC and its association with HCC prognosis. The results of bioinformatics analysis showed that RFC4 was overexpressed in HCC tissues compared with noncancerous tissues and significantly correlated with poor prognosis for HCC. Through immunohistochemistry, the association between RFC4 expression and clinical-pathological features of HCC patients was evaluated. Western blots were applied to investigate relative protein expression. Then in vivo and in vitro experiments were performed to explore the function of RFC4 in HCC tumor cells. The present results suggest that high level expression of RFC4 is associated with tumor size. In addition, RFC4 knockdown suppressed the cell proliferation and sphere formation of hepatoma cells in vitro. Moreover, silencing of RFC4 significantly decreased the growth of tumors in a xenograft tumor model. In conclusion, our study indicates that RFC4 is a potential prognostic predictor associated with poor outcomes for HCC patients. Furthermore, knocking down RFC4 could significantly inhibit tumor progression both in vitro and in vivo. Therefore, the present study can shed new light on the understanding of molecular mechanisms of HCC and may provide molecular targets and diagnostic biomarkers for the treatment of HCC.

Identification of a five-gene signature in association with overall survival for hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is considered to be a malignant tumor with a high incidence and a high mortality. Accurate prognostic models are urgently needed. The present study was aimed at screening the critical genes for prognosis of HCC.

Methods

The GSE25097, GSE14520, GSE36376 and GSE76427 datasets were obtained from Gene Expression Omnibus (GEO). We used GEO2R to screen differentially expressed genes (DEGs). A protein-protein interaction network of the DEGs was constructed by Cytoscape in order to find hub genes by module analysis. The Metascape was performed to discover biological functions and pathway enrichment of DEGs. MCODE components were calculated to construct a module complex of DEGs. Then, gene set enrichment analysis (GSEA) was used for gene enrichment analysis. ONCOMINE was employed to assess the mRNA expression levels of key genes in HCC, and the survival analysis was conducted using the array from The Cancer Genome Atlas (TCGA) of HCC. Then, the LASSO Cox regression model was performed to establish and identify the prognostic gene signature. We validated the prognostic value of the gene signature in the TCGA cohort.

Results

We screened out 10 hub genes which were all up-regulated in HCC tissue. They mainly enrich in mitotic cell cycle process. The GSEA results showed that these data sets had good enrichment score and significance in the cell cycle pathway. Each candidate gene may be an indicator of prognostic factors in the development of HCC. However, hub genes expression was weekly associated with overall survival in HCC patients. LASSO Cox regression analysis validated a five-gene signature (including CDC20, CCNB2, NCAPG, ASPM and NUSAP1). These results suggest that five-gene signature model may provide clues for clinical prognostic biomarker of HCC.

Identification of Core Genes Related to Progression and Prognosis of Hepatocellular Carcinoma and Small-Molecule Drug Predication

Background

Hepatocellular carcinoma (HCC) is one of the most leading causes of cancer death with a poor prognosis. However, the underlying molecular mechanisms are largely unclear, and effective treatment for it is limited. Using an integrated bioinformatics method, the present study aimed to identify the key candidate prognostic genes that are involved in HCC development and identify small-molecule drugs with treatment potential.

Methods and Results

In this study, by using three expression profile datasets from Gene Expression Omnibus database, 1,704 differentially expressed genes were identified, including 671 upregulated and 1,033 downregulated genes. Then, weighted co-expression network analysis revealed nine modules are related with pathological stage; turquoise module was the most associated module. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway analyses (KEGG) indicated that these genes were enriched in cell division, cell cycle, and metabolic related pathways. Furthermore, by analyzing the turquoise module, 22 genes were identified as hub genes. Based on HCC data from gene expression profiling interactive analysis (GEPIA) database, nine genes associated with progression and prognosis of HCC were screened, including ANLN, BIRC5, BUB1B, CDC20, CDCA5, CDK1, NCAPG, NEK2, and TOP2A. According to the Human Protein Atlas and the Oncomine database, these genes were highly upregulated in HCC tumor samples. Moreover, multivariate Cox regression analysis showed that the risk score based on the gene expression signature of these nine genes was an independent prognostic factor for overall survival and disease-free survival in HCC patients. In addition, the candidate small-molecule drugs for HCC were identified by the CMap database.

Conclusion

In conclusion, the nine key gene signatures related to HCC progression and prognosis were identified and validated. The cell cycle pathway was the core pathway enriched with these key genes. Moreover, several candidate molecule drugs were identified, providing insights into novel therapeutic approaches for HCC.

MYCT1 Inhibits the Adhesion and Migration of Laryngeal Cancer Cells Potentially Through Repressing Collagen VI

MYCT1, a target of c-Myc, inhibits laryngeal cancer cell migration, but the underlying mechanism remains unclear. In the study, we detected differentially expressed genes (DEGs) from laryngeal cancer cells transfected by MYCT1 using RNA-seq (GSE123275). DEGs from head and neck squamous cell carcinoma (HNSCC) were first screened by comparison of transcription data from the Gene Expression Omnibus (GSE6631) and the Cancer Genome Atlas (TCGA) datasets using weighted gene co-expression network analysis (WGCNA). GO and KEGG pathway analysis explained the functions of the DEGs. The DEGs overlapped between GSE6631and TCGA datasets were then compared with ours to find the key DEGs downstream of MYCT1 related to the adhesion and migration of laryngeal cancer cells. qRT-PCR and Western blot were applied to validate gene expression at mRNA and protein levels, respectively. Finally, the cell adhesion, migration, and wound healing assays were to check cell adhesion and migration abilities, respectively. As results, 39 overlapping genes were enriched in the GSE6631 and TCGA datasets, and most of them revealed adhesion function. Thirteen of 39 genes including COL6 members COL6A1, COL6A2, and COL6A3 were overlapped in GSE6631, TCGA, and GSE123275 datasets. Similar to our RNA-seq results, we confirmed that COL6 is a target of MYCT1 in laryngeal cancer cells. We also found that MYCT1 inhibited the adhesion and migration of laryngeal cancer cells via COL6. These indicate that COL6 is a potential target of MYCT1 and participates the adhesion and migration of laryngeal cancer cells, which provides an important clue for further study on how MYCT1 regulating COL6 in laryngeal cancer progression.

Identification of DNA repair-related genes predicting pathogenesis and prognosis for liver cancer

Background

Liver cancer (LC) is one of the most fatal cancers throughout the world. More efficient and sensitive gene signatures that could accurately predict survival in LC patients are vitally needed to promote a better individualized and effective treatment.

Material/methods

422 LC and adjacent normal tissues with both RNA-Seq and clinical data in TCGA were embedded in our study. Gene set enrichment analysis (GSEA) was applied to identify genes and hallmark gene sets that are more valuable for liver cancer therapy. Cox regression analysis was used to identify genes related to overall survival (OS) and build the prediction model. cBioPortal database was used to examine the alterations of the panel mRNA signature. ROC curves and Kaplan–Meier curves were used to validate the prediction model. Besides, the expression of the genes in the model were validated using quantitative real-time PCR in clinical tissue specimens.

Results

The panel of DNA repair-related mRNA signature consisted of seven mRNAs: RFC4 (replication factor C subunit 4), ZWINT (ZW10 interacting kinetochore protein), UPF3B (UPF3B regulator of nonsense mediated mRNA decay), NCBP2 (nuclear cap binding protein subunit 2), ADA (adenosine deaminase), SF3A3 (splicing factor 3a subunit 3) and GTF2H1 (general transcription factor IIH subunit 1). On-line analysis of cBioPortal database found that the expression of the panel mRNA has a wide variation ranging from 7 to 10%. All the mRNAs were significantly upregulated in LC tissues compared to normal tissues (P < 0.05). The risk model is closely related to the OS of LC patients. The hazard ratio (HR) is 2.184 [95% CI (confidence interval) 1.523–3.132] and log-rank P-value < 0.0001. For clinical specimen validation, we found that all of the genes in the model upregulated in liver cancer tissues versus normal liver tissues, which was consistent with the results predicted.

Conclusions

Our study demonstrated a mRNA signature including seven mRNA for prognosis prediction of LC. This panel gene signature provides a new criterion for accurate diagnosis and therapeutic target of LC.

The Clinical Significance and Potential Molecular Mechanism of PTTG1 in Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is the major histological type of esophageal cancers worldwide. Transcription factor PTTG1 was seen highly expressed in a variety of tumors and was related to the degree of tumor differentiation, invasion, and metastasis. However, the clinical significance of PTTG1 had yet to be verified, and the mechanism of abnormal PTTG1 expression in ESCC was not clear. In this study, the comprehensive analysis and evaluation of PTTG1 expression in ESCC were completed by synthesizing in-house immunohistochemistry (IHC), clinical sample tissue RNA-seq (in-house RNA-seq), public high-throughput data, and literature data. We also explored the possible signaling pathways and target genes of PTTG1 in ESCC by combining the target genes of PTTG1 (displayed by ChIP-seq), differentially expressed genes (DEGs) of ESCC, and PTTG1-related genes, revealing the potential molecular mechanism of PTTG1 in ESCC. In the present study, PTTG1 protein and mRNA expression levels in ESCC tissues were all significantly higher than in non-cancerous tissues. The pool standard mean difference (SMD) of the overall PTTG1 expression was 1.17 (95% CI: 0.72-1.62, P < 0.01), and the area under curve (AUC) of the summary receiver operating characteristic (SROC) was 0.86 (95% CI: 0.83-0.89). By combining the target genes displayed by ChIP-seq of PTTG1, DEGs of ESCC, and PTTG1-related genes, it was observed that PTTG1 may interact with these genes through chemokines and cytokine signaling pathways. By constructing a protein-protein interaction (PPI) network and combining ChIP-seq data, we obtained four PTTG1 potential target genes, SPTAN1, SLC25A17, IKBKB, and ERH. The gene expression of PTTG1 had a strong positive correlation with SLC25A17 and ERH, which suggested that PTTG1 might positively regulate the expression of these two genes. In summary, the high expression of PTTG1 may play an important role in the formation of ESCC. These roles may be completed by PTTG1 regulating the downstream target genes SLC25A17 and ERH.

Identification of Hub Genes Associated With Hepatocellular Carcinoma Using Robust Rank Aggregation Combined With Weighted Gene Co-expression Network Analysis

Background

Bioinformatics provides a valuable tool to explore the molecular mechanisms underlying pathogenesis of hepatocellular carcinoma (HCC). To improve prognosis of patients, identification of robust biomarkers associated with the pathogenic pathways of HCC remains an urgent research priority.

Methods

We employed the Robust Rank Aggregation method to integrate nine qualified HCC datasets from the Gene Expression Omnibus. A robust set of differentially expressed genes (DEGs) between tumor and normal tissue samples were screened. Weighted gene co-expression network analysis was applied to cluster DEGs and the key modules related to clinical traits identified. Based on network topology analysis, novel risk genes derived from key modules were mined and biological verification performed. The potential functions of these risk genes were further explored with the aid of miRNA–mRNA regulatory networks. Finally, the prognostic ability of these genes was assessed by constructing a clinical prediction model.

Results

Two key modules showed significant association with clinical traits. In combination with protein–protein interaction analysis, 29 hub genes were identified. Among these genes, 19 from one module showed a pattern of upregulation in HCC and were associated with the tumor node metastasis stage, and 10 from the other module displayed the opposite trend. Survival analyses indicated that all these genes were significantly related to patient prognosis. Based on the miRNA-mRNA regulatory network, 29 genes strongly linked to tumor activity were identified. Notably, five of the novel risk genes, ABAT, DAO, PCK2, SLC27A2, and HAO1, have rarely been reported in previous studies. Gene set enrichment analysis for each gene revealed regulatory roles in proliferation and prognosis of HCC. Least absolute shrinkage and selection operator regression analysis further validated DAO, PCK2, and HAO1 as prognostic factors in an external HCC dataset.

Conclusion

Analysis of multiple datasets combined with global network information presents a successful approach to uncover the complex biological mechanisms of HCC. More importantly, this novel integrated strategy facilitates identification of risk hub genes as candidate biomarkers for HCC, which could effectively guide clinical treatments.

A Key mRNA-miRNA-lncRNA Competing Endogenous RNA Triple Sub-network Linked to Diagnosis and Prognosis of Hepatocellular Carcinoma

Growing evidence has illustrated critical roles of competing endogenous RNA (ceRNA) regulatory network in human cancers including hepatocellular carcinoma. In this study, we aimed to find promising diagnostic and prognostic biomarkers for patients with hepatocellular carcinoma. Three novel unfavorable prognosis-associated genes (CELSR3, GPSM2, and CHEK1) was first identified. We also demonstrated that these genes were significantly upregulated in hepatocellular carcinoma cell lines and tissues. Next, 154 potential miRNAs of CELSR3, GPSM2, and CHEK1 were predicted. CHEK1-hsa-mir-195-5p/hsa-mir-497-5p and GPSM2-hsa-mir-122-5p axes were defined as two key pathways in carcinogenesis of hepatocellular carcinoma by combination of in silico analysis and experimental validation. Subsequently, lncRNAs binding to hsa-mir-195-5p, hsa-mir-497-5p, and hsa-mir-122-5p were predicted via starBase and miRNet databases. After performing expression analysis and survival analysis for these predicted lncRNAs, we showed that nine lncRNAs (SNHG1, SNHG12, LINC00511, HCG18, FGD5-AS1, CERS6-AS1, NUTM2A-AS1, SNHG16, and ASB16-AS1) were markedly increased in hepatocellular carcinoma and their upregulation indicated poor prognosis. Moreover, a similar mRNA-miRNA-lncRNA analysis for six “known” genes (CLEC3B, DNASE1L3, PTTG1, KIF2C, XPO5, and UBE2S) was performed. Subsequently, a comprehensive mRNA-miRNA-lncRNA triple ceRNA network linked to prognosis of patients with hepatocellular carcinoma was established. Moreover, all RNAs in this network exhibited significantly diagnostic values for patients with hepatocellular carcinoma. In summary, the current study constructed a mRNA-miRNA-lncRNA ceRNA network associated with diagnosis and prognosis of hepatocellular carcinoma.