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

Transcriptomic Analysis and Finding of Potential Key mRNA Expression Profile in Human Cumulus Cells During in Vitro Culture and Different Passages Based on Integrated Bioinformatics Analysis

This study leveraged microarray datasets to investigate differentially expressed genes (DEGs) in cumulus cells and their relevance in predicting the successful implantation of embryos in human in-vitro fertilization procedures. The microarray data were obtained from the GEO database, encompassing samples of cumulus cells during in vitro culture and different passages. To ensure data consistency, inter-batch normalization was performed, and Principal Component Analysis (PCA) was applied to assess the impact of normalization on sample group clustering. The integrated dataset included samples from cumulus cells during in vitro culture, comprising 17,662 genes. Utilizing the "limma" software package, 1906 DEGs were identified, with 437 genes downregulated and 589 genes upregulated in the cumulus cells of infertility cases, while 748 genes were upregulated, and 1317 genes were downregulated in cumulus cells of successful implantation cases. Functional enrichment analysis utilized Gene Ontology, Metascape, and Kyoto Encyclopedia of Genes and Genomes pathway enrichment tools. Biological processes and molecular functions were enriched, including protein targeting, mRNA processing, and molecular binding among the identified DEGs. Furthermore, target prediction and functional enrichment analysis of microRNAs (miRNAs) revealed 25 key genes and 13 relevant miRNAs were identified. Notably, hsa-miR-149, hsa-miR23b, hsa-miR-877, hsa-miR593, hsa-miR-18a, hsa-miR25, hsa-miR185, mmu-miR-207, hsa-miR425, hsa-miR214, hsa-miR-129, hsa-miR-629, and hsa-miR-194 emerged as the most prominent miRNAs with potential regulatory roles in successful embryo implantation. This comprehensive analysis provides valuable insights into the molecular mechanisms underlying embryo implantation, offering potential targets for further research and therapeutic interventions in assisted reproductive technologies.

RNA-Interference-Mediated miR-122-Based Gene Regulation in Colon Cancer, a Structural In Silico Analysis

The role of microRNA 122 (miR-122) in colorectal cancer (CRC) has not been widely investigated. In the current study, we aimed to identify the prominent gene and protein interactors of miR122 in CRC. Based on their binding affinity, these targets were chosen as candidate genes for the creation of miR122-mRNA duplexes. Following this, we examined the miRNA-mediated silencing mechanism using the gene-silencing complex protein Argonaute (AGO). Public databases, STRING, and GeneMANIA were utilized to identify major proteins and genes interacting with miR-122. DAVID, PANTHER, UniProt, FunRich, miRwalk, and KEGG were used for functional annotation, pathway enrichment, binding affinity analysis, and expression of genes in different stages of cancer. Three-dimensional duplexes of hub genes and miR-122 were created using the RNA composer, followed by molecular interaction analysis using molecular docking with the AGO protein. We analyzed, classified, and scrutinized 93 miR-122 interactors using various bioinformatic approaches. A total of 14 hub genes were categorized as major interactors of miR-122. The study confirmed the role of various experimentally documented miR-122 interactors such as MTDH (Q86UE4), AKT1 (P31749), PTPN1 (P18031), MYC (P01106), GSK3B (P49841), RHOA (P61586), and PIK3CG (P48736) and put forth several novel interactors, with AKT3 (Q9Y243), NCOR2 (Q9Y618), PIK3R2 (O00459), SMAD4 (P61586), and TGFBR1 (P36897). Double-stranded RNA duplexes of the strongest interactors were found to exhibit higher binding affinity with AGO. In conclusions, the study has explored the role of miR-122 in CRC and has identified a closely related group of genes influencing the prognosis of CRC in multiple ways. Further, these genes prove to be targets of gene silencing through RNA interference and might serve as effective therapeutic targets in understanding and treating CRC.

Identification of hub genes and transcription factor regulatory network for heart failure using RNA-seq data and robust rank aggregation analysis

Background

Heart failure (HF) is the end stage of various cardiovascular diseases with a high mortality rate. Novel diagnostic and therapeutic biomarkers for HF are urgently required. Our research aims to identify HF-related hub genes and regulatory networks using bioinformatics and validation assays.

Methods

Using four RNA-seq datasets in the Gene Expression Omnibus (GEO) database, we screened differentially expressed genes (DEGs) of HF using Removal of Unwanted Variation from RNA-seq data (RUVSeq) and the robust rank aggregation (RRA) method. Then, hub genes were recognized using the STRING database and Cytoscape software with cytoHubba plug-in. Furthermore, reliable hub genes were validated by the GEO microarray datasets and quantitative reverse transcription polymerase chain reaction (qRT-PCR) using heart tissues from patients with HF and non-failing donors (NFDs). In addition, R packages “clusterProfiler” and “GSVA” were utilized for enrichment analysis. Moreover, the transcription factor (TF)–DEG regulatory network was constructed by Cytoscape and verified in a microarray dataset.

Results

