研究用 MPRA 与 CROP-seq 系统解析 727 个 NDD 非编码风险变异的功能

神经发育障碍（NDD）在遗传上高度多基因。常见变异方面，大规模孤独症GWAS已定位出五个全基因组显著位点，并提示不同亚型间的多基因架构存在异质性^[3]；编码区研究也通过大样本外显子测序锁定102个风险基因，其功能多集中于基因表达调控与神经元通讯^[2]，新发蛋白截断变异的作用亦借助人群参照队列得到更精确的界定^[4]。然而，关联信号大多落在非编码区，功能解释严重滞后：等位基因特异性开放染色质研究提示，此类变异可能通过改变转录因子结合、经染色质远程接触作用于远端基因，但究竟哪些变异真正具备调节活性、靶基因是谁，缺乏系统检验^[5]。领域内也早有呼吁，应在相关的神经发育模型中对风险位点做功能刻画^[1]。

这篇预印本正是对这一难题的回应。McAfee等人用大规模并行报告基因检测（MPRA）在人神经祖细胞中检测了727个NDD关联的非编码变异，鉴定出38个具有调节活性的表达调节变异（emVars），并预测其改变与神经发育信号及细胞代谢相关的转录因子结合基序；再结合表达数量性状位点（eQTL）与activity-by-contact（ABC）模型，把emVars映射到候选靶基因；最后用CROP-seq在人多能干细胞分化的神经元中鉴定内源靶基因并刻画其下游转录效应。结果显示，CROP-seq优先支持ABC模型的预测，同时揭示了额外的调节关系，且下游效应汇聚于涉及神经元连接、突触功能与线粒体生物学的发育调控程序。

这项工作把非编码风险变异从统计关联推进到“变异—调节活性—内源靶基因—共享转录程序”的联系链条，其揭示的汇聚通路也与外显子测序研究中风险基因集中于基因表达调控和神经元通讯的功能图景相呼应^[2]，并为非编码变异的解读提供了从体外报告基因筛选、三维基因组建模到内源CRISPR验证的组合路径。该团队在此方向已有积累：第一作者2023年曾系统考察等位基因调节活性^[6]，团队2026年又发表了大规模基因调节活性评估^[7]。

当然，727个变异中仅38个在神经祖细胞条件下检出活性，其余究竟是确实无功能，还是受限于检测灵敏度与所测细胞状态，仍待分辨；emVars与靶基因的关系在更成熟的神经元乃至体内发育情境中能否重现，这些汇聚的转录程序又如何转化为疾病风险，都是悬而未决的后续问题。