Identification of Potential Biomarkers in Synovial Tissue of Rheumatoid Arthritis via Integrated Microarray Analysis. Identification of Biomarkers of Rheumatoid Arthritis

2021 ◽  
Author(s):  
Li Tao ◽  
ChaoLiang Xiong ◽  
Li Xue
Keyword(s):  
Rheumatoid Arthritis ◽  
Gene Expression ◽  
Pathway Analysis ◽  
Expression Profiles ◽  
Diagnostic Value ◽  
Functional Enrichment ◽  
Venn Diagram ◽  
Ppi Network ◽  
Hub Genes ◽  
Gene Markers

Abstract Background: Rheumatoid arthritis (RA) is a chronic autoimmune disease characterized by synovitis and subsequent destruction of cartilage and bone. This study aimed to explore RA-related gene markers and the underlying molecular mechanism.Material and Methods: The expression profiles of GSE77298, GSE55235 and GSE12021 were obtained from the Gene Expression Omnibus database. Then, the differential gene expression analysis was conducted between GSE77298 and GSE55235 datasets. Limma package and a Venn diagram were utilized to screen the overlapping differentially expressed genes (DEGs), and Functional enrichment and pathway analysis were performed by using DAVID database. Subsequently, a protein-protein interaction (PPI) network was established, and candidate hub genes were recognized by using STRING and Cytoscape software. Finally, another dataset (GSE12021) was used for the validation of diagnostic value of the candidate hub genes and to identify real hub genes by using receiver operating characteristic (ROC) curves.Results: A total of 385 DEGs were detected, which include 19 downregulated genes and 366 upregulated genes. GO and KEGG pathway analysis showed that DEGs was mainly enriched in various immune and inflammatory response-related functions and pathways. The PPI network was composed of 374 nodes and 767 edges. A total of 8 real hub genes (HLA-DRA, HLA-DRB1, LCK, VAV1, HLA-DPA1, HLA-DPB1, C3AR1 and CD3D) which displayed an excellent diagnostic value for RA were identified.Conclusion: these findings may provide novel and reliable biomarkers for RA, which have some interesting implications for early diagnosis, prognosis and targeted therapy.

scholarly journals Identification of Potentially Therapeutic Target Genes of Hepatocellular Carcinoma

2020 ◽  
Vol 17 (3) ◽  
pp. 1053
Author(s):  
Chengzhang Li ◽  
Jiucheng Xu
Keyword(s):  
Gene Expression ◽  
Target Genes ◽  
Expression Profiles ◽  
Differential Expression Analysis ◽  
Gene Expression Profiles ◽  
Functional Enrichment ◽  
Ppi Network ◽  
Hub Genes ◽  
Gene Differential Expression ◽  
Oncomine Database

Background: Hepatocellular carcinoma (HCC) is a major threat to public health. However, few effective therapeutic strategies exist. We aimed to identify potentially therapeutic target genes of HCC by analyzing three gene expression profiles. Methods: The gene expression profiles were analyzed with GEO2R, an interactive web tool for gene differential expression analysis, to identify common differentially expressed genes (DEGs). Functional enrichment analyses were then conducted followed by a protein-protein interaction (PPI) network construction with the common DEGs. The PPI network was employed to identify hub genes, and the expression level of the hub genes was validated via data mining the Oncomine database. Survival analysis was carried out to assess the prognosis of hub genes in HCC patients. Results: A total of 51 common up-regulated DEGs and 201 down-regulated DEGs were obtained after gene differential expression analysis of the profiles. Functional enrichment analyses indicated that these common DEGs are linked to a series of cancer events. We finally identified 10 hub genes, six of which (OIP5, ASPM, NUSAP1, UBE2C, CCNA2, and KIF20A) are reported as novel HCC hub genes. Data mining the Oncomine database validated that the hub genes have a significant high level of expression in HCC samples compared normal samples (t-test, p < 0.05). Survival analysis indicated that overexpression of the hub genes is associated with a significant reduction (p < 0.05) in survival time in HCC patients. Conclusions: We identified six novel HCC hub genes that might be therapeutic targets for the development of drugs for some HCC patients.

scholarly journals Network analyses of circular RNAs expression profiles and their hub genes in pancreatic ductal adenocarcinoma

