No. 2 Department of Orthopaedics, Hebei Province Cangzhou Hospital of Integrated Traditional and Western Medicine (Cangzhou No. 2 Hospital), Cangzhou, Hebei 061001, P.R. China.
Department of Orthopaedics, Tianjin Medical University General Hospital, Tianjin 300052, P.R. China.
Mol Med Rep. 2017 Dec;16(6):8665-8672. doi: 10.3892/mmr.2017.7741. Epub 2017 Oct 9.
The present study aimed to reveal the potential genes associated with the pathogenesis of intervertebral disc degeneration (IDD) by analyzing microarray data using bioinformatics. Gene expression profiles of two regions of the intervertebral disc were compared between patients with IDD and controls. GSE70362 containing two groups of gene expression profiles, 16 nucleus pulposus (NP) samples from patients with IDD and 8 from controls, and 16 annulus fibrosus (AF) samples from patients with IDD and 8 from controls, was downloaded from the Gene Expression Omnibus database. A total of 93 and 114 differentially expressed genes (DEGs) were identified in NP and AF samples, respectively, using a limma software package for the R programming environment. Gene Ontology (GO) function enrichment analysis was performed to identify the associated biological functions of DEGs in IDD, which indicated that the DEGs may be involved in various processes, including cell adhesion, biological adhesion and extracellular matrix organization. Pathway enrichment analysis using the Kyoto Encyclopedia of Genes and Genomes (KEGG) demonstrated that the identified DEGs were potentially involved in focal adhesion and the p53 signaling pathway. Further analysis revealed that there were 35 common DEGs observed between the two regions (NP and AF), which may be further regulated by 6 clusters of microRNAs (miRNAs) retrieved with WebGestalt. The genes in the DEG‑miRNA regulatory network were annotated using GO function and KEGG pathway enrichment analysis, among which extracellular matrix organization was the most significant disrupted biological process and focal adhesion was the most significant dysregulated pathway. In addition, the result of protein‑protein interaction network modules demonstrated the involvement of inflammatory cytokine interferon signaling in IDD. These findings may not only advance the understanding of the pathogenesis of IDD, but also identify novel potential biomarkers for this disease.
本研究旨在通过生物信息学分析微阵列数据,揭示与椎间盘退行性变(IDD)发病机制相关的潜在基因。比较了 IDD 患者和对照组椎间盘两个区域的基因表达谱。从基因表达综合数据库(GEO)下载了包含两组基因表达谱的数据 GSE70362,其中一组为 16 个椎间盘核(NP)样本,分别来自 IDD 患者和对照组的 8 个样本;另一组为 16 个纤维环(AF)样本,分别来自 IDD 患者和对照组的 8 个样本。使用 R 编程语言中的 limma 软件包,分别在 NP 和 AF 样本中鉴定了 93 个和 114 个差异表达基因(DEG)。基因本体论(GO)功能富集分析用于鉴定 IDD 中 DEG 的相关生物学功能,结果表明 DEG 可能参与各种过程,包括细胞黏附、生物黏附和细胞外基质组织。京都基因与基因组百科全书(KEGG)通路富集分析表明,所鉴定的 DEG 可能与粘着斑和 p53 信号通路有关。进一步分析显示,在这两个区域(NP 和 AF)之间观察到 35 个共同的 DEG,这些基因可能进一步受到从 WebGestalt 检索到的 6 个 miRNA 簇的调节。使用 GO 功能和 KEGG 通路富集分析对 DEG-miRNA 调控网络中的基因进行注释,其中细胞外基质组织是最显著的失调生物学过程,粘着斑是最显著的失调途径。此外,蛋白质-蛋白质相互作用网络模块的结果表明,炎症细胞因子干扰素信号在 IDD 中发挥作用。这些发现不仅可以深入了解 IDD 的发病机制,还可以为该疾病鉴定新的潜在生物标志物。