《iScience》：Genomic structural variations distinguishing archaic and modern humans are associated with neurodevelopment
编辑推荐：
理解已灭绝的古人类与现代人之间神经发育的差异是进化生物学的核心目标之一。结构变异（structural variants, SVs）是表型变化的关键基因组驱动因素。本研究通过筛选来自多种非人灵长类（non-human primate, NHP）谱系的固定SV事
理解已灭绝的古人类与现代人之间神经发育的差异是进化生物学的核心目标之一。结构变异（structural variants, SVs）是表型变化的关键基因组驱动因素。本研究通过筛选来自多种非人灵长类（non-human primate, NHP）谱系的固定SV事件，鉴定了一组存在于古人类基因组中的SV，称为古人类结构变异（archaic hominin structural variants, AHSVs）。AHSVs显著富集于调控大脑发育和功能的基因中，其中许多基因与神经精神疾病相关。研究发现若干AHSVs与大脑顺式调控元件（cis-regulatory elements, CREs）重叠。此外，在三个神经相关基因（PDE10A、RAPGEF4和DCTN4）中，现代人特异性染色质环锚与古人类共有的缺失（deletions, DELs）重合，提示AHSVs可能影响了三维染色质架构。总体而言，本研究拓展了我们对现代人脑适应性进化的理解。
论文解读：古人类结构变异与神经发育的关联研究
一、研究背景与科学问题
理解使我们成为独特现代人类的遗传机制是进化生物学中最关键的问题之一，而探索现代人脑功能复杂性演化的分子机制是这一探究中不可或缺的部分。尼安德特人和丹尼索瓦人作为与现代人亲缘关系最近的已灭绝古人类，其基因组为研究现代人演化过程中的基因组及表型变化提供了独特窗口。约60–70万年前，尼安德特-丹尼索瓦祖先从非洲扩张至欧亚大陆，随后分化为东部和西部亚群——丹尼索瓦人和尼安德特人。约4–5万年前，现代人走出非洲并基本取代了二者。有研究认为，这种取代在一定程度上与古人类和现代人之间脑功能复杂性的差异有关。
以往研究多聚焦于单核苷酸变异（single-nucleotide variant, SNV）水平，例如现代人特有的单氨基酸替换影响KIF18A、KNL1和SPAG5等与脑发育中神经祖细胞有丝分裂活性相关的蛋白。然而，结构变异（structural variants, SVs）作为长度超过50碱基对（bp）的大规模基因组变异，相比SNV往往对表型有更显著的影响，可影响蛋白质功能、顺式调控元件（cis-regulatory elements, CREs）或改变三维基因组结构。由于古DNA样本稀缺且基因组DNA序列碎片化等技术限制，获取古人类完整SV图谱十分困难，针对古人类与现代人差异SV的研究仍相对匮乏，这限制了对SVs如何促成现代人与古人类表型差异（尤其是神经发育过程和脑功能复杂性）的理解。
二、研究内容与核心发现
研究人员利用Mao等人先前发表的、基于长读长测序数据鉴定的非人灵长类（non-human primate, NHP）谱系固定SV数据集（共1,338,997个事件），通过基因分型方法对四个高质量古人类基因组——阿尔泰山丹尼索瓦洞穴的30×覆盖度丹尼索瓦人基因组、52×覆盖度Altai尼安德特人基因组、克罗地亚Vindija洞穴的30×覆盖度尼安德特人基因组以及阿尔泰山Chagyrskaya洞穴的27×覆盖度尼安德特人基因组——进行SV鉴定。最终筛选出951个高置信度SVs，称为古人类结构变异（archaic hominin structural variants, AHSVs）。这些AHSVs约占0.47 Mb，平均片段大小为493 bp，范围从50 bp到8.5 kb，在基因组中均匀分布。约71.4%的AHSVs与重复元件重叠，其中SINE类的Alu家族和LINE类的L1家族占主导，提示转座元件可能是AHSVs形成的主要驱动力。
三、关键技术与方法
研究人员使用了以下主要技术方法：利用Paragraph（v2.3）这一基于图的SV基因分型工具对古人类样本进行NHP SV基因分型；使用CrossMap（v0.7.3）进行基因组坐标转换（GRCh38至GRCh37）；采用Ensembl Variant Effect Predictor（VEP）（v101.0）进行SV功能注释；使用BEDTools（v2.31.1）进行基因组区间重叠分析；利用clusterProfiler（v4.12.6）进行基因本体（Gene Ontology, GO）和疾病本体语义与富集（Disease Ontology Semantic and Enrichment, DOSE）分析。数据来源包括四个古人类基因组BAM文件、1000 Genomes Project的现代人测序数据（HG00098、HG00155、HG00139）、ENCODE候选顺式调控元件（cCREs）数据库以及已发表的脑组织H3K27ac ChIP-seq数据。
四、研究结果
AHSVs的鉴定与基因组特征：通过Paragraph对四个古人类样本进行基因分型，共鉴定出5,381个变异，包括4,492个插入（insertions, INSs）和889个缺失（deletions, DELs）。经筛选获得951个高置信度AHSVs。三个尼安德特人分别携带654、552和605个SVs，丹尼索瓦人携带225个。39个SVs为四个古人类样本共有，77个仅在三个尼安德特人中共享。AHSVs的纯合变异频率高于被排除的变异集。AHSVs在长度分布上富集于较短片段（DELs约250 bp，INSs约50 bp），且与群体水平人类SVs不同，未表现出端粒邻近偏好。
AHSVs的功能注释：VEP注释显示，35.6%的AHSVs为内含子变异，19.9%为特征延长（feature elongation），16.3%为特征截短（feature truncation），编码区变异仅占0.1%。约5.7%的AHSVs（N=54）与ENCODE cCREs重叠，其中远端增强子样信号占比最大。有4个AHSVs被注释为编码序列变异，涉及ENOSF1（与先天性角化不良相关）、C6orf15（与低级别胶质瘤和重度抑郁障碍相关）、KBTBD7（与脓毒性脑损伤相关）和FMO1（与心血管疾病相关）基因。
AHSV相关基因的富集分析：对464个AHSV相关基因进行GO富集分析，发现显著富集的生物学过程主要与神经元发育和活性相关，如突触信号传导、跨突触信号传导、神经元分化中的细胞形态发生和神经肌肉接头发育等。相比之下，全部NHP SV相关基因（N=21,311）的GO富集结果未呈现明显的神经相关模式，而主要富集于基础生理和代谢过程（如蛋白质分解代谢过程和DNA损伤应答），表明AHSVs作为NHP SVs的一个子集，对神经相关基因产生了强烈而集中的影响。DOSE分析显示，AHSV相关基因显著富集于神经发育和精神疾病相关表型，如精神健康发育障碍、癫痫和自闭症谱系障碍。
AHSVs对脑CREs的潜在调控影响：研究人员利用已发表的包含人类、黑猩猩和猕猴多个脑区的60,702个H3K27ac富集区域数据，发现39个AHSVs与39个脑CREs重叠。一个典型例子是位于CACNB4基因内的H3K27ac富集区域与三个尼安德特人样本中鉴定到的DEL重叠，该DEL与所有NHP物种共享，且在前额叶皮层（prefrontal cortex, PFC）中重叠区域的H3K27ac信号显著降低。CACNB4编码电压门控钙通道β亚基家族成员，与癫痫、精神分裂症和阿尔茨海默病等疾病相关。此外，研究人员还发现3个与黑猩猩共享的AHSVs与人类和黑猩猩间差异富集的CREs重叠，18个与猕猴共享的AHSVs与人类和猕猴间差异富集的CREs重叠。一个例子是ZNF831基因座内的H3K27ac富集区域与Altai和Chagyrskaya中检测到的DEL重叠，该区域中现代人的H3K27ac信号显著高于黑猩猩和猕猴。
AHSVs对染色质环锚的潜在影响：研究人员利用已发表的人类和黑猩猩神经祖细胞（neural progenitor cells, NPCs）中的染色质环数据，发现三个值得关注的例子：PDE10A、RAPGEF4和DCTN4基因区域中，每个环有一个锚在人类和黑猩猩间保守，而另一个锚在物种间位置不同；在人类中，这些人类特异性锚区域与尼安德特人和黑猩猩共享的DELs重叠。利用诱导多能干细胞（iPSCs）来源的神经干细胞（neural stem cells, NSCs）的ChIP-seq数据，研究人员发现这些区域中黑猩猩的CTCF信号和H3K27ac信号较人类降低，与DNA序列的DEL一致，可能导致这些区域在古人类和黑猩猩中无法作为环锚发挥功能。PDE10A主要在纹状体表达，在运动调控中起关键作用，其突变与多动障碍相关；RAPGEF4对静息膜电位和内向钠电流的成熟至关重要，参与灵长类前额叶皮层发育关键期的神经元成熟，并与自闭症等神经发育障碍相关；DCTN4与胶质瘤和阿尔茨海默病相关。
五、讨论与结论
研究人员指出，AHSVs的特征在于其存在于NHP和古人类中，且与现代人不同，这使其成为理解现代人适应性进化的宝贵资源。与进化距离较远的NHP物种共享的SVs可能源于平行进化、趋同进化或不完全谱系分选等因素。本研究发现的大多数与AHSVs直接重叠的基因在脑发育和功能相关生物学过程中显著富集，这一特征在全部NHP SVs中并不明显，提示AHSVs可能在塑造现代人与古人类神经相关表型差异中发挥了重要作用。在三维尺度上，研究人员观察到古人类中三个影响神经相关基因（PDE10A、RAPGEF4和DCTN4）的DELs与黑猩猩共享，并与现代人特异性环锚重叠，这可能促成人类、古人类和黑猩猩之间环锚的差异，进而影响基因表达调控模式的差异。
研究结论：本研究通过比较基因组分析，鉴定出951个区分古人类与现代人的AHSVs，这些SVs显著富集于脑发育和功能相关基因，其中许多与神经精神疾病相关。部分AHSVs与脑CREs重叠，且在PDE10A、RAPGEF4和DCTN4三个神经相关基因中，现代人特异性染色质环锚与古人类共有的DELs重合，提示AHSVs可能影响了三维染色质架构。本研究拓展了对现代人脑适应性进化分子机制的理解，为未来古人类SV图谱的扩展、理解这些重要变异在人类进化中的作用以及阐明它们如何塑造独特的现代人提供了重要资源。
