scholarly journals Exploring the Key Genes and Pathways in the Formation of Corneal Scar Using Bioinformatics Analysis

2020 ◽  
Vol 2020 ◽  
pp. 1-10
Author(s):  
Binfeng Liu ◽  
Ang Li ◽  
Hongbo Wang ◽  
Jialin Wang ◽  
Gongwei Zhai ◽  
...  
Keyword(s):  
Gene Expression ◽  
Differentially Expressed Genes ◽  
Protein Interaction ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Signal Pathways ◽  
Hub Genes ◽  
Protein Protein Interaction ◽  
Ppi Networks ◽  
Corneal Scar

The Corneal wound healing results in the formation of opaque corneal scar. In fact, millions of people around the world suffer from corneal scars, leading to loss of vision. This study aimed to identify the key changes of gene expression in the formation of opaque corneal scar and provided potential biomarker candidates for clinical treatment and drug target discovery. We downloaded Gene expression dataset GSE6676 from NCBI-GEO, and analyzed the Differentially Expressed Genes (DEGs), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway analyses, and protein-protein interaction (PPI) network. A total of 1377 differentially expressed genes were identified and the result of Functional enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) identification and protein-protein interaction (PPI) networks were performed. In total, 7 hub genes IL6 (interleukin-6), MMP9 (matrix metallopeptidase 9), CXCL10 (C-X-C motif chemokine ligand 10), MAPK8 (mitogen-activated protein kinase 8), TLR4 (toll-like receptor 4), HGF (hepatocyte growth factor), EDN1 (endothelin 1) were selected. In conclusion, the DEGS, Hub genes and signal pathways identified in this study can help us understand the molecular mechanism of corneal scar formation and provide candidate targets for the diagnosis and treatment of corneal scar.

scholarly journals Exploration and validation of related hub gene expression during SARS-CoV-2 infection of human bronchial organoids

2021 ◽  
Vol 15 (1) ◽  
Author(s):  
Ke-Ying Fang ◽  
Wen-Chao Cao ◽  
Tian-Ao Xie ◽  
Jie Lv ◽  
Jia-Xin Chen ◽  
...  
Keyword(s):  
Gene Expression ◽  
Differentially Expressed Genes ◽  
Protein Interaction ◽  
Social Order ◽  
Expression System ◽  
Differentially Expressed ◽  
Dimensional Structure ◽  
Hub Genes ◽  
R Software ◽  
Protein Protein Interaction

Abstract Background In the novel coronavirus pandemic, the high infection rate and high mortality have seriously affected people’s health and social order. To better explore the infection mechanism and treatment, the three-dimensional structure of human bronchus has been employed in a better in-depth study on severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Methods We downloaded a separate microarray from the Integrated Gene Expression System (GEO) on a human bronchial organoids sample to identify differentially expressed genes (DEGS) and analyzed it with R software. After processing with R software, Gene Ontology (GO) and Kyoto PBMCs of Genes and Genomes (KEGG) were analyzed, while a protein–protein interaction (PPI) network was constructed to show the interactions and influence relationships between these differential genes. Finally, the selected highly connected genes, which are called hub genes, were verified in CytoHubba plug-in. Results In this study, a total of 966 differentially expressed genes, including 490 upregulated genes and 476 downregulated genes were used. Analysis of GO and KEGG revealed that these differentially expressed genes were significantly enriched in pathways related to immune response and cytokines. We construct protein-protein interaction network and identify 10 hub genes, including IL6, MMP9, IL1B, CXCL8, ICAM1, FGF2, EGF, CXCL10, CCL2, CCL5, CXCL1, and FN1. Finally, with the help of GSE150728, we verified that CXCl1, CXCL8, CXCL10, CCL5, EGF differently expressed before and after SARS-CoV-2 infection in clinical patients. Conclusions In this study, we used mRNA expression data from GSE150819 to preliminarily confirm the feasibility of hBO as an in vitro model to further study the pathogenesis and potential treatment of COVID-19. Moreover, based on the mRNA differentiated expression of this model, we found that CXCL8, CXCL10, and EGF are hub genes in the process of SARS-COV-2 infection, and we emphasized their key roles in SARS-CoV-2 infection. And we also suggested that further study of these hub genes may be beneficial to treatment, prognostic prediction of COVID-19.

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.

scholarly journals Comprehensive analysis of an lncRNA-miRNA-mRNA competing endogenous RNA network in pulpitis

PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e7135 ◽  
Author(s):  
Fangcao Lei ◽  
Han Zhang ◽  
Xiaoli Xie
Keyword(s):  
Differentially Expressed Genes ◽  
Protein Interaction ◽  
Dental Pulp ◽  
Expression Profiles ◽  
Differentially Expressed ◽  
Hub Genes ◽  
Regulatory Relationship ◽  
Protein Protein Interaction ◽  
Qrt Pcr ◽  
Ppi Networks

Background Pulpitis is a common inflammatory disease that affects dental pulp. It is important to understand the molecular signals of inflammation and repair associated with this process. Increasing evidence has revealed that long noncoding RNAs (lncRNAs), via competitively sponging microRNAs (miRNAs), can act as competing endogenous RNAs (ceRNAs) to regulate inflammation and reparative responses. The aim of this study was to elucidate the potential roles of lncRNA, miRNA and messenger RNA (mRNA) ceRNA networks in pulpitis tissues compared to normal control tissues. Methods The oligo and limma packages were used to identify differentially expressed lncRNAs and mRNAs (DElncRNAs and DEmRNAs, respectively) based on expression profiles in two datasets, GSE92681 and GSE77459, from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were further analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Protein–protein interaction (PPI) networks and modules were established to screen hub genes using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the Molecular Complex Detection (MCODE) plugin for Cytoscape, respectively. Furthermore, an lncRNA-miRNA-mRNA-hub genes regulatory network was constructed to investigate mechanisms related to the progression and prognosis of pulpitis. Then, quantitative real-time polymerase chain reaction (qRT-PCR) was applied to verify critical lncRNAs that may significantly affect the pathogenesis in inflamed and normal human dental pulp. Results A total of 644 upregulated and 264 downregulated differentially expressed genes (DEGs) in pulpitis samples were identified from the GSE77459 dataset, while 8 up- and 19 downregulated probes associated with lncRNA were identified from the GSE92681 dataset. Protein–protein interaction (PPI) based on STRING analysis revealed a network of DEGs containing 4,929 edges and 623 nodes. Upon combined analysis of the constructed PPI network and the MCODE results, 10 hub genes, including IL6, IL8, PTPRC, IL1B, TLR2, ITGAM, CCL2, PIK3CG, ICAM1, and PIK3CD, were detected in the network. Next, a ceRNA regulatory relationship consisting of one lncRNA (PVT1), one miRNA (hsa-miR-455-5p) and two mRNAs (SOCS3 and PLXNC1) was established. Then, we constructed the network in which the regulatory relationship between ceRNA and hub genes was summarized. Finally, our qRT-PCR results confirmed significantly higher levels of PVT1 transcript in inflamed pulp than in normal pulp tissues (p = 0.03). Conclusion Our study identified a novel lncRNA-mediated ceRNA regulatory mechanisms in the pathogenesis of pulpitis.

scholarly journals Key Genes and Prognostic Analysis in HER2+ Breast Cancer

2021 ◽  
Vol 20 ◽  
pp. 153303382098329
Author(s):  
Yujie Weng ◽  
Wei Liang ◽  
Yucheng Ji ◽  
Zhongxian Li ◽  
Rong Jia ◽  
...  
Keyword(s):  
Breast Cancer ◽  
Gene Expression ◽  
Differentially Expressed Genes ◽  
Protein Interaction ◽  
Differentially Expressed ◽  
Protein Kinase Activity ◽  
Protein Protein Interaction ◽  
Risk Of Death ◽  
Her2 Breast Cancer ◽  
Key Genes

Human epidermal growth factor 2 (HER2)+ breast cancer is considered the most dangerous type of breast cancers. Herein, we used bioinformatics methods to identify potential key genes in HER2+ breast cancer to enable its diagnosis, treatment, and prognosis prediction. Datasets of HER2+ breast cancer and normal tissue samples retrieved from Gene Expression Omnibus and The Cancer Genome Atlas databases were subjected to analysis for differentially expressed genes using R software. The identified differentially expressed genes were subjected to gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses followed by construction of protein-protein interaction networks using the STRING database to identify key genes. The genes were further validated via survival and differential gene expression analyses. We identified 97 upregulated and 106 downregulated genes that were primarily associated with processes such as mitosis, protein kinase activity, cell cycle, and the p53 signaling pathway. Visualization of the protein-protein interaction network identified 10 key genes ( CCNA2, CDK1, CDC20, CCNB1, DLGAP5, AURKA, BUB1B, RRM2, TPX2, and MAD2L1), all of which were upregulated. Survival analysis using PROGgeneV2 showed that CDC20, CCNA2, DLGAP5, RRM2, and TPX2 are prognosis-related key genes in HER2+ breast cancer. A nomogram showed that high expression of RRM2, DLGAP5, and TPX2 was positively associated with the risk of death. TPX2, which has not previously been reported in HER2+ breast cancer, was associated with breast cancer development, progression, and prognosis and is therefore a potential key gene. It is hoped that this study can provide a new method for the diagnosis and treatment of HER2 + breast cancer.

