目的:通过生物信息学方法分析骨关节炎(OA)患者与正常人血清微小RNAs(miRNAs)表达差异。方法:从基因表达数据库(GEO)和ArrayExpress数据库中查找OA血清相关的芯片数据，并下载GSE105027数据集。筛选样本包括8例OA患者血清和12例正常(NC)人样本，其中女性患者5例，健康女性7例。用R语言Limma包分别筛选出男/女性OA和男/女性NC之间的差异表达miRNAs(DEMs)，设定阈值为 P<0.05和|log2FC|>1，用ggplot2包绘制火山图和韦恩图，得到2个DEMs：hsa-miR-33b-3p和hsa-miR-4284。用STRING和Cytoscape软件构件蛋白互作网络图(PPI)网络，依据LncBase和MNDR数据库筛选相关的长链非编码RNA(lncRNAs)，绘制内源竞争RNA(ceRNA网络互作图，Cytohubba筛选出关键基因，结合NONCODE和lncRNA SNP2数据库筛选与骨关节炎相关的lncRNAs，得出关键调控轴。 结果:共筛选出男、女性共有的DEMs 2个：hsa-miR-33b-3p和hsa-miR-4284，均为下调。结合ceRNA网络及数据库综合分析出lncRNA KCNQ1OT1-miR-33b-3p-PIK3CA轴。结论:通过生物信息学分析结果显示OA和NC的血清miRNAs表达差异集中在细胞凋亡、炎症和代谢，而lncRNA KCNQ1OT1-miR-33b-3p-PIK3CA轴在OA发病机制中尤为很重要。

作者：许炎；杨威；邢欢欢；裴洪

来源：中华实验外科杂志 2023 年 40卷 7期