A total of 201 robust DEGs were identified in patients with HF and NFDs. STRING and Cytoscape analysis recognized six hub genes, among which ASPN, COL1A1, and FMOD were confirmed as reliable hub genes through microarray datasets and qRT-PCR validation. Functional analysis showed that the DEGs and hub genes were enriched in T-cell-mediated immune response and myocardial glucose metabolism, which were closely associated with myocardial fibrosis. In addition, the TF–DEG regulatory network was constructed, and 13 significant TF–DEG pairs were finally identified.

Conclusion

Our study integrated different RNA-seq datasets using RUVSeq and the RRA method and identified ASPN, COL1A1, and FMOD as potential diagnostic biomarkers for HF. The results provide new insights into the underlying mechanisms and effective treatments of HF.

Modular network inference between miRNA-mRNA expression profiles using weighted co-expression network analysis

Connecting transcriptional and post-transcriptional regulatory networks solves an important puzzle in the elucidation of gene regulatory mechanisms. To decipher the complexity of these connections, we build co-expression network modules for mRNA as well as miRNA expression profiles of breast cancer data. We construct gene and miRNA co-expression modules using the weighted gene co-expression network analysis (WGCNA) method and establish the significance of these modules (Genes/miRNAs) for cancer phenotype. This work also infers an interaction network between the genes of the turquoise module from mRNA expression data and hubs of the turquoise module from miRNA expression data. A pathway enrichment analysis using a miRsystem web tool for miRNA hubs and some of their targets, reveal their enrichment in several important pathways associated with the progression of cancer.

A Minimal Subset of Seven Genes Associated with Tumor Hepatocyte Differentiation Predicts a Poor Prognosis in Human Hepatocellular Carcinoma

Simple Summary

Liver cancer is one of the most commonly diagnosed cancers worldwide and the fourth leading cause of cancer-related deaths. Hepatocellular carcinoma (HCC) accounts for at least 80% of all malignant liver primary tumors. A better characterization of molecular mechanisms underlying HCC onset and progression may lead to discover new therapeutic targets and biomarkers. In this study, we performed an integrative transcriptomics analysis to evaluate the clinical relevance of genes associated with hepatocyte differentiation in human HCC. The HepaRG cell line model was used to define a gene expression signature reflecting the status of tumor hepatocyte differentiation. This signature was able to stratify HCC patients into clinically relevant molecular subtypes. Then, a minimal subset of seven differentiation-associated genes was identified to predict a poor prognosis in several cancer datasets.

Abstract

Hepatocellular carcinoma (HCC) is a deadly cancer worldwide as a result of a frequent late diagnosis which limits the therapeutic options. Tumor progression in HCC is closely correlated with the dedifferentiation of hepatocytes, the main parenchymal cells in the liver. Here, we hypothesized that the expression level of genes reflecting the differentiation status of tumor hepatocytes could be clinically relevant in defining subsets of patients with different clinical outcomes. To test this hypothesis, an integrative transcriptomics approach was used to stratify a cohort of 139 HCC patients based on a gene expression signature established in vitro in the HepaRG cell line using well-controlled culture conditions recapitulating tumor hepatocyte differentiation. The HepaRG model was first validated by identifying a robust gene expression signature associated with hepatocyte differentiation and liver metabolism. In addition, the signature was able to distinguish specific developmental stages in mice. More importantly, the signature identified a subset of human HCC associated with a poor prognosis and cancer stem cell features. By using an independent HCC dataset (TCGA consortium), a minimal subset of seven differentiation-related genes was shown to predict a reduced overall survival, not only in patients with HCC but also in other types of cancers (e.g., kidney, pancreas, skin). In conclusion, the study identified a minimal subset of seven genes reflecting the differentiation status of tumor hepatocytes and clinically relevant for predicting the prognosis of HCC patients.

Identification of Five Hub Genes as Key Prognostic Biomarkers in Liver Cancer via Integrated Bioinformatics Analysis

Liver cancer is one of the most common cancers and the top leading cause of cancer death globally. However, the molecular mechanisms of liver tumorigenesis and progression remain unclear. In the current study, we investigated the hub genes and the potential molecular pathways through which these genes contribute to liver cancer onset and development. The weighted gene co-expression network analysis (WCGNA) was performed on the main data attained from the GEO (Gene Expression Omnibus) database. The Cancer Genome Atlas (TCGA) dataset was used to evaluate the association between prognosis and these hub genes. The expression of genes from the black module was found to be significantly related to liver cancer. Based on the results of protein-protein interaction, gene co-expression network, and survival analyses, DNA topoisomerase II alpha (TOP2A), ribonucleotide reductase regulatory subunit M2 (RRM2), never in mitosis-related kinase 2 (NEK2), cyclin-dependent kinase 1 (CDK1), and cyclin B1 (CCNB1) were identified as the hub genes. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses showed that the differentially expressed genes (DEGs) were enriched in the immune-associated pathways. These hub genes were further screened and validated using statistical and functional analyses. Additionally, the TOP2A, RRM2, NEK2, CDK1, and CCNB1 proteins were overexpressed in tumor liver tissues as compared to normal liver tissues according to the Human Protein Atlas database and previous studies. Our results suggest the potential use of TOP2A, RRM2, NEK2, CDK1, and CCNB1 as prognostic biomarkers in liver cancer.

