Bioinformatics analysis of constructing a HCV-related hepatocellular carcinoma miRNA-mRNA regulation network.
Ontology highlight
ABSTRACT: Abstract
Hepatocellular carcinoma (HCC) is one of the tumors with a higher mortality rate globally, which significantly threatens people's health. Hepatitis C virus (HCV) infection is a major driving factor of HCC. This study aims to determine the key microRNA (miRNA), hub genes, and related pathways, construct potential miRNA-mRNA regulatory networks, and clarify the new molecular mechanism of HCV-related HCC. In this study, 16 differentially expressed miRNAs (DE miRNAs) were identified. The prediction of potential transcription factors and target genes not only found that SP1 and ERG1 may potentially regulate most of the screened DE miRNAs, but it also obtained 2923 and 1782 predicted target genes for the up-regulation and down-regulation of DE miRNAs, respectively. Subsequently, the introduction of differentially expressed genes dataset GSE62232 for target gene verification yielded 98 and 147 potential up-regulation and down-regulation target genes. The gene ontology (GO) and Kyoto encyclopedia of genes and genomes pathway enrichment analysis showed that they were mainly enriched in the cell cycle process, that is, subsequently, 20 hub genes were screened out through the protein-protein interaction network, and related genes were further evaluated using the GEPIA database. Based on the above analysis, the miRNA-hub gene regulatory network was constructed. In short, this research's hub genes and miRNAs closely related to HCV-related HCC were screened and identified through bioinformatics analysis and then built their connection. These results are expected to find potential therapeutic targets for HCV-related HCC.
SUBMITTER: Hao R
PROVIDER: S-EPMC8376384 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Aims: We aimed to explore the crucial miRNA-mRNA axis through bioinformatics analysis and provide evidences for the development of pathophysiological mechanisms and new therapies for HBV-related HCC. Methods: MiRNA (GSE76903) and mRNA (GSE77509) dataset were used to screen differentially expressed miRNAs (DE-miRNAs) and differentially expressed mRNAs (DE-mRNAs) using R software. Overlapping genes between DE-mRNAs and target genes of DE-miRNAs were identified as candidate genes. Hub genes were obtained via cytohubba analysis. The expression at protein and mRNA levels and prognostic value of hub genes were evaluated based on The Cancer Genome Atlas (TCGA) data. Key miRNA-mRNA axes were constructed according to predicted miRNA-mRNA pairs. MiRNA expression and prognostic role were respectively identified using starBase v3.0 and Kaplan-Meier plotter database. Real-time PCR was performed to verify the expression of crucial miRNAs and mRNAs. Coexpression of crucial miRNA and mRNA were analyzed using starBase v3.0. Results: CDK1, CCNB1, CKS2 and CCNE1 were screened as hub genes, which were significantly upregulated at protein and mRNA levels. These up-regulated hub genes were also significantly associated with poor prognosis. Hsa-mir-195-5p/CDK1, hsa-mir-5589-3p/CCNB1 and hsa-let-7c-3p/CKS2 were screened as critical miRNA-mRNA axes. Critical miRNAs were decreased in HCC, which indicates unfavourable prognosis. QPCR results showed that crucial miRNAs were decreased, whereas critical mRNAs were increased in HBV-related HCC. A reverse relationship between miRNA and mRNA in crucial axis was further verified. Conclusion: This study identified several miRNA-mRNA axes in HBV-related HCC. Hsa-mir-195-5p/CDK1, hsa-mir-5589-3p/CCNB1 and hsa-let-7c-3p/CKS2 might serve as potential prognostic biomarkers and therapeutic targets for HBV-related HCC.
| S-EPMC7757140 | biostudies-literature
Project description:Liver hepatocellular carcinoma (LIHC) is one of the major causes of cancer-related death worldwide with increasing incidences, however there are very few studies about the underlying mechanisms and pathways in the development of LIHC. We obtained LIHC samples from The Cancer Genome Atlas (TCGA) to screen differentially expressed mRNAs, lncRNAs, miRNAs and driver mutations. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, Gene ontology enrichment analyses and protein-protein interaction (PPI) network were performed. Moreover, we constructed a competing endogenous lncRNAs-miRNAs-mRNAs network. Finally, cox proportional hazards regression analysis was used to identify important prognostic differentially expressed genes. Total of 1284 mRNAs, 123 lncRNAs, 47 miRNAs were identified within different tissues of LIHC patients. GO analysis indicated that upregulated and downregulated differentially expressed mRNAs (DEmRNAs) were mainly associated with cell division, DNA replication, mitotic sister chromatid segregation and complement activation respectively. Meanwhile, KEGG terms revealed that upregulated and downregulated DEmRNAs were primarily involved in DNA replication, Metabolic pathways, cell cycle and Metabolic pathways, chemical carcinogenesis, retinol metabolism pathway respectively. Among the DERNAs, 542 lncRNAs-miRNAs-mRNAs pairs were predicted to construct a ceRNA regulatory network including 35 DElncRNAs, 26 DEmiRNAs and 112 DEmRNAs. In the Kaplan-Meier analysis, total of 43 mRNAs, 14 lncRNAs and 3 miRNAs were screened out to be significantly correlated with overall survival of LIHC. The mutation signatures were analyzed and its correlation with immune infiltrates were evaluated using the TIMER in LIHC. Among the mutation genes, TTN mutation is often associated with poor immune infiltration and a worse prognosis in LIHC. This work conducted a novel lncRNAs-miRNAs-mRNAs network and mutation signatures for finding potential molecular mechanisms underlying the development of LIHC. The biomarkers also can be used for predicting prognosis of LIHC.
| S-EPMC8051809 | biostudies-literature
Project description:AimsThis study aimed to construct a miRNA-mRNA network to elucidate the molecular mechanism of exosome function in metastatic HCC.MethodsWe explored the Gene Expression Omnibus (GEO) database and then analyzed the RNAs of 50 samples to obtain differentially expressed miRNAs (DEMs) and mRNAs (DEGs) involved in the progression of metastatic HCC. Next, a miRNA-mRNA network related to exosomes in metastatic HCC was constructed on the basis of the identified DEMs and DEGs. Finally, the function of the miRNA-mRNA network was explored by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Immunohistochemistry was performed to validate expression of NUCKS1 in HCC specimens. Based on the immunohistochemistry, the score of the NUCKS1 expression was calculated, and the patients were divided into high- and low-expression patients, and the differences in survival between the two groups were compared.ResultsThrough our analysis, 149 DEMs and 60 DEGs were identified. In addition, a miRNA-mRNA network, including 23 miRNAs and 14 mRNAs, was constructed. Low expression of NUCKS1 was validated in the majority of HCCs compared with their matched adjacent cirrhosis specimens (P < 0.001), which was consistent with our result of differential expression analyses. HCC patients with low expression of NUCKS1 had shorter overall survival than those with high NUCKS1 expression (P = 0.0441).ConclusionsThe novel miRNA-mRNA network will provide new insights into the underlying molecular mechanisms of exosomes in metastatic HCC. NUCKS1 might serve a potential therapeutic target to restrain the development of HCC.
| S-EPMC10123261 | biostudies-literature
Project description:Immune cells in the tumor microenvironment play a crucial role in regulating tumor progression. The circular RNA (circRNA) regulatory network involved in immune cell infiltration in hepatocellular carcinoma (HCC) remains largely unknown. In this study, the "estimate the proportion of immune and cancer cells" (EPIC) application is used to evaluate the fractions of immune cells, cancer-associated fibroblasts, and endothelial cells in HCC from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Patients with a high macrophage fraction have better overall survival, and macrophage fraction is an independent prognostic factor for HCC. Next, the common differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs), and circRNAs (DEcircRNAs) between paired tumor and non-tumor tissues are screened out from the TCGA and/or GEO databases. Through spearman correlation analysis, the macrophage-related DEmRNAs are identified to construct a circRNA-miRNA-mRNA regulatory network, which includes 6 DEcircRNAs, 7 DEmiRNAs, and 45 DEmRNAs. Functional enrichment analysis reveals that these DEmRNAs are mainly involved in immune-related processes. Furthermore, six hub DEmRNAs are identified to establish a hub circRNA regulatory network. Among the DEmRNAs in the network, PRC1 is identified as the most influential node. PRC1 high expression is correlated with poor prognosis and low macrophage infiltration in HCC. Taken together, we identify a certain circRNA regulatory network related to macrophage infiltration and provide novel insight into the mechanism of study and therapeutic targets for HCC.
| S-EPMC7500212 | biostudies-literature
Project description:Increasing studies have reported that circular RNAs (circRNAs) play critical roles in tumorigenesis and cancer progression. However, the underlying regulatory mechanisms of circRNA-related competing endogenous RNA (ceRNA) in liver hepatocellular carcinoma (LIHC) are still unclear. In the present study, we discovered dysregulated circRNAs through Gene Expression Omnibus (GEO) analysis and validated the expression of the top seven circRNAs with upregulated expression by qRT-PCR and Sanger sequencing. Then, the Cancer-Specific CircRNA Database (CSCD) was used to predict the downstream miRNAs of seven circRNAs, and expression and survival analyses through The Cancer Genome Atlas (TCGA) were performed to identify the key miRNA in LIHC. Thereafter, the hsa_circ_0017264-hsa-miR-195-5p subnetwork was successfully constructed. Subsequently, we predicted downstream target genes of hsa-miR-195-5p with TargetScan, miRDB, and mirtarbase and overlapped them with differentially expressed mRNAs to obtain 21 target genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed to predict the biological and functional roles of these target genes. Finally, with Pearson correlation and prognostic value analysis, a survival-related hsa_circ_0017264-hsa-miR-195-5p-CHEK1/CDC25A/FOXK1 axis was established. Gene set enrichment analysis (GSEA) was performed to determine the function of CHEK1/CDC25A/FOXK1 in the ceRNA network. Moreover, immune infiltration analysis revealed that the ceRNA network was markedly associated with the levels of multiple immune cell infiltrates, immune cell biomarkers and immune checkpoints. Overall, the hsa_circ_0017264-hsa-miR-195-5p-CHEK1/CDC25A/FOXK1 network might provide novel insights into the potential mechanisms underlying LIHC onset and progression.
| S-EPMC8924452 | biostudies-literature
Project description:Long non-coding RNAs (lncRNAs) are ncRNA transcripts >200 nucleotides that are important genetic regulators. LncRNAs can directly regulate mRNA through a lncRNA-mRNA regulatory mode and can also regulate mRNA through competitive binding to micro (mi)RNA, which is generally known as the competitive endogenous RNA (ceRNA) network. The present study evaluated the functional roles and regulatory networks of lncRNAs in chronic glomerulonephritis (CGN). The proliferative ability of mouse glomerular mesangial cells (GMCs) induced by different concentrations of lipopolysaccharide (LPS) was assessed using the Cell Counting Kit-8 assay, and RNA sequencing (RNA-seq) was performed to identify differentially expressed lncRNAs in LPS-induced GMCs. Based on the sequencing results, six lncRNAs were selected for validation using reverse transcription-quantitative PCR (RT-qPCR). Furthermore, the lncRNA-mRNA regulatory network and the lncRNA-miRNA-mRNA ceRNA network were constructed to assess the role and mechanism of CGN-related lncRNAs. To elucidate the biological functions of lncRNAs, Gene Ontology (GO) biological process term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on all mRNAs involved in the lncRNA-mRNA regulatory network and in the ceRNA network. A total of 1,532 differentially expressed lncRNAs, including 594 upregulated lncRNAs and 938 downregulated lncRNAs, were identified using RNA-seq. The results of RT-qPCR validation were consistent with RNA-seq results. An lncRNA-mRNA regulatory network, including 236 lncRNAs and 556 mRNAs, and a ceRNA network, including 6 lncRNAs, 18 miRNAs and 419 mRNAs, were successfully constructed. The GO biological process term enrichment and KEGG pathway enrichment analyses demonstrated that those lncRNAs were often related to inflammatory response and substance metabolism. The present study identified key CGN-related lncRNAs in LPS-induced GMCs, and further demonstrated a global view of the lncRNA-mRNA regulatory network and ceRNA network involved in CGN. These results offered novel insights into the roles of lncRNAs in the pathogenesis of CGN and identified potential diagnostic biomarkers.
| S-EPMC10375445 | biostudies-literature
Project description:This study was to identify important circRNA-miRNA-mRNA (ceRNAs) regulatory mechanisms in hepatocellular carcinoma (HCC). The circRNA dataset GSE97332 and miRNA dataset GSE57555 were used for analyses. Functional enrichment analysis for miRNA and target gene was conducted using cluster Profiler. Survival analysis was conducted through R package Survival. The ceRNAs and drug-gene interaction networks were constructed. The ceRNAs network contained five miRNAs including hsa-miR-25-3p, hsa-miR-3692-5p, hsa-miR-4270, hsa-miR-331-3p, and hsa-miR-125a-3p. Among the network, hsa-miR-25-3p targeted the most genes, hsa-miR-3692-5p and hsa-miR-4270 were targeted by more circRNAs than other miRNAs, hsa-circ-0034326 and hsa-circ-0011950 interacted with three miRNAs. Furthermore, target genes, including NRAS, ITGA5, SLC7A1, SEC14L2, SLC12A5, and SMAD2 were obtained in drug-gene interaction network. Survival analysis showed NRAS, ITGA5, SLC7A1, SEC14L2, SLC12A5, and SMAD2 were significantly associated with prognosis of HCC. NRAS, ITGA5, and SMAD2 were significantly enriched in proteoglycans in cancer. Moreover, hsa-circ-0034326 and hsa-circ-0011950 might function as ceRNAs to play key roles in HCC. Furthermore, miR-25-3p, miR-3692-5p, and miR-4270 might be significant for HCC development. NRAS, ITGA5, SEC14L2, SLC12A5, and SMAD2 might be prognostic factors for HCC patients via proteoglycans in cancer pathway. Taken together, the findings will provide novel insight into pathogenesis, selection of therapeutic targets and prognostic factors for HCC.
| S-EPMC8687197 | biostudies-literature
Project description:There is increasing evidence that miRNAs play an important role in the prognosis of HCC. There is currently a lack of acknowledged models that accurately predict patient prognosis. The aim of this study is to create a miRNA-based model to precisely forecast a patient's prognosis and a miRNA-mRNA network to investigate the function of a targeted mRNA. TCGA miRNA dataset and survival data of HCC patients were downloaded for differential analysis. The outcomes of variance analysis were subjected to univariate and multivariate Cox regression analyses and LASSO analysis. We constructed and visualized prognosis-related models and subsequently used violin plots to probe the function of miRNAs in tumor cells. We predicted the target mRNAs added those to the String database, built PPI protein interaction networks, and screened those mRNA using Cytoscape. The hub mRNA was subjected to GO and KEGG analysis to determine its biological role. Six of them were associated with prognosis: hsa-miR-139-3p, hsa-miR-139-5p, hsa-miR-101-3p, hsa-miR-30d-5p, hsa-miR-5003-3p, and hsa-miR-6844. The prognostic model was highly predictive and consistently performs, with the C index exceeding 0.7 after 1, 3, and 5 years. The model estimated significant differences in the Kaplan-Meier plotter and the model could predict patient prognosis independently of clinical indicators. A relatively stable miRNA prognostic model for HCC patients was constructed, and the model was highly accurate in predicting patients with good stability over 5 years. The miRNA-mRNA network was constructed to explore the function of mRNA.
| S-EPMC10904818 | biostudies-literature
Project description:Background: Hepatocellular carcinoma (HCC) accounts for the majority of liver cancer, with the incidence and mortality rates increasing every year. Despite the improvement of clinical management, substantial challenges remain due to its high recurrence rates and short survival period. This study aimed to identify potential diagnostic and prognostic biomarkers in HCC through bioinformatic analysis. Methods: Datasets from GEO and TCGA databases were used for the bioinformatic analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were carried out by WebGestalt website and clusterProfiler package of R. The STRING database and Cytoscape software were used to establish the protein-protein interaction (PPI) network. The GEPIA website was used to perform expression analyses of the genes. The miRDB, miRWalk, and TargetScan were employed to predict miRNAs and the expression levels of the predicted miRNAs were explored via OncomiR database. LncRNAs were predicted in the StarBase and LncBase while circRNA prediction was performed by the circBank. ROC curve analysis and Kaplan-Meier (KM) survival analysis were performed to evaluate the diagnostic and prognostic value of the gene expression, respectively. Results: A total of 327 upregulated and 422 downregulated overlapping DEGs were identified between HCC tissues and noncancerous liver tissues. The PPI network was constructed with 89 nodes and 178 edges and eight hub genes were selected to predict upstream miRNAs and ceRNAs. A lncRNA/circRNA-miRNA-mRNA network was successfully constructed based on the ceRNA hypothesis, including five lncRNAs (DLGAP1-AS1, GAS5, LINC00665, TYMSOS, and ZFAS1), six circRNAs (hsa_circ_0003209, hsa_circ_0008128, hsa_circ_0020396, hsa_circ_0030051, hsa_circ_0034049, and hsa_circ_0082333), eight miRNAs (hsa-miR-150-5p, hsa-miR-19b-3p, hsa-miR-23b-3p, hsa-miR-26a-5p, hsa-miR-651-5p, hsa-miR-10a-5p, hsa-miR-214-5p and hsa-miR-486-5p), and five mRNAs (CDC6, GINS1, MCM4, MCM6, and MCM7). The ceRNA network can promote HCC progression via cell cycle, DNA replication, and other pathways. Clinical diagnostic and survival analyses demonstrated that the ZFAS1/hsa-miR-150-5p/GINS1 ceRNA regulatory axis had a high diagnostic and prognostic value. Conclusion: These results revealed that cell cycle and DNA replication pathway could be potential pathways to participate in HCC development. The ceRNA network is expected to provide potential biomarkers and therapeutic targets for HCC management, especially the ZFAS1/hsa-miR-150-5p/GINS1 regulatory axis.
| S-EPMC8977626 | biostudies-literature
Project description:BACKGROUND As all we know, gastric cancer (GC) is a highly aggressive disease. Recently, circular RNA (circRNA) was found to play a vital role in regulation of GC. Some circRNAs could regulate messenger RNA (mRNA) expression by functioning as a microRNA (miRNA) sponge. Nevertheless, the circRNA-miRNA-mRNA regulatory network involved GC rarely has been explored and researched. MATERIAL AND METHODS All the differentially expressed circRNAs, miRNAs, and mRNA were derived from Gene Expression Omnibus (GEO) microarray data (GSE78092, GSE89143, GSE93415, and GSE54129). GC level 3 miRNA-sequencing data and clinical information were downloaded from The Cancer Genome Atlas (TCGA) database. Furthermore, a circRNA-miRNA-mRNA regulatory network was constructed by Cytoscape (version 3.6.1). Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway revealed the functions and signaling pathways associated with these target genes. Hub genes of protein-protein interaction (PPI) network were identified by STRING database and cytoHubba. RESULTS The regulatory network consists of 3 circRNAs, 22 miRNAs, and 128 mRNAs. Only 3 miRNAs of the network were consistent with the expression of TCGA and were associated with some clinical features. The results of the functional analysis of 128 mRNAs showed that GO analysis and KEGG pathways of inclusion criteria were 49 and 24, respectively. PPI network and Cytoscape showed that the top 10 hub genes were MYC, CTGF, TGFBR2, TGFBR1, SERPINE1, KRAS, ZEB1, THBS1, CDK6, and TNS1; 4 of which were verified by GEPIA based on TCGA. Highly expressed SERPINE1 had a poor OS (over survival) and DFS (disease-free survival), and TGFBR1 expression increased along with the increase of clinical stages. CONCLUSIONS This study looked at a circRNA-miRNA-mRNA regulatory network associated with GC and explored the potential functions of mRNA in the network, then identified a new molecular marker for prediction, prognosis, and therapeutic targets for clinical patients.
| S-EPMC6880644 | biostudies-literature