利用加权基因共表达网络分析和机器学习筛选神经损伤的生物标志物

Screening Biomarkers for Nerve Injury Using Weighted Gene Co-Expression Network Analysis and Machine Learning.

作者信息Shuming Cao, Chengyue Yu, Nana Wang, Jianhua Xu, Weiguo Xu
PMID41733593
发布时间2026-02
DOI10.1002/brb3.71279

实验完整度

包含生物信息学分析(WGCNA、机器学习)、独立数据集验证、免疫浸润分析、单细胞测序以及体外细胞实验(PC12细胞敲低Fcgr2b验证功能),但主要验证集中在生信层面,体内实验缺失。

主要模型

PC12细胞系 大鼠神经损伤与假手术样本 背根神经节单细胞数据

重点核对

训练集合并GSE236754、GSE202111、GSE172133,使用sva去除批次效应 验证集包括GSE102721、GSE175760等独立数据集 WGCNA软阈值β=5 差异表达基因筛选阈值|Log2FC|>0.5且p<0.05 PC12细胞用100 ng/mL NGF刺激0、1、3、5天

摘要

背景:神经损伤触发复杂的分子反应,涉及免疫激活和神经元损伤,然而关键调控基因及其机制仍知之甚少。在此,我们整合了多转录组数据集和机器学习,以识别和验证神经损伤的新型生物标志物,并阐明其功能作用。方法:神经损伤和假手术样本的RNA-seq数据和单细胞转录组数据来源于基因表达综合(GEO)数据库。使用加权基因共表达网络分析(WGCNA)、差异表达分析和三种机器学习算法来识别与神经损伤相关的hub基因。在独立数据集中验证了这些hub基因的表达模式和诊断价值。使用CIBERSORT算法分析了这些基因与免疫细胞浸润的相关性。最后,利用单细胞RNA测序(scRNA-seq)数据研究了hub基因在神经细胞中的细胞特异性表达模式。结果:通过WGCNA和三种机器学习方法鉴定了七个神经损伤相关基因,其中Atf3、Bin2、Fcgr2b和Ucn在验证队列中表现出稳健的诊断性能(AUC > 0.7)。功能富集表明这些基因参与神经炎症、神经元命运决定和JAK-STAT/NF-κB信号通路。免疫浸润分析将其表达与M2巨噬细胞极化和CD4+ T细胞耗竭相关联,而scRNA-seq突出了细胞特异性模式:Atf3和Ucn在神经元中富集,而Fcgr2b和Bin2在巨噬细胞/NK细胞中占主导。此外,Fcgr2b促进了PC12细胞中神经突的生长。结论:我们的研究揭示了Atf3、Bin2、Fcgr2b和Ucn是关键的神经损伤生物标志物,在神经免疫串扰中具有双重作用,为神经修复的治疗靶向提供了新的见解。此外,Fcgr2b可能参与神经损伤后的神经突生长。

实验结论

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

研究问题
神经损伤中关键的调控基因及其在神经免疫串扰中的功能机制是什么?
核心机制
Atf3、Bin2、Fcgr2b和Ucn可能通过JAK-STAT和NF-κB信号通路参与神经炎症和免疫调节,与M2巨噬细胞极化正相关,与CD4+ T细胞耗竭负相关。
主要证据
在训练集和验证集中,Atf3、Bin2、Fcgr2b和Ucn的AUC均大于0.7;GSEA显示富集于神经炎症、神经元命运决定及JAK-STAT/NF-κB通路;scRNA-seq显示Atf3和Ucn在神经元中高表达,Fcgr2b和Bin2在巨噬细胞/NK细胞中高表达;体外敲低Fcgr2b抑制NGF诱导的PC12细胞神经突生长及NF200表达。
研究意义
本研究为神经损伤的诊断和治疗提供了潜在生物标志物,特别是Fcgr2b可能在神经损伤后的神经突生长中发挥作用,为神经修复的治疗靶点提供了新见解。

研究路径

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

1

数据获取与预处理

整合多数据集并去除批次效应,为后续分析提供可靠数据。

从GEO数据库下载RNA-seq数据集GSE236754、GSE202111、GSE172133作为训练集,GSE102721、GSE175760作为验证集,利用sva包去除批次效应。

