基于生物信息学分析鉴定和验证骨关节炎中线粒体功能障碍的核心基因

Identification and Verification of Hub Mitochondrial Dysfunction Genes in Osteoarthritis Based on Bioinformatics Analysis.

作者信息Hui Niu, Xingxing Deng, Qian Zhang, Yijun Zhao, Jinfeng Wen, Wenyu Li, Huan Liu, Xiong Guo, Cuiyan Wu
PMID38292759
发布时间2024-01-23
DOI10.1155/2024/6822664

实验完整度

中

包含生物信息学分析(差异表达、WGCNA、LASSO)和体外实验验证(qRT-PCR和Western blot),但缺乏体内模型或功能实验。

主要模型

人骨关节炎软骨细胞 正常软骨细胞

重点核对

GEO数据集GSE55235、GSE82107和GSE114007的样本类型和数量 DEMDFGs筛选阈值(adjusted P<0.05,|log2FC|>1) WGCNA软阈值β=4,MM>0.8,GS>0.5 qRT-PCR和Western blot验证的细胞来源和抗体稀释度

摘要

目的:年龄相关的线粒体功能障碍及其相关的氧化应激可能促进骨关节炎的发展。本研究旨在鉴定骨关节炎患者中线粒体功能障碍相关的核心基因,帮助预测OA风险,并揭示OA进展的机制。方法:从GEO(GSE55235、GSE82107和GSE114007)和GeneCard数据库下载OA表达数据和线粒体功能障碍基因。筛选OA与对照样本之间差异表达的线粒体功能障碍基因(DEMDFGs)。对DEMDFGs进行基因本体(GO)和京都基因与基因组百科全书(KEGG)通路分析。通过WGCNA和LASSO回归分析确定核心基因。ROC曲线显示了每个核心基因的诊断效能。构建并验证了预测OA风险的列线图模型。通过外部数据集、qRT-PCR和蛋白质印迹验证了核心基因在OA和正常软骨细胞中的表达。结果:共鉴定出31个DEMDFGs,其中15个基因上调,16个基因下调。GO功能富集分析显示,DEMDFGs富集在与能量代谢和细胞呼吸相关的生物过程中。通过加权基因共表达网络分析,我们确定了四个不同的共表达模块,其中蓝色模块与OA的相关性最强。DEMDFGs与该模块的交集产生了八个候选基因。在对数据进行LASSO分析后,确定了四个核心基因(ACADL、CYBA、SLC19A2和UCP2)作为OA的潜在生物标志物。这四个基因的表达水平在GSE114007数据集中得到外部验证。并且这四个基因在OA和正常软骨细胞中的生物学差异表达已得到验证。此外,通过ROC曲线分析,这四个核心基因具有良好的敏感性和特异性,用这四个基因构建的风险模型显示出良好的性能。总之,我们的研究可能提供新的线粒体功能障碍核心基因,对理解OA的病理、诊断和治疗具有潜在的临床应用价值。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
骨关节炎中与线粒体功能障碍相关的核心基因及其诊断价值如何?
核心机制
线粒体功能障碍导致ATP产生减少、氧化应激增加和脂质代谢紊乱,核心基因ACADL、CYBA、SLC19A2和UCP2可能通过影响这些过程促进OA进展。
主要证据
基于GEO数据集和GeneCards数据库,通过WGCNA和LASSO筛选出四个核心基因,外部数据集GSE114007和体外软骨细胞qRT-PCR及Western blot验证了其差异表达,ROC曲线显示良好的诊断效能。
研究意义
为理解OA的病理机制、诊断和治疗提供新的潜在生物标志物和靶点。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

数据获取与处理

获取OA基因表达数据并识别差异表达的线粒体功能障碍基因

从GEO数据库下载GSE55235、GSE82107和GSE114007数据集,合并GSE55235和GSE82107并去除批次效应,从GeneCards获取MDFGs,筛选DEMDFGs。

2

核心基因筛选

通过WGCNA和LASSO回归确定与OA相关的核心线粒体功能障碍基因

对合并数据进行WGCNA分析,选择软阈值β=4,确定蓝色模块与OA显著相关,将模块中候选基因与DEMDFGs取交集,再用LASSO回归进一步筛选。

3

诊断价值评估和风险模型构建

评估核心基因的诊断效能并构建预测OA风险的列线图模型

绘制ROC曲线计算各核心基因的AUC,使用rms包构建列线图模型,并用校准曲线评估预测性能。

4

通路富集和PPI网络分析

分析核心基因参与的通路和蛋白相互作用

使用GSEA进行KEGG富集分析,使用GeneMANIA构建PPI网络。

5

外部数据集验证

验证核心基因在独立数据集中的表达和预测能力

在GSE114007数据集中验证四个核心基因的表达差异和ROC曲线。

6

体外实验验证

验证核心基因在OA和正常软骨细胞中的mRNA和蛋白表达

从OA患者和正常供体获取软骨样本,分离培养软骨细胞,进行qRT-PCR和Western blot检测。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
胰蛋白酶HyClone--
II型胶原酶Gibco--
DMEM/F12培养基HyClone--
胎牛血清Gbico--
青霉素和链霉素Sigma--
SteadyPure通用RNA提取试剂盒IIAccurateAG21022
StarScript III逆转录预混液GenStar--
CFX96实时荧光定量PCR检测系统Bio-Rad--
SYBR Green预混液Pro Taq HS qPCR试剂盒AccurateAG11701
RIPA裂解液Epizyme--
PVDF膜----
ACADL抗体ImmunowayYT6498
CYBA抗体ImmunowayYN1733
SLC19A2抗体Biossbs-10738R
UCP2抗体ImmunowayYT4813
GAPDH抗体Proteintech10494-1-AP
山羊抗兔IgG二抗ProteintechSA00001-2
ECL化学发光试剂Epizyme--
化学发光成像系统Bio-Rad--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
数据获取与处理
GEO数据集GSE55235、GSE82107和GSE114007的样本数量,批次效应去除方法
阅读提示:Methods 2.1
差异表达基因筛选
DEMDFGs筛选阈值(adjusted P<0.05,|log2FC|>1)
阅读提示:Methods 2.2
WGCNA分析
软阈值β=4,SD>0.75,MM>0.8,GS>0.5
阅读提示:Methods 2.3
诊断模型构建
列线图模型的构建和验证方法
阅读提示:Methods 2.4
qRT-PCR验证
样本来源(3个OA患者和3个正常供体),qRT-PCR实验条件(引物序列、内参GAPDH、生物学重复和技术重复次数)
阅读提示:Methods 2.7, 2.8, Table 1
Western blot验证
抗体稀释度(ACADL 1:1000, CYBA 1:1000, SLC19A2 1:1000, UCP2 1:1000, GAPDH 1:10000),蛋白提取和检测方法
阅读提示:Methods 2.9