Methylene blue

点击蓝字 关注我 发现更多精彩内容 ❄ 引言 ❄ 生信朋友们上午好呀，今天解读发表于Aging Cell的TMFI多组学研究。 骨骼肌不仅是运动器官，更是维持全身代谢稳态的核心组织。然而，传统研究往往将肌内脂肪浸润视为单纯的影像学现象或衰老伴随的被动改变，缺乏系统性、因果性的多组学解析。本研究跳出“单一表型—单一疾病”的传统分析框架，以大腿肌间脂肪浸润（TMFI）为切入点，构建了一条从人群健康关联 → 生活方式中介 → 全基因组关联 → 单细胞定位 → 多组织转录调控 → 蛋白代谢标志物验证的完整证据链，系统论证了TMFI作为跨系统健康风险标志物的生物学基础。 对生信研究者而言，这套框架的核心启示在于：不满足于“发现相关”，而是通过PRS分群、SMR/TWAS共定位、蛋白代谢中介分析等多维手段，将表型与底层分子机制串联为可验证的因果链条，有效规避了单纯相关性分析难以支撑生物学结论的尴尬。该研究的分析思路可直接迁移至其他影像学表型或组织特异性性状的深度挖掘中。 还在为课题注水，选题焦虑，发文卡壳发愁吗，点击下方名片⬇️，关注走进生信，后台滴滴，带你复现高分思路，探索更多科研领域~ ✦ 研究亮点✦ 🌍从“局部影像指标”升级为“跨系统健康哨兵” 传统上，TMFI多被置于骨关节炎或老年肌少症的狭窄语境下讨论。本研究首次在55,120人的规模上系统证明：TMFI不仅与全因死亡显著相关，还同时关联12大系统、数十种疾病的发病风险，且呈明确剂量反应关系。这意味着，大腿肌肉的脂肪浸润程度具备了“系统性健康晴雨表”的潜力。 🌍PRS＋孟德尔化＋中介分析三阶递进 研究通过GWAS鉴定出79个独立TMFI风险位点，并构建PRS验证了遗传易感性与多系统疾病之间的因果关系。进一步，SMR与TWAS共定位分析在9种组织中交叉验证了TMFI相关基因的表达效应，而蛋白代谢中介分析则揭示了生活方式与遗传风险如何通过循环分子影响TMFI。这种“遗传工具→分子中介→临床终点”的递进逻辑，极大提升了结论的证据等级。 🌍从“靶点发现”到“老药新用”无缝衔接 研究者将GWAS定位基因映射至DGIdb数据库，筛选出15种可作用于≥3个TMFI相关基因的已有药物（如双嘧达莫、己酮可可碱、亚甲蓝等），并同时警示了多种化疗药物可能加重TMFI的不良效应。这一设计使基础发现直通临床转化，提供了明确的干预候选清单。 🌍一套流程可复制至任意组织影像表型 从UKB影像数据提取 → 全表型组关联 → 单细胞转录组整合 → 多组织eQTL共定位 → 蛋白代谢深度学习建模，本研究搭建了一条标准化的影像表型深度解析管线。对于生信从业者而言，只要手中拥有可量化的影像指标（如肝脏脂肪分数、冠脉钙化积分、脑区体积等），即可几乎无缝迁移这套分析框架，大幅提升课题的系统性和发表潜力。 ☁ 文献解析 ☁ ❄ 研究背景 ❄ 为什么关注大腿肌间脂肪？ 骨骼肌占人体总重约40%，其质量与功能直接决定代谢稳态。TMFI指的是脂肪在 thigh 肌纤维间及肌细胞内的异常沉积，MRI是当前最敏感的无创定量手段。既往研究已零星报道TMFI与糖尿病、心衰、骨关节炎等单病种的关联，但始终缺乏一幅全局性的“健康影响地图”。 ✦ 研究方法&数据来源 ✦ 研究方法 表型关联与中介：偏相关＋Cox回归评估TMFI与衰老/炎症指标及12类系统性疾病风险的关联；限制性立方样条确认剂量反应；中介分析量化生活方式通过TMFI影响健康结局的比例。 GWAS与多组学整合：46,347人GWAS鉴定79个独立SNP；MAGMA/FUMA/GWAS catalog多角度基因注释；构建PRS在独立队列验证遗传风险与临床结局关联，并检验PRS与生活方式的相乘/相加交互。 单细胞锚定：scPagwas整合GWAS与vastus lateralis单细胞转录组，计算trait-relevant score，精准定位关键肌细胞亚群及其KEGG通路。 多组织共定位：整合GTEx v8九种组织eQTL，通过SMR＋HEIDI检验与TWAS交叉验证，筛选跨组织一致关联TMFI的基因。 蛋白代谢全景与深度学习：Olink（2923蛋白）和NMR（251代谢物）关联分析；筛选显著分子检验其中介效应；1D-CNN构建TMFI蛋白/代谢特征评分，验证与远期健康结局的关联。 数据来源 本研究所有数据均来源于UK Biobank。TMFI影像数据来自55,120名参与者的颈-膝Dixon MRI扫描，自动分割四区域脂肪浸润率。基因型数据源于46,347名欧洲裔参与者的Affymetrix Axiom芯片，质控后保留约930万个SNP；另设362,286人的独立非GWAS欧洲裔队列用于PRS验证。单细胞转录组数据取自33,975个vastus lateralis细胞（Lovrić等，2022）。eQTL数据整合自GTEx v8的9种组织（骨骼肌、脂肪、肝、肾、脾、肺、全血等）。蛋白组学采用Olink平台检测54,306份血浆中的2923种蛋白；代谢组学采用NMR平台检测约31,991份样本中的251种代谢物。药物靶点注释依托DGIdb数据库，涵盖901种已知药物。 ✿ 研究结果 ✿ 01 TMFI全景关联分析 四个部位分别做偏相关分析，保留各自的空间异质性，不直接取均值。TMFI与全因死亡及12类系统性疾病之间呈剂量依赖的正相关。最高四分位组40岁后的预期寿命明显缩短，早发疾病组的TMFI水平也比晚发组更高。寿命差距和时序顺序这两条线索放在一起看，因果方向就更清楚了，不是单纯的相关性。 图1 大腿肌肉脂肪浸润与多系统健康结果相关 干货启发：这种“先分部位独立算、确认一致后再汇总”的策略值得借鉴——遇到多部位或多亚型数据时，先分后合既能验证一致性，也能发现异常亚型，避免一上来就平均丢掉信息。 02 生活方式中介分析 五种生活方式与TMFI的RCS拟合结果显示，睡眠呈U型关联（最佳时长7.52小时），其余四项近似线性。在此基础上逐一检验TMFI的中介比例，范围在6.7%到71.7%之间，把“生活方式→TMFI→疾病”三段独立关联串成了一条因果链，而不是仅仅停留在相关性层面。 图2 大腿肌肉脂肪浸润介导生活方式因素对健康结果的影响 干货启发：中介分析的价值在于把两个独立的相关性（A→C和B→C）变成三段式因果链（A→B→C）。如果你的表型已有关联证据，中介分析是提升故事完整性的低成本高回报操作。 03 GWAS与PRS验证 GWAS结果中，79个独立SNP映射到47个基因组区域，大部分为首次报道。三种基因注释方法分别进行后再取交集，得到45个高置信候选基因。GWAS目录中的富集分析显示这些基因主要关联体型、脂肪分布和炎症相关表型。后续在36万独立样本中构建PRS，验证了遗传风险与多系统疾病的显著正相关关系。药物靶点映射筛选出15种已有药物，可作为后续重定位研究的候选。整个流程从发现到注释再到验证和转化，逻辑上是连续的。 图3 股四头肌脂肪浸润的遗传位点和定位基因 干货启发：三种注释取交集可大幅降低假阳性，PRS独立验证是审稿人最看重的证据之一。若样本量不够分队列，可用交叉验证或公开GWAS汇总数据做外部验证，但独立队列权重最高。 04 单细胞定位 四种分型在分子层面上的特征是各不相同的。免疫原性型的一个显著特点是血清中ADAbs滴度升高；代谢型则主要表现为粪便中的短链脂肪酸水平下降，同时Escherichia/Shigella属的丰度上升；表观型的标志物集中在血清HDAC8水平和Foxp3甲基化程度；屏障缺陷型的特征是肠道通透性增加，粪便中ZO-1和Occludin减少。这些标志物都可以通过血清或粪便样本采集，临床操作上是可行的，不需要进行侵入性的肠活检。由于每种分型背后的核心调控环路不同，对应的干预方向也需要相应调整。这意味着分型本身可以作为临床个体化治疗方案设计的一个参考依据。 图4 与股四头肌脂肪浸润相关的细胞类型和富集途径 干货启发：GWAS信号是混合细胞群体的平均效应，scPagwas可帮你拆解到单细胞精度。只要你有公开单细胞数据（大多疾病都能找到），这套定位方法可直接迁移，快速锚定“主犯”细胞类型。 05 跨组织转录共定位 在9种组织中同时做了SMR和TWAS，最后只有4个基因在两个方法、所有组织里方向一致。跨方法、跨组织的这种严格筛选能大幅降低假阳性风险，意味着保留下来的4个基因在遗传调控层面的一致性相对较高，可以优先考虑做后续的生物学验证。 图5  多个组织的基因表达影响大腿肌肉脂肪浸润 干货启发：不只要看单个组织的eQTL结果，更要在多组织中交叉验证——如果一个基因在骨骼肌、脂肪、肝脏中方向一致，其效应可能更稳健，后续验证成功率更高。 06 蛋白代谢全景与深度学习建模 不同分型对应的干预方向，需要根据核心调控通路的差异来匹配。代谢加重表观异常的类型，可以考虑丁酸盐配合HDAC抑制剂。免疫原性和屏障缺陷并存的情况，用粪菌移植加维得利珠单抗可能更合适。单纯免疫原性为主的，倾向于用低剂量免疫调节剂联合生物制剂，同时配合血药浓度监测。屏障功能缺陷为主且没有明显免疫原性特征的类型，则适合选择作用于肠道局部的生物制剂，辅以膳食纤维和益生菌来改善屏障稳定性。这些对应关系是建立在前面几层分子分型基础上的，可以作为后续临床适配和干预方案设计的参考框架。 图6 与股部肌肉脂肪浸润相关的循环蛋白和代谢物 干货启发：如果你的课题涉及循环组学数据，建议加上“特征评分建模+独立验证”这一步。它不仅验证了标志物的临床价值，还能把“需影像设备才能测”的表型变成“抽管血就能估”，转化价值翻倍。 ✦ 结论 ✦ 该研究首次系统证明大腿肌间脂肪浸润是全因死亡及多系统疾病的独立风险因子，并揭示生活方式可通过TMFI间接致病；通过GWAS-单细胞-SMR/TWAS-蛋白代谢五位一体分析，绘制了TMFI从遗传变异到循环分子的完整机制图谱，为精准干预提供了可操作靶点与候选药物。 可扩展方向： 纵向轨迹分析：本研究仅用单时间点TMFI，若引入重复影像数据，可构建TMFI变化速率与健康结局的动态关联模型 跨表型孟德尔随机化：以TMFI相关SNP为工具，系统评估TMFI与各系统疾病间的双向因果流 多模态深度学习：将TMFI影像特征与基因组、蛋白组整合，构建多模态早筛模型 其他组织脂肪浸润：肝脏、胰腺、心肌、骨髓等部位的异位脂肪沉积，均可套用本研究全管线框架 性别/年龄分层机制：本研究发现TMFI效应在<60岁及男性中更强，针对性开展分层多组学分析具有重要临床意义 声明：本人所有内容均为个人生信学习思路与经验分享，仅作学术交流和个人复盘使用，不属于专业权威科普。若有专业错误，欢迎评论指出! 如果也想用GWAS+多组学+单细胞思路进行生信分析，复现高分思路，走进生信助力每一个科研梦想! END ❄ ✦ ❄ ✦ ❄ 往/期/回/顾 IF 8.2！中科院一区TOP抢发｜单细胞两图一表，完整呈现CLL 13年克隆演化与趋同耐药 IF25.5！多祖源GWAS+PRS-CS+外部验证：一个PRS串起AD风险、病理与认知