人皮下脂肪组织转录组连接线粒体转录与肥胖代谢异质性：MHO与MUO在Reactome通路及模块水平较单基因水平更可重复

《Phycology》：Human Subcutaneous Adipose Transcriptomes Link Mitochondrial Transcription to Metabolic Heterogeneity in Obesity

　　

背景：代谢健康型肥胖（metabolically healthy obesity, MHO）与代谢不健康型肥胖（metabolically unhealthy obesity, MUO）可能占据脂肪组织功能障碍连续谱中的不同位置，但稳定的分子区分仍不确定。方法：研究人员整合了三个公共人类皮下脂肪组织（subcutaneous adipose tissue, SAT）转录组数据集，同时保持参与者和平台独立性。GSE244118用于发现，GSE152991被限制于14名与发现队列不重叠的参与者，GSE55200提供跨平台支持证据。研究人员进行了差异表达分析、相关性校正的Reactome检验、方向感知的跨数据集整合、发现集定义的模块评分和脂肪组织去卷积。结果：在发现数据集中，50个基因满足代谢不健康型肥胖与代谢健康型肥胖比较的常规阈值，但没有任何基因通过预设的fold-change-aware检验或严格的跨数据集重复。在1570个共享Reactome通路中，409个在signed-Z整合中达到错误发现率（false-discovery rate, FDR）< 0.05，387个在所有数据集中方向一致。代谢不健康型肥胖持续显示较低的线粒体翻译、呼吸链和能量代谢程序，并伴随细胞外基质（extracellular matrix, ECM）重塑的继发性增加。在仅限非发现数据集且采用rank和z-mean评分时，发现集定义的线粒体模块在代谢不健康型肥胖中较低。仅14人保留集未建立显著的模块差异，细胞外基质模块综合不显著。有序横断面模式和组成敏感性描述了组织水平关联。结论：因此，代谢性肥胖异质性在通路和模块水平比单基因水平更可重复。未解决的临床混杂和bulk组织组成限制了功能和因果解释。

研究背景：肥胖具有生物学异质性，MHO与MUO是定义依赖且非固定诊断的代谢状态。MHO并非一致良性，仍可能与长期心代谢风险相关，并可能向代谢功能障碍转变。皮下脂肪组织（subcutaneous adipose tissue, SAT）是异质性的核心，其安全储存脂质的能力、组织重塑、脂肪细胞更新、内分泌、血管、免疫和结构功能均参与代谢健康。当扩张适应不良时，缺氧、炎症、纤维化、血管生成受损和线粒体功能障碍可共同出现，促成系统性胰岛素抵抗。人类研究支持线粒体与基质重塑成分，但bulk SAT转录组平均了多种细胞信号，MHO/MUO队列小、表型定义和平台不同，跨研究重复性受限。方法学建议分别预处理、评估研究独立性并整合可比统计量，而非直接合并表达矩阵。研究人员因此假设代谢性肥胖异质性作为协调的线粒体和ECM程序，比稳定单基因效应更可重复。

研究人员整合三个公共人类SAT转录组数据集，保持参与者和平台独立，并按单基因、Reactome通路、样本模块和估计细胞组成层级评估可重复性；主要比较为MUO?MHO。结论是：代谢性肥胖异质性在通路和模块水平比单基因水平更可重复，MUO显示较低的线粒体转录、呼吸链和能量代谢程序，并伴随ECM重塑增加；其意义在于提供组织水平的分层线索，但小样本、未解决临床混杂和bulk组成限制因果与临床解释。

关键技术方法：研究人员对三个公共人类SAT转录组数据集进行二次分析并保持参与者和平台独立。GSE244118（15 MHL、19 MHO、19 MUO）用于发现和模块定义；GSE152991排除30名与GSE244118重叠者及1名不适用MHL后，以6 MHO和8 MUO作为非重叠RNA测序保留集；GSE55200（7 MHL、8 MHO、8 MUO）提供跨平台支持。各矩阵分别预处理，以edgeR/limma进行差异表达和TREAT检验，以CAMERA进行相关性校正Reactome通路检验，以方向感知signed-Z整合跨数据集证据，以发现集定义模块并用rank和z-mean评分，用dtangle和CLR进行脂肪组织去卷积及组成敏感性分析。

研究结果：
队列和分析角色：分析包含90个分别分析的SAT表达样本。GSE244118用于发现和模块定义；GSE152991最终6 MHO和8 MUO构成参与者非重叠RNA测序保留集；GSE55200含7 MHL、8 MHO、8 MUO，作为跨平台支持证据。每个表达矩阵单独分析，未跨平台合并，跨数据集证据仅在基因方向、通路统计或标准化模块效应水平合并。

