第4章 PVS1(无功能变异)证据详解与复现验证

PVS1是ACMG/AMP指南中力度最强的致病性证据之一。说白了,就是认定某个变异完全摧毁了基因功能——要么蛋白截短,要么移码乱套,要么关键剪接位点直接报废。但力度越强,责任越大,误判的代价也越高。这一章我们就把PVS1的底裤扒干净,讲讲怎么才算“无功能”,以及如何可靠地、可复现地拿到这条证据。

4.1 到底什么样的变异配得上PVS1?

ACMG官方给的框架很清晰:无功能变异(loss-of-function variant)是指那些预期会导致基因产物完全缺失或功能完全丧失的变异。

我个人习惯把它拆成三层核查:

  • 变异类型:无义、移码、起始密码子丢失、经典±1或±2剪接位点变异、单/多外显子缺失。
  • 基因约束:该基因对LOF变异是否敏感?LOF是已知致病机制吗?
  • 例外排除:有没有可能绕过无功能机制?比如NMD逃逸、选择性剪接、基因冗余。

我在项目里遇到过不止一次,明明是经典剪接位点+2T>C,初看PVS1妥妥的,结果RNA测序发现一个隐蔽的隐秘剪接位点被激活,产生了仍保留部分功能的截短蛋白。嗯,这时候只能降级到PS3甚至PM4了。

核心原则: 拿到PVS1之前,先问自己三个问题——类型对不对?基因受不受LOF?有没有旁路?

4.2 PVS1证据的逐层拆解

我们按变异类型逐一梳理。表格里列出了最常见的几类,我把判断要点和常见雷区都列上。

变异类型 PVS1适用条件 常见陷阱
无义突变(nonsense) 位于最后一个外显子-50 bp之前,预期触发NMD。 若位于倒数第二外显子的50 bp内,可能逃逸NMD,产生C端截短蛋白,功能可能残留。
移码突变(frameshift) 插入/缺失长度不是3的倍数,导致阅读框改变,提前出现终止密码子。 同样存在NMD逃逸风险,特别当终止密码子位于最后一个外显子时。
经典剪接位点(±1/2) 内含子-外显子边界的GT/AG破坏。 可能存在隐秘剪接位点,需RNA验证。
起始密码子丢失 ATG变为其他密码子,无法启动翻译。 下游是否存在替代起始位点?若有,可能产生N端截短但仍部分有功能的蛋白。
整外显子/全基因缺失 通过qPCR、MLPA或NGS覆盖度确认为纯合/杂合缺失。 注意区分假缺失(多态性导致比对失败)。
曾经踩过的坑: 有个无义突变位于倒数第二个外缘子靠近3′端,预测NMD逃逸,但我们RNA实验确实检测到少量全长转录本。最终只给了PM4(蛋白质功能改变证据),没有上PVS1。永远不要只看基因组坐标推断NMD。

4.3 PVS1的知识图谱:评估决策流

下面这张图我整理成了工具卡,每次做PVS1评估时直接对标。

PVS1证据评估决策流程 PVS1证据评估决策流程 变异类型:无义/移码/剪接/缺失 基因约束:LOF为已知致病机制 位置核查:是否触发NMD? 例外排除:NMD逃逸/隐秘剪接/替代起始 PVS1 是否适用? 不适用,降级至PM4/PS3 适用,PVS1成立 ⚠️ 如果RNA实验证实功能缺失,即使基因组预测有旁路,仍可考虑PVS1 但必须记录实验细节和对照数据
使用建议: 我一般把这个流程贴在实验室的共享wiki里,每次评估前走一遍,防止漏掉例外检查。

4.4 复现验证的硬门槛

ACMG解读最怕什么?不同实验室对同一变异给不同证据。PVS1尤其容易打架。我整理了几条必须严格复现的验证步骤。

  1. 变异注释一致性:统一使用GRCh38/hg38参考基因组,VEP或SnpEff版本锁定。代码示例:
vep --cache --everything --assembly GRCh38 \
    -i variants.vcf -o variants_vep.txt --fork 4
  1. NMD预测标准化:使用NM_转录本编号,取MANE Select。不能随意选一个异构体。预测NMD时按≤50 bp规则判断。
# 获取最后一个外显子坐标
grep "NM_001271208.2" gencode.v38.annotation.gtf | \
awk '$3=="exon" {print $1,$4,$5,$7}' | sort -k2 -n | tail -n1
  1. 剪接预测工具结果归档:SpliceAI、MaxEntScan、dbscSNV至少用两个,截断值统一。我曾经因为SpliceAI和MaxEntScan结果矛盾,最后加做了mini-gene实验才定论。
# SpliceAI 输出示例(阈值0.2)
Gene    Variant        DS_AG  DS_AL  DS_DG  DS_DL
BRCA2   chr13:32332545 0.89   0.02   0.03   0.01
  1. 群体频率核查:gnomAD v4.0,过滤高通量测序批次效应。如果变异在东亚人群中频率>0.001,立刻警惕——很可能不是真正的LOF。
曾经有一例: 某个无义变异在gnomAD v2.1中频率0.0005,但在v3.1中因为测序深度提升,发现是测序错误。还好我们复查了原始BAM文件,避免了一次假阳性PVS1发布。

4.5 实操避坑清单

下面这些是我在多年评估工作中总结出来最容易出问题的地方,写下来供你参考:

  • 转录本选择差异:不同实验室用不同NM号,导致NMD预测结果相反。我建议团队统一使用MANE Select + ClinGen推荐转录本。
  • 剪接位点标注错误:有些数据库会把内含子+1G>A标成“剪接位点变异”,但实际距离外显子边界超过2bp,不应给PVS1。
  • 起始密码子下游的替代起始:有个基因ATG下游90 bp处存在另一个in-frame ATG,我们通过Western blot发现表达量只下降了40%,最终只给了PVS1_strong(降级)。
  • NMD逃逸的实验验证:纯干预测不够,有条件就做cycloheximide处理后的qPCR,看转录本丰度是否变化。
给我自己的检查清单: 每个PVS1候选变异,我都会在共享文档里贴出:①变异类型截图②NMD预测结果③剪接预测汇总④RNA实验有无。这样同事review时一目了然。