IF31.7｜人类皮下脂肪组织单细胞DNA甲基化与三维基因组图谱

脂肪细胞低甲基化、髓系细胞高甲基化，且脂肪细胞短程互作富集，共同塑造腹部肥胖遗传风险。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

腹部肥胖是全球范围内日益严重的健康问题，其风险部分由组织及细胞水平的低度炎症介导。皮下脂肪组织（SAT）是缓冲肥胖的关键脂肪库，高度异质，包含脂肪细胞、脂肪干细胞和祖细胞（ASPC）、血管周细胞、内皮细胞、髓系细胞、淋巴样细胞和肥大细胞等多种类型。单核RNA测序（snRNA-seq）已揭示SAT的细胞类型特异性基因表达，但DNA甲基化和染色质构象等表观基因组层面的细胞类型差异仍不清楚。已有研究表明，基因调控区的动态甲基化影响基因表达，而染色质构象组织为区室、结构域和环，对基因调控至关重要。然而，人类SAT中细胞类型水平的表观基因组景观尚未被系统描绘，这阻碍了对腹部肥胖遗传风险的精准评估。

科学问题

卡在哪里： 此前研究主要依赖bulk甲基化或转录组数据，无法区分SAT内部不同细胞类型的表观基因组异质性。单细胞甲基化与染色质构象联合分析在脑等组织已有应用，但人类SAT这一关键代谢组织的细胞类型水平DNA甲基化、三维基因组结构及其与腹部肥胖遗传风险的关联仍属空白。

本文要回答： 本文旨在利用snm3C-seq技术，绘制人类SAT的单细胞DNA甲基化、染色质构象和基因表达图谱，鉴定细胞类型特异的差异甲基化区域和三维基因组特征，并评估这些表观基因组位点对腹部肥胖及相关炎症标志物遗传风险的贡献。

技术路线

作者主要用了：snm3C-seq、snRNA-seq、甲基化通路基因表达分析、TF motif富集、GWAS富集与分区PRS。

作者从芬兰女性SAT活检组织中分离细胞核，利用snm3C-seq同时获取单细胞DNA甲基化和染色质构象信息，并对同一样本进行snRNA-seq。通过整合两种模态，鉴定出七个主要细胞类型及一个过渡态细胞群。随后，作者系统分析了细胞类型水平的差异甲基化区域（DMR）、转录因子结合基序、染色质区室、结构域和环，并评估了甲基化通路基因（如DNMT3A和TET1）的表达模式。最后，将细胞类型特异的表观基因组区域与UK Biobank的腹部肥胖、BMI、CRP和MASLD等GWAS数据整合，进行分区多基因风险评分和GWAS变异富集分析，以揭示表观基因组对疾病遗传风险的贡献。

图 1　Fig. 1。图1展示了研究设计：从芬兰女性SAT活检中分离细胞核，分别进行snm3C-seq（同时测DNA甲基化和染色质构象）和snRNA-seq。后续分析包括细胞类型鉴定、差异甲基化区域、染色质构象动态、甲基化通路基因表达与脂肪生成基因的共表达模式，以及GWAS变异富集。

核心亮点

亮点 1｜脂肪细胞低甲基化、髓系细胞高甲基化

全基因组鉴定出705,063个差异甲基化区域（DMR），覆盖5.39%的基因组。脂肪细胞和ASPC中分别有56.3%和50.6%的DMR呈低甲基化状态，而髓系细胞中73.0%的DMR呈高甲基化状态。47.3%的DMR在脂肪细胞/ASPC与髓系细胞之间呈现相反的甲基化模式。这种细胞类型间的甲基化差异主要由DNMT3A（髓系高表达）和TET1（脂肪细胞高表达）驱动，提示二者在SAT细胞类型特异性甲基化模式中发挥关键调控作用。

图 2　Fig. 2。图2a显示基于5 kb bin mCG、100 kb bin染色质构象及联合分析的UMAP图，细胞按SAT细胞类型着色。图2b的Sankey图显示三种注释高度一致，但存在一个过渡态细胞群（63个核），其mCG注释为血管周细胞，染色质构象注释为脂肪细胞。图2c-e展示基因体mCG低甲基化与RNA表达的相关性及跨模态整合的一致性。

亮点 2｜脂肪细胞富集短程染色质互作

脂肪细胞中短程（100 kb–2 Mb）与长程（10–100 Mb）染色质互作的比例显著高于其他细胞类型（-log10P > 91），且与ASPC相比差异更大（-log10P > 197）。脂肪细胞还拥有更多、更短的结构域（平均4,120个/细胞，中位长度679 kb），以及更多染色质环（8,852个环峰，中位环长230 kb）。这些局部三维基因组特征与脂肪细胞较高的转录活性相关（UMI总数增加1.5倍），提示短程互作可能支持脂肪生成相关基因的精细调控。

图 3　Fig. 3。图3a点图显示PPAR信号通路基因（ACSL1、ADIPOQ、LPL、PCK1、PLIN1、PLIN4）在脂肪细胞和过渡态细胞中基因体mCG低甲基化且高表达，其他细胞类型则相反。图3b堆叠条形图显示各细胞类型中DMR的甲基化状态比例：脂肪细胞和ASPC以低甲基化为主，髓系细胞以高甲基化为主。

亮点 3｜染色质区室反映细胞谱系关系

基于100 kb分辨率的区室分析显示，内皮细胞和髓系细胞具有更频繁的区室切换，而脂肪细胞、ASPC和血管周细胞则呈现更长的同型区室片段。差异构象区域（44.3%的基因组）中，活性A区室富集至少1.36倍。层次聚类将内皮细胞和髓系细胞归为一支，与基于DMR的谱系树一致，支持髓系祖细胞可产生血管内皮细胞的假说。

图 4　Fig. 4。图4a显示各细胞类型接触频率随基因组距离的分布，脂肪细胞短程互作比例更高。图4b箱线图展示短程/长程互作比，脂肪细胞显著高于其他类型（除过渡态）。图4e展示12号染色体上各细胞类型的区室模式，内皮和髓系细胞区室切换更频繁。图4f显示差异构象区域在A区室中富集。

亮点 4｜脂肪细胞A区室和低甲基化区域富集腹部肥胖遗传风险

利用UK Biobank数据，所有细胞类型A区室构建的WHRadjBMI多基因风险评分（PRS）均显著（P < 0.05）且富集（Pperm < 0.05），解释了全基因组PRS至少80%的方差，而B区室无富集。脂肪细胞低甲基化区域和A区室同时富集WHRadjBMI GWAS变异，髓系A区室则富集CRP GWAS变异。这表明脂肪细胞活性染色质和低甲基化区域是腹部肥胖遗传风险的主要载体，而髓系A区室与炎症风险相关。

图 5　Fig. 5。图5a示意DNA甲基化和去甲基化通路。图5b点图显示TET1在脂肪细胞高表达，DNMT3A在髓系细胞高表达。图5c堆叠条形图显示各细胞类型低甲基化和高甲基化DMR比例，与TET1/DNMT3A表达模式一致。图5d UMAP显示细胞平均全局mCG水平，脂肪细胞和ASPC较低，髓系细胞较高。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： snm3C-seq、snRNA-seq、甲基化通路基因表达分析、TF motif富集、GWAS富集与分区PRS
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 关注脂肪组织生物学、单细胞表观基因组学或代谢性疾病遗传机制的研究者。

带走一句

单细胞多组学揭示脂肪细胞与髓系细胞的表观基因组对立模式，为解析复杂性状的细胞类型特异性遗传调控提供了新范式。