饱和突变诱变能否对意义不明的 RAD51D 变异进行分类？

存活读出把变异转化为功能评分

这项研究要回答的问题是：实验室测量能否区分有害的 RAD51D 错义变异与无害变异。RAD51D 是参与同源重组（修复 DNA 双链断裂的通路）的抑癌基因，其致病变异与乳腺癌和卵巢癌相关。患者中发现的错义变异大多仍属意义不明变异，临床医生因此缺少对其分类的功能依据。

作者采用多重变异效应检测。他们构建的文库覆盖了 6,888 个可能的 RAD51D 单氨基酸变异中的 97.6%，每个细胞最多携带一份稳定整合的变异拷贝，宿主是经 CRISPR 敲除内源 RAD51D 的 U2OS 细胞系。细胞在 250 nM 奥拉帕利（一种 PARP 抑制剂，可杀死无法通过同源重组修复损伤的细胞）中处理两个传代。选择前后用深度扩增子测序读取变异频率，并为每个变异计算功能评分；评分为负表示该变异在选择中被耗竭，即未能恢复修复功能。

这一规模的意义在于，它把逐个变异的生化问题变成了覆盖几乎整个蛋白的排序测量。评分本身不是临床结论：它是在一个细胞系、一种药物选择条件下得到的功能读出，作者将其作为分类框架中的一项证据，而非独立诊断。

评分阈值把变异分为三类功能状态

为了让连续评分可用，作者用对照变异设定阈值。不改变氨基酸的同义变异定义中性端，截短蛋白的无义变异定义功能缺失端。评分在 −0.170 至 −0.446 之间的变异被判为中间型，低于 −0.446 的判为功能缺失，前提是同时通过局部错误符号率（local false sign rate）的显著性阈值。

对照变异的表现符合预期。305 位残基之前的终止获得型变异平均评分为 −0.989，同义变异平均为 0.023，两组在 494 个变异上被完全分开，精确率-召回率曲线下面积为 1.00。306 位残基及之后的截短变异评分与同义变异相当，作者将其解释为该 C 端区域大多无结构，在本检测中并非必需。在错义变异中，74.5%（4,593 个）为中性，平均评分 −0.07；16.8%（1,038 个）为功能缺失，平均 −0.78；8.6%（536 个）为中间型，平均 −0.31。

这些类别边界是后续临床比较所依赖的单位。被判为中间型的变异与功能缺失变异并不等同，论文在统计有多少意义不明变异可被重新分类时对二者区别对待。

已知致病与良性变异落在图谱两端

对图谱最强的检验，是与已有临床分类的变异是否一致。作者从 ClinVar 整理了 211 个蛋白层面的唯一变异，包括 155 个良性或可能良性条目和 36 个致病或可能致病条目，另有 28 个提交结果相互冲突的条目，作者将其中的 1 个判为致病、27 个判为良性。全部 37 个致病或可能致病变异都评为功能缺失，全部 174 个良性或可能良性变异都评为中性。

个别案例也支持同一图景。常见多态性 R165Q 评为中性，与其在 gnomAD 中 14.1% 的次要等位基因频率一致。此前被列为致病或可能致病的两个错义变异 S207L 和 S207W 都被判为功能缺失，其中 S207L 较温和，评分为 −0.503，S207W 更严重，为 −1.703。

作者还把评分换算为 ClinGen 的 OddsPath 比值，用来量化功能检测区分已知致病与已知良性变异的能力。功能异常变异的比值为 174，功能正常变异为 0.027（即 1/37），分别对应支持致病和支持良性的强证据等级。正是这一换算框架，使该检测能在临床分类中贡献分值，而不只是描述生化现象。

多数意义不明变异被评为功能正常

最具临床意义的检验对象，是 ClinVar 中列为意义不明或结果冲突的 789 个 RAD51D 错义变异。作者排除了 49 个预测影响剪接的变异，以及 68 个中间型或未通过显著性阈值的变异，剩下 672 个可被赋予方向的变异。其中 87% 评为功能正常，作者指出，在没有冲突证据的情况下，这足以将它们推向可能良性。其余 85 个变异（12.6%）评为功能缺失，可获得支持致病分类的分值。

这一分布是论文对分类问题的实际回答：该检测并不能重新分类每一个意义不明变异，但能为其中大多数给出一个明确方向的强推动。作者还指出，85 个功能缺失变异中近三分之二未出现在 gnomAD 中；在已出现的变异中，累积次要等位基因频率为 2.03 × 10−4。他们将该数字作为下限估计，认为 RAD51D 错义功能缺失变异可能使约 1/20,000 的个体面临遗传性乳腺癌和卵巢癌风险。

这一风险数字是从人群频率推断而来，并非实测的外显率，且适用于错义功能缺失这一子集，而非全部 RAD51D 变异。重新分类的表述也依赖 ClinGen 的分值体系，因此该检测提供的是证据，其权重由该框架设定。

