基底放射状胶质细胞如何影响新皮层大小
这篇论文要回答的是:人类特有的基因调控发生在皮层发育的哪个环节。作者的起点是一条已有充分证据的关联:基底放射状胶质细胞(bRG)在脑室下区分裂、不与脑室直接接触,其数量和增殖潜力在不同物种间与神经元数量和新皮层大小同步变化。bRG 的形态复杂度、谱系行为和代谢,都被提出可能参与这一关联。
此前未被刻画的是基因组层面:哪些顺式调控元件(CRE,即控制基因在何时、何处表达的 DNA 序列)在人类 bRG 中获得了活性。由于无法获取大猿胎儿组织,作者用大猩猩、黑猩猩和人类的诱导多能干细胞构建大脑类器官(CeO),并从中分离出顶端放射状胶质细胞(aRG)、bRG、基底中间祖细胞和神经元。他们在类器官发育第 8 周分析染色质可及性和转录,并以妊娠第 13–14 周的人类胎儿皮层作为参照,因为在这一阶段四类细胞均已存在。
人类 bRG 中 599 个可及性升高的元件
核心测量是一个比较结果,而不是绝对计数。只有当某个 CRE 在人类 bRG 中的可及性显著高于黑猩猩 bRG 和大猩猩 bRG 时,才被认定为人 bRG 候选元件,阈值为 log2 倍数变化至少 1.5、p ≤ 0.05。正是这种双重比较,使该集合代表人类升高而非仅仅灵长类间可变;这些区域还以人类胎儿 bRG 的峰作为参照。
由此得到 599 个人类 bRG 候选 CRE。论文展示了这些元件在三个物种重复样本中的可及性 z 分数,k-means 聚类把 599 个元件按类器官样本和细胞类型归为三种主要可及性模式。其中两个元件属于人类加速区,即在人类谱系上序列变化异常快的区域。作者利用放射状胶质细胞和基底中间祖细胞的细胞类型特异性 Hi-C 接触数据,把元件与基因关联,得到 589 个与这 599 个元件存在物理接触的独特靶基因。
生物学信号首先出现在靶基因集合上:基因本体分析显示,这些基因富集于前脑发育、细胞间接触调控和成熟神经细胞形态调控。这一结果描述的是与人类活性升高元件接触的基因,尚不涉及实测的细胞形态。
基因表达指向同一形态主题
研究的转录组部分提出的是另一个问题,分母也不同。比较黑猩猩与人类 bRG,作者找到 939 个差异表达基因,标准为 log2 倍数变化至少 1、校正后 p ≤ 0.05。这些基因在人类类器官 bRG 中的表达模式与人类胎儿 bRG 高度相关,支持在这一阶段用类器官替代胎儿情境。
在人类 bRG 中上调的基因主要富集于调控神经细胞形态的过程。“神经元或细胞突起调控”和“轴突发生”等条目只出现在人类 bRG 上调基因中,在 aRG 和基底中间祖细胞中未检测到,作者据此认为这是 bRG 特有的表型。此前被报道与基底祖细胞形态复杂度和增殖潜力相关的 PALMD,在人类 bRG 中显著高于黑猩猩 bRG,但在人类 aRG 或中间祖细胞中并未升高。
为在基因集层面汇总这些结果,作者用三个基因本体分组构建了形态调控通路活性评分:细胞形态发生、细胞突起形态发生和细胞突起组织。在黑猩猩类器官中,评分从 aRG 到 bRG 和中间祖细胞升高,再到神经元进一步升高。人类细胞也呈现同样的总体上升,但人类 bRG 的评分几乎与人类神经元相同,高于人类 aRG 和中间祖细胞,也高于黑猩猩祖细胞。该评分是表达的加权汇总,并非对细胞形状的直接测量,因此反映的是通路层面的活性,而不是形态本身。
人类 bRG 似乎复用了神经元程序
推动形态富集和通路评分的许多基因,传统上与神经元形态发生相关。作者把 465 个人类 bRG 上调基因投射到已发表的人类胎儿、猕猴胎儿和小鼠胚胎皮层数据集中,检验祖先神经元程序是否在人类 bRG 中活跃。
在人类胎儿细胞类型中,这些基因分为两类:一类在基底中间祖细胞和神经元中表达更高,另一类在放射状胶质细胞中相对更高。在妊娠第 11 至 18 周的人类胎儿发育阶段中,这些基因稳定出现在放射状胶质细胞以及兴奋性和抑制性神经元中。在猕猴胎儿皮层中,同一基因集主要在神经祖细胞和兴奋性神经元中表达,并在类 bRG 群体中富集。作者把这种双重模式解释为“征用”:人类 bRG 在形态调控元件和基因上获得活性,而这些元件和基因本就嵌入祖先神经元程序,而非发明了一套全新程序。这一解释是作者对表达重叠的解读,并非已证明的因果序列。
FAM107A 缺失使 bRG 形态简化
功能实验从相关性转向扰动。作者用 CRISPR/Cas9 在人类大脑类器官中敲除 FAM107A 和 CNGA3,并分析脑室下区中电转后 GFP 阳性、SOX2 阳性、EOMES 阴性的 bRG。每个细胞经 Z 轴层扫分割并判定形态类型,同时量化突起数量。
敲除 FAM107A 显著减少了全部突起数量,以及直接从胞体长出的初级突起数量,但不影响突起长度、分支或树状化。敲除 CNGA3 未显著改变任何所测形态参数。两个基因在 aRG 中也有表达,但敲除未影响 aRG 的典型形态,支持效应具有 bRG 特异性。因此两个基因表现并不相同:在该实验中,FAM107A 对突起数量有可测量的作用,而 CNGA3 的贡献未被这些形态读数捕捉到。
作者还报告,敲除 FAM107A 减少了脑室下区/皮层板中 HuCD 阳性、KI67 阴性的神经元;在黑猩猩类器官中过表达人类 FAM107A,产生的表型与人类敲除相似。他们把后者解释为 FAM107A 水平需要在两个方向上严格调控,这与既往关于形态调控蛋白的工作一致。
类器官证据及其发育上限
这些比较依赖类器官,而非大猿胎儿组织。作者把这一点同时视为必要性和优势:类器官含有 bRG,并在细胞形态、组织结构、整体基因表达等方面与妊娠早期神经发育平行,其表观基因组景观也大体保守。选择第 8 周类器官阶段,是因为四类细胞均已存在,且此时类器官与胎儿大脑在细胞和转录组层面高度相似。
作者指出,与皮层演化最相关的局限是类器官中外侧脑室下区的扩张和形成有限,而这一结构在人类胎儿大脑妊娠中期会显著扩张。他们提出,改进培养基(例如加入白血病抑制因子)可能提高 bRG 增殖并支持形成有结构的外侧脑室下区。他们还指出,物种间发育时序差异可能影响物种特异表型出现的时间,但同时认为神经发生中期的关键里程碑大体保守,且大猩猩、黑猩猩和人类的中期神经发生类器官发育看起来高度同步。
另一条边界在于“征用”的证据性质:与神经元表达的重叠是跨数据集的关联,论文中的因果主张来自对两个基因的敲除实验,而非直接改写候选元件本身。
调控层面的增益而非新基因
这篇论文的独特主张在于人类 bRG 差异发生的层级。最强的信号不是一组人类独有的基因,而是 599 个在人类 bRG 中相对另外两种大猿可及性升高的 CRE,以及黑猩猩与人类 bRG 之间 939 个差异表达基因,其中上调部分富集于形态相关条目。作者把这一结果与既往关于共享基因差异表达和影响放射状胶质细胞行为的人类加速区的例子并列,并称其贡献在于把这一概念扩展到数百个在大猿间 CRE 活性存在差异的基因组区域。
对评估这一结果的读者而言,重点在于两条独立线索的汇合:染色质可及性和转录都指向细胞突起与形态调控,而扰动数据为 FAM107A 一个基因赋予了在人类 bRG 中可测量的形态作用。CNGA3 仍由表达和增强子关联支持,但未被所报告的形态读数支持。最稳妥的解读是:人类特有的形态调控特征把人类 bRG 与其他大猿区分开来,并可能是这些祖细胞形态复杂化的候选机制,其中因果证据集中在 FAM107A 扰动上,元件层面的因果性尚未被检验。
正确理解物种比较
论文中的两组比较回答不同问题,不应合并。599 个元件的集合是相对黑猩猩和大猩猩共同定义的,因此捕捉的是所采样三个物种中人类谱系特有的变化。939 个基因的集合是黑猩猩与人类的对比,作者通过检查人类上调基因在大猩猩 bRG 中是否也更低来判断方向,结论是该模式反映人类上调而非黑猩猩下调。因此,一个相对黑猩猩升高的人类基因,若没有大猩猩比较支持,并不自动意味着相对大猩猩也属人类特异。
实际后果是:元件层面的主张更严格,而基因层面和通路评分的主张描述的是人类与黑猩猩的差异,大猩猩仅用作方向性核对。把本工作与其他人类获得性增强子研究相比较的读者还应注意细胞类型和阶段框架:测量来自类器官发育第 8 周分离的 bRG,参照妊娠第 13–14 周的人类胎儿皮层,而非成体组织或妊娠中期外侧脑室下区。
本文引用的文献
人类基底放射状胶质细胞的形态调控特征与新皮层演化
对大猩猩、黑猩猩和人类大脑类器官的比较染色质与转录分析,鉴定出 599 个人类活性升高的 bRG 候选 CRE,关联到 589 个富集于细胞形态过程的靶基因;在人类类器官中敲除 FAM107A 减少了 bRG 的突起数量。
