《International Journal of Molecular Sciences》：Single-Cell Transcriptomic and Regulatory Network Analyses Identify a Putative HNF4A–OXR1 Association Linked to Podocyte Mitochondrial Dysfunction in Diabetic Kidney Disease
编辑推荐：
足细胞损伤是糖尿病肾病（DKD）的核心病理事件，但足细胞功能障碍背后的线粒体调控机制仍未完全阐明。研究人员整合了GSE195460中人肾活检组织的单核RNA测序数据，并结合足细胞特异性差异表达分析、功能富集分析、线粒体基因注释、共表达网络分析、机器学习特征选择
足细胞损伤是糖尿病肾病（DKD）的核心病理事件，但足细胞功能障碍背后的线粒体调控机制仍未完全阐明。研究人员整合了GSE195460中人肾活检组织的单核RNA测序数据，并结合足细胞特异性差异表达分析、功能富集分析、线粒体基因注释、共表达网络分析、机器学习特征选择、拟时序分析、计算机模拟敲除、转录因子富集分析、基序分析以及分子对接；研究结果在GSE209781中进行了外部验证。足细胞表现出显著的相关疾病转录组扰动，共鉴定出677个探索性细胞水平差异表达基因，并富集了细胞骨架、黏附、发育、氧化磷酸化（OXPHOS）及代谢相关通路。共鉴定出31个线粒体相关差异表达基因，其中OXR1基于一致的表达下调、网络中心性、机器学习筛选及外部验证被优先确定为候选调控枢纽基因。HNF4A被预测为与OXR1及其他线粒体靶基因相关联的上游调控因子。计算机模拟OXR1扰动使预测的损伤状态概率从15%增加至82%，而启动子基序分析进一步为HNF4A–OXR1的潜在调控关系提供了计算层面的支持。上述发现突出了OXR1作为DKD足细胞线粒体功能障碍相关候选基因的重要性，并提示了HNF4A–OXR1的潜在调控关系。
论文解读
一、研究背景与科学问题
糖尿病肾病（DKD）是糖尿病的重要微血管并发症，也是全球终末期肾病的主要病因。足细胞损伤和丢失是DKD的关键病理特征，直接导致蛋白尿和进行性肾功能障碍。尽管目前肾素–血管紧张素–醛固酮系统阻断剂和钠–葡萄糖协同转运蛋白2抑制剂等疗法已改善肾脏结局，但仍有大量患者面临疾病进展风险。线粒体功能障碍已被认为是DKD足细胞损伤的重要驱动因素——足细胞具有高能量需求，氧化磷酸化、脂肪酸氧化、氧化还原平衡及线粒体稳态的紊乱可能损害其结构和功能完整性。然而，驱动这些线粒体改变的上游调控机制仍不清楚。此外，传统的批量转录组研究可能掩盖肾脏异质微环境中细胞类型特异性的信号，限制了对足细胞内在调控网络的识别。因此，本研究旨在利用单细胞转录组学结合多层级计算分析，系统刻画DKD足细胞特异性的线粒体转录改变，并识别潜在的调控基因和调控网络。
二、研究设计与主要结论
研究人员整合了人肾活检组织的单核RNA测序数据（GSE195460作为发现队列），结合线粒体基因注释、机器学习特征选择、共表达网络分析、拟时序分析、计算机模拟敲除、转录因子富集、基序分析及分子对接等多层次计算策略，并在GSE209781和GSE271188两个独立队列中进行外部验证。研究发现足细胞是DKD中转录组改变最显著的肾细胞群体之一，OXR1被多维度一致地优先确定为候选线粒体调控枢纽基因，HNF4A被预测为其潜在上游调控因子。该研究为理解DKD足细胞线粒体功能障碍的调控机制提供了网络层面的框架，并提出了可供后续实验验证的HNF4A–OXR1调控假说。
三、主要关键技术方法
研究人员使用了以下关键技术方法：整合GSE195460人肾单核RNA测序数据，经质量控制（线粒体转录本阈值5%）、SCTransform归一化和Harmony批次校正后，基于典型谱系标志物注释八种主要肾细胞群体并提取足细胞；采用Wilcoxon秩和检验进行探索性细胞水平差异表达分析，并以DESeq2进行供体水平伪批量分析；利用MitoCarta 3.0线粒体基因清单鉴定线粒体相关差异表达基因；使用clusterProfiler进行GO/KEGG富集分析和GSEA；采用随机森林（RF）、支持向量机–递归特征消除（SVM-RFE）和LASSO三种机器学习方法进行特征选择；基于Pearson相关构建共表达网络并计算度、介数中心性、接近中心性和PageRank等拓扑指标；使用ChEA3进行转录因子富集分析；采用Slingshot进行拟时序轨迹推断；使用scTenifoldKnk进行计算机模拟OXR1扰动；利用JASPAR数据库进行HNF4A基序筛选；使用CB-Dock2（AutoDock Vina评分）进行PGE1–EP2分子对接。外部验证使用GSE209781（3例DKD和3例参考供体）和GSE271188（7例对照和4例DKD供体）队列。
四、研究结果
4.1 单细胞转录组谱与质量控制敏感性分析
整合分析鉴定出八种主要肾细胞群体。转录组距离分析显示DKD相关差异在肾细胞群体中广泛存在，足细胞是改变最显著的谱系之一。质量控制敏感性分析显示，在5%线粒体转录本阈值下，OXR1在对照和DKD足细胞中的平均表达分别为1.14和0.69（log2FC = ?0.724；p = 0.021），10%阈值下结果一致（log2FC = ?0.691；p = 0.028）。
4.2 DKD足细胞的差异表达与功能富集
细胞水平差异表达分析鉴定出677个DKD相关基因（337个上调、340个下调）。患者水平伪批量分析鉴定出18个显著基因，OXR1在患者水平也呈现显著下调（log2FC = ?0.735；BH-FDR = 0.0086）。KEGG富集显示IgSF CAM信号和细胞骨架相关通路富集；GO富集显示肌动蛋白丝相关过程调控和细胞–基质黏附富集。GSEA显示DKD足细胞中氧化磷酸化通路（NES = ?2.17）和代谢通路基因集（NES = ?1.97）均呈负富集。
4.3 DKD足细胞中的线粒体基因表达与功能模块评分
与MitoCarta 3.0清单交叉后鉴定出31个线粒体相关差异表达基因（10个上调、21个下调）。上调基因包括MT-CO2（平均log2FC = 1.529）、MT-ND4L（0.545）和MT-ND4（0.534）；下调基因包括ALDH6A1（?2.037）、GLDC（?1.568）、ACSM2B（?1.393）、ACSM2A（?1.041）和OXR1（?0.717）。模块评分分析显示DKD足细胞中OXPHOS评分显著升高（0.3689 vs. 0.1769；p = 8.88 × 10?17），而TCA循环评分（?0.0247 vs. 0.0173；p = 6.92 × 10?5）和脂肪酸氧化评分（?0.0341 vs. 0.0078；p = 1.28 × 10?3）均显著降低。GO分析显示线粒体相关差异表达基因富集于酰基辅酶A代谢过程（富集倍数 = 47.94）。
4.4 线粒体相关差异表达基因的机器学习与网络分析
三种机器学习方法的交集鉴定出8个共识候选基因：ACSM2A、ACSM2B、FHIT、OXR1、MT-CO2、MT-CO3、MT-ND4和MT-ND4L。MT-CO2在随机森林和SVM-RFE中均排名第一。OXR1在随机森林中排名第四，在SVM-RFE中排名第八，并被LASSO保留。网络分析显示OXR1的度为16，介数中心性为50，接近中心性为0.1881，PageRank为0.0361，在31个线粒体相关差异表达基因中具有最高的介数和接近中心性。
4.5 独立单细胞队列中的OXR1表达
在GSE209781中，DKD足细胞中OXR1平均表达为0.504，参考足细胞为0.713，DKD组约低29%。在GSE271188中，排除足细胞数少于10个的供体后，DKD组OXR1平均表达为0.633，对照组为1.119（p = 0.0238）。
4.6 转录因子富集与预测调控网络分析
ChEA3分析将HNF1B排为第一（综合平均秩分 = 9.33），HNF4A排为第六（18.67）。预测的转录因子–靶基因网络将HNF4A与OXR1、ACSM2A、ACSM2B、FHIT及其他足细胞差异表达基因相关联。
4.7 足细胞轨迹与OXR1扰动的预测损伤相关效应
Slingshot分析纳入1574个足细胞核，分布于5个亚簇（C0–C4），推断出两条转录谱系。OXR1表达与两条谱系的推断拟时序均呈弱负相关（谱系1：ρ = ?0.0829；谱系2：ρ = ?0.0658）。DKD足细胞中OXR1平均表达在两条谱系中均低于对照。计算机模拟OXR1扰动后，预测的损伤状态概率从基线15%增加至82%。调控距离分析鉴定出6个FDR < 0.05的基因：OXR1（d = 2.77 × 10?4）、CERS6、ZDHHC6、SRGAP2B、CBLB和MAGI2。
4.8 预测的HNF4A基序与PGE1–EP2相互作用
基序分析在OXR1上游调控区鉴定出推定的HNF4A结合序列CAAAGTCA，匹配HNF4A矩阵MA0114.5（基序得分 = 12.11742，相对得分 = 0.9098152）。分子对接将PGE1定位于人EP2受体的预测结合口袋内，建模的配体–受体接触约为2.8–3.2 ?。
五、讨论与结论
讨论部分指出，OXR1的优先排序并非仅基于差异表达，而是由线粒体基因筛选、共识机器学习特征选择、网络中心性分析和外部单细胞验证等多个独立分析层面共同支持，患者水平伪批量分析进一步强化了这一证据。然而，有限的生物学重复和研究的计算性质要求谨慎解读。拟时序分析推断的两条谱系结构对根簇选择敏感，反映的是推断的转录状态关系而非时间上的DKD进展。计算机模拟OXR1扰动提示了网络层面的潜在后果，但仍属假说生成性质，不能确立因果关系。HNF4A作为候选上游调控因子的证据来自转录因子富集、调控网络连接和启动子基序分析的多重汇聚，但尚不能确立直接转录调控关系。PGE1–EP2对接结果应被解读为结构兼容性的假说生成证据，而非连续信号通路的证据。研究局限性包括：主要基于公开数据集的计算机分析、单细胞测序无法直接评估蛋白表达和线粒体呼吸功能、计算机模拟敲除不能替代实验性功能研究、肾切除来源参考组织的潜在混杂、外部队列样本量小、线粒体过滤可能排除部分应激细胞群体，以及个体间异质性和数据集技术效应可能影响结果普适性。
结论部分翻译：OXR1被确定为与DKD足细胞线粒体功能障碍相关的候选基因，HNF4A被预测为潜在上游调控因子。计算扰动和启动子基序分析为OXR1的生物学相关性及HNF4A–OXR1的潜在调控关系提供了假说生成层面的支持。需要进一步的实验研究来确定OXR1在DKD足细胞中的功能作用和调控机制。