2019 ◽  
Author(s):  
Yanyan Tang ◽  
Ping Zhang
Keyword(s):  
Gene Expression ◽  
Pancreatic Ductal Adenocarcinoma ◽  
Signaling Pathway ◽  
Pathway Analysis ◽  
Kegg Pathway ◽  
Expression Profiles ◽  
Ductal Adenocarcinoma ◽  
Ppi Network ◽  
Hub Genes ◽  
Kegg Pathway Analysis

Abstract Pancreatic ductal adenocarcinoma (PDAC) is one of the most common malignant tumor in digestive system. CircRNAs involve in lots of biological processes through interacting with miRNAs and their targeted mRNA. We obtained the circRNA gene expression profiles from Gene Expression Omnibus (GEO) and identified differentially expressed genes (DEGs) between PDAC samples and paracancerous tissues. Bioinformatics analyses, including GO analysis, KEGG pathway analysis and PPI network analysis, were conducted for further investigation. We also constructed circRNA‑microRNA-mRNA co-expression network. A total 291 differentially expressed circRNAs were screened out. The GO enrichment analysis revealed that up-regulated DEGs were mainly involved metabolic process, biological regulation, and gene expression, and down-regulated DEGs were involved in cell communication, single-organism process, and signal transduction. The KEGG pathway analysis, the upregulated circRNAs were enriched cGMP-PKG signaling pathway, and HTLV-I infection, while the downregulated circRNAs were enriched in protein processing in endoplasmic reticulum, insulin signaling pathway, regulation of actin cytoskeleton, etc. Four genes were identified from PPI network as both hub genes and module genes, and their circRNA‑miRNA-mRNA regulatory network also be constructed. Our study indicated possible involvement of dysregulated circRNAs in the development of PDAC and promoted our understanding of the underlying molecular mechanisms.

scholarly journals Bioinformatics analysis of differentially expressed genes and pathways in the development of cervical cancer

BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Baojie Wu ◽  
Shuyi Xi
Keyword(s):  
Gene Expression ◽  
Cervical Cancer ◽  
Differentially Expressed Genes ◽  
Expression Profiles ◽  
Gene Expression Profiles ◽  
Molecular Targets ◽  
Enrichment Analysis ◽  
Ppi Network ◽  
Hub Genes ◽  
Kegg Pathways

Abstract Background This study aimed to explore and identify key genes and signaling pathways that contribute to the progression of cervical cancer to improve prognosis. Methods Three gene expression profiles (GSE63514, GSE64217 and GSE138080) were screened and downloaded from the Gene Expression Omnibus database (GEO). Differentially expressed genes (DEGs) were screened using the GEO2R and Venn diagram tools. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. Gene set enrichment analysis (GSEA) was performed to analyze the three gene expression profiles. Moreover, a protein–protein interaction (PPI) network of the DEGs was constructed, and functional enrichment analysis was performed. On this basis, hub genes from critical PPI subnetworks were explored with Cytoscape software. The expression of these genes in tumors was verified, and survival analysis of potential prognostic genes from critical subnetworks was conducted. Functional annotation, multiple gene comparison and dimensionality reduction in candidate genes indicated the clinical significance of potential targets. Results A total of 476 DEGs were screened: 253 upregulated genes and 223 downregulated genes. DEGs were enriched in 22 biological processes, 16 cellular components and 9 molecular functions in precancerous lesions and cervical cancer. DEGs were mainly enriched in 10 KEGG pathways. Through intersection analysis and data mining, 3 key KEGG pathways and related core genes were revealed by GSEA. Moreover, a PPI network of 476 DEGs was constructed, hub genes from 12 critical subnetworks were explored, and a total of 14 potential molecular targets were obtained. Conclusions These findings promote the understanding of the molecular mechanism of and clinically related molecular targets for cervical cancer.

scholarly journals LASSO and Bioinformatics Analysis in the Identification of Key Genes for Prognostic Genes of Gynecologic Cancer

2021 ◽  
Vol 11 (11) ◽  
pp. 1177
Author(s):  
Shao-Hua Yu ◽  
Jia-Hua Cai ◽  
De-Lun Chen ◽  
Szu-Han Liao ◽  
Yi-Zhen Lin ◽  
...  
Keyword(s):  
Gene Expression ◽  
Endometrial Cancer ◽  
Pathway Analysis ◽  
Bioinformatics Analysis ◽  
Expression Profiles ◽  
Gynecologic Cancer ◽  
Functional Enrichment ◽  
Lasso Regression ◽  
Potential Biomarker ◽  
Study Gene Expression

