哈佛与牛津开发基因致病性预测模型，现已成功预测 3600 万个致命基因突变

直接从基因层面预测疾病，这一直是近现代医学研究的主要方向之一。

然而，全体人类的基因变异体数量远超现有的探测技术，甚至仅仅是不同个体的蛋白质区编码也会展现出巨大的差异性。

因此，超过 98% 的基因变异给人体带来的影响依旧是未知且无法预测的。

但最近，来自哈佛医学院和牛津大学的科学家合作开发了一种 AI 模型，成功预测了 3219 个疾病基因中超过 3600 万个变体的致病性，并将超过 25 万个未知变体进行了归类。

这项研究现已登上 Nature。

"从进化中预测致病性"

其实，现在临床上已有用于预测基因变异影响的模型。

但这些模型往往是在经过标注的临床数据集上进行有监督学习，一旦进入现实场景，标签偏差、标签稀疏以及噪音就会造成其准确率的下降，并不能作为基因变异体分类的可靠依据。

而这次的研究团队提出了一个叫做 EVE（Evolutionary model of Variant Effect）的模型。

这是一个仅根据进化序列训练的无监督生成模型。

模型预测变异基因的致病性主要分为两步：

第一步，使用变型自动编码器 VAE 来学习蛋白质的氨基酸序列分布。

学习了多个领域的复杂高维分布之后，模型就捕捉到了进化过程中的自然序列约束，包括各种位置之间的复杂依赖关系。

再从得到的近似后验分布（Approximate Posterior Distribution）中取样，评估每个单一氨基酸变体相对于野生型的相对可能性。

这种相对可能性被称为"进化指数"，与临床标签进行比较后发现，区分致病性和良性标签的数值在不同的蛋白质中是一致的，这说明无监督的方法能够有效推断致病性。

第二步，在所有单一氨基酸变体的进化指数分布上拟合了一个双组分（two-component）的全局-局部高斯混合模型。

这一步的输出是在区间 [0，1] 内定义的连续致病性值，0 代表良性，1 代表致病性。

然后将 EVE 模型运用于 ClinVar 数据库中的 3219 个人类基因上，得到的结果图中的平均曲线面积（AUC）为 0.91，说明 EVE 模型对绝大多数的基因变异都能做到具有临床意义的预测：

优于已知模型，与实验预测效果一致

研究团队也将 EVE 模型与已知的模型进行了对比，可以看到，在预先确定已知的已标注临床数据的预测上，其效果优于同类计算模型：

那么这样一个 AI 计算模型与用于预测致病性的经典方法 -- 深度突变扫描实验（Deep Mutational Scan Experiment）相比效果又如何呢？

对比实验后可以看到，EVE 模型在临床预测方面的总体表现与经典方法效果基本一致：

而当从 ClinVar 数据库中选择一组数量规模更大，但高质量标注较小的数据时，EVE 模型的表现甚至更好：

哈佛 & 牛津合作出品

这篇论文有三位共同一作，其中 Jonathan Frazer 和 Mafalda Dias 都来自哈佛大学的系统生物学，他们同时也是 Marks Group 实验室中的一员。

而 Pascal Notin 则是来自牛津大学的计算机科学专业的博士生，主要研究领域包括贝叶斯深度学习、生成模型、因果推理和计算生物学的交叉领域。

论文链接：

https://www.nature.com/articles/s41586-021-04043-8