项目背景

某三甲医院儿科遗传诊断课题组,常年接收疑难罕见病 trio 家系(先证者+父母)的全外显子组测序数据,累积队列已超过 120 个家系。此前分析全部外送,每轮候选变异筛选与返修周期 2-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 天。