The aim of this study is to identify potential biomarkers for early diagnosis of gynecologic cancer in order to improve survival. Cervical cancer (CC) and endometrial cancer (EC) are the most common malignant tumors of gynecologic cancer among women in the world. As the underlying molecular mechanisms in both cervical and endometrial cancer remain unclear, a comprehensive and systematic bioinformatics analysis is required. In our study, gene expression profiles of GSE9750, GES7803, GES63514, GES17025, GES115810, and GES36389 downloaded from Gene Expression Omnibus (GEO) were utilized to analyze differential gene expression between cancer and normal tissues. A total of 78 differentially expressed genes (DEGs) common to CC and EC were identified to perform the functional enrichment analyses, including gene ontology and pathway analysis. KEGG pathway analysis of 78 DEGs indicated that three main types of pathway participate in the mechanism of gynecologic cancer such as drug metabolism, signal transduction, and tumorigenesis and development. Furthermore, 20 diagnostic signatures were confirmed using the least absolute shrink and selection operator (LASSO) regression with 10-fold cross validation. Finally, we used the GEPIA2 online tool to verify the expression of 20 genes selected by the LASSO regression model. Among them, the expression of PAMR1 and SLC24A3 in tumor tissues was downregulated significantly compared to the normal tissue, and found to be statistically significant in survival rates between the CC and EC of patients (p < 0.05). The two genes have their function: (1.) PAMR1 is a tumor suppressor gene, and many studies have proven that overexpression of the gene markedly suppresses cell growth, especially in breast cancer and polycystic ovary syndrome; (2.) SLC24A3 is a sodium–calcium regulator of cells, and high SLC24A3 levels are associated with poor prognosis. In our study, the gene signatures can be used to predict CC and EC prognosis, which could provide novel clinical evidence to serve as a potential biomarker for future diagnosis and treatment.

scholarly journals Integrative transcriptome data mining for identification of core lncRNAs in breast cancer

PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7821 ◽  
Author(s):  
Xiaoming Zhang ◽  
Jing Zhuang ◽  
Lijuan Liu ◽  
Zhengguo He ◽  
Cun Liu ◽  
...  
Keyword(s):  
Breast Cancer ◽  
Gene Expression ◽  
Prognostic Value ◽  
Early Stage ◽  
Expression Profiles ◽  
Diagnostic Value ◽  
The Cancer Genome Atlas ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Data Set

Background Cumulative evidence suggests that long non-coding RNAs (lncRNAs) play an important role in tumorigenesis. This study aims to identify lncRNAs that can serve as new biomarkers for breast cancer diagnosis or screening. Methods First, the linear fitting method was used to identify differentially expressed genes from the breast cancer RNA expression profiles in The Cancer Genome Atlas (TCGA). Next, the diagnostic value of all differentially expressed lncRNAs was evaluated using a receiver operating characteristic (ROC) curve. Then, the top ten lncRNAs with the highest diagnostic value were selected as core genes for clinical characteristics and prognosis analysis. Furthermore, core lncRNA-mRNA co-expression networks based on weighted gene co-expression network analysis (WGCNA) were constructed, and functional enrichment analysis was performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID). The differential expression level and diagnostic value of core lncRNAs were further evaluated by using independent data set from Gene Expression Omnibus (GEO). Finally, the expression status and prognostic value of core lncRNAs in various tumors were analyzed based on Gene Expression Profiling Interactive Analysis (GEPIA). Results Seven core lncRNAs (LINC00478, PGM5-AS1, AL035610.1, MIR143HG, RP11-175K6.1, AC005550.4, and MIR497HG) have good single-factor diagnostic value for breast cancer. AC093850.2 has a prognostic value for breast cancer. AC005550.4 and MIR497HG can better distinguish breast cancer patients in early-stage from the advanced-stage. Low expression of MAGI2-AS3, LINC00478, AL035610.1, MIR143HG, and MIR145 may be associated with lymph node metastasis in breast cancer. Conclusion Our study provides candidate biomarkers for the diagnosis and prognosis of breast cancer, as well as a bioinformatics basis for the further elucidation of the molecular pathological mechanism of breast cancer.