Mechanism-Centric Approaches for Biomarker Detection and Precision Therapeutics in Cancer

Biomarker discovery is at the heart of personalized treatment planning and cancer precision therapeutics, encompassing disease classification and prognosis, prediction of treatment response, and therapeutic targeting. However, many biomarkers represent passenger rather than driver alterations, limiting their utilization as functional units for therapeutic targeting. We suggest that identification of driver biomarkers through mechanism-centric approaches, which take into account upstream and downstream regulatory mechanisms, is fundamental to the discovery of functionally meaningful markers. Here, we examine computational approaches that identify mechanism-centric biomarkers elucidated from gene co-expression networks, regulatory networks (e.g., transcriptional regulation), protein-protein interaction (PPI) networks, and molecular pathways. We discuss their objectives, advantages over gene-centric approaches, and known limitations. Future directions highlight the importance of input and model interpretability, method and data integration, and the role of recently introduced technological advantages, such as single-cell sequencing, which are central for effective biomarker discovery and time-cautious precision therapeutics.

A gene module identification algorithm and its applications to identify gene modules and key genes of hepatocellular carcinoma

To further improve the effect of gene modules identification, combining the Newman algorithm in community detection and K-means algorithm framework, a new method of gene module identification, GCNA-Kpca algorithm, was proposed. The core idea of the algorithm was to build a gene co-expression network (GCN) based on gene expression data firstly; Then the Newman algorithm was used to initially identify gene modules based on the topology of GCN, and the number of clusters and clustering centers were determined; Finally the number of clusters and clustering centers were input into the K-means algorithm framework, and the secondary clustering was performed based on the gene expression profile to obtain the final gene modules. The algorithm took into account the role of modularity in the clustering process, and could find the optimal membership module for each gene through multiple iterations. Experimental results showed that the algorithm proposed in this paper had the best performance in error rate, biological significance and CNN classification indicators (Precision, Recall and F-score). The gene module obtained by GCNA-Kpca was used for the task of key gene identification, and these key genes had the highest prognostic significance. Moreover, GCNA-Kpca algorithm was used to identify 10 key genes in hepatocellular carcinoma (HCC): CDC20, CCNB1, EIF4A3, H2AFX, NOP56, RFC4, NOP58, AURKA, PCNA, and FEN1. According to the validation, it was reasonable to speculate that these 10 key genes could be biomarkers for HCC. And NOP56 and NOP58 are key genes for HCC that we discovered for the first time.

Identification of hub genes and biological pathways in hepatocellular carcinoma by integrated bioinformatics analysis

Background

Hepatocellular carcinoma (HCC), the main type of liver cancer in human, is one of the most prevalent and deadly malignancies in the world. The present study aimed to identify hub genes and key biological pathways by integrated bioinformatics analysis.

Methods

A bioinformatics pipeline based on gene co-expression network (GCN) analysis was built to analyze the gene expression profile of HCC. Firstly, differentially expressed genes (DEGs) were identified and a GCN was constructed with Pearson correlation analysis. Then, the gene modules were identified with 3 different community detection algorithms, and the correlation analysis between gene modules and clinical indicators was performed. Moreover, we used the Search Tool for the Retrieval of Interacting Genes (STRING) database to construct a protein protein interaction (PPI) network of the key gene module, and we identified the hub genes using nine topology analysis algorithms based on this PPI network. Further, we used the Oncomine analysis, survival analysis, GEO data set and random forest algorithm to verify the important roles of hub genes in HCC. Lastly, we explored the methylation changes of hub genes using another GEO data (GSE73003).

Results

Firstly, among the expression profiles, 4,130 up-regulated genes and 471 down-regulated genes were identified. Next, the multi-level algorithm which had the highest modularity divided the GCN into nine gene modules. Also, a key gene module (m1) was identified. The biological processes of GO enrichment of m1 mainly included the processes of mitosis and meiosis and the functions of catalytic and exodeoxyribonuclease activity. Besides, these genes were enriched in the cell cycle and mitotic pathway. Furthermore, we identified 11 hub genes, MCM3, TRMT6, AURKA, CDC20, TOP2A, ECT2, TK1, MCM2, FEN1, NCAPD2 and KPNA2 which played key roles in HCC. The results of multiple verification methods indicated that the 11 hub genes had highly diagnostic efficiencies to distinguish tumors from normal tissues. Lastly, the methylation changes of gene CDC20, TOP2A, TK1, FEN1 in HCC samples had statistical significance (P-value < 0.05).

Conclusion

MCM3, TRMT6, AURKA, CDC20, TOP2A, ECT2, TK1, MCM2, FEN1, NCAPD2 and KPNA2 could be potential biomarkers or therapeutic targets for HCC. Meanwhile, the metabolic pathway, the cell cycle and mitotic pathway might played vital roles in the progression of HCC.