发现阶段基因水平差异：GSE244118中，MHO?MHL得1969个差异表达基因（differentially expressed genes, DEGs），MUO?MHL得2958个，MUO?MHO仅50个常规阈值DEGs。glmTreat/TREAT检验确认386个MHO?MHL和1207个MUO?MHL基因，但MUO?MHO无基因。表达模式以肥胖共享变化为主，MHO保留/MUO特异7个，MHO保护候选0个，进展性失调37个。结论：发现阶段大效应MHO?MUO差异稀少，未通过预设TREAT。

候选基因重复有限：GSE152991保留集中无基因满足全基因组DEG定义；34个候选可明确映射，14个方向一致，效应量一致性弱（Spearman ρ=0.049，p=0.783），2个有名义方向支持，1个通过候选集FDR，0个严格重复。GSE55200中41个候选可映射，33个方向一致（精确二项p=5.61×10^?5），但幅度不相关（Spearman ρ=0.040，p=0.806），3个名义支持，无候选集或严格全基因组重复。结论：单基因水平重复有限且不一致。

线粒体和ECM通路的跨数据集趋同：CAMERA竞争性基因集检验显示，1570个共享Reactome通路中409个meta-analysis FDR<0.05，387个方向一致。MUO下调线粒体翻译、线粒体核糖体相关质量控制、呼吸电子传递、有氧呼吸、嵴形成、TCA循环和支链氨基酸分解，上调纤连蛋白基质形成和ECM蛋白聚糖通路。14条代表性通路均FDR<0.05且方向一致，留一数据集分析保留预期方向。结论：可重复区分集中在协调的线粒体和基质重塑程序，而非稳定单基因。

模块证据和有序横断面模式：仅用GSE244118定义线粒体生物能量学模块（R-HSA-5368287、R-HSA-1428517）和ECM重塑模块（R-HSA-1474244），去除共享基因后分别含351和281个基因。非发现数据集综合中，线粒体模块在MUO较低：rank评分Hedges’ g=?1.004（95% CI ?1.735至?0.274；FDR=0.0140），z-mean评分g=?0.955（95% CI ?1.679至?0.230；FDR=0.0197），I²=0%；14人GSE152991单独未达FDR<0.05。ECM模块方向为正但未达预设阈值。GSE244118和GSE55200均显示MHL?MHO?MUO有序横断面模式，MHO处于中间。结论：支持组织水平有序状态，但不证明个体纵向进展。

细胞组成敏感性：dtangle估计显示MUO的估计脂肪细胞丰度低于MHO（随机效应Hedges’ g=?0.601，95% CI ?1.073至?0.129，FDR=0.0378，I²=0%）；巨噬细胞和血管细胞估计异质且未稳定。调整CLR转换组成主成分后，GSE244118模块对比减弱，非发现数据集方向未被一致消除。结论：仅较低估计脂肪细胞丰度有稳定跨数据集支持；组成调整是敏感性分析而非中介分析，bulk不能区分细胞丰度与细胞内表达变化。

讨论部分总结：线粒体转录程序跨数据集趋同，但不等同于直接线粒体功能测量或生物标志物验证；其广度涉及翻译、呼吸链、能量代谢等多程序。ECM重塑为次要、情境依赖特征，广泛ECM模块未达预设阈值。组成不是必然混杂，脂肪细胞减少、脂肪生成受损或免疫募集可能属于生物过程。有序横断面模式提示MHO可能为中间状态，但干预研究显示脂肪转录变化兼具可逆与持久成分，如DiOGenes中SCD、FADS1、FADS2变化及LEP持续下调，减重手术后IGF1、DUSP1、GPX3、GLUL降低可保留两年。研究优势包括重叠审计、分别预处理、相关性校正通路检验、发现集定义模块、非发现数据集综合、双评分、去卷积和组成调整；局限包括保留集仅14人、代谢健康定义不同、横断面bulk、临床链接不可用、年龄性别BMI等残余混杂、仅两个非发现集合及14通路集由现有数据选择。未来需完全独立临床定义队列、连续代谢表型、单核/空间/蛋白质组/功能线粒体测量。

结论翻译：对人类SAT转录组的分别分析识别出与代谢性肥胖异质性相关的线粒体转录程序。通路水平趋同和仅局限于非发现数据集的线粒体模块综合，提供了比单基因重复更强的支持；而广泛的ECM模块综合未达到预设显著性阈值。有序横断面表达模式和组成敏感性描述的是组织水平关联，但未建立纵向进展、细胞内在线粒体功能障碍或所识别转录模式的时间持续性。小样本和未解决的临床混杂排除了因果、诊断或治疗选择主张。