scholarly journals Identification of Potentially Functional Circular RNAs hsa_circ_0070934 and hsa_circ_0004315 as Prognostic Factor of Hepatocellular Carcinoma by Integrated Bioinformatics Analysis

2021 ◽  
Author(s):  
Pejman Morovat ◽  
Saman Morovat ◽  
Arash M. Ashrafi ◽  
Shahram Teimourian
Keyword(s):  
Gene Expression ◽  
Hepatocellular Carcinoma ◽  
Prognostic Factor ◽  
Sequence Data ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Circular Rnas ◽  
Ppi Network ◽  
Hub Genes ◽  
Cerna Network

Abstract Hepatocellular carcinoma (HCC) is one of the most prevalent cancers worldwide, which has a high mortality rate and poor treatment outcomes with yet unknown molecular basis. It seems that gene expression plays a pivotal role in the pathogenesis of the disease. Circular RNAs (circRNAs) can interact with microRNAs (miRNAs) to regulate gene expression in various malignancies by acting as competitive endogenous RNAs (ceRNAs). However, the potential pathogenesis roles of the ceRNA network among circRNA/miRNA/mRNA in HCC are unclear. In this study, first, the HCC circRNA expression data were obtained from three Gene Expression Omnibus microarray datasets (GSE164803, GSE94508, GSE97332), and the differentially expressed circRNAs (DECs) were identified using R limma package. Also, the liver hepatocellular carcinoma (LIHC) miRNA and mRNA sequence data were retrieved from TCGA, and differentially expressed miRNAs (DEMIs) and mRNAs (DEGs) were determined using the R DESeq2 package. Second, CSCD website was used to uncover the binding sites of miRNAs on DECs. The DECs' potential target miRNAs were revealed by conducting an intersection between predicted miRNAs from CSCD and downregulated DEMIs. Third, some related genes were uncovered by intersecting targeted genes predicted by miRWalk and targetscan online tools with upregulated DEGs. The ceRNA network was then built using the Cytoscape software. The functional enrichment and the overall survival time of these potential targeted genes were analyzed, and a PPI network was constructed in the STRING database. Network visualization was performed by Cytoscape, and ten hub genes were detected using the CytoHubba plugin tool. Four DECs (hsa_circ_0000520, hsa_circ_0008616, hsa_circ_0070934, hsa_circ_0004315) were obtained and six miRNAs (hsa-miR-542-5p, hsa-miR-326, hsa-miR-511-5p, hsa-miR-195-5p, hsa-miR-214-3p, and hsa-miR-424-5p) which are regulated by the above DECs were identified. Then 543 overlapped genes regulated by six miRNAs mentioned above were predicted. Functional enrichment analysis showed that these genes are mostly associated with cancer regulation functions. Ten hub genes (TTK،AURKB, KIF20A، KIF23، CEP55، CDC6، DTL، NCAPG، CENPF، PLK4) have been screened from the PPI network of the 204 survival-related genes. KIF20A, NCAPG, TTK, PLK4, and CDC6 were selected for the highest significant p-values. In the end, a circRNA-miRNA-mRNA regulatory axis was established for five final selected hub genes. This study implies the potential pathogenesis of the obtained network and proposes that the two DECs (has_circ_0070934 and has_circ_0004315) may be important prognostic factor for HCC.

scholarly journals Using Weighted Gene Co-Expression Network Analysis to Identify Key Genes Related to the Onset of Preeclampsia

2020 ◽  
Author(s):  
Xinyang Shen ◽  
Zhijian Wang ◽  
Zhirui Zeng ◽  
Zhenqin Huang ◽  
Xiaowei Huang ◽  
...  
Keyword(s):  
Enrichment Analysis ◽  
Diagnostic Value ◽  
Functional Enrichment ◽  
Immune Response Gene ◽  
Control Group ◽  
Ppi Network ◽  
Hub Genes ◽  
Ontology Term ◽  
Gene Data ◽  
Core Genes

