科学问题
人类大脑皮层的扩张主要依靠外放射状胶质细胞等祖细胞驱动,但这些胶质亚型各自的基因调控程序长期不清楚。此前单细胞染色质构象测序的分辨率停留在 5 kb 左右,很难看清远端增强子与启动子之间的精细连接,大量精神疾病风险变异因此找不到对应的靶基因。本研究要回答的是:能否在纯化的胶质亚型中建立全基因组高分辨率三维调控网络,并据此解释非编码变异与人类加速区如何影响皮层发育。
背景与领域脉络
人类大脑的演化伴随着皮层的显著扩张,外放射状胶质细胞是其中特异扩增的祖细胞。既往研究用单细胞转录组与染色质开放性描绘了皮层细胞多样性,但染色质构象捕获技术的分辨率限制了对 2 kb 尺度启动子—增强子环的解析。这一空白既妨碍理解胶质谱系特化,也让精神疾病相关的非编码风险变异长期无法定位靶基因。在纯化细胞亚型水平建立高分辨率三维调控网络,因此成为破解皮层发育与发育异常的关键一步。
方法与技术路线
团队用荧光激活细胞分选纯化妊娠中期人胎脑皮层中的四类胶质细胞。先按 SOX2 与 HOPX 的组合,从 6 例孕 15–18 周样本中分出心室放射状胶质细胞(vRG)与外放射状胶质细胞(oRG)。再按 PU.1 与 OLIG2 标记,从 9 例孕 22–24 周样本中分出小胶质细胞(MG)与少突胶质前体细胞(OPC)。每类细胞每种检测保留 2–4 个重复。随后平行开展四类检测:RNA 测序、ATAC 测序、全基因组重亚硫酸盐测序与 H3K4me3 介导的 PLAC 测序。功能验证则借助小鼠胚胎转基因报告、shRNA 敲减单细胞测序与机器学习模型完成。
主要结果
四类胶质细胞中各注释出六万到七万个候选调控元件
团队在 vRG、oRG、OPC 和 MG 中分别鉴定出 69,141、72,450、65,295 与 60,958 个候选顺式调控元件(cCRE),其中六成以上落在启动子之外。与已发表的大规模并行报告实验比对,四类细胞中 44.5%–51.1% 的 cCRE 显示真实增强子活性,说明这套注释在功能上站得住。
2 kb 分辨率把启动子互作图谱做厚了 4 倍
H3K4me3-PLAC 测序在四类细胞中分别检出 136,389、144,089、140,297 与 135,123 个显著互作,全部定位到 2 kb 精度。这个数量是既往 5 kb 分辨率神经发育图谱的 4 倍。约 80% 的互作发生在启动子与远端元件之间,平均跨距 188–233 kb。含 cCRE 的远端片段与启动子的接触强度,明显高于不含 cCRE 的片段,提示调控元件本身也在塑造基因组三维架构。
小鼠体内实验确认了这些元件的增强子身份
团队挑选 18 个放射状胶质细胞与中间前体细胞元件做胚胎期转基因报告实验,14 个(77.8%)在脑内显示增强子活性,11 个阳性元件中有 10 个同时是 vRG 与 oRG 的 cCRE。对照的神经元元件仅 3/11 阳性,与取材时点处于祖细胞丰度高峰(E12.5)相符。1,233 个 VISTA 神经增强子在 vRG、oRG 与 OPC 的互作 cCRE 中显著富集,并能跨越 800 kb 连接到神经系统基因 PTPRG 的启动子。
染色质构象比 DNA 甲基化更能分开 vRG 与 oRG
vRG 与 oRG 之间检出 7,941 个差异开放区域,差异甲基化区域只有 756 个,两者相差十倍。26% 的 vRG 与 63% 的 oRG 差异表达基因可由差异开放区域(启动子处或远端互作)解释,而甲基化只能解释 3.8% 与 5.6%。这一结果与此前单细胞多组学仅凭染色质构象就能区分放射状胶质细胞亚型的观察相互印证。
LHX2 维持 oRG 的自我更新并压制星形胶质化
足迹分析显示,ASCL1、TCF4、TCF12 的结合位点偏向 vRG,LHX2 的结合位点偏向 oRG。预测结合差异与转录因子表达差异显著相关(PCC=0.54,P=0.009)。LHX2 在 oRG 中有 126 个特异预测结合位点,结合染色质互作串联出 87 个高表达靶基因,覆盖 18.4%(16 个)的 oRG 差异表达基因。shRNA 敲减 LHX2 后,oRG 向星形胶质细胞方向偏移,G1 期细胞比例升高,自我更新能力下降。
风险变异与人类加速区都落在 oRG 的远端元件上
机器学习模型在四类胶质的 cCRE 中优先排选出 112 个精神分裂症风险变异,并在体内证实 rs4449074 的风险等位基因会削弱 vRG 增强子活性。oRG 特异 cCRE 显著富集人类加速区,其中部分元件的调控活性与黑猩猩直系同源序列存在人类特异的差异,并与神经发育关键基因形成立体互作。
创新点
相比只关注神经元、或分辨率停留在 5 kb 的既往工作,本研究首次在纯化胶质亚型上把三维互作图谱推到 2 kb。更关键的是它给出了一条方法学结论:高度接近的 vRG 与 oRG,用染色质构象区分远比用 DNA 甲基化灵敏。作者还把人类加速区直接挂到 oRG 的远端调控环上,为从三维基因组理解皮层演化提供了可复用的路径。
意义与影响
这是目前最完整的人类胎脑胶质细胞三维调控资源之一。它把非编码风险变异从“染色体上的位置”推进到“远端靶基因”,为精神分裂症等疾病的机制解释提供参照。同时,它从亚型特异的构象差异出发,为灵长类皮层扩张这一演化难题提供了可检验的分子线索。
局限与待验证
样本来自中期胎脑,供体数量有限且孕周不连续,无法覆盖完整的发育动态。功能验证主要依赖小鼠胚胎报告系统,而小鼠本身几乎不存在外放射状胶质细胞,人属特异的调控元件与互作难以在天然细胞环境中复现。此外增强子活性只抽验了少量元件,全基因组范围的因果靶基因仍有待独立验证。
对我们的相关性
做单细胞与空间组学的人可以直接把这套 2 kb 启动子互作数据当作评估预测算法的真值基准;做 AI 与计算生物学的人可以借鉴“多组学加三维构象先验”筛选非编码变异的范式。对我们自己的研究来说,胶质亚型特异的 cCRE 与互作网络可用于注释临床 GWAS 位点,把风险变异落到远端靶基因上。
术语速查
全文 PDF 已归档 Google Drive
04-学术文献/Paper/;原文版权归 Nature 及作者所有