种系结构变异检测中的隐性偏差

Hidden biases in germline structural variant detection

作者信息Michael M Khayat, Sayed Mohammad Ebrahim Sahraeian, Samantha Zarate, Andrew Carroll, Huixiao Hong, Bohu Pan, Leming Shi, Richard A Gibbs, Marghoob Mohiyuddin, Yuanting Zheng, Fritz J Sedlazeck
PMID34930391
发布时间2021-12
DOI10.1186/s13059-021-02558-x

实验完整度

中

包含多个样本、重复和中心的测序实验,以及长读长验证和基因分型,但主要是分析性研究,缺乏功能性验证。

主要模型

LCL5-8淋巴母细胞系 中国家庭四重奏样本

重点核对

使用四种比对方法(BWA-MEM、Isaac、Stampy、Bowtie2)进行比对 测序中心1-3的测序覆盖度分别为24.37x、26.33x、36.33x,插入片段大小分别为321.55bp、318.53bp、399.12bp 使用MetaSV和Parliament2进行SV调用 使用PacBio长读长数据进行验证

摘要

背景:基因组结构变异(SV)是许多生物体中基因型和表型变化的重要决定因素。然而,从下一代测序数据中检测SV仍然具有挑战性。结果:在本研究中,一个中国家庭四重奏的DNA在三个不同的测序中心以三次重复进行测序。利用不同的分析流程生成了总共288个衍生数据集,并进行比较以识别分析变异性的来源。比对方法对变异性的贡献最大,其次是测序中心和重复。有趣的是,仅由一个中心或重复支持的SV通常是真阳性,分别有47.02%和45.44%与长读长SV调用集重叠。这与中心和重复中SV调用的总体较高假阴性率(15.72%)一致。最后,我们观察到SV调用变异性在基因分型方法中仍然存在,表明底层测序和制备方法的影响。结论:本研究首次提供了下一代测序中SV鉴定变异性来源的详细见解,并强调了在大队列中SV调用所面临的剩余挑战。我们进一步就如何减少SV调用变异性以及比对方法的选择提出了建议。补充信息:在线版本包含补充材料,可在10.1186/s13059-021-02558-x获取。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
下一代测序中结构变异检测的变异性来源是什么?
核心机制
比对方法的选择是变异性的最大来源,测序中心和重复的贡献次之;仅单中心或单重复支持的SV中约一半是真阳性,表明存在较高的假阴性率。
主要证据
基于中国家庭四重奏的288个SV调用集,与PacBio长读长SV调用集比较,以及SVTyper基因分型验证。
研究意义
为降低SV检测变异性提供建议,强调多种比对方法联合使用可提高灵敏度和精确度,促进SV调用的常规应用。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

样本制备与测序

获得多个测序中心和重复的高质量测序数据,用于评估变异来源。

从中国家庭四重奏的淋巴母细胞系提取DNA,在三个测序中心以三次重复制备文库并进行Illumina X10测序。

2

短读长SV调用

使用多种比对方法和SV调用器生成SV调用集,用于评估变异性。

使用BWA-MEM、Isaac、Stampy、Bowtie2四种比对方法映射到GRCh38参考基因组,经去重复和碱基质量校准后,使用MetaSV和Parliament2进行SV调用,共生成288个SV调用集。

3

长读长SV验证

建立高质量SV调用集作为参考,以评估短读长SV调用的准确性。

对LCL5-8样本进行PacBio测序,使用minimap2比对和Sniffles进行SV调用,与短读长SV调用集比较。

4

变异来源分层分析

独立评估不同因素(比对方法、测序中心、重复、去重复/碱基质量校准)对SV调用变异性的贡献。

使用SURVIVOR对SV调用集进行合并和分层,分别比较不同因素引起的SV差异,并评估其与PacBio调用集的重叠。

5

基因分型验证

独立验证SV检测的变异性是否在基因分型层面也存在。

使用SVTyper对中心特异性和比对方法特异性的SV进行基因分型,评估各样本支持情况。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
RPMI 1640培养基Gibco31870-082
胎牛血清Gibco10091-148
血液和细胞培养DNA Maxi试剂盒QIAGEN13362
Illumina TruSeq DNA nano文库制备试剂盒Illumina--
Covaris LE220超声破碎仪Covaris--
Qubit 3.0荧光计ThermoFisher Scientific--
Quant-iT dsDNA HS检测试剂盒ThermoFisher ScientificQ32854
Agilent 2100生物分析仪Agilent--
SMRTbell Express模板制备试剂盒2.0Pacific Biosciences--
g-TUBEs剪切管Covaris--
BluePippin切胶仪Sage Science--
AMPure PB磁珠Pacific Biosciences--
Illumina X10测序平台Illumina--
PacBio Sequel测序平台Pacific Biosciences--

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
样本制备与测序
细胞培养条件(培养基、血清、传代次数)、DNA提取试剂盒、文库制备试剂盒、片段化参数、测序平台和读长
阅读提示:Methods: Ethics, consent, and sample ascertainment; Next-generation sequencing
短读长SV调用
比对工具及版本、参考基因组、去重复和碱基质量校准步骤、SV调用工具及版本、合并参数
阅读提示:Methods: Short-read SV identification; Mapper comparisons
长读长SV验证
PacBio文库制备试剂盒和流程、测序平台、比对工具(minimap2)版本、SV调用工具(Sniffles)版本及参数
阅读提示:Methods: Generation and comparisons to long-read-based SV calls
变异来源分层分析
SURVIVOR合并参数、分层策略(如何固定其他因素)、比较标准(重叠距离、类型、大小)
阅读提示:Methods: Mapper comparisons; Center comparisons; Replicate comparisons; Dedup/recalibration comparisons
基因分型验证
SVTyper版本和参数、输入VCF和bam文件、过滤标准(至少一个reads支持)、SURVIVOR合并参数
阅读提示:Methods: SV genotyping