Extracting Biological Significant Subnetworks from Protein-Protein Interactions Induced by Differentially Expressed Genes of HIV-1 Vpr Variants

2015 ◽  
Vol 4 (4) ◽  
pp. 35-51 ◽  
Author(s):  
Bandana Barman ◽  
Anirban Mukhopadhyay
Keyword(s):  
Differentially Expressed Genes ◽  
Protein Interaction ◽  
Protein Interactions ◽  
Protein Interaction Network ◽  
Interaction Network ◽  
Differentially Expressed ◽  
Wild Type ◽  
Protein Protein Interaction ◽  
Ppi Networks ◽  
Hiv 1

Identification of protein interaction network is very important to find the cell signaling pathway for a particular disease. The authors have found the differentially expressed genes between two sample groups of HIV-1. Samples are wild type HIV-1 Vpr and HIV-1 mutant Vpr. They did statistical t-test and found false discovery rate (FDR) to identify the genes increased in expression (up-regulated) or decreased in expression (down-regulated). In the test, the authors have computed q-values of test to identify minimum FDR which occurs. As a result they found 172 differentially expressed genes between their sample wild type HIV-1 Vpr and HIV-1 mutant Vpr, R80A. They found 68 up-regulated genes and 104 down-regulated genes. From the 172 differentially expressed genes the authors found protein-protein interaction network with string-db and then clustered (subnetworks) the PPI networks with cytoscape3.0. Lastly, the authors studied significance of subnetworks with performing gene ontology and also studied the KEGG pathway of those subnetworks.

scholarly journals Identification of key genes, pathways and potential therapeutic agents for IgA nephropathy using an integrated bioinformatics analysis

2020 ◽  
Vol 21 (2) ◽  
pp. 147032032091963
Author(s):  
Xiaoxue Chen ◽  
Mindan Sun
Keyword(s):  
Differentially Expressed Genes ◽  
Immunoglobulin A ◽  
Therapeutic Targets ◽  
Gene Expression Omnibus ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Signal Pathways ◽  
Immunoglobulin A Nephropathy ◽  
Hub Genes ◽  
Limma Package

Purpose: This study aims to identify immunoglobulin-A-nephropathy-related genes based on microarray data and to investigate novel potential gene targets for immunoglobulin-A-nephropathy treatment. Methods: Immunoglobulin-A-nephropathy chip data was obtained from the Gene Expression Omnibus database, which included 10 immunoglobulin-A-nephropathy and 22 normal samples. We used the limma package of R software to screen differentially expressed genes in immunoglobulin-A-nephropathy and normal glomerular compartment tissues. Functional enrichment (including cellular components, molecular functions, biological processes) and signal pathways were performed for the differentially expressed genes. The online analysis database (STRING) was used to construct the protein-protein interaction networks of differentially expressed genes, and Cytoscape software was used to identify the hub genes of the signal pathway. In addition, we used the Connectivity Map database to predict possible drugs for the treatment of immunoglobulin-A-nephropathy. Results: A total of 348 differentially expressed genes were screened including 107 up-regulated and 241 down-regulated genes. Functional analysis showed that up-regulated differentially expressed genes were mainly concentrated on leukocyte migration, and the down-regulated differentially expressed genes were significantly enriched in alpha-amino acid metabolic process. A total of six hub genes were obtained: JUN, C3AR1, FN1, AGT, FOS, and SUCNR1. The small-molecule drugs thapsigargin, ciclopirox and ikarugamycin were predicted therapeutic targets against immunoglobulin-A-nephropathy. Conclusion: Differentially expressed genes and hub genes can contribute to understanding the molecular mechanism of immunoglobulin-A-nephropathy and providing potential therapeutic targets and drugs for the diagnosis and treatment of immunoglobulin-A-nephropathy.

scholarly journals Hub Genes Related to the Pathogenesis and Progression of Colorectal Cancer and Adenoma by Integrative Bioinformatics Approaches

2021 ◽  
Author(s):  
Yuxuan HUANG ◽  
Ge CUI
Keyword(s):  
Colorectal Cancer ◽  
Cell Cycle ◽  
Differentially Expressed Genes ◽  
Protein Interaction ◽  
Biological Networks ◽  
Interaction Network ◽  
Differentially Expressed ◽  
Venn Diagram ◽  
Hub Genes ◽  
Protein Protein Interaction

