整合机器学习与多层转录组分析将JAK2和ANXA5锁定为连接阻塞性睡眠呼吸暂停与氧化应激的关键基因，并在患者脂肪组织、间歇缺氧小鼠和CPAP治疗后样本中得到验证
核心概要
该研究整合limma差异分析、WGCNA、GeneCards氧化应激基因集、PPI网络与LASSO、随机森林、SVM-RFE三种机器学习方法，从阻塞性睡眠呼吸暂停（OSA）脂肪组织转录组中筛出57个共有差异基因并锁定JAK2与ANXA5两个核心基因，随后用单细胞测序、scTenifoldKnk虚拟敲除、免疫浸润解卷积以及RT-qPCR和Western blot临床样本验证，发现JAK2在OSA中显著上调、ANXA5显著下调，二者主要富集于单核细胞，且CPAP治疗后JAK2下降、ANXA5回升。
Schematic overview of the study design.
PubMed深度剖析
研究在OSA脂肪组织中鉴定出JAK2上调、ANXA5下调这一对方向相反的核心基因，并在独立数据集GSE135917中复现了相同趋势，ROC分析显示JAK2的AUC在训练集为0.838、验证集为0.790，ANXA5分别为0.775和0.919。 以往OSA与氧化应激的关联研究多聚焦单一分子通路，缺少多层转录组整合；该工作把差异表达、共表达模块与外部氧化应激基因集交叉，再用三种机器学习算法取交集，把候选范围收敛到两个基因。 基于GSE38792（10例OSA、8例对照）训练与GSE135917（34例OSA、8例对照）验证的转录组数据，配合ROC曲线与95%置信区间报告，属于计算层面的重复性证据。
单细胞分析显示JAK2与ANXA5在OSA样本的单核细胞中高表达，并识别出八种单核细胞亚型，其中Mono-inflammatory亚群在OSA患者中扩张、Mono-non-classical亚群相对减少，Mono-ISG亚群也发生改变。 该研究把两个核心基因的表达定位到具体免疫细胞亚群，并给出OSA中单核细胞亚型构成重塑的描述，而不仅是组织水平的基因表达差异。 基于人单细胞数据集GSE214865，经线粒体比例、基因数等质控后得到16个聚类并注释为14种主要细胞群，亚型身份由CD14、FCGR3A、CX3CR1、ISG15、IL1B等标记基因的点图与小提琴图支持。
scTenifoldKnk虚拟敲除JAK2或ANXA5后，S100P、ALOX5AP、PROK2、PADI4等免疫与炎症相关基因出现强扰动（log2FC > 6），敲除ANXA5还使HLA-DPB1、HLA-DPA1、HLA-DQA1、HLA-DRA等HLA基因显著失调，通路富集指向抗原加工与呈递、干扰素反应和炎症调控。 该研究用计算模拟的基因功能缺失预测下游调控网络，把ANXA5与抗原呈递过程联系起来，提出其此前未被充分认识的免疫调节功能。 属于基于单细胞数据的机器学习推断结果，作者以柱状图、火山图、MA图和点图多形式呈现，并报告富集显著性（如p < 1e-15、−log10(p) > 20），但为计算预测而非实验敲除。
临床与动物验证显示OSA患者脂肪组织中JAK2 mRNA和蛋白显著升高、ANXA5显著降低，间歇缺氧小鼠脂肪组织呈现同样方向，10例OSA患者经CPAP治疗后JAK2表达下降、ANXA5上调，且与非OSA肥胖及代谢肥胖组相比该失调在OSA患者中更突出。 该研究把生物信息学筛选结果推进到人脂肪组织、小鼠模型和治疗前后配对样本三个层面的实验验证，形成从计算到实验的闭环。 采用RT-qPCR与Western blot，样本来自上海新华医院耳鼻咽喉头颈外科40例受试者（20例OSA、20例对照）及每组5只的C57BL/6间歇缺氧小鼠模型，CPAP前后比较基于10例患者。
启示与展望
该研究为关注OSA分子机制、生物标志物和免疫代谢的研究者提供了可复用的分析路径：从公共转录组数据出发，经WGCNA与机器学习收敛候选基因，再用单细胞、虚拟敲除和实验验证逐层确认。其结论适用于脂肪组织（训练数据为内脏脂肪，验证样本为软腭脂肪）和单核细胞这一细胞背景，作者提出JAK2与ANXA5的联合检测有望作为PSG之外的辅助诊断思路，并提示JAK2可能作为CPAP不耐受患者的辅助治疗探索方向。
读者仍需关注几个开放问题：核心转录组队列样本量有限（训练集18例、验证集42例），临床验证为单中心40例，CPAP前后比较仅10例，这些规模对结论稳定性的影响有待更大规模、多中心前瞻性队列检验；虚拟敲除是计算推断，其预测的抗原呈递与干扰素通路变化尚需真实基因敲除或干预实验确认；研究聚焦脂肪组织与单核细胞，其他器官与细胞类型中的表现仍属未知；此外，本文为全文解析，图1至图16的具体数值与统计细节未在文本中逐项呈现，若需核对效应量应回到原文图表。