长期自主运动部分逆转肥胖相关棕色脂肪组织（BAT）转录组重塑：免疫炎症及收缩/Ca2+相关转录特征

《Metabolism Open》：Long-term voluntary exercise partially reverses obesity-associated transcriptomic remodeling in brown adipose tissue

背景：尽管运动可能抵消肥胖诱导的棕色脂肪组织（brown adipose tissue，BAT）功能和分子重塑引起的一些变化，但长期自主运动（voluntary exercise，VE）在多大程度上逆转肥胖相关BAT转录变化仍不清楚。

方法：将雄性C57BL/6J小鼠分为对照组、高脂饮食（high-fat diet，HFD）组或HFD加自主运动（HFD + VE）组，为期四个月（每组n = 4）。使用Clariom S Mouse Array分析肩胛间BAT转录组谱。采用eBayes ANOVA探索差异表达，并选择在较高表达组中绝对倍数变化>2、名义P值<0.05且平均log₂信号强度≥8的转录簇用于探索性分析。另外检查错误发现率（false discovery rate，FDR）校正P值以考虑多重检验。使用DAVID进行基因本体论（Gene Ontology，GO）富集分析。

结果：使用预先设定的探索性标准（在较高表达组中绝对倍数变化>2、名义P<0.05且平均log₂信号强度≥8），对照组与HFD组之间有579个转录簇存在差异。在HFD与HFD + VE比较中，89个转录簇符合这些探索性标准；其中74个在VE后表达降低，15个在VE后表达升高。在全基因组多重检验校正后，89个转录簇中没有1个在FDR<0.05时保持显著。然而，GO富集分析仍识别出协调的免疫/炎症和肌肉收缩相关生物学过程，其中若干过程在GO条目检验的FDR校正后仍显著。

结论：长期自主运动与转录组变化相关，这些变化与HFD相关改变的一个子集相反，尤其是免疫和炎症特征。还识别出一个受限的收缩/Ca²⁺相关转录特征，尽管仅凭转录组数据不能确定Ca²⁺信号改变。这些探索性发现强调需要进一步研究运动相关BAT重塑背后的细胞和信号机制。

论文解读：长期自主运动部分逆转肥胖相关棕色脂肪组织转录组重塑

研究背景与问题
肥胖导致脂肪组织结构和功能重塑，是2型糖尿病和心血管疾病等代谢紊乱的重要危险因素。运动可改善肥胖相关代谢功能障碍，既往研究多关注骨骼肌和白色脂肪组织（white adipose tissue，WAT），但越来越多证据提示棕色脂肪组织（brown adipose tissue，BAT）也参与运动诱导的代谢适应。BAT线粒体含量高、产热能力强，受交感刺激尤其冷暴露激活，成人代谢活跃BAT活性与肥胖程度负相关。肥胖会改变BAT形态和代谢，包括脂质蓄积、线粒体功能与氧化代谢受损及炎症重塑，损害其代谢特性。运动可改变BAT代谢、血管化和炎症状态，但长期自主运动（voluntary exercise，VE）能否将肥胖相关转录变化逆转向非肥胖状态尚不清楚。因此，研究人员旨在识别高脂饮食（high-fat diet，HFD）改变、随后被运动反向改变的基因和生物学过程，假设长期VE可对抗BAT中部分HFD相关转录改变，尤其是代谢和炎症过程。该研究作者为Takuya Sakurai、Eriko Nozaki和Tetsuya Iwawa。

研究人员开展的研究、结论及意义
研究人员将6周龄雄性C57BL/6J小鼠随机分为对照组、HFD组和HFD + VE组，干预四个月，每组n = 4，分析肩胛间BAT转录组。结果显示，慢性HFD诱导BAT广泛转录重塑，涉及免疫/炎症、细胞黏附、脂质代谢、棕色脂肪细胞分化和产热功能；长期VE与部分HFD相关改变方向相反的转录变化相关，尤其是免疫和炎症特征，另识别出一个受限的收缩/Ca²⁺相关转录特征。尽管无个体转录簇在全基因组错误发现率（false discovery rate，FDR）校正后显著，但探索性基因本体论（Gene Ontology，GO）分析发现免疫/炎症和肌肉收缩相关过程在GO条目FDR校正后仍富集。该研究提示运动相关BAT重塑可能涉及免疫微环境改善及潜在Ca²⁺相关过程，但需进一步细胞和信号机制研究。

主要关键技术方法
样本来源为Sankyo Labo Service Corporation, Tokyo, Japan的6周龄雄性C57BL/6J小鼠，随机分为C、HFD和HFD + VE组，干预四个月，每组n = 4。采集肩胛间BAT，提取总RNA，要求RNA完整性值等效（RINe）≥7.0。使用Clariom S Mouse Array进行全基因组表达谱分析，在Transcriptome Analysis Console（TAC）中采用信号空间转换-稳健多阵列平均法（SST-RMA）归一化，以eBayes ANOVA检验差异表达。预设探索标准为绝对倍数变化>2、名义P<0.05、较高表达组平均log₂信号强度≥8，并检查FDR校正P值。使用DAVID进行GO生物学过程富集分析；比较三组识别HFD反应性且VE后反向变化的转录簇；检查肌源性和骨骼肌相关基因表达。

研究结果
3.1 高脂饮食诱导BAT转录组变化
使用Clariom S Mouse Array，按预设标准，C与HFD组间579个转录簇不同，其中348个上调、231个下调。最强烈改变包括Mmp12、Gpnmb、Lep、Lgals3、Atp6v0d2、Cyp2e1、Tnnt3、Cyp2b10、Mup12和Mfsd2a等。GO富集分别针对上调和下调转录簇。整体上，HFD相关重塑涉及免疫/炎症以及细胞黏附、脂质代谢、棕色脂肪细胞分化和产热功能相关过程。

