IF50.5｜人类下丘脑时空细胞图谱HYPOMAP

433,369个单核转录组+空间转录组，首次系统解析人类下丘脑细胞构成与肥胖遗传关联。

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

研究背景

下丘脑是协调体温、睡眠、能量稳态、生殖和昼夜节律等基本生理功能的核心脑区。过去对下丘脑细胞组成和神经环路的理解主要来自啮齿类研究。人类遗传学研究已揭示瘦素-黑皮质素通路在食欲调控中的关键作用，POMC和AgRP神经元是该通路的中心节点。目前获批的肥胖和糖尿病药物如司美格鲁肽、替西帕肽和司美诺肽均靶向下丘脑，但支持其分子机制的证据多来自小鼠。单细胞测序技术的成熟为在人类中直接绘制下丘脑细胞图谱提供了可能。

科学问题

卡在哪里： 此前缺乏高分辨率的人类下丘脑单细胞转录组图谱，且未与空间信息整合。人类与小鼠下丘脑神经元类型在转录组层面的保守性尚不清楚，尤其是POMC神经元亚型和关键GPCR的表达差异。此外，肥胖GWAS信号富集于哪些下丘脑细胞类型、由哪些效应基因驱动，仍缺乏系统分析。

本文要回答： 构建人类下丘脑的单细胞时空转录组图谱，解析细胞类型多样性、空间分布、跨物种保守性，并鉴定与BMI遗传关联相关的细胞类型和效应基因。

技术路线

作者主要用了：snRNA-seq、空间转录组、cell2location、CELLECT、MAGMA、罕见变异负荷分析。

作者对8例正常BMI供体的下丘脑进行snRNA-seq，并整合公共全脑数据，共获得433,369个细胞核，构建HYPOMAP。通过多级层次聚类定义452个细胞类型（C4水平）。利用Visium空间转录组对7例供体的9个切片进行测序，并用cell2location将单细胞簇映射到空间位置，实现区域注释。进一步与小鼠下丘脑图谱进行跨物种比较，重点关注POMC/AgRP神经元和GPCR表达。最后，将HYPOMAP与BMI GWAS数据整合，用CELLECT和MAGMA鉴定富集细胞类型，并在UK Biobank外显子组数据中验证效应基因的罕见变异与BMI的关联。

图 1　Fig. 1。图1展示HYPOMAP构建流程和整体细胞图谱。a为数据集和分析结果示意图；b为整合后UMAP，按主要细胞类群着色，显示神经元、少突胶质细胞、星形-室管膜细胞等；c聚焦C2-36 Mid-3 GLU-3簇，该簇包含PVN的AVP神经元和乳头体亚簇，插图为全局UMAP中该簇位置；d为该簇的树状图和亚簇结构，热图显示平均表达。

核心亮点

亮点 1｜HYPOMAP覆盖43万细胞

整合8例正常BMI供体的snRNA-seq数据（311,964个核）与公共全脑数据（121,405个核），共433,369个细胞核，平均每个核检测到4,541个UMI和2,040个基因。UMAP显示主要细胞类群：神经元（166,475）、少突胶质细胞（175,109）、星形-室管膜细胞（63,111）和其他非神经元细胞（28,674）。多级层次聚类产生452个C4簇，涵盖下丘脑主要核团如ARC、VMH、PVN、SCN等。

图 2　Fig. 2。图2为人类下丘脑细胞类型的层次共识树（圆形树状图）。最高层C0分为神经元、少突胶质细胞、星形-室管膜细胞和其他非神经元细胞。C1-C3进一步定义细胞类型结构，神经元在C4有452个簇。边缘标注簇名，热图环显示两个研究来源的细胞相对贡献。该图展示细胞类型分类体系的层级关系。

亮点 2｜空间映射揭示区域异质性

