《Molecular and Cellular Probes》：Pooled Multiplex Amplicon Massive Parallel Sequencing for High-Throughput Confirmation of Whole-Genome Sequencing Potentially Pathogenic Variants in Patients with Parkinsonism
编辑推荐：
•开发并验证一种基于混合多重扩增子大规模平行测序（maMPS）的快速、高通量全基因组测序（WGS）变异验证策略，仅需微量DNA投入。
•针对86例帕金森综合征患者及对照中的307个变异进行靶向评估，将DNA样本池化为16个DNA样本混合物，实现经济高效的正交确
•开发并验证一种基于混合多重扩增子大规模平行测序（maMPS）的快速、高通量全基因组测序（WGS）变异验证策略，仅需微量DNA投入。
•针对86例帕金森综合征患者及对照中的307个变异进行靶向评估，将DNA样本池化为16个DNA样本混合物，实现经济高效的正交确认。
•验证准确性高，成功确认96%以上的可评估单核苷酸变异（SNVs）和小片段插入缺失（indels，≤12 bp），与Sanger测序表现出优异的一致性。
•混合样本中观察到的变异等位基因频率比值与样本池组成所衍生的理论预期比例密切相关。
•为大规模WGS临床研究中正交变异验证提供了一种稳健、可扩展且经济的Sanger测序替代方案。
研究背景与问题：高通量分子遗传学测序按人类基因组覆盖范围可分为靶向测序、全外显子组测序（WES）与全基因组测序（WGS）。WGS覆盖编码区、非编码区、调控元件及结构变异，每人可检出约300万至500万个单核苷酸变异（SNVs）与小片段插入/缺失（indels）以及数千个结构变异，其中数十个可能具有潜在临床影响。这些变异可能需依据实验室政策、变异类型与质量指标进行正交确认，以确保所检变异真实存在，而非初始测序技术产生的假象。目前Sanger测序仍是确认NGS变异的金标准，但其仅适用于每样本少量变异的验证；若将Sanger测序用于86个样本全部约500万变异的确认，则极为费时、昂贵，并消耗大量DNA与试剂。同时，NGS流程为追求高灵敏度可能提高假阳性率，而现行专业指南已建议将确认限定于测序质量置信度较低或结果意外的变异，不再要求普遍确认。因此，亟需一种低DNA消耗、低成本、可在同一运行中并行验证大量变异的替代验证策略，这正是本研究开展的出发点。
研究内容与结论：研究人员利用捷克共和国Hornacko地区隔离人群（该人群遗传异质性低、帕金森综合征家族负担高且患病率增加）的86例独特DNA样本（其中21例来自Hornacko地区患者、44例来自Hornacko地区对照、21例为Hornacko地区以外多症状患者）开展研究，将经帕金森综合征相关表型（178个基因）、人群等位基因频率阈值（<0.001，gnomAD v3.1.2）、ACMG标准及ClinVar分类多重过滤后筛选出的307个变异（138个SNVs、126个小片段indels、43个≥12 bp的较大indels，含14个线粒体（mt）DNA变异）设计进290个独特扩增子，将归一化后的86份DNA池化为16个混合物，在Ion Torrent平台完成多重扩增子大规模平行测序（maMPS）。研究结论为：该方法可准确可靠地评估264个SNPs与小片段indels中的254个（超过96%），成功确认96%以上的可评估SNVs与小片段indels（≤12 bp），且混合样本中观察到的变异等位基因频率比值与理论预期比例密切相关，为WGS临床研究中变异验证提供了稳健、可扩展且经济的Sanger测序替代方案，对帕金森综合征及类似多基因疾病的遗传学研究与诊断具有重要应用价值。
关键技术方法：研究首先采用Illumina TruSeq DNA PCR-free建库试剂盒在NovaSeq 6000平台完成WGS初筛，以Varsome Clinical 13.1.6进行数据分析与变异注释；通过Ion AmpliSeq Custom Panels及AmpliSeq Designer工具设计引物，构建含290个独特扩增子的多重扩增子面板；86份DNA归一化至10 ng/μL后池化为16个混合物（前10个各含5份样本、其余6个各含6份样本）并稀释至0.7 ng/μL；经Ion Chef系统自动化建库与模板制备后，在Ion S5系统（530B芯片）上测序，数据由Torrent Suite Software 5.12.2处理并比对至GRCh38.p2.mask1参考基因组；选取19个代表性变异以Sanger测序进行正交方法学验证；采用R 4.5.2软件进行精确双侧二项检验，并以Benjamini–Hochberg程序控制错误发现率（FDR）。样本队列来源于捷克Hornacko隔离人群的86例DNA样本。
研究结果：maMPS评估总体结果显示，86个样本共评估693个变异观测，243个独特变异（347次观测）获确认，其中A1规则189个、A2规则24个、A1/A2规则19个、A3规则（同聚物区校正）29个；64个独特变异（346次观测）未获确认或不可评估，包括B规则（推定WGS假象）18个独特变异/252次观测、C1（未覆盖）5个、C2（同聚物区信号均一化）5个、C3（≥12 bp indel）36个独特变异/81次观测。与Sanger测序一致性结果显示，19个经Sanger测序的代表性变异中，maMPS按A1–A3规则确认的9个变异全部获Sanger确认；7个B规则变异中6个Sanger未检出、1个（CAT基因缺失）获Sanger确认；2个低深度C1变异经Sanger检出；1个C3变异两种方法均未确认，表明maMPS确认变异与Sanger具有良好一致性。预期与观察比值比较结果显示，278个（不含同聚物区与mtDNA变异）经分析并确认的变异中，33个（11.8%）经FDR校正后与理论预期比例显著偏离，246个（88.2%）与预期一致；显著偏离者中7个（占显著结果21.2%）高于预期、26个（78.6%）低于预期；indels在显著偏离变异中占比过高，提示测序方法对indel等位基因频率的估计准确性低于SNVs。特殊区域结果显示，全部14个mtDNA变异均在多样本中可靠确认，包括2个异质性比例分别为12%与13%者；5个位于同聚物区的变异无法可靠评估；43个>12 bp大片段indel中36个频繁无法检测，仅7个例外获确认。
讨论部分指出，NGS已成为临床遗传学检测的常见标准，各专业机构现建议实验室依据检测内容、测序平台、变异类型与质量指标制定各自的定向确认政策，不再要求普遍确认；Sanger测序存在等位基因脱扣导致假阴性、同聚物区及高相似度区域（如假基因）引物设计困难等局限。本研究方法的创新性不在于扩增子测序本身，而在于确认工作流设计：个体DNA按预设混合使86份样本缩减为16个测序文库，因每池组成已知，可计算每个WGS变异的预期等位基因比例，确认依赖三个互补特征——独立测序平台检测、出现于预期DNA混合物、与理论预期等位基因比例一致，该设计可在DNA投入有限的条件下，于同一测序运行中可靠、独立地验证低至1/12预期等位基因比例的杂合态变异。预期与观察比值的显著偏离宜解读为技术变异性的指示因素，可能源于加样与定量变异、覆盖不均与位点脱落、变异检出的系统偏倚，以及相对二项模型的过度离散；未来研究可采用beta-binomial模型、技术重复或位点特异性归一化以更好区分真实混合偏差与随机测序噪声。
研究结论：在排除该法认为不可靠或不可评估的≥12 bp indel变异后，研究可准确可靠地评估264个SNPs与小片段indels中的254个（超过96%）。该方法表现出稳健性能，若引入更长读长，还可扩展用于验证更大规模的基因组结构变异，从而提高其在帕金森综合征及类似多基因疾病遗传学研究与诊断中的适用性。
