基因测序解读天书-基因测序分析师的年终总结

基因测序，解读天书——基因测序分析师的年终总结测序仪的荧光信号在实验室的暗室里闪烁了365天，我面前的分析队列里还躺着最后3份2026年的临床全外显子组测序样本。今年实验室累计完成12760例测序分析任务，其中罕见病诊断样本4217例，肿瘤伴随诊断样本5832例，健康风险筛查样本2021例，科研合作样本690例，整体分析准确率达到99.92%，比去年提升了0.07个百分点。这串数字背后是312例罕见病患儿的明确诊断，是2146例肿瘤患者获得精准用药方案的依据，是127个家庭避免了再发遗传病的风险，也是我作为基因测序分析师一年来最厚重的工作注脚。今年第一季度处理的一例脊髓性肌萎缩症（SMA）疑似病例让我印象格外深刻。患儿只有10个月大，已经出现了明显的肌张力减退和运动发育迟缓，外院临床高度怀疑SMA，但常规PCR检测只发现SMN1基因外显子7的杂合缺失，按照常规诊断标准无法确诊，临床医生把样本送到我们实验室做补充测序分析。拿到下机数据的第一时间，我先核对了测序深度：SMN1基因区域的平均测序深度达到287×，覆盖度100%，不存在大片段缺失漏检的可能。接下来比对SMN1和SMN2的同源序列时，我发现常规分析流程默认过滤的一个位于外显子7剪接位点的单核苷酸变异（SNV），恰恰出现在患儿仅存的那条SMN1等位基因上。这个变异的频率在人群数据库里只有0.0012%，之前没有被正式列入致病变异目录，但功能预测显示它会完全破坏剪接供体位点，导致SMN1蛋白完全失去功能。我没有直接按照流程标注为“意义未明变异”，而是调取了患儿父母的测序数据进行家系验证，发现父亲携带该位点杂合变异，母亲携带SMN1外显子7杂合缺失，完美符合常染色体隐性遗传模式。我又联系了实验室的功能验证团队，通过迷你基因实验确认了该变异确实会导致外显子跳读，最终给临床出具了明确的阳性诊断报告。患儿在确诊后第三周就用上了靶向治疗药物，上个月随访的时候已经能完成抬头的动作，这也是我今年收到的最有分量的临床反馈。罕见病样本的分析难点永远在于变异的致病性判定。今年我们实验室更新了ACMG（美国医学遗传学与基因组学学院）变异分类指南的本地化细则，我作为规则制定小组的成员，参与了127个中国人群特有变异的致病性证据梳理工作。之前通用的数据库里，东亚人群的变异数据只占不到15%，很多在中国人群里相对常见的良性变异会被误判为疑似致病变异，给临床诊断带来不必要的困扰。比如我们今年统计到的位于GJB2基因的c.109G>A变异，在东亚人群中的携带率达到1.87%，之前的数据库里标注为疑似致病变异，但我们收集到的32个纯合变异携带者里，有28个没有任何听力损失的表型，最终结合功能实验证据，我们把这个变异的致病性等级下调为“可能良性”，仅当合并另一个明确致病变异时才考虑致病性。今年全年我们通过本地化的变异判定规则，一共纠正了117例之前可能出现的误诊，让39个家庭避免了错误的生育指导，这比单纯提升分析速度更有价值。肿瘤伴随诊断的分析工作又是另一种逻辑。和罕见病的生殖系变异分析不同，肿瘤样本的分析需要区分体细胞变异和生殖系变异，还要计算肿瘤突变负荷（TMB）、微卫星不稳定（MSI）这些和免疫治疗相关的指标，每一个数值都直接关系到患者的治疗方案选择。今年3月我们承接了某三甲医院的晚期非小细胞肺癌队列研究样本，其中有一例患者的TMB值在第一次分析时只有8.2个突变/Mb，刚好低于免疫治疗获益的阈值（10个突变/Mb），但临床医生反馈患者的临床表型更符合高TMB的特征。我重新调取了样本的白细胞对照数据，发现该患者的造血系统存在克隆性造血突变，常规分析流程把这些克隆性造血来源的变异当成了生殖系变异过滤掉了，导致TMB计算值偏低。我调整了过滤规则，把等位基因频率在5%-30%之间的造血系统相关变异单独拎出来进行人工校验，最终重新计算的TMB值达到14.7个突变/Mb，符合免疫治疗的适应症。患者后续接受了PD-1抑制剂单药治疗，6个月后复查CT显示肿瘤缩小了62%，达到了部分缓解的标准。这件事也推动我们实验室更新了肿瘤TMB计算的标准流程，增加了克隆性造血变异的校正步骤，今年全年通过这个校正流程，一共为124例患者修正了TMB计算结果，其中37例患者因此调整了治疗方案。今年我们实验室引入了纳米孔长读长测序平台，我作为第一批接受培训的分析师，花了三个月的时间吃透了长读长数据的分析方法。短读长测序对于大片段的结构变异、重复区域的变异检测始终存在短板，而长读长测序刚好可以弥补这个缺陷。今年7月我们遇到了一例疑似脆性X综合征的患儿，常规短读长测序只检测到FMR1基因的CGG重复数为58次，属于灰区范围，无法确诊。我们用纳米孔长读长测序对该样本进行重测，直接读取到了完整的CGG重复区域，发现患儿的CGG重复数实际达到了212次，属于全突变范围，而且还存在重复区域的甲基化异常，直接明确了诊断。更重要的是，我们通过长读长测序还发现患儿母亲的FMR1基因存在CGG重复的前突变（重复数78次），而且是嵌合型，之前的短读长测序完全没有检测到这个嵌合变异。今年下半年我们用长读长测序一共解决了42例短读长测序无法确诊的疑难病例，其中包括17例动态突变类疾病、12例复杂结构变异导致的罕见病、13例肿瘤复杂融合基因的鉴定。现在我已经牵头编写了实验室的长读长测序数据分析SOP，培训了8名新的分析师掌握这项技术，明年我们计划把长读长测序的临床应用范围进一步扩大到产前诊断领域。分析流程的自动化和质控体系的升级也是今年的重点工作。去年我们的分析流程人工干预率还有12.3%，今年我牵头开发了3个自动化校验脚本，针对变异的家系共分离验证、常见假阳性变异的自动过滤、临床表型与变异的匹配度自动评分这三个场景，实现了自动化处理。今年全年的分析人工干预率已经下降到了4.7%，分析周期从原来的平均7个工作日缩短到了4.2个工作日，对于加急的产前诊断样本，我们甚至可以做到24小时内出具初步分析报告。质控方面，今年我们参加了12次国家卫健委临床检验中心组织的室间质评，所有项目的合格率都是100%，其中罕见病变异识别、肿瘤体细胞变异检测两个项目的得分在全国所有实验室里排名前5%。我还整理了今年所有出现过的分析偏差案例，一共23例，每一例都做了根因分析，其中17例是样本制备阶段的问题，6例是分析流程的参数设置问题，针对这些问题我们一共更新了11项质控规则，明年的目标是把分析准确率进一步提升到99.95%以上。当然今年的工作也不是没有遗憾。有一例DMD（杜氏肌营养不良症）的胎儿产前诊断样本，我们在分析的时候发现胎儿携带母亲来源的DMD基因外显子45-52杂合缺失，按照常规判断，男性胎儿应该会发病，所以临床医生根据我们的报告建议孕妇终止了妊娠。但后来我们对引产组织进行验证的时候发现，胎儿存在非常高比例的体细胞嵌合，正常的DMD基因占比达到了78%，这种情况下患儿出生后很可能只是轻症携带者，甚至不会出现明显的临床表型。虽然我们的分析流程本身没有问题，产前诊断的局限性也提前告知了家属，但这件事还是让我难受了很久。之后我查阅了大量的文献，开发了一个针对产前样本的嵌合变异定量分析流程，对于男性胎儿的X染色体连锁遗传病变异，只要检测到嵌合信号，就会自动启动深度分析和验证，目前这个流程已经在实验室试运行了3个月，已经发现了2例之前可能漏检的嵌合病例。还有一个遗憾是，今年我们的科研合作产出没有达到预期。年初我们计划和高校合作发表2篇关于中国人群特有致病变异的研究论文，但因为临床样本分析任务太重，大部分时间都扑在日常检测上，科研数据的整理和分析进度滞后了很多，目前只有1篇论文在投稿阶段，另一篇还在数据整理阶段。明年我打算每周抽出半天的时间专门做科研相关的工作，把这两年积累的12000多例中国人群的测序数据好好挖掘一下，尤其是那些意义未明的变异，很多都可能是中国人群特有的致病变异，把这些变异的致病性搞清楚，能给以后的临床诊断带来很大的帮助。实验室的测序仪还在不停运转，下机队列里已经排上了2027年的第一批样本。基因这本天书，我们现在能读懂的内容还不到10%，但每多读懂一个字符，就可能多救一个患者，多帮一个家庭。明年的目标很明确：把长读长测序的分析流程磨得更成