[发明专利]一种全基因组测序的数据质量评估方法及其应用在审
申请号: | 202011593896.0 | 申请日: | 2020-12-29 |
公开(公告)号: | CN112599189A | 公开(公告)日: | 2021-04-02 |
发明(设计)人: | 张静波;李小雨;王伟伟;伍启熹;王建伟;刘倩;唐宇 | 申请(专利权)人: | 北京优迅医学检验实验室有限公司 |
主分类号: | G16B20/10 | 分类号: | G16B20/10;G16B30/10;G16B40/00 |
代理公司: | 北京路浩知识产权代理有限公司 11002 | 代理人: | 商秀玲 |
地址: | 100195 北京市海淀区*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 基因组 数据 质量 评估 方法 及其 应用 | ||
1.一种全基因组测序的数据质量评估方法,其特征在于,包括:
获得待测样本的全基因组测序结果,将所述全基因组测序结果比对到人类参考基因组,对比对结果数据进行预处理后计算相应的观察值;
进行观察值趋势拟合后,通过纵向波动指标和横向波动指标进行质量评估;
所述纵向波动指标为基于观察值的变异系数反应观察值在趋势拟合线上下波动幅度的标准;
所述横向波动指标为基于小片段信号占目标区域的比例反应观察值碎片化分布情况的标准。
2.根据权利要求1所述的数据质量评估方法,其特征在于,所述纵向波动指标具体为:
对于总窗口数为N的全基因组测序结果中任一窗口的残差x,纵向波动指标VI计算方法如下:
和/或,
所述横向波动指标具体为:
对于整体趋势拟合后,异常区域占全部检测区域的比例Rglobal,细节趋势拟合后,异常区域占全部检测区域的比例Rdetail,横向波动指标HI的计算方法如下:
HI=Rdetail-Rglobal。
3.根据权利要求1所述的数据质量评估方法,其特征在于,所述全基因组测序结果为基于二代测序平台的低深度测序结果,或二代测序平台基于深度进行拷贝数变异检测的全基因组测序结果。
4.根据权利要求3所述的数据质量评估方法,其特征在于,所述全基因组测序结果为NIPT、NIPTPlus、CNVSeq或PGS中的一种或多种。
5.根据权利要求1所述的数据质量评估方法,其特征在于,所述预处理包括如下步骤:
将所述比对结果数据划分为多个窗口,对各个窗口中的数据进行数据偏好校正和归一化处理。
6.根据权利要求5所述的数据质量评估方法,其特征在于,所述划分为多个窗口为通过等长窗口划分、等数据量窗口划分、单碱基或单独靶向目标区域进行窗口划分;和/或,
所述数据偏好校正为GC校正、Mappability校正或PCA校正中的一种或多种。
7.根据权利要求5或6所述的数据质量评估方法,其特征在于,在所述划分为多个窗口后,还包括:使用重叠窗口进行序列数平滑处理。
8.根据权利要求1所述的数据质量评估方法,其特征在于,所述观察值趋势拟合采用隐马尔可夫算法、环状二元分割算法或Fused Lasso算法中的一种或多种算法。
9.一种测序方法,其特征在于,使用权利要求1-7任一项所述方法进行质量控制。
10.权利要求1-8任一项所述数据质量评估方法在排除低质量测序样品中的应用;所述应用优选为:
设置纵向波动指标和横向波动指标的阈值,通过权利要求1-8任一项所述数据质量评估方法进行待测样品的测序结果质量评估后,将纵向波动指标和横向波动指标高于阈值的样品作为低质量测序样品排除。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京优迅医学检验实验室有限公司,未经北京优迅医学检验实验室有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202011593896.0/1.html,转载请声明来源钻瓜专利网。
- 数据显示系统、数据中继设备、数据中继方法、数据系统、接收设备和数据读取方法
- 数据记录方法、数据记录装置、数据记录媒体、数据重播方法和数据重播装置
- 数据发送方法、数据发送系统、数据发送装置以及数据结构
- 数据显示系统、数据中继设备、数据中继方法及数据系统
- 数据嵌入装置、数据嵌入方法、数据提取装置及数据提取方法
- 数据管理装置、数据编辑装置、数据阅览装置、数据管理方法、数据编辑方法以及数据阅览方法
- 数据发送和数据接收设备、数据发送和数据接收方法
- 数据发送装置、数据接收装置、数据收发系统、数据发送方法、数据接收方法和数据收发方法
- 数据发送方法、数据再现方法、数据发送装置及数据再现装置
- 数据发送方法、数据再现方法、数据发送装置及数据再现装置