3.2 自主运动重塑HFD喂养小鼠BAT中的生物学过程
HFD + VE组平均日跑距离为4067 ± 711 m/day（mean ± SD，n = 4）。HFD与HFD + VE比较中，89个转录簇符合标准，74个在VE后降低、15个升高。74个降低转录簇主要富集免疫和炎症过程；GO条目FDR校正后显著者包括炎症反应、受体介导内吞、先天免疫反应、通过主要组织相容性复合体II类（MHC class II）对外源性肽抗原的加工和呈递及免疫系统过程。15个升高转录簇主要富集肌肉收缩和肌节过程；FDR校正后显著者包括肌节组织、肌肉收缩、骨骼肌收缩、心肌收缩、横纹肌收缩调节和骨骼肌细丝组装。个体转录簇层面无全基因组FDR<0.05显著者，故为探索性候选，但GO分析支持协调富集。

进一步识别C vs HFD和HFD vs HFD + VE中均符合标准且方向相反的HFD反应性转录簇，共66个：60个HFD升高/VE降低，6个HFD降低/VE升高。6个HFD降低/VE升高转录簇富集肌肉收缩相关过程，骨骼肌收缩FDR = 3.71 × 10^?4，横纹肌收缩调节FDR = 0.017，肌肉收缩调节FDR = 0.0339；肌节组织名义P = 0.00884、FDR = 0.0663。这6个转录簇为Acta1、Tnnt3、Tnni2、Casq1、Myh8和Myot，代表收缩机器和Ca²⁺处理相关受限基因集，而非广泛肌源性分化。60个HFD升高/VE降低转录簇强富集免疫/炎症过程；FDR校正后显著者包括炎症反应（FDR = 2.21 × 10^?5）、先天免疫反应（FDR = 1.21 × 10^?4）、受体介导内吞（FDR = 2.38 × 10^?4）、通过MHC class II对外源性肽抗原的加工和呈递（FDR = 0.00381）和免疫系统过程（FDR = 0.0102），另包括通过MHC class I的抗原加工和呈递、中性粒细胞趋化、血管生成正调控、中性粒细胞稳态和淀粉样β清除。这些结果为探索性通路层面证据，而非转录本水平恢复证明。

讨论部分总结
讨论认为，慢性HFD诱导BAT广泛转录重塑，长期VE与部分HFD相关改变相反的转录变化相关。无个体转录簇达全基因组FDR显著，但若干GO过程在GO条目FDR校正后显著，应解释为功能相关转录本协调变化的探索性证据。经典棕色脂肪细胞与骨骼肌细胞均源自表达Myf5的祖细胞，肩胛间BAT邻近骨骼肌；但经典肌源性调节因子（myogenic differentiation 1、myogenin、Myf5、myogenic factor 6、myocyte enhancer factor 2C）表达很低或三组无显著差异，代表性骨骼肌相关基因（creatine kinase (muscle)、myoglobin、titin、nebulin、desmin）也无显著改变，故不支持大量骨骼肌污染。6个反向变化基因与收缩机器和Ca²⁺处理相关；原文提到Fu等报道有氧运动增强HFD诱导肥胖小鼠BAT中血管内皮生长因子（VEGF）信号，并提出VEGF受体激活、磷脂酶Cγ（PLCγ）、IP3/Ca²⁺信号、钙调磷酸酶/活化T细胞核因子（NFAT）和环氧合酶2（COX2）级联。本研究基因不同，但Ca²⁺处理和Ca²⁺依赖收缩机器的反向变化提供独立转录观察，然而这些基因变化本身不证明Ca²⁺信号改变，全组织转录组也不能确定细胞来源。免疫方面，肥胖伴慢性低度脂肪组织炎症和巨噬细胞募集激活，BAT也易感；VE降低表达转录簇与炎症、先天免疫、抗原呈递和白细胞过程相关，与既往HFD增加BAT免疫/炎症基因、有氧运动抗炎的报道一致。但全组织分析不能区分棕色脂肪细胞、免疫细胞和基质细胞的贡献。HFD还抑制细胞黏附和BAT代谢功能，下调GO包括同型细胞-细胞黏附、细胞黏附、脂肪酸代谢过程、棕色脂肪细胞分化、产热和脂肪酸β-氧化；细胞黏附相关包括protocadherin gamma簇成员转录簇，其在BAT中作用尚不清楚。限制包括缺少对照饮食加VE组、RNA耗尽导致无独立qPCR验证、个体转录本发现为探索性；运动范式不同：Fu等为HFD诱导肥胖后4周强迫跑台，本研究为HFD持续期间长期自主跑轮。原文结论为：慢性HFD与BAT广泛转录重塑相关，长期VE与部分HFD相关改变相反的转录变化相关，尤其是免疫相关过程；探索性GO分析识别出免疫/炎症和肌肉收缩相关过程，受限六基因收缩/Ca²⁺相关特征可能是有趣组成部分，但转录组数据本身不能证明Ca²⁺信号改变；需整合细胞类型特异分析和功能评估以明确机制。

研究结论（译）
长期自主运动与转录组变化相关，这些变化与HFD相关改变的一个子集相反，特别是免疫和炎症特征。还识别出一个受限的收缩/Ca²⁺相关转录特征，尽管仅凭转录组数据不能确定Ca²⁺信号改变。这些探索性发现强调需要进一步研究运动相关BAT重塑背后的细胞和信号机制。