Abstract Background: Preeclampsia is a form of hypertension in pregnancy, which induced by complicated factors. However, the pathogenesis of the disease is unclear. The present study was aimed to discover the critical biomarkers associated with the occurrence and development of preeclampsia. Methods：Gene data profile GSE75010 was downloaded from the Gene Expression Omnibus (GEO) database and used as discovery cohort to establish a WGCNA network determining significant modules which associated with clinical traits. Subsequently, functional enrichment analysis, pathway analysis and protein-protein interaction (PPI) network construction were performed on the core genes in significant modules to identify hub genes. Then, gene data profile GSE25906 was used as verified cohort to determine their diagnostic value of hub genes. The protein expression levels of these hub genes in preeclampsia and control placental tissues were verified using immunohistochemistry method. Finally, GSEA was performed to analyze their enrichment pathways. Results: Total 33 co-expression modules were identified after the establishment of WGCNA, of which 4 gene modules were identified as significant modules because they were related to multiple (>3) clinical traits. Total 75 core genes in significant modules were analyzed, and results showed that they were mainly enriched in adaptive immune response (Gene Ontology term) and platelet activation (Kyoto Encyclopedia of Genes and Genomes term). Finally, a total of 5 genes including TYROBP, PLEK, LCP2, HCK, ITGAM were identified as hub genes which scored high in PPI network and had high diagnostic value. Furthermore, the protein level of these 5 genes in placental tissues of preeclampsia was lower than that of the control group. Moreover, these 5 genes were all enriched in 17 pathways, including autoimmunity pathway. Conclusions：These 5 genes (TYROBP, PLEK, LCP2, HCK, ITGAM) may be closely related to the pathogenesis of preeclampsia, which may also help the diagnosis and therapy of preeclampsia.

scholarly journals Identification and validation of hub genes of synovial tissue for patients with osteoarthritis and rheumatoid arthritis

Hereditas ◽  
2021 ◽  
Vol 158 (1) ◽  
Author(s):  
Yanzhi Ge ◽  
Zuxiang Chen ◽  
Yanbin Fu ◽  
Xiujuan Xiao ◽  
Haipeng Xu ◽  
...  
Keyword(s):  
Rheumatoid Arthritis ◽  
Gene Expression ◽  
Signaling Pathway ◽  
Quantitative Polymerase Chain Reaction ◽  
Functional Enrichment ◽  
Proximal Promoter ◽  
Integrated Analysis ◽  
Data Sets ◽  
Hub Genes ◽  
Multiple Gene

Abstract Background Osteoarthritis (OA) and rheumatoid arthritis (RA) were two major joint diseases with similar clinical phenotypes. This study aimed to determine the mechanistic similarities and differences between OA and RA by integrated analysis of multiple gene expression data sets. Methods Microarray data sets of OA and RA were obtained from the Gene Expression Omnibus (GEO). By integrating multiple gene data sets, specific differentially expressed genes (DEGs) were identified. The Gene Ontology (GO) functional annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and protein–protein interaction (PPI) network analysis of DEGs were conducted to determine hub genes and pathways. The “Cell Type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT)” algorithm was employed to evaluate the immune infiltration cells (IICs) profiles in OA and RA. Moreover, mouse models of RA and OA were established, and selected hub genes were verified in synovial tissues with quantitative polymerase chain reaction (qPCR). Results A total of 1116 DEGs were identified between OA and RA. GO functional enrichment analysis showed that DEGs were enriched in regulation of cell morphogenesis involved in differentiation, positive regulation of neuron differentiation, nuclear speck, RNA polymerase II transcription factor complex, protein serine/threonine kinase activity and proximal promoter sequence-specific DNA binding. KEGG pathway analysis showed that DEGs were enriched in EGFR tyrosine kinase inhibitor resistance, ubiquitin mediated proteolysis, FoxO signaling pathway and TGF-beta signaling pathway. Immune cell infiltration analysis identified 9 IICs with significantly different distributions between OA and RA samples. qPCR results showed that the expression levels of the hub genes (RPS6, RPS14, RPS25, RPL11, RPL27, SNRPE, EEF2 and RPL19) were significantly increased in OA samples compared to their counterparts in RA samples (P < 0.05). Conclusion This large-scale gene analyses provided new insights for disease-associated genes, molecular mechanisms as well as IICs profiles in OA and RA, which may offer a new direction for distinguishing diagnosis and treatment between OA and RA.

scholarly journals Identification of potential biomarkers with colorectal cancer based on bioinformatics analysis and machine learning

2021 ◽  
Vol 18 (6) ◽  
pp. 8997-9015
Author(s):  
Ahmed Hammad ◽  
◽  
Mohamed Elshaer ◽  
Xiuwen Tang ◽  
◽  
...  
Keyword(s):  
Gene Expression ◽  
Colorectal Cancer ◽  
Machine Learning ◽  
Roc Curve ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Ppi Network ◽  
Hub Genes ◽  
Survival Analyses ◽  
Potential Biomarkers

