《Frontiers in Immunology》：Integrative multi-omics analysis prioritizes compartment-specific candidate targets of gut microbiota metabolites in diabetic kidney disease
编辑推荐：
摘要
目的：肠-肾轴已成为一个关键研究领域。然而，肠道菌群代谢物在糖尿病肾病（diabetic kidney disease, DKD）肾小球与肾小管间质区室中的区室特异性调控机制仍不完全清楚。
方法：研究整合了gutMGene、SEA、STP、GEO
摘要
目的：肠-肾轴已成为一个关键研究领域。然而，肠道菌群代谢物在糖尿病肾病（diabetic kidney disease, DKD）肾小球与肾小管间质区室中的区室特异性调控机制仍不完全清楚。
方法：研究整合了gutMGene、SEA、STP、GEO、Nephroseq v5及KIT数据库，采用多组学与实验方法，包括差异表达分析、加权基因共表达网络分析（Weighted Gene Co-expression Network Analysis, WGCNA）、机器学习、孟德尔随机化（Mendelian randomization, MR）、基因集富集分析（Gene Set Enrichment Analysis, GSEA）、免疫浸润分析、基因集变异分析（gene set variation analysis, GSVA）、临床相关性分析、单细胞谱分析、分子对接、分子动力学模拟、CCK-8实验及RT-qPCR。研究人员针对每个区室构建了“微生物-底物-代谢物-靶点”（“Microbiota-Substrate-Metabolite-Target”, M-S-M-T）网络。
结果：机器学习鉴定出6个肾小球核心基因（IGFBP6、PLA2G4A、CTSK、HTR2B、PDGFRA、MMP7）和5个肾小管间质核心基因（CA2、HSD11B2、NQO2、MMP7、CYP24A1）。孟德尔随机化、单细胞谱分析及临床队列为这些基因与估算肾小球滤过率（estimated glomerular filtration rate, eGFR）及DKD的关联提供了支持性的遗传学和临床证据。免疫浸润分析和GSVA提示，在两种区室中，M2巨噬细胞的估计比例均较高，且肥大细胞的估计组成由活化态向静息态转变，这些变化与核心基因表达相关。通过“M-S-M-T”网络、药物相似性与毒性筛选，咖啡酸和柚皮素查尔酮被优先确定为共享候选代谢物。分子对接与分子动力学模拟提示其与HTR2B和MMP7具有有利的计算机模拟结合特性。两种代谢物均可减轻高糖（high glucose, HG）诱导的HK-2细胞和足细胞中核心基因的失调。
结论：本研究提出了DKD肾小球和肾小管间质中区室特异性的“M-S-M-T”调控网络。靶向HTR2B和MMP7的共享代谢物咖啡酸和柚皮素查尔酮在体外产生了初步的保护性转录反应。这些发现具有假设生成性质，在得出任何治疗学推论之前，需要在真实世界的微生物组、代谢组及体内研究中加以验证。
论文解读文章
**研究背景与科学问题**
糖尿病肾病（diabetic kidney disease, DKD）是终末期肾病（end-stage renal disease, ESRD）的主要病因。尽管近年来钠-葡萄糖共转运蛋白2（Sodium-Glucose Cotransporter 2, SGLT2）抑制剂、胰高血糖素样肽-1受体激动剂（Glucagon-Like Peptide-1 Receptor Agonists, GLP-1 RAs）及非甾体类盐皮质激素受体拮抗剂（non-steroidal Mineralocorticoid Receptor Antagonists, ns-MRAs）等药物在临床试验中显示出独立于降糖作用的心肾保护效应，但其精确机制尚未完全阐明。这提示器官间串扰（inter-organ crosstalk）可能是减缓DKD进展的关键治疗轴心。肠-肾轴（gut-kidney axis）作为DKD发病机制中重要的器官间通讯机制，已受到广泛关注。临床研究一致发现DKD患者粪便样本中存在显著的肠道菌群失调，且这些改变与脂质代谢、血糖稳态及肾功能等临床指标密切相关。功能上，微生物代谢物的异常可驱动疾病进程。然而，既往基于通用数据库（如GeneCards）的靶点识别方法无法充分解析DKD中关键病理事件究竟起源于肾小球还是肾小管间质区室。肾脏的细胞异质性对靶点解析构成了挑战。为此，研究人员利用区室特异性数据集和高级生物信息学方法，构建了分区室解析的调控网络，旨在阐明肠道菌群代谢物在DKD肾小球与肾小管间质中差异化的分子靶点及调控机制。
**研究内容与结论**
该研究首先从gutMGene数据库中筛选出248种肠道菌群代谢物，运用SwissTargetPrediction（STP）和Similarity Ensemble Approach（SEA）数据库预测其靶点并取交集，获得556个共享靶基因。利用来自基因表达综合数据库（GEO）的多个转录组数据集，分别对肾小球和肾小管间质进行差异表达分析和加权基因共表达网络分析（WGCNA），并通过机器学习（构建了基于12种算法的137个模型）结合SHAP分析，鉴定出肾小球区室的6个核心基因（IGFBP6、PLA2G4A、CTSK、HTR2B、PDGFRA、MMP7）和肾小管间质区室的5个核心基因（CA2、HSD11B2、NQO2、MMP7、CYP24A1）。双样本孟德尔随机化（MR）分析利用eQTL、pQTL和mQTL遗传工具，为核心基因与估算肾小球滤过率（eGFR）及DKD相关肾脏结局的关联提供了支持性遗传证据。单细胞分析揭示了核心基因在DKD各肾脏细胞类型中的特异性表达变化。临床队列（Nephroseq v5）相关性分析显示，肾小球核心基因与eGFR呈负相关，而肾小管间质中的CA2、HSD11B2、NQO2与eGFR呈正相关，MMP7和CYP24A1与eGFR呈负相关。免疫浸润分析和基因集变异分析（GSVA）一致提示DKD两种区室中M2巨噬细胞估计比例升高，肥大细胞估计组成由活化态向静息态转变，且与核心基因表达相关。基于核心基因构建了区室特异性的“微生物-底物-代谢物-靶点”（“Microbiota-Substrate-Metabolite-Target”, M-S-M-T）网络。通过药物相似性（Lipinski五规则）和毒性预测筛选，咖啡酸（caffeic acid）和柚皮素查尔酮（naringenin chalcone）被确定为两个区室共享的候选代谢物。分子对接和分子动力学模拟提示这两种代谢物与HTR2B和MMP7具有稳定的结合构象及合理的结合自由能。体外细胞实验证实，在高糖（HG）诱导的HK-2细胞和足细胞中，咖啡酸和柚皮素查尔酮可逆转HG引起的核心基因转录失调，并通过HTR2B和MMP-7依赖性途径恢复Bcl-2表达，发挥初步的抗凋亡保护效应。该论文发表在《Frontiers in Immunology》。
**主要关键技术方法**
该研究整合了多种计算与实验技术。数据来源方面，利用了GEO数据库中的10个转录组数据集（肾小球：训练集GSE104948、GSE96804、GSE47183；验证集GSE99339；测试集GSE30528、GSE142025。肾小管间质：训练集GSE104954、GSE30529；验证集GSE99325；测试集GSE47184），以及gutMGene、SEA、STP、Nephroseq v5、KIT等数据库。核心方法包括：limma差异表达分析、WGCNA共表达网络构建、基于12种算法的137个机器学习模型与SHAP特征重要性分析、基于eQTL/pQTL/mQTL的双样本孟德尔随机化(MR)、单基因GSEA、CIBERSORT免疫浸润反卷积、GSVA通路活性评分、单细胞RNA测序数据分析、分子对接、分子动力学模拟、MM/GBSA结合自由能计算、CCK-8细胞活力检测及RT-qPCR基因表达验证。体外实验使用HK-2人肾小管上皮细胞和条件永生化的足细胞（HPC）。
**研究结果分述**
**3.1 与肠道菌群代谢物相关的DKD肾小球和肾小管间质靶点的逐步筛选：** 通过差异表达分析、WGCNA及与代谢物靶基因的逐步取交集，分别获得肾小球区室9个和肾小管间质区室15个关键基因。PPI网络及GO/KEGG富集分析提示这些基因富集于GnRH信号通路、Fc gamma R介导的吞噬作用、TRP通道炎症介质调节、血清素突触及肾素-血管紧张素系统等通路，提示两区室存在共享病理机制。
**3.2 具有跨队列泛化能力的核心基因机器学习筛选：** 通过多模型策略和独立测试集验证，确定肾小球核心基因（IGFBP6, CTSK, PLA2G4A, HTR2B, PDGFRA, MMP7）和肾小管间质核心基因（MMP7, HSD11B2, CYP24A1, CA2, NQO2）。SHAP分析显示IGFBP6和CTSK是肾小球模型的主要贡献者，MMP7是肾小管间质模型的主要贡献者。
**3.3 来自MR的支持性遗传证据：** MR分析表明NQO2、CTSK、PLA2G4A、HTR2B的遗传预测表达或水平与eGFR降低相关；IGFBP6与慢性肾脏病及肾衰竭风险相关；HSD11B2甲基化保护eGFR，CYP24A1高表达对DKD具有保护作用；PDGFRA表现出分期特异性效应。
**3.4 核心基因的免疫相关特征的计算推断：** 单基因GSEA显示核心基因高表达与免疫和纤维化通路正相关，与代谢通路负相关。免疫浸润分析显示DKD中M2巨噬细胞和静息肥大细胞浸润分数升高，活化肥大细胞降低。GSVA进一步显示核心基因表达与M2巨噬细胞极化及肥大细胞脱颗粒相关通路富集分数共变。
**3.5 核心基因的单细胞分析及与肾功能的相关性：** 单细胞数据显示核心基因在DKD肾脏各细胞类型（如系膜细胞、壁层上皮细胞、成纤维细胞、近曲小管、髓袢、集合管等）中存在区室和细胞类型特异性的表达改变。临床队列相关性分析验证了核心基因表达与eGFR的显著相关关系。
**3.6 “M-S-M-T”假设生成框架：** 构建了肾小球和肾小管间质特异性的“M-S-M-T”网络，分别涉及79种和116种肠道菌群，25种和52种底物，33种和75种代谢物，为肠-肾交互提供了预测性框架。
**3.7 药物相似性与毒性分析：** 通过Lipinski规则和ADMETlab 3.0毒性预测，从两区室网络中筛选出共享的且药物相似性良好、无预测肾毒性/致癌性的代谢物：咖啡酸和柚皮素查尔酮。
**3.8 分子对接：** 分子对接显示柚皮素查尔酮与HTR2B（-8.602 kcal/mol）和MMP7具有较好的亲和力，咖啡酸与HTR2B结合亲和力为-6.897 kcal/mol，分别与关键残基形成氢键和π-堆积相互作用。
**3.9 分子动力学模拟：** MM/GBSA计算显示HTR2B-柚皮素查尔酮和HTR2B-咖啡酸复合物的结合自由能（ΔGbind）分别为-108.19±11.67 kJ/mol和-115.37±7.78 kJ/mol，显著稳定于MMP7-柚皮素查尔酮复合物（-40.54±10.91 kJ/mol）。RMSD、RMSF、SASA及回旋半径分析表明三种复合物在100 ns模拟中总体构象稳定。
**3.10 柚皮素查尔酮和咖啡酸的体外治疗效果及靶点依赖性：** 在HG诱导的HK-2细胞和足细胞中，两种代谢物均能逆转HG引起的核心基因mRNA水平的失调（HG导致肾小管间质核心基因下调，而代谢物使其恢复；HG导致足细胞核心基因下调，而代谢物使其恢复）。Western blot证实两种代谢物可恢复HG抑制的MMP-7和HTR2B蛋白表达以及Bcl-2抗凋亡蛋白水平。使用选择性HTR2B抑制剂RS-127445和MMP-7抑制剂marimastat可消除代谢物的保护效应，表明其作用依赖于这些靶点。
**讨论与结论**
讨论部分指出，本研究构建的肾小球和肾小管间质核心基因谱具有生物学合理性：肾小球核心基因反映了机械感觉-蛋白酶激活-炎症-纤维化级联，而肾小管间质核心基因则与重吸收和酸碱平衡生理学密切相关。MMP7作为两区室共享节点，参与基底膜重塑。免疫相关分析提示肠道菌群代谢物可能作为系统性信号分子扰动核心基因，从而促进DKD的促纤维化免疫微环境。与已报道的临床多组学数据比较，该网络能部分复现临床DKD中胆汁酸和芳香族氨基酸代谢物的改变，支持其框架合理性。咖啡酸和柚皮素查尔酮显示出初步的保护性转录反应，但其口服生物利用度、化学稳定性及临床转化尚存在障碍。作者也明确指出了研究的局限性：MR分析基于欧洲人群血液来源的遗传工具，缺乏组织特异性；M-S-M-T网络为基于数据库的预测，缺乏本研究自身的微生物组和代谢组学实测数据；免疫浸润结果基于反卷积推断，属相关性而非因果性；体外实验仅限于急性高糖暴露的两种细胞系。
研究结论总结为：本研究鉴定出DKD肾小球和肾小管间质特异性核心基因，通过MR、免疫浸润分析、单细胞转录组学和临床队列探索了其在DKD中的作用与机制，构建了相应的区室特异性“M-S-M-T”网络，并通过分子对接、分子动力学模拟和体外实验，将咖啡酸和柚皮素查尔酮确定为值得进一步验证的共享候选代谢物。这些发现为DKD中的肠-肾轴建立了假设生成框架，并为临床微生物组和代谢组学研究定义了具体的验证优先事项。