Abstract Aims: To utilize the bioinformatics to analyze the differentially expressed genes (DEGs), interaction proteins, perform gene enrichment analysis, protein-protein interaction network (PPI) and map the hub genes between colorectal cancer(CRC) and colorectal adenocarcinomas(CA).Methods: We analyzed a microarray dataset (GSE32323 and GSE4183) from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) in tumor tissues and non-cancerous tissues were identified using the dplyr and Venn diagram packages of the R Studio software. Functional annotation of the DEGs was performed using the Gene Ontology (GO) website. Pathway enrichment (KEGG) used the WebGestalt to analyze the data and R Studio to generate the graph. We constructed a protein–protein interaction (PPI) network of DEGs using STRING and Cytoscape software was used for visualization. Survival analysis of the hub genes and was performed using the online platform GEPIA to determine the prognostic value of the expression of hub genes in cell lines from CRC patients. The expression of molecules with prognostic values was validated on the UALCAN database. The expression of hub genes was examined using the Human Protein Atlas. Results: Applying the GEO2R analysis and R studio, we identified a total of 471 upregulated and 278 downregulated DEGs. By using the online database WebGestalt, we identified the most relevant biological networks involving DEGs with statistically significant differences in expression were mainly associated with biological processes involved in the cell proliferation, cell cycle transition, cell homeostasis and indicated the role of each DEGs in cell cycle regulation pathways. We found 10 hub genes with prognostic values were overexpressed in the CRC and CA samples.Conclusion: we found out ten hub genes and three core genes closely associated with the pathogenesis and prognosis of CRC and CA, which is of great significance for colorectal tumor early detection and prognosis evaluation.

scholarly journals Weighted gene co-expression network analysis to identify key modules and hub genes associated with paucigranulocytic asthma

2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Min Li ◽  
Wenye Zhu ◽  
Chu Wang ◽  
Yuanyuan Zheng ◽  
Shibo Sun ◽  
...  
Keyword(s):  
Gene Expression ◽  
Network Analysis ◽  
Differentially Expressed Genes ◽  
Immune Status ◽  
Expression Patterns ◽  
Gene Expression Pattern ◽  
Gene Set Enrichment Analysis ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Hub Genes

Abstract Background Asthma is a heterogeneous disease that can be divided into four inflammatory phenotypes: eosinophilic asthma (EA), neutrophilic asthma (NA), mixed granulocytic asthma (MGA), and paucigranulocytic asthma (PGA). While research has mainly focused on EA and NA, the understanding of PGA is limited. In this study, we aimed to identify underlying mechanisms and hub genes of PGA. Methods Based on the dataset from Gene Expression Omnibus(GEO), weighted gene coexpression network analysis (WGCNA), differentially expressed genes (DEGs) analysis and protein–protein interaction (PPI) network analysis were conducted to construct a gene network and to identify key gene modules and hub genes. Functional enrichment analyses were performed to investigate the biological process, pathways and immune status of PGA. The hub genes were validated in a separate dataset. Results Compared to non-PGA, PGA had a different gene expression pattern, in which 449 genes were differentially expressed. One gene module significantly associated with PGA was identified. Intersection between the differentially expressed genes (DEGs) and the genes from the module that were most relevant to PGA were mainly enriched in inflammation and immune response regulation. The single sample Gene Set Enrichment Analysis (ssGSEA) suggested a decreased immune infiltration and function in PGA. Finally six hub genes of PGA were identified, including ADCY2, CXCL1, FPRL1, GPR109B, GPR109A and ADCY3, which were validated in a separate dataset of GSE137268. Conclusions Our study characterized distinct gene expression patterns, biological processes and immune status of PGA and identified hub genes, which may improve the understanding of underlying mechanism and provide potential therapeutic targets for PGA.

Integrated computational analysis revealed hub genes in lung adenocarcinoma

2021 ◽  
Author(s):  
Muhammad Jamal ◽  
Abdul Saboor Khan ◽  
Hina Iqbal Bangash ◽  
Tian Xie ◽  
Tianbao Song ◽  
...  
Keyword(s):  
Gene Expression ◽  
Lung Cancer ◽  
Survival Analysis ◽  
Differentially Expressed Genes ◽  
Functional Enrichment ◽  
Differentially Expressed ◽  
Normal Lung ◽  
Hub Genes ◽  
Altered Expression ◽  
Study Gene Expression

