[发明专利]基于全外显子测序的非靶向区域基因型填充方法

说明书

本发明公开了基于全外显子测序的非靶向区域基因型填充方法、系统、设备和计算机可读存储介质，方法其包括：获取目标队列的全外显子测序数据、参考全基因组测序数据集；对参考全基因组测序数据集中的位点进行过滤，输出参考全基因组测序数据集的SNP位点信息；基于SNP位点信息和全外显子测序数据，计算每个样本SNP位点的GL，并得到每个样本的GL；合并每个样本的GL，得到全部样本的GL；对参考全基因组测序数据集进行片段化分割处理得到基因组片段信息；基于全部样本的GL、SNP位点信息、基因组片段信息，利用机器学习算法对目标队列中单个样本非靶向区域的基因型进行估算，得到单个样本非靶向区域的基因型估算结果。

技术领域

本发明涉及基因预测的技术领域，更具体地，涉及一种基于全外显子测序的非靶向区域基因型填充方法及其系统。

背景技术

随着高通量测序技术(Next generation sequencing,NGS)的快速发展，越来越多的GWAS试图采用测序技术进行基因型鉴定。全基因组测序技术(Whole genomesequencing,WGS)能够有效地覆盖大多数的遗传突变位点，但缺点在于价格昂贵，不能应用到大规模队列样本上。另外，如果为了成本而降低测序深度，则无法准确鉴定低频变异的基因型状态，该低频变异往往位于蛋白编码区域，具有重要的生物学功能。

全外显子测序技术(Whole Exome sequencing,WES)已经广泛用于研究蛋白质编码变异在遗传疾病中的作用。然而，基因组中超过98％的序列为非编码序列，存在有很多能够显著影响基因表达的重要调控元件。大量研究表明非编码变异与疾病的形成和发展具有密切的关系，因此，如何准确测量全基因组范围内的基因变异信息成为制约科研人员探究疾病机理的关键。WES会针对基因外显子区域设计探针，靶向捕获该区域后进行测序，在数据通量一定的情况下可以显著提高外显子区域的测序深度；然而，该技术的缺点是无法监测位于基因组中绝大部分区域的非编码变异，而这些非编码变异往往具有更显著的遗传效力。研究表明，在WES测序数据不仅包含了设计探针部位的序列，还有一部分位于探针区域外的序列被捕获，进而被测序仪检测到。相较于探针区域，非探针区域/非靶向区域(Off-target sequence)的测序深度往往偏低，无法利用传统的基因组鉴定算法进行准确识别。随着大型参考基因组的逐步公开，新的基因型鉴定算法可以利用参考基因型panel对低覆盖度的区域进行基因型推断。

然而，大型参考基因组panel具有很强的人群特异性，目前主流的数据库多数以欧洲或者混合人群为主，由于人种差异导致单体型的不同，上述panel无法契合不同国家/种族/地域的特定人群的低覆盖度基因型鉴定需求。

发明内容

本发明旨在至少解决现有技术中存在的技术问题之一。为此，本发明提供一种基于全外显子测序的非靶向区域基因型填充方法及其系统；本发明方法基于全外显子测序技术和参考全基因组测序数据集，计算目标队列中非靶向区域的基因型可能性，得到非靶向区域的基因型估算结果，整合低深度测序数据和大型参考基因组panel，高效准确鉴定WES非靶向区域遗传变异位点的基因型；从深层次挖掘隐含在测序数据背后的生命规律，解决相关的生命科学问题。

本申请第一方面公开一种基于全外显子测序的非靶向区域基因型填充方法，包括：

获取目标队列的全外显子测序数据、参考全基因组测序数据集；

对所述参考全基因组测序数据集中的位点进行过滤，输出参考全基因组测序数据集的SNP位点信息；所述SNP位点信息包括染色体号、基因组坐标及等位基因、已去除非SNP的位点的人群基因型信息；

基于所述参考全基因组测序数据集的SNP位点信息和所述目标队列的全外显子测序数据，计算目标队列中每个样本各个SNP位点的基因型可能性结果，并得到每个样本中的基因型可能性结果；

合并所述每个样本中的基因型可能性结果，得到目标队列全部样本中基因型可能性结果；