《Frontiers in Endocrinology》：Tissue-resolved transcriptomic and proteomic profiling reveals organ-specific molecular signatures associated with inter-animal variation in glycemic response to metformin in obese diabetic Mongolian gerbils
编辑推荐：
背景：2型糖尿病（Type 2 diabetes mellitus，T2DM）影响全球数亿人。二甲双胍仍是肥胖T2DM的一线治疗药物；然而，30%–40%的患者单药治疗无法达到充分的血糖控制，HbA1c降幅在个体间存在2至3倍差异，且这种血糖应答个体差异背后的
背景：2型糖尿病（Type 2 diabetes mellitus，T2DM）影响全球数亿人。二甲双胍仍是肥胖T2DM的一线治疗药物；然而，30%–40%的患者单药治疗无法达到充分的血糖控制，HbA1c降幅在个体间存在2至3倍差异，且这种血糖应答个体差异背后的多器官协调机制仍不清楚。
方法：研究人员对肥胖糖尿病蒙古沙鼠给予二甲双胍（100 mg/kg/天，口服灌胃；n = 4）或溶媒对照（n = 3）处理5周。研究人员对肝脏、骨骼肌和腹股沟白色脂肪组织（inguinal white adipose tissue，iWAT）进行了匹配的转录组学（RNA-seq）和数据非依赖采集蛋白质组学（data-independent acquisition，DIA）分析。跨组学分析锚定血糖参数，关键肝脏发现通过实时定量聚合酶链反应（quantitative real-time polymerase chain reaction，qRT-PCR）验证。肠道微生物群通过宏基因组测序进行谱系分析。
结果：二甲双胍降低了2小时负荷后血糖（2-hour post-load glucose，2hPG）并改善了肝脏脂肪变性。在肝脏中，整合的转录组-蛋白质组分析揭示线粒体底物载体（Slc25a47、Slc25a29、Slc25a30）和脂质处理因子（如Apoa5、Apoc2）的一致性下调，并与处理动物中的系统性血糖负担相关。在骨骼肌中，通路水平整合支持代谢重塑程序的激活，包括PPAR信号通路，该通路与处理动物中血糖改善的程度相关。在iWAT中，免疫和趋化因子通路的转录组富集，得到蛋白质组趋势的支持，提示在处理动物中，较高的残余炎症活性与较弱的血糖应答相关。肠道宏基因组分析未检测到微生物多样性或群落结构的显著改变。
结论：肝脏线粒体-脂质重塑和肌肉代谢激活与血糖改善一致，而脂肪组织免疫特征可能促成血糖应答的个体差异，为二甲双胍的系统性药理学提供了一个多器官、数据支持的框架。
论文解读
一、研究背景与问题
2型糖尿病（Type 2 diabetes mellitus，T2DM）影响全球5.89亿成年人（2024年），预计到2050年将增至8.53亿。约五分之四的患者同时合并肥胖，肥胖加剧胰岛素抵抗和代谢功能障碍。二甲双胍仍是肥胖T2DM的一线治疗药物，并拥有美国糖尿病协会（ADA）最强的推荐等级；然而，其获益并不均衡。约30%–40%的患者单用二甲双胍无法达到充分的血糖控制，且HbA1c降幅在个体间存在2至3倍差异。既往多数研究孤立地考察肝脏、骨骼肌和脂肪组织。肝脏通过AMPK依赖和非AMPK依赖通路抑制糖异生并重塑脂质代谢；肌肉通过摄取和氧化葡萄糖调控全身葡萄糖处置；脂肪组织作为免疫代谢枢纽，脂肪因子和浸润免疫细胞调节胰岛素敏感性。单器官研究虽已识别出AMPK等重要靶点，但无法捕捉二甲双胍如何在器官间重新分配代谢和炎症信号。因此，需要一种组织分辨的多组学方法。研究人员选用自发肥胖糖尿病蒙古沙鼠模型，该模型无需饮食干预即可发展出肥胖、胰岛素抵抗、高血糖和肝脏脂肪变性，为研究药物反应提供了遗传学上一致的平台。
二、研究设计与主要发现
研究人员对肥胖糖尿病蒙古沙鼠给予二甲双胍（100 mg/kg/天，口服灌胃；n = 4）或溶媒对照（n = 3）处理35天，对肝脏、骨骼肌和腹股沟白色脂肪组织（inguinal white adipose tissue，iWAT）进行匹配的RNA测序（RNA-seq）和数据非依赖采集质谱（DIA-MS）蛋白质组学分析，并通过宏基因组测序分析肠道微生物群。研究建立了组织分辨的转录组和蛋白质组分析框架，将匹配的转录组-蛋白质组特征锚定到个体动物水平的表现型，为解析复杂代谢疾病中血糖应答的个体差异提供了可扩展的策略。
三、关键技术方法
研究人员采用的主要技术方法包括：基于Illumina NovaSeq 6000平台的RNA测序（paired-end 150 bp），以蒙古沙鼠参考基因组MunDraft-v1.0（Ensembl Release 111）为比对参考；基于Orbitrap Astral质谱仪的数据非依赖采集质谱（DIA-MS）蛋白质组学，以UniProt蒙古沙鼠蛋白质组数据库（Taxonomy ID: 10047）为检索库，DIA-NN（v1.8）处理原始数据；定量实时聚合酶链反应（qRT-PCR）和Western blot验证关键发现；基于Illumina NovaSeq平台的粪便宏基因组测序，经Kraken2比对NCBI NR数据库进行物种分类；跨组学整合分析采用九象限图、KEGG通路富集（clusterProfiler）和模块评分等方法。样本来源于首都医科大学伦理委员会批准的肥胖糖尿病蒙古沙鼠模型（批准号AEEI-2021-052），动物为8–10月龄雄性，入组体重88–122 g。
四、研究结果
二甲双胍改善肥胖糖尿病蒙古沙鼠的高血糖并缓解胰腺-肝脏病理
二甲双胍治疗显著降低第5周2小时负荷后血糖（2hPG）（5.03 ± 0.28 vs 6.90 ± 0.36 mmol/L，Welch's t检验P = 0.0138，Cohen's d = ?3.23），空腹血糖（FPG）呈非显著性降低（4.75 ± 0.42 vs 6.03 ± 0.27 mmol/L，P = 0.0518）。体重变化在组间相当（?13.5 ± 3.5 vs ?14.3 ± 4.6 g，P = 0.89），表明血糖改善并非由体重下降所致。组织病理学分析显示，二甲双胍部分恢复了对照组中受损的胰岛结构，免疫荧光染色显示处理后胰岛内胰岛素阳性区域扩大。在肝脏中，二甲双胍减少脂滴积累并改善肝脏形态，定量分析确认肝脏脂滴面积显著减少（1.20 ± 0.26% vs 5.44 ± 0.91%，P = 0.0346，降幅77.9%）。
三个代谢组织的转录组景观及肝脏PC1与血糖负担的关联
主成分分析（PCA）显示肝脏、骨骼肌和iWAT的转录组分布（PC1/PC2解释方差：肝脏35.5%/22.5%，肌肉31.2%/21.8%，iWAT 42.1%/36.3%）。在肝脏中，PC1与第5周2hPG呈负相关（Spearman ρ = ?0.964，精确置换P = 0.0028），与FPG亦呈负相关（ρ = ?0.919，P = 0.0063），且不受技术变量或RNA完整性影响。肝脏PC1与血糖改善幅度（ΔFPG/Δ2hPG）无显著相关，表明其更接近终点疾病严重程度而非应答程度。
肝脏：RNA-蛋白质一致性与一致的通路转变
在肝脏中，RNA与蛋白质log2倍数变化的总体相关性中等（Spearman ρ = 0.161），但29个在两层均显著的基因显著超过随机预期（超几何P = 0.0027）。应用严格候选标准（方向一致，两层|log2FC| ≥ 0.585）得到13个肝脏基因，集中于两个功能模块：线粒体载体家族成员（Slc25a47、Slc25a29、Slc25a30）和关键脂质处理基因（PLIN2、APOA5、APOC2、LCAT），以及Lepr。通路富集显示PPAR信号、脂肪酸代谢和过氧化物酶体被抑制（FDR约10??–10?3），而胰岛素信号及相关通路上调（FDR约10?2）。六个聚焦代谢通路（氧化磷酸化、产热、过氧化物酶体、三羧酸循环、脂肪酸降解和PPAR信号）在两层一致下调。qRT-PCR确认了六个候选基因的下调，Western blot确认了Slc25a47蛋白的降低。与已建立的二甲双胍响应基因的基准比较确认了预期方向性：Pck1（log2FC = ?1.04）、Cpt1a（?2.22）、Acadm（?1.26）等均下调，而AMPKα1/α2转录本不变，与AMPK的翻译后激活一致。
骨骼肌：跨组学一致性与疗效关联的模块活性
骨骼肌基因水平整合显示重叠有限，RNA-蛋白质log2倍数变化相关性弱（ρ = 0.037）。在5,882个匹配基因产物中，仅Qdpr和Apob在RNA（FDR < 0.05）和蛋白质（名义P < 0.05）两层均显著，仅Apob满足严格双显著性和效应量标准。通路水平分析更具信息性：13个KEGG通路在转录组显著富集并具有一致的蛋白质组趋势，上调通路包括PPAR信号、胆汁分泌、视黄醇代谢、类固醇激素生物合成和细胞色素P450药物代谢。二甲双胍在肝脏中抑制相同通路，表明其激活肌肉分解代谢同时抑制肝脏脂质处理。在处理动物中，PPAR信号模块与血糖负担呈负相关（Spearman ρ = ?1.00，精确置换P = 0.083，n = 4），为探索性关联。qRT-PCR确认Myef2和Apob在肌肉中上调，Western blot确认AMPK调节亚基Prkag2蛋白增加。
iWAT：免疫/微环境特征与蛋白质组趋势支持及与血糖应答个体差异的关联
脂肪组织呈现不同模式：178个差异丰度蛋白但仅22个差异表达基因（FDR < 0.05），无探索性跨层候选基因，RNA-蛋白质log2倍数变化相关性弱（ρ = 0.074），应答主要为翻译后水平。转录组富集涉及免疫程序，包括趋化因子信号、Th17细胞分化和细胞因子-细胞因子受体相互作用。基于秩的蛋白质组趋势检验显示趋化因子信号具有最强的跨层支持（趋势P = 0.0467）。在处理动物中，整合的趋化因子模块评分与Δ2hPG呈正相关（Spearman ρ = 1.00，精确置换P = 0.083，n = 4），为探索性、假设生成的模式。qRT-PCR确认Ccr7、Ccl19、Cd40和Cd40lg在iWAT中上调。Western blot确认Sirt3蛋白增加，但因Sirt3为线粒体去乙酰化酶，研究人员将其解释为线粒体功能的翻译后调控而非免疫证据。
肠道微生物组分析未发现二甲双胍处理的显著结构变化
粪便宏基因组测序未发现二甲双胍对肠道微生物组的显著影响。α多样性指数（ACE、Chao1、Shannon、Simpson）组间无显著差异（所有P > 0.05），β多样性排序（PCA、NMDS）未显示处理相关聚类（PERMANOVA，P > 0.05）。LEfSe仅识别出溶媒组富集的低丰度、大多未表征谱系，无已知代谢相关分类群在组间差异。
数据支持的多器官二甲双胍应答模型
跨三个代谢组织，二甲双胍留下三个可识别的特征。在肝脏中，以线粒体载体（Slc25a家族）和脂质处理基因（Plin2、Apo系列）下调为核心的一致候选基因程序，与全局转录组PC1一致并紧密追踪系统性血糖负担的降低。在骨骼肌中，代谢模块（特别是PPAR信号）的协调激活可能与治疗疗效相关，模块活性越高，处理动物中血糖改善越大。在iWAT中，脂肪微环境表现出免疫/趋化因子特征的上调。与肌肉不同，该炎症模块活性较高与处理动物中血糖应答减弱相关。多组织模块分析进一步支持器官间协调：肝脏候选基因模块与肌肉PPAR信号模块在蛋白质层呈负相关（ρ = ?0.75，P = 0.052），形成肝脏下调/肌肉上调的代谢镜像；肌肉PPAR模块与iWAT趋化因子模块呈负相关（ρ = ?0.75，P = 0.052）。研究人员将13个探索性跨层候选基因和六个组织水平通路模块分为四类：仅被二甲双胍改变但无可检测的血糖变化相关性（LEPR、PLIN2、ECI2及肝脏蛋白层模块）；显示探索性血糖变化相关性但未被二甲双胍显著改变（肌肉RNA和蛋白层模块及iWAT RNA层模块）；既被二甲双胍改变又显示探索性血糖变化相关性（11个特征，包括SLC25A29、SLC25A30、APOA5、APOC2、APOB、LCAT及肝脏RNA层模块）；两者均无（iWAT蛋白层模块）。
五、讨论与结论
讨论部分指出，该研究通过同时分析肝脏、肌肉和脂肪组织，证明二甲双胍协调跨组织的器官特异性程序。肝脏应答以Slc25a47及相关线粒体溶质载体为核心，Slc25a47编码线粒体NAD+转运体，其抑制激活AMPKα，识别Slc25a47为与二甲双胍处理相关的候选肝脏分子。肌肉应答通过PPAR信号和胆汁酸代谢的通路水平激活实现，PPAR激活促进脂肪酸氧化和线粒体生物发生，类似运动诱导的代谢适应。脂肪组织免疫激活占主导，蛋白质组变化远超转录组变化，应答主要为翻译后水平。趋化因子模块评分与残余葡萄糖正相关，与脂肪组织作为免疫代谢枢纽的视角一致。肝脏Lepr的一致下调提供了系统性背景，瘦素抵抗是该模型的标志，Lepr下调可能反映改善的代谢稳态下的适应性有机体水平应答。肠道微生物组未检测到显著变化，为描述性阴性发现。研究局限性包括：蛋白质层使用名义P值未进行多重检验校正、单终点设计无法区分药物原发作用与继发适应、组内关联基于4只处理动物为探索性、排除雌性动物及使用沙鼠模型限制普适性。结论指出，二甲双胍通过不同的组织特异性程序发挥作用：抑制肝脏线粒体载体和脂质代谢，开启肌肉PPAR信号和分解代谢程序，触发脂肪组织免疫通路。每个特征均与血糖应答的个体差异相关；这些组织程序是否介导二甲双胍疗效仍需功能验证。组织分辨的多组学分析指向候选通路模块，在考虑作为预测性生物标志物或治疗靶点之前需在独立队列中验证。