Abstract Background Lung cancer (LUCA) is the leading cause of cancer-related morbidities and mortalities globally. Despite the recent advancements in lung cancer research, understanding of the molecular mechanism underlying LUCA tumorigenesis and prognosis remains suboptimal. This study aims to identify the candidate biomarkers and therapeutic genes in lung cancer. Methods In this study, gene expression profiles of GSE30219, GSE33532, GSE32863 and GSE43458 were downloaded from GEO. The differentially expressed genes (DEGs) in LUAD tissue and normal lung tissue with a p-value < 0.05 and a |log fold change (FC)| >1.0 were identified by GEO2R. For functional enrichment analysis of these DEGs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed with KOBAS and DAVID tools. Next, the candidate hub genes were filtered out with Cytoscape using CytoHubba plugin. These hub genes were validated by (the Cancer Genome Atlas) TCGA-based gene expression analysis, protein-protein network interaction (PPI) analysis, survival analysis. Moreover, the expression of these genes in cancer and normal tissue was assessed in the Human Protein Atlas (HPA) database. In addition, miRNA network of the hub genes was constructed. Finally, DGIdb database was used to check the drug-targeting potentials of the hub genes. Results a total of 332 overlapping differentially expressed genes (DEGs) including 73 upregulated and 259 downregulated, respectively were identified. GO analysis revealed that the DEGs were principally regulating various cancer-associated functions and pathways. The module analysis revealed 55 hub genes in 4 modules. The survival analysis through Kaplan-Meier (KM) plotter indicated that the altered expression of these genes resulted in the poor overall survival (OS) of LUCA patients. Moreover, these genes show a differential expression on both protein and mRNA level in cancer patient compared to the normal. In addition, in addition, 6 potential microRNAs (miRNAs) interacting with hub genes were identified. Finally, a list of 117 therapeutic small molecules was tabulated that could facilitate LUCA treatment. Conclusions the findings of this study may help in the development of novel and reliable biomarkers for diagnosis, prognosis and therapeutic intervention for LUAD.

scholarly journals Aberrantly hydroxymethylated differentially expressed genes and the associated protein pathways in osteoarthritis

PeerJ ◽  
2019 ◽  
Vol 7 ◽  
pp. e6425 ◽  
Author(s):  
Yang Fang ◽  
Pingping Wang ◽  
Lin Xia ◽  
Suwen Bai ◽  
Yonggang Shen ◽  
...  
Keyword(s):  
Gene Expression ◽  
Gene Ontology ◽  
Differentially Expressed Genes ◽  
Kegg Pathway ◽  
Interaction Network ◽  
Differentially Expressed ◽  
Hub Genes ◽  
Protein Protein Interaction ◽  
Highly Expressed Genes ◽  
Protein Protein Interaction Network

Background The elderly population is at risk of osteoarthritis (OA), a common, multifactorial, degenerative joint disease. Environmental, genetic, and epigenetic (such as DNA hydroxymethylation) factors may be involved in the etiology, development, and pathogenesis of OA. Here, comprehensive bioinformatic analyses were used to identify aberrantly hydroxymethylated differentially expressed genes and pathways in osteoarthritis to determine the underlying molecular mechanisms of osteoarthritis and susceptibility-related genes for osteoarthritis inheritance. Methods Gene expression microarray data, mRNA expression profile data, and a whole genome 5hmC dataset were obtained from the Gene Expression Omnibus repository. Differentially expressed genes with abnormal hydroxymethylation were identified by MATCH function. Gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the genes differentially expressed in OA were performed using Metascape and the KOBAS online tool, respectively. The protein–protein interaction network was built using STRING and visualized in Cytoscape, and the modular analysis of the network was performed using the Molecular Complex Detection app. Results In total, 104 hyperhydroxymethylated highly expressed genes and 14 hypohydroxymethylated genes with low expression were identified. Gene ontology analyses indicated that the biological functions of hyperhydroxymethylated highly expressed genes included skeletal system development, ossification, and bone development; KEGG pathway analysis showed enrichment in protein digestion and absorption, extracellular matrix–receptor interaction, and focal adhesion. The top 10 hub genes in the protein–protein interaction network were COL1A1, COL1A2, COL2A1, COL3A1, COL5A1, COL5A2, COL6A1, COL8A1, COL11A1, and COL24A1. All the aforementioned results are consistent with changes observed in OA. Conclusion After comprehensive bioinformatics analysis, we found aberrantly hydroxymethylated differentially expressed genes and pathways in OA. The top 10 hub genes may be useful hydroxymethylation analysis biomarkers to provide more accurate OA diagnoses and target genes for treatment of OA.

Sign in / Sign up

Export Citation Format

Share Document