Identifying hub genes and immune infiltration of osteoarthritis using comprehensive bioinformatics analysis.
Ontology highlight
ABSTRACT: Background
Osteoarthritis (OA) is the most common chronic degenerative joint disorder globally that is characterized by synovitis, cartilage degeneration, joint space stenosis, and sub-cartilage bone hyperplasia. However, the pathophysiologic mechanisms of OA have not been thoroughly investigated.Methods
In this study, we conducted various bioinformatics analyses to identify hub biomarkers and immune infiltration in OA. The gene expression profiles of synovial tissues from 29 healthy controls and 36 OA samples were obtained from the gene expression omnibus database to identify differentially expressed genes (DEGs). The CIBERSORT algorithm was used to explore the association between immune infiltration and arthritis.Results
Eighteen hub DEGs were identified as critical biomarkers for OA. Through gene ontology and pathway enrichment analyses, it was found that these DEGs were primarily involved in PI3K-Akt signaling pathway and Rap1 signaling pathway. Furthermore, immune infiltration analysis revealed differences in immune infiltration between patients with OA and healthy controls. The hub gene ZNF160 was closely related to immune cells, especially mast cell activation in OA.Conclusion
Overall, this study presented a novel method to identify hub DEGs and their correlation with immune infiltration, which may provide novel insights into the diagnosis and treatment of patients with OA.
SUBMITTER: Wu ZY
PROVIDER: S-EPMC8527722 | biostudies-literature | 2021 Oct
REPOSITORIES: biostudies-literature
Publications
Identifying hub genes and immune infiltration of osteoarthritis using comprehensive bioinformatics analysis.
Wu Zheng-Yuan ZY Du Gang G Lin Yi-Cai YC
Journal of orthopaedic surgery and research 20211020 1
<h4>Background</h4>Osteoarthritis (OA) is the most common chronic degenerative joint disorder globally that is characterized by synovitis, cartilage degeneration, joint space stenosis, and sub-cartilage bone hyperplasia. However, the pathophysiologic mechanisms of OA have not been thoroughly investigated.<h4>Methods</h4>In this study, we conducted various bioinformatics analyses to identify hub biomarkers and immune infiltration in OA. The gene expression profiles of synovial tissues from 29 hea ...[more]
Similar Datasets
Project description:BackgroundEarly-onset schizophrenia (EOS) occurs between the ages of 13 and 17 years, and neurobiological factors leading to cognitive deficits and psychotic symptoms with varying degrees of positive and negative symptoms. Numerous studies have demonstrated a broad link between immune dysregulation and the central nervous system in EOS, and its pathogenesis involves immune dysfunction, but the exact biological mechanisms have not been elucidated. This study employs immune infiltration analysis and bioinformatics to unveil the pathogenic mechanisms of EOS and identify potential diagnostic biomarkers, aiming for more precise clinical interventions.MethodsIn this study, we recruited 26 EOS patients and 27 healthy controls (HCs), and microarray data were collected. Crossover genes were identified using weighted gene co-expression network analysis (WGCNA) and differential expression genes (DEGs) analysis. These genes were subjected to genome enrichment analysis (GSEA) and gene ontology (GO) analysis. Hub genes were identified through protein-protein interactions (PPIs) and the GeneMANIA database. The diagnostic potential of immune-associated hub genes was evaluated using ROC analysis. Immune infiltration in EOS was analyzed with CIBERSORT. Regulatory miRNAs for the hub genes were predicted using miRNet, and the correlation between mRNAs and miRNAs was analyzed and validated in clinical samples.ResultsBy WGCNA and DEGs analysis, 330 relevant genes were screened in EOS patients compared to HCs. Functional enrichment analysis using Metascape showed significant enrichment in immune system pathways. Subsequently, a PPI network was constructed to select the top 10 potential hub genes, and functional analysis was performed by GeneMANIA, resulting in the identification of four immune-related genes. In addition, significant differences were observed among the four immune cell types in the two groups of samples. ROC analysis showed clinical relevance of the immune-related hub genes, and the AUC of all genes was greater than 0.7. A miRNA-mRNA regulatory network was constructed from miRNA data, and three miRNAs were found to be significantly associated with the immune-related hub genes.ConclusionOur findings demonstrated that CCL3, IL1B, CXCL8, CXCL10 and miR-34a-5p may be biomarkers that play crucial roles in the underlying mechanisms of EOS immune-related pathways. These findings contribute to the understanding of EOS pathophysiology and may help identify new diagnostic and therapeutic targets.
| S-EPMC11748289 | biostudies-literature
Project description:Androgenetic alopecia (AGA) is a type of non-scarring hair loss. Current drugs for AGA are accompanied by adverse reactions and a high recurrence rate. Thus, the discovery of diagnostic biomarkers and therapeutic targets for AGA remains imperatively warranted. The GSE90594 dataset, which contained scalp skin biopsies from 14 male AGA cases and healthy volunteers, was used to identify the differentially expressed genes (DEGs). Functional enrichment analysis was subsequently performed. Next, the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database combined with the cytoHubba plugin of Cytoscape were used to obtain the key genes of AGA. Thereafter, the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm was performed to evaluate the relative abundance of immune cells between male AGA patients and healthy controls. The correlation between key genes and infiltrating immune cells was analyzed to obtain the significant immune-cell related genes (IRGs), then intersected with the DEGs between immortalized balding and non-balding human dermal papilla cells (DPCs) of the GSE93766 dataset as well as the DEGs obtained by the GSE90594 dataset, thus obtaining the hub genes of AGA. Finally, the hub genes were validated using GSE36169, which contained expression profiling of tissues biopsied from haired and bald scalps of five individuals with AGA. A total of 234 DEGs were obtained from the GSE90594 dataset, which were mainly enriched in the extracellular matrix (ECM)-related pathways and immune-related activities. The STRING database and ten algorithms in the cytoHubba plugin of Cytoscape disclosed 21 key DEGs. The results of the CIBERSORT algorithm revealed the relative abundances of 20 kinds of immune cells between diseased and healthy individuals, and yielded 15 IRGs involved in the pathogenesis of AGA. Next, the intersection analysis identified four hub genes of AGA, comprising COL1A2, PCOLCE, ITGAX, and LOX. The GSE36169 dataset validated the expression pattern of hub genes in the haired scalp of AGA patients. We discovered that the hub genes identified are closely linked with the causative factors of AGA, which could be used as the viable diagnostic and therapeutic target in the clinical applications.
| S-EPMC9761178 | biostudies-literature
Project description:Ulcerative colitis (UC) is a chronic inflammatory disease of the intestine, whose pathogenesis is not fully understood. Given that immune infiltration plays a key role in UC progression, our study aimed to assess the level of immune cells in UC intestinal mucosal tissues and identify potential immune-related genes. The GSE65114 UC dataset was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between healthy and UC tissues were identified using the "limma" package in R, while their Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were determined with the clusterProfiler package. Protein-protein interaction network analysis and visualization were performed with STRING and Cytoscape. Immune cell infiltration was calculated with CIBERSORT. The relationship between hub genes and immune-infiltrated cells in UC was determined by Pearson correlation. A total of 206 DEGs were identified, of which 174 were upregulated and 32 downregulated. GO and KEGG functional classification indicated DEG enrichment in immune response pathways, including Toll-like receptor signaling, IL-17 signaling, and immune system process and chemokine signaling. 13 hub genes were identified. Infiltration matrix analysis of immune cells showed abundant plasma cells, memory B cells, resting CD4 memory T cells, γδ T cells, M0 and M1 macrophages, and neutrophils in UC intestinal tissues. Correlation analysis revealed 13 hub genes associated with immune-infiltrated cells in UC. 13 hub genes associated with immune-infiltrated cells in UC were identified; they included CXCL13, CXCL10, CXCL9, CXCL8, CCL19, CTLA4, CCR1, CD69, CD163, IL7R, PECAM1, TLR8 and TLR2. These genes could potentially serve as markers for the diagnosis and treatment of UC.
| S-EPMC10101977 | biostudies-literature
Project description:ObjectivesAllergic rhinitis (AR) refers to a form of respiratory inflammation that mainly affects the sinonasal mucosa. The purpose of this study was to explore the level of immune cell infiltration and the pathogenesis of AR.MethodsWe performed a comprehensive analysis of two gene expression profiles (GSE50223 and GSE50101, a total of 30 patients with AR and 31 healthy controls). CIBERSORT was used to evaluate the immune cell infiltration levels. Weighted gene coexpression network analysis was applied to explore potential genes or gene modules related to immune status, and enrichment analyses including gene ontology, Kyoto Encyclopedia of Genes and Genomes, gene set enrichment analysis, and gene set variation analysis, were performed to analyze the potential mechanisms in AR. A protein-protein interaction network was constructed to investigate the hub genes, and consensus clustering was conducted to identify the molecular subtypes of AR.ResultsCompared to the healthy controls, patients with AR had high abundance levels and proportions of CD4+ memory-activated T cells. One hundred and eight immune-related differentially expressed genes were identified. Enrichment analysis suggested that AR was mainly related to leukocyte cell-cell adhesion, cytokine-cytokine receptor interaction, T-cell activation, and T-cell receptor signaling pathway. Ten hub genes, including TYROBP, CSF1R, TLR8, FCER1G, SPI1, ITGAM, CYBB, FCGR2A, CCR1, and HCK, which were related to immune response, might be crucial to the pathogenesis of AR. Three molecular subtypes with significantly different immune statuses were identified.ConclusionThis study improves our understanding of the molecular mechanisms in AR via comprehensive strategies and provides potential diagnostic biomarkers and therapeutic targets of AR.
| S-EPMC10696276 | biostudies-literature
Project description:Osteoarthritis (OA) is a progressive cartilage degradation disease, concomitant with synovitis, osteophyte formation, and subchondral bone sclerosis. Over 37% of the elderly population is affected by OA, and the number of cases is increasing as the global population ages. Therefore, the objective of this study was to identify and analyze the hub genes of OA combining with comprehensive bioinformatics analysis tools to provide theoretical basis in further OA effective therapies. Two sample sets of GSE46750 contained 12 pairs OA synovial membrane and normal samples harvested from patients as well as GSE98918 including 12 OA and non-OA patients were downloaded from the Gene Expression Omnibus database (GEO) database. Differentially expressed genes (DEGs) were identified using Gene Expression Omnibus 2R (GEO2R), followed by functional enrichment analysis, protein-protein interaction networks construction. The hub genes were identified and evaluated. An OA rat model was constructed, hematoxylin and eosin staining, safranin O/fast green staining, cytokines concentrations of serum were used to verify the model. The hub genes expression level in the knee OA samples were verified using RT-qPCR. The top 20 significantly up-regulated and down-regulated DEGs were screened out from the two datasets, respectively. The top 18 GO terms and 10 KEGG pathways were enriched. Eight hub genes were identified, namely MS4A6A, C1QB, C1QC, CD74, CSF1R, HLA-DPA1, HLA-DRA and ITGB2. Among them, the hub genes were all up-regulated in in vivo OA rat model, compared with healthy controls. The eight hub genes identified (MS4A6A, C1QB, C1QC, CD74, CSF1R, HLA-DPA1, HLA-DRA and ITGB2) were shown to be associated with OA. These genes can serve as disease markers to discriminate OA patients from healthy controls.
| S-EPMC10739719 | biostudies-literature
Project description:Background Osteoarthritis (OA) is the most common type of arthritis. OA can cause joint pain, stiffness, and loss of function. The pathogenesis of OA is not completely clear. Moreover, there is no effective treatment, and clinical management is limited to symptomatic relief or joint surgery. This study utilized bioinformatics to analyze normal and OA articular cartilage samples to find biomarkers and therapeutic targets for OA. Methods The GSE169077 gene chip dataset was downloaded from the public gene chip data platform of the National Biotechnology Information Center. The dataset included 6 samples of OA tissues and 5 samples of healthy cartilage tissues. Differentially expressed genes (DEGs) were screened using the R language “limma” function package under the threshold of log2[fold change (FC)] ≥2 and a P value <0.05. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) signal pathways of the target genes were enriched and analyzed using the database for annotation, visualization, and integrated discovery (DAVID), and a protein-protein interaction (PPI) network was further constructed using the search tool for the retrieval of interacting genes/proteins (STRING) database. The coexpression relationship of the genes in the module was visualized and screened with Cytoscape. Results A total of 27 DEGs were identified, including 9 downregulated genes and 18 upregulated genes. GO signal pathway enrichment analysis showed involvement in hypoxic response, fibrous collagen trimer, and extracellular matrix structural components. KEGG analysis demonstrated associations with protein digestion and absorption, extracellular matrix receptor interaction, and the peroxisome proliferator-activated receptor signal pathway, among several other pathways. A PPI network was obtained through STRING analysis, and the results were imported into Cytoscape software. The 27 DEGs were sequenced by the cytoHubba plug-in by various calculation methods, and 5 hub genes (COL1A1, COL1A2, POSTN, BMP1, and MMP13) were finally selected. These genes were analyzed by PPI again and annotated with GO and KEGG in different colors. Conclusions Bioinformatics technology effectively identified differential genes in the knee cartilage tissue of healthy controls and patients with OA, providing opportunities to further explore the mechanism and treatment of OA on a transcriptional level.
| S-EPMC9929772 | biostudies-literature
Project description:Heart failure (HF) and osteoarthritis (OA) are medical conditions that can significantly impact daily activities. Evidence has shown that HF and OA may share some pathogenic mechanisms. However, the underlying genomic mechanisms remain unclear. This study aimed to explore the underlying molecular mechanism and identify diagnostic biomarkers for HF and OA. With the cutoff criteria of fold change (FC) > 1.3 and P < .05, 920, 1500, 2195, and 2164 differentially expressed genes (DEGs) were identified in GSE57338, GSE116250, GSE114007, and GSE169077, respectively. After making the intersection of DEGs, we obtained 90 upregulated DEGs and 51 downregulated DEGs in HF datasets and 115 upregulated DEGs and 75 downregulated DEGs in OA datasets. Afterward, we conducted genome ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, protein-protein interaction (PPI) networks, and hub genes screening based on DEGs. Then, 4 common DEGs (fibroblast activation protein alpha [FAP], secreted frizzled-related protein 4 (SFRP4), Thy-1 cell surface antigen (THY1), matrix remodeling associated 5 [MXRA5]) between HF and OA were screened and validated in GSE5406 and GSE113825 datasets, based on which we established the support vector machine (SVM) models. The combined area under the receiver operating characteristic curve (AUC) of THY1, FAP, SFRP4, and MXRA5 in the HF training and test sets reached 0.949 and 0.928. While in the OA training set and test set, the combined AUC of THY1, FAP, SFRP4, and MXRA5 reached 1 and 1, respectively. The analysis of immune cells in HF revealed high levels of dendritic cell (DC), B cells, natural killer T cell (NKT), Type 1 regulatory T cell (Tr1), cytotoxic T cell (Tc), exhausted T cell (Tex), and mucosal-associated invariant T cell (MAIT), while displaying lower levels of monocytes, macrophages, NK, CD4 + T, gamma delta T (γδ T), T helper type 1 (Th1), T helper type 2 (Th2), and effector memory T cell (Tem). Moreover, the 4 common DEGs were positively correlated with DCs and B cells and negatively correlated with γδ T. In OA patients, the abundance of monocyte, macrophage, CD4 + naïve, and natural T regulatory cell (nTreg) was higher, while the infiltration of CD8 + T, γδ T, CD8 + naïve, and MAIT was lower. The expression of THY1 and FAP was significantly correlated with macrophage, CD8 + T, nTreg, and CD8 + naïve. SFRP4 was correlated with monocyte, CD8 + T, γδ T, CD4 + naïve, nTreg, CD8 + naïve and MAIT. MXRA5 was correlated with macrophage, CD8 + T, nTreg and CD8 + naïve. FAP, THY1, MXRA5, and SFRP4 may be diagnostic biomarkers for both HF and OA, and their correlation with immune cell infiltrations suggests shared immune pathogenesis.
| S-EPMC10313258 | biostudies-literature
Project description:BackgroundPsoriasis is a chronic, prolonged, and recurrent skin inflammatory disease. However, the pathogenesis of psoriasis is not completely clear, thus we aimed to explore potential molecular basis of it.MethodsTwo datasets were downloaded from the Gene Expression Omnibus database. After identifying the differentially expressed genes of psoriasis skin lesion samples and healthy controls, three kinds of analyses, namely functional annotation, protein-protein interaction (PPI) network, and immune infiltration analyses, were performed.ResultsA total of 152 up-regulated genes and 38 down-regulated genes were selected for subsequent analyses. Evaluation of the PPI network identified the most important module containing 13 hub genes. Gene ontology analysis showed that the hub genes have a significant enrichment effect on positive regulation of cell migration, defense response to the other organism and epithelial cell differentiation. KEGG signaling pathway analysis showed that the hub genes were significantly enriched in chemokine signaling, Toll-like receptor signaling pathway, and IL-17 signaling pathway. Compared with the normal control sample, naive B cells, CD8+ T cells, activated memory CD4+ T cells, follicular helper T cells, gamma delta T cells, resting NK cells, monocytes, M0 macrophages, M1 macrophages, activated dendritic cells and neutrophils infiltrated more, while memory B cells, naive CD4+ T cells, regulatory T cells (Tregs), activated NK cells, resting mast cells, and eosinophils infiltrated less.ConclusionTo conclude, the hub genes and pathways identified from psoriasis lesions and normal controls along with the immune infiltration profile may provide new insights into the study of psoriasis.
| S-EPMC7886814 | biostudies-literature
Project description:Osteoarthritis (OA) is a common cause of morbidity and disability worldwide. However, the pathogenesis of OA is unclear. Therefore, this study was conducted to characterize the pathogenesis and implicated genes of OA. The gene expression profiles of GSE82107 and GSE55235 were downloaded from the Gene Expression Omnibus database. Altogether, 173 differentially expressed genes including 68 upregulated genes and 105 downregulated genes in patients with OA were selected based on the criteria of ∣log fold-change | >1 and an adjusted p value < 0.05. Protein-protein interaction network analysis showed that FN1, COL1A1, IGF1, SPP1, TIMP1, BGN, COL5A1, MMP13, CLU, and SDC1 are the top ten genes most closely related to OA. Quantitative reverse transcription-polymerase chain reaction showed that the expression levels of COL1A1, COL5A1, TIMP1, MMP13, and SDC1 were significantly increased in OA. This study provides clues for the molecular mechanism and specific biomarkers of OA.
| S-EPMC7428874 | biostudies-literature
Project description:Osteoarthritis (OA) is one of the most common causes of disability and its development is associated with numerous factors. A major challenge in the treatment of OA is the lack of early diagnosis. In the present study, a bioinformatics method was employed to filter key genes that may be responsible for the pathogenesis of OA. From the Gene Expression Omnibus database, the datasets GSE55457, GSE12021 and GSE55325 were downloaded, which comprised 59 samples. Of these, 30 samples were from patients diagnosed with osteoarthritis and 29 were normal. Differentially expressed genes (DEGs) were obtained by downloading and analyzing the original data using bioinformatics. The Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathways were analyzed using the Database for Annotation, Visualization and Integrated Discovery online database. Protein-protein interaction network analysis was performed using the Search Tool for the Retrieval of Interacting Genes/proteins online database. BSCL2 lipid droplet biogenesis associated, seipin, FOS-like 2, activator protein-1 transcription factor subunit (FOSL2), cyclin-dependent kinase inhibitor 1A (CDKN1A) and kinectin 1 (KTN1) genes were identified as key genes by using Cytoscape software. Functional enrichment revealed that the DEGs were mainly accumulated in the ErbB, MAPK and PI3K-Akt pathways. Reverse transcription-quantitative PCR analysis confirmed a significant reduction in the expression levels of FOSL2, CDKN1A and KTN1 in OA samples. These genes have the potential to become novel diagnostic and therapeutic targets for OA.
| S-EPMC7903481 | biostudies-literature