利用Visium对9个下丘脑切片进行空间转录组分析，通过cell2location将单细胞簇映射到空间位置。VMH的5个C3簇在空间上呈现不同分布。对空间spot进行Leiden聚类得到23个区域簇，可区分ARC、VMH、SCN、PVN、SON、LTN、ME和TMN等核团。非神经元细胞也显示空间异质性：星形胶质细胞分为3个空间分离的亚群，室管膜细胞和伸长细胞在第三脑室壁和正中隆起有不同标记基因表达。

图 3　Fig. 3。图3展示空间转录组结果。a为9个下丘脑切片的参考图谱，从前到后排列；b为SLC17A6、SLC32A1、TBX3、FEZF1和SIM1的空间表达，分别标记VMH、ARC和PVN；c-e为C2-35 Mid-2 GLU-2的5个C3分支，c为UMAP高亮，d为分支着色，e为空间映射，显示VMH内不同亚簇的空间分布差异。

亮点 3｜POMC神经元跨物种差异

人类POMC神经元分为3个C4簇：C4-373 POMC/PRDM12（高表达LEPR，68%共表达）、C4-374 POMC/CALCR、C4-375。空间上C4-373靠近正中隆起，C4-374靠近第三脑室。与小鼠比较，C4-373对应小鼠Anxa2.Pomc亚簇，C4-374对应Ttr.Pomc，但C4-375与Glipr1.Pomc相关性低。关键差异：GLP1R在人类POMC/PRDM12簇中表达，而小鼠中Glp1r和Lepr分别在不同POMC亚群表达。人类POMC/CALCR簇表达CALCR，对应小鼠Glp1r+ POMC神经元。

图 4　Fig. 4。图4比较人类与小鼠神经元保守性，并聚焦POMC/AgRP神经元。a为全局保守性树状图，内热图显示跨物种关系类型（1:1、1:N、M:1），外热图显示匹配小鼠簇的相关性；b为UMAP高亮POMC+和AGRP+簇；c为POMC簇空间映射；d为POMC簇与小鼠对应簇的相关性；e-g为POMC簇中物种富集GPCR的表达；h为AgRP神经元中GPCR表达。

亮点 4｜BMI GWAS富集291个神经元簇

将HYPOMAP与BMI GWAS（N=806,834）整合，用CELLECT和MAGMA分析，发现452个C4簇中291个显著富集BMI GWAS信号（Bonferroni校正P<0.05/452）。富集簇主要位于中下丘脑，尤其Mid-1（96%）、Mid-2（64%）和Mid-3（86%）。最显著簇C4-77（P=1.43×10^-26）表达PTGFR和ONECUT1/2/3，映射到DMH。其他富集簇包括HMX3+、NPSR1+、GAL+和SST+神经元。

图 5　Fig. 5。图5展示BMI GWAS富集和效应基因验证。a为452个C4簇的BMI GWAS富集显著性排序，虚线为Bonferroni阈值，291个簇显著富集；b-g为UK Biobank中6个效应基因（BSN、CALCR、CORO1A、MC4R、PCSK1、POMC）的罕见变异与BMI关联的森林图，显示效应量和置信区间。

亮点 5｜效应基因与罕见变异关联

通过CELLECT-GENES鉴定426个效应基因（95%分位数细胞类型特异性且位于MAGMA基因关联前1000）。在UK Biobank外显子组数据（n=419,692）中进行罕见变异负荷分析（MAF<0.1%，PTV或CADD≥25错义变异），发现6个基因与BMI显著相关（P<0.05/426）：MC4R、PCSK1、POMC、CALCR、BSN和CORO1A。其中BSN近期已被报道与肥胖相关，CORO1A是全新发现（β=0.98±0.215，P=5.6×10^-6）。

生信可带走

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

• 方法栈： snRNA-seq、空间转录组、cell2location、CELLECT、MAGMA、罕见变异负荷分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事下丘脑、代谢调控、单细胞空间组学或肥胖遗传学研究的科研人员。

带走一句

单细胞图谱与空间转录组整合，再叠加GWAS和罕见变异分析，能系统鉴定疾病相关细胞类型和候选基因，为靶点发现提供可复用的框架。