多维注释多类别变异致病性预测方法

随着我国医疗健康水平的发展，许多既往临床疑难未诊断病例得到了明确的诊断，但是仍有相当一部分的病例没有得到及时准确的诊断，给病人、家庭和社会带来了沉重的负担。 全基因组测序（WGS）和全外显子测序（WES）技术的普及和机器学习相关技术的快速发展给解决上述问题带来了新的途径。高通量测序会产生大量的基因变异，外显子区域仅单核苷酸变异即超过7800万个，但其中绝大部分并不致病，只有约2%有对应临床分类，所以如何高效地在剩余变异中筛选出真正的致病变异，是一个重要且复杂的课题。 机器学习 方法可以通过利用生物数据中的模式来预测未注释变异的致病性，从而填补这个变异解释的空白。已有许多基于机器学习的计算方法已被开发用于预测变异致病性，但是已有算法中的大部分只能专注于其中某一种生物学机制或某一种变异，在临床诊断中有较大局限性，限制了疾病特别是罕见病的诊断率，难以开发或应用针对潜在遗传因素的临床治疗。该算法是一个预测外显区域内除同义变异外其他类型变异的广泛适用模型，增强了致病性预测领域多类别变异整合预测的能力。将变异预测的结果与临床表型等结合能够在临床诊断、筛查疾病等方面具有重要价值。 

由于导致遗传性疾病的不同种类致病变异通过多种生物学机制发挥作用，且算法需要面对多种不同的变异类型进行致病性的预测，研究团队开发了 MAGPIE 算法。MAGPIE算法主要分为三个阶段：多维特征注释，带有分离特征选择的自动特征工程，以及模型分布训练。 

第一步将来自于ClinVar的经过验证的变异与来自gnomAD的良性罕见变异进行多维特征注释，特征类别包括：1）表观基因组学、2）功能效应、3）剪接效应、4）基于人群的特征、5）生化特性、6）保守性。这里加入gnomAD中罕见良性变异是为了消除原有ClinVar中数据中可能存在的bias，即数据中人群等位基因频率与变异致病性高度相关。 第二步，研究团队使用Oba等人提出的基于贝叶斯PCA的缺失值估计方法对经过注释的数据进行缺失值填补并通过一系列条件进行过滤筛选得到候选训练数据集，应用特征工程将特征数量扩充至3,000个以上。 为了避免过拟合问题，研究团队将所有特征分为两个类别，即核心集和附加集。前者包括与致病性直接相关的特征，如基于人口统计学的特征、保守性和功能效应。这些特征在训练阶段被全部保留下来。后者包括其他特征和自动特征工程生成的特征。研究团队将使用这些特征进行不超过50轮的训练，将特征减少到200个或更少。随后，在每一轮之后，研究人员评估每个特征的重要性，并丢弃相对重要性得分低于0.001的特征。研究团队定义了五轮训练，每轮使用5折交叉验证，对于Light GBM中各参数进行分步调节，并在保证准确性和分类性能指标没有下降的条件下降低模型的过拟合问题，增强其泛化能力。 梯度提升方法的应用使得MAGPIE成为一个可解释的机器学习模型。为了更好地理解模型学习到的表示，研究团队计算了不同特征模态的重要性，并以网络视图的形式呈现。MAGPIE中的特征重要性表示每个特征对模型的贡献程度，并反映了它们的信息增益。在分类过程中，MAGPIE集成了六种特征模态。其中，重要的特征组是功能效应，包括功能丧失得分（LoF_score）、人类基因损伤指数（GDI）等指标。之前研究结论为功能丧失变异更有可能导致疾病，与之相符的是研究人员发现MAGPIE预测的致病变异的LoF得分低于良性变异。GDI是每个编码蛋白质的人类基因的累积变异损伤，高度受损基因的变异更不可能导致疾病。研究团队观察到GDI特征被赋予了很高的权重，表明GDI较低的基因倾向于具有较高的MAGPIE预测分数，并更可能是致病的，与先前发表的研究一致。 

研究团队分别在三个不同的数据集，分别是ClinVar与gnomAD混合数据集independent test set、SwissProt与gnomAD混合数据集orthogonal validation set、由ACMG指导的临床验证数据集ACMG-guided dataset中对MAGPIE的性能进行了评估，其中independent test set为类别平衡数据集，即阳性样本与阴性样本数基本一致；orthogonal validation set中阳性样本占比超过90%，ACMG-guided dataset中阳性样本占比不足20%。研究人员发现MAGPIE可以在多种类别分布的数据集中稳定的实现0.95以上的AUC，0.88以上的AUPRC，以及0.9以上的准确率（accuracy）。MAGPIE同样在各数据集中准确预测所有变异类型，能够在移码变异（fr[x]ameshift）、提前终止变异（stopgain）等变异类型中有均超过0.85的准确率。 

为了检测MAGPIE在更符合真实临床数据下的性能表现，研究团队在上述三个数据集上筛选出罕见变异（AF<0.01）数据集并测试MAGPIE的性能。一般来讲，在罕见变异数据集上各方法的性能会略差于未经筛选的数据集，但相比于其他方法，MAGPIE在上述三个罕见变异数据集中拥有小的性能削弱，并能够在各罕见变异数据集中达到0.9以上的AUC、0.85以上的AUPRC，与0.85以上的准确率. 

此外，研究团队为了进一步验证算法的性能，将MAGPIE应用于四个基因（ATP7B、CFTR、FBN1和LMNA），这些基因中包含大量已知经过实验验证的导致不同孟德尔遗传病的致病变异。对于上述基因，与其他方法的默认参数相比，MAGPIE能够识别出多已知的致病变异，准确率平均在85%以上。例如，对于FBN1基因，该基因变异导致马凡综合征（Marfan syndrome），并包含1621个致病变异。MAGPIE成功预测出了1415（87%）个候选变异为致病。对于与囊性纤维化相关的CFTR基因，MAGPIE预测出了164（95%）个变异为致病。与MutationTaster、VEST4和其他工具相比，MAGPIE分数集中在1附近，这表明MAGPIE对致病变异的预测有更高的置信度。