I.

The molecular mechanisms of osteoarthritis, the most common chronic disease, remain unexplained. This study aimed to use bioinformatic methods to identify the key biomarkers and immune infiltration in osteoarthritis. Gene expression profiles (GSE55235, GSE55457, GSE77298, and GSE82107) were selected from the Gene Expression Omnibus database. A protein-protein interaction network was created, and functional enrichment analysis and genomic enrichment analysis were performed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) databases. Immune cell infiltration between osteoarthritic tissues and control tissues was analyzed using the CIBERSORT method. Identify immune patterns using the ConsensusClusterPlus package in R software using a consistent clustering approach. Molecular biological investigations were performed to discover the important genes in cartilage cells. A total of 105 differentially expressed genes were identified. Differentially expressed genes were enriched in immunological response, chemokine-mediated signaling pathway, and inflammatory response revealed by the analysis of GO and KEGG databases. Two distinct immune patterns (ClusterA and ClusterB) were identified using the ConsensusClusterPlus. Cluster A patients had significantly lower resting dendritic cells, M2 macrophages, resting mast cells, activated natural killer cells and regulatory T cells than Cluster B patients. The expression levels of TCA1, TLR7, MMP9, CXCL10, CXCL13, HLA-DRA, and ADIPOQSPP1 were significantly higher in the IL-1β-induced group than in the osteoarthritis group in an in vitro qPCR experiment. Explaining the differences in immune infiltration between osteoarthritic tissues and normal tissues will contribute to the understanding of the development of osteoarthritis.

Background: Intervertebral disc degeneration (IDD), characterized by diverse pathological changes, causes low back pain (LBP). However, prophylactic and delaying treatments for IDD are limited. The aim of our study was to investigate the gene network and biomarkers of IDD and suggest potential therapeutic targets.Methods: Differentially expressed genes (DEGs) associated with IDD were identified by analyzing the mRNA, miRNA, and lncRNA expression profiles of IDD cases from the Gene Expression Omnibus (GEO). The protein–protein interaction (PPI) network, Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis as well as miRNA–lncRNA–mRNA networks were conducted. Moreover, we obtained 71 hub genes and performed a comprehensive analysis including GO, KEGG, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), Disease Ontology (DO), methylation analysis, receiver operating characteristic (ROC) curve analysis, immune infiltration analysis, and potential drug identification. We finally used qRT-PCR to verify 13 significant DEGs in normal and degenerative nucleus pulposus cells (NPCs).Results: We identified 305 DEGs closely related to IDD. The GO and KEGG analyses indicated that changes in IDD are significantly associated with enrichment of the inflammatory and immune response. GSEA analysis suggested that cell activation involved in the inflammatory immune response amide biosynthetic process was the key for the development of IDD. The GSVA suggested that DNA repair, oxidative phosphorylation, peroxisome, IL-6-JAK-STAT3 signaling, and apoptosis were crucial in the development of IDD. Among the 71 hub genes, the methylation levels of 11 genes were increased in IDD. A total of twenty genes showed a high functional similarity and diagnostic value in IDD. The result of the immune cell infiltration analysis indicated that seven genes were closely related to active natural killer cells. The most relevant targeted hub genes for potential drug or molecular compounds were MET and PIK3CD. Also, qRT-PCR results showed that ARHGAP27, C15orf39, DEPDC1, DHRSX, MGAM, SLC11A1, SMC4, and LINC00887 were significantly downregulated in degenerative NPCs; H19, LINC00685, mir-185-5p, and mir-4306 were upregulated in degenerative NPCs; and the expression level of mir-663a did not change significantly in normal and degenerative NPCs.Conclusion: Our findings may provide new insights into the functional characteristics and mechanism of IDD and aid the development of IDD therapeutics.