中国科学院北京基因组研究所（国家生物信息中心）研究员陈华团队在一项新研究中构建了比较群体基因组学新算法（CEGA）。相关结果近日发表于《基因组生物学》期刊。

随着基因组测序技术的发展，物种和群体水平基因组数据呈指数增长。这些数据为从基因组水平鉴定和解析自然选择机制提供了前所未有的机遇。但目前的分析方法仍面临技术瓶颈，其中一个关键问题是如何高效、准确地检测作用于非编码区的自然选择效应，同时高效、高性能地分析多物种大样本数据。

陈华团队在多物种联合等位基因频谱理论以及哈德森-克雷特曼-阿瓜德生物统计学检验（HKA）的框架上开发了CEGA，新算法整合微进化过程与宏观进化过程模型，有效刻画自然选择和群体历史在非编码区形成的遗传多态性“印记”，可高效、准确地检测作用于非编码区上的正向选择及平衡选择信号。CEGA同时分析物种间的分歧位点和物种内的多态位点信息，当两物种分化时间比较短时，多态位点蕴含的信息有助于准确地推断分化时间、有效群体大小等信息，从而有利于区分自然选择效应与群体历史干扰，因此该方法在不同物种分化时间尺度上具有更广泛的适用性。

研究者通过仿真分析表明，对于不同的选择强度以及物种分化时间，CEGA检测正选择及平衡选择的效果均优于现有方法。尤其对于选择强度较弱或者物种分化时间比较短的情景，CEGA的优势更为明显。除了用于检测自然选择外，研究者往往希望提供对自然选择发生过程的深入认识。 鉴于此，CEGA还基于群体遗传学模型提供了对自然选择强度等关键参数的推断。

研究团队将CEGA应用于已发表9个现代人类及9个黑猩猩的群体基因组数据，进行了编码区、非编码区两个层面上的比较分析，鉴定了在人类基因组中受自然选择作用而快速进化基因，并发现这些基因的功能显著富集在与大脑容量、大脑皮层的总面积以及大脑皮层的厚度等相关表型和分子通路。此外，在与免疫反应和病原体抵抗相关的区域（如主要组织相容性复合体MHC）存在显著的平衡选择信号。以上仿真分析以及人与黑猩猩基因组真实数据分析的结果表明，CEGA是一种有效的算法工具，可用于大规模群体基因组测序数据的高效分析。

相关论文信息：

https://doi.org/10.1186/s13059-023-03068-8