损伤集中在 DNA 结合区与 ATPase 核心

把评分映射到结构上，可以看出替换在哪些位置伤害最大。埋藏残基上的变异平均比表面暴露残基上的变异更有害，平均评分分别为 −0.272 和 −0.113。但筛选也发现，RAD51C 或 XRCC2 界面附近的表面暴露残基上存在高度有害的变异，这是简单的“埋藏对暴露”规则会漏掉的。

若干残基几乎不容忍任何替换。在 Walker A 基序中，保守的苏氨酸 T114 只容忍丝氨酸替换。在 Walker B 的天冬氨酸 D206 处，19 种替换中有 18 种产生强功能缺失，平均评分 −1.22，化学性质相近的谷氨酸较温和，为 −0.52。单链 DNA 骨架附近 263–267 位的残基强烈排斥谷氨酸和天冬氨酸替换，作者将其归因于与带负电骨架的排斥，而甘氨酸 265 完全不容忍任何替换。

这些模式识别出有害变异的热点，也赋予图谱机制层面的解读：在本检测中，DNA 结合环和 ATPase 核心是氨基酸改变最严重损害同源重组的区域。

RAD51D 似乎抑制 BCDX2 的 ATP 水解

生化实验要回答的是 RAD51D 在其复合物中究竟做什么。RAD51D 与 RAD51B、RAD51C 和 XRCC2 组成异源四聚体 BCDX2，同时还与 XRCC3、RAD51C 和 XRCC2 组成第二个复合物 X3CDX2。只有 RAD51B–RAD51C 这一对水解 ATP；RAD51D–XRCC2 这一对通过一种此前未知的机制调节该活性，这正是论文要阐明的。

使用纯化的亚复合物，作者发现把 RAD51D–XRCC2 加入 RAD51B–RAD51C 会降低 ATP 水解。单链 DNA 结合缺陷的变异无法产生这种降低，水解水平与单独的 RAD51B–RAD51C 相当。因此 BCDX2 内部的 RAD51D–RAD51C 界面似乎是调节 ATPase 活性所必需的，作者据此提出一个假设：RAD51D 的主要功能反而是减慢 BCDX2 的 ATP 水解，从而为 RAD51 核蛋白丝组装留出足够的时间和空间。

这是作为假设提出的，而非已证实的体内因果链。但它解释了为什么在结构表面图上看似无害的界面变异仍可能有害：它们破坏的是调节，而不是折叠。

一个等位基因解耦 BCDX2 与 X3CDX2 活性

筛选还发现了一个把两种功能分开的变异。RAD51D-L4H 在 BCDX2 内未见明显生化缺陷，也仍能形成 X3CDX2 复合物，但它在 X3CDX2 内结合 DNA 更弱，并且无法在 RPA 包被的单链 DNA 上促进 RAD51 介导的链交换。作者将其描述为一种功能分离等位基因，可能解耦 BCDX2 与 X3CDX2 的活性。

在链交换实验中，野生型 X3CDX2 能在 RPA 包被的 DNA 上促进 RAD51 介导的交换，而受试的 X3CDX2 突变复合物不能，这与细胞检测中观察到的同源重组缺陷一致。在检查的 7 个癌症相关有害变异中，有 6 个改变了 RAD51B–RAD51C 的 ATPase 活性调节，方式或是削弱 RAD51D 与 RAD51C 的相互作用，或是降低 BCDX2 的 DNA 结合；其中 6 个还影响了 X3CDX2 的形成。

要点在于，RAD51D 至少通过两条途径支持同源重组：调节 BCDX2 的 ATPase 活性，以及支持 X3CDX2 的链交换活性。因此一个变异的功能评分可能反映其中任一条途径，这也是为什么在混合存活读出之外，生化随访同样重要。

把图谱当作强证据而非最终结论

论文中最具决策意义的比较，是该检测两个方向的证据。功能缺失评分带来 4 分支持致病，中性评分带来 4 分支持良性，但中间型评分和未通过显著性阈值的变异被排除在意义不明变异的统计之外，而不是被强行归入某一类。这一排除是读者应留意的边界：87% 和 12.6% 描述的是 672 个可解读变异，而不是全部 789 个意义不明条目。

与计算预测的比较进一步凸显了重点。在被 AlphaMissense 高置信度判为致病的变异中，有 35% 在本检测中功能正常，作者将其解读为该预测工具可能存在较高的假阳性率。对 RAD51D 而言，实测功能评分因此携带了基于序列的预测所没有的信息，论文自身的表述也把最大权重放在实验评分给出的功能缺失和中性判定上。

因此实际优先级是：在图谱给出明确方向判定时使用它，把中间型和未达显著性的变异视为未解决，并把生化机制理解为界面变异和 ATPase 核心变异为何如此表现的解释，而不是另一项独立的临床检测。

本文引用的文献