2

WGCNA和差异表达分析

识别与神经损伤相关的共表达模块和差异表达基因。

使用WGCNA构建基因共表达网络,筛选出与神经损伤显著相关的turquoise模块(1126个基因);同时利用limma包进行差异表达分析,得到71个DEGs;取交集获得42个候选基因。

3

机器学习筛选hub基因

通过多种算法缩小候选基因范围,筛选具有诊断潜力的核心基因。

使用LASSO、SVM-RFE和随机森林三种机器学习算法对候选基因进行筛选,取交集获得7个hub基因,进一步验证后确定Atf3、Bin2、Fcgr2b和Ucn为关键基因。

4

验证hub基因表达与诊断价值

在独立数据集中验证关键基因的表达差异和诊断价值。

在训练集和验证集GSE102721、GSE175760中比较hub基因在神经损伤和假手术组的表达,并绘制ROC曲线计算AUC。

5

功能富集和免疫浸润分析

探究hub基因的生物学功能和免疫细胞浸润相关性。

使用GSEA对四个关键基因进行GO和KEGG富集分析;使用CIBERSORT计算22种免疫细胞比例,并分析关键基因表达与免疫浸润的相关性。

6

单细胞表达模式分析

在单细胞水平上研究关键基因的细胞类型特异性表达。

对GSE255014数据集进行scRNA-seq分析,聚类和注释细胞类型,并可视化关键基因的表达。

7

Fcgr2b功能验证实验

验证Fcgr2b对神经突生长的影响。

在PC12细胞中敲低Fcgr2b,用NGF刺激后观察神经突长度变化,并通过RT-qPCR和Western blot检测Fcgr2b和NF200的表达。

研究方法

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

产品清单

实验环节名称品牌货号
RPMI-1640培养基IMMOCELLPM150110
胎牛血清ProcellPM164210
马血清IMMOCELLIMC-108
青霉素-链霉素ProcellP1400
Lipofectamine 2000转染试剂Invitrogen--
神经生长因子Thermo450-34-20UG
TriQuick总RNA提取试剂SolarbioR1100
Evo M-MLV逆转录预混液Accurate Biotechnology (Hunan) Co., Ltd.AG11706-S
2×SuperStar通用SYBR预混液CWBIOCW3360
PCR仪(SLAN-96S)Shanghai Hongshi Medical Technology Co., LtdSLAN-96S
RIPA裂解液SolarbioR0010
BCA蛋白定量试剂盒CWBIOCW0014S
β-肌动蛋白抗体Proteintech--
Fcgr2b抗体InvitrogenMA5-47232
NF200抗体Proteintech60331-1-Ig
辣根过氧化物酶标记的山羊抗兔IgG(H+L)(亲和纯化)Beijing Zhongshan Jinqiao Biological Technology Co., Ltd.ZB-2301
辣根过氧化物酶标记的山羊抗兔IgG(H+L)(亲和纯化)Beijing Zhongshan Jinqiao Biological Technology Co., Ltd.ZB-2305
全自动化学发光图像分析系统ClinxChemi6000

关键环节

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

环节核对要点
数据获取与预处理
训练集组成(GSE236754, GSE202111, GSE172133)和验证集(GSE102721, GSE175760);批次效应去除方法(sva包,n.sv=2)
阅读提示:Materials and Methods 2.1 数据获取
WGCNA和差异表达分析
WGCNA软阈值β=5;差异表达基因筛选阈值|Log2FC|>0.5且p<0.05
阅读提示:Materials and Methods 2.2 和 2.3
机器学习筛选hub基因
LASSO、SVM-RFE、随机森林的具体参数(如lambda.min、mtry=2、ntree=500等)
阅读提示:Materials and Methods 2.5 多种机器学习算法
Fcgr2b功能验证实验
PC12细胞来源(IMMOCELL, IM-R014);NGF浓度(100 ng/mL)和刺激时间(0、1、3、5天);siRNA序列和转染试剂(Lipofectamine 2000);检测方法(RT-qPCR、Western blot)
阅读提示:Materials and Methods 2.9 细胞培养与转染、2.10 RT-qPCR、2.11 Western blot;Results 3.8