<abstract> <p>Colorectal cancer (CRC) is one of the most common malignancies worldwide. Biomarker discovery is critical to improve CRC diagnosis, however, machine learning offers a new platform to study the etiology of CRC for this purpose. Therefore, the current study aimed to perform an integrated bioinformatics and machine learning analyses to explore novel biomarkers for CRC prognosis. In this study, we acquired gene expression microarray data from Gene Expression Omnibus (GEO) database. The microarray expressions GSE103512 dataset was downloaded and integrated. Subsequently, differentially expressed genes (DEGs) were identified and functionally analyzed via Gene Ontology (GO) and Kyoto Enrichment of Genes and Genomes (KEGG). Furthermore, protein protein interaction (PPI) network analysis was conducted using the STRING database and Cytoscape software to identify hub genes; however, the hub genes were subjected to Support Vector Machine (SVM), Receiver operating characteristic curve (ROC) and survival analyses to explore their diagnostic values. Meanwhile, TCGA transcriptomics data in Gene Expression Profiling Interactive Analysis (GEPIA) database and the pathology data presented by in the human protein atlas (HPA) database were used to verify our transcriptomic analyses. A total of 105 DEGs were identified in this study. Functional enrichment analysis showed that these genes were significantly enriched in biological processes related to cancer progression. Thereafter, PPI network explored a total of 10 significant hub genes. The ROC curve was used to predict the potential application of biomarkers in CRC diagnosis, with an area under ROC curve (AUC) of these genes exceeding 0.92 suggesting that this risk classifier can discriminate between CRC patients and normal controls. Moreover, the prognostic values of these hub genes were confirmed by survival analyses using different CRC patient cohorts. Our results demonstrated that these 10 differentially expressed hub genes could be used as potential biomarkers for CRC diagnosis.</p> </abstract>

scholarly journals Comprehensive analysis of DNA methylation and gene expression profiles in cholangiocarcinoma

2019 ◽  
Vol 19 (1) ◽  
Author(s):  
Cheng Zhang ◽  
Bingye Zhang ◽  
Di Meng ◽  
Chunlin Ge
Keyword(s):  
Gene Expression ◽  
Dna Methylation ◽  
Differentially Expressed Genes ◽  
Expression Profiles ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Hub Genes ◽  
Ppi Networks ◽  
Differentially Methylated Genes ◽  
Tumorigenesis And Progression

Abstract Background The incidence of cholangiocarcinoma (CCA) has risen in recent years, and it has become a significant health burden worldwide. However, the mechanisms underlying tumorigenesis and progression of this disease remain largely unknown. An increasing number of studies have demonstrated crucial biological functions of epigenetic modifications, especially DNA methylation, in CCA. The present study aimed to identify and analyze methylation-regulated differentially expressed genes (MeDEGs) involved in CCA tumorigenesis and progression by bioinformatics analysis. Methods The gene expression profiling dataset (GSE119336) and gene methylation profiling dataset (GSE38860) were obtained from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) were identified using the limma packages of R and GEO2R, respectively. The MeDEGs were obtained by overlapping the DEGs and DMGs. Functional enrichment analyses of these genes were then carried out. Protein–protein interaction (PPI) networks were constructed using STRING and visualized in Cytoscape to determine hub genes. Finally, the results were verified based on The Cancer Genome Atlas (TCGA) database. Results We identified 98 hypermethylated, downregulated genes and 93 hypomethylated, upregulated genes after overlapping the DEGs and DMGs. These genes were mainly enriched in the biological processes of the cell cycle, nuclear division, xenobiotic metabolism, drug catabolism, and negative regulation of proteolysis. The top nine hub genes of the PPI network were F2, AHSG, RRM2, AURKB, CCNA2, TOP2A, BIRC5, PLK1, and ASPM. Moreover, the expression and methylation status of the hub genes were significantly altered in TCGA. Conclusions Our study identified novel methylation-regulated differentially expressed genes (MeDEGs) and explored their related pathways and functions in CCA, which may provide novel insights into a further understanding of methylation-mediated regulatory mechanisms in CCA.

Sign in / Sign up

Export Citation Format

Share Document