项目背景
某三甲医院儿科遗传诊断课题组,常年接收疑难罕见病 trio 家系(先证者+父母)的全外显子组测序数据,累积队列已超过 120 个家系。此前分析全部外送,每轮候选变异筛选与返修周期 2-4 周,遇到需要补做家系成员时排期更久。
分析流程与瓶颈
- 外包交付的是 VCF + PDF 报告,课题组无法按自己的候选基因列表灵活重筛
- 积累的家系数据分散在外包公司,联合分析(120 家系 joint call)无法开展
- 住院疑难病例等不起 4 周周期,临床医生多次提出加速需求
解决方案
课题组按 WES 算力计算器的建议配置了 16 核 / 64GB / 8TB(NVMe+HDD 分层) 的本地服务器,部署标准化流程:
| 环节 | 工具链 | 单家系耗时 |
|---|---|---|
| FASTQ→BAM(BWA-MEM2+BQSR) | GATK Best Practices | 约 4 小时 |
| 单家系变异检测 | HaplotypeCaller + GenotypeGVCFs | 约 2 小时 |
| 注释与筛选 | VEP + ClinVar/HGMD 本地库 | 20 分钟 |
| trio 筛选(AD/AR/X 连锁) | 自定义 R 脚本 | 10 分钟 |
| 120 家系 joint call(队列) | GenomicsDBImport | 周末批处理 36 小时 |
最终产出
本地化后新家系从数据到候选变异列表缩短到 7 天内(含测序周期后本地分析仅 1 天),疑难病例可按临床表型实时调整筛选策略。120 家系队列联合分析支撑了课题组一项遗传异质性研究,鉴定出两个既往被漏检的复合杂合变异。
配置建议
临床级 WES 本地化的甜点配置是 64GB 内存 + 16 核:日常 trio 足够,队列 joint call 靠批处理消化。测序量再翻倍(如加做 WGS)时再考虑 128GB 档。医院内网部署注意预留 8TB 以上存储与双盘冗余。
常见问题
WES trio 分析本地化需要多大内存?
GATK HaplotypeCaller 单样本 16GB 可跑,但 BQSR+_joint-call 多样本时 32-64GB 更稳。临床级 WES 本地化建议 64GB 起步。
本地分析和送外包比快多少?
外包一轮注释返修通常 2-4 周;本地化后新增家系当天入库、当天出 joint call,疑难病例迭代一轮缩短到 1-2 天。