CN110211630A - 致病性单亲二倍体的筛查装置及存储介质和处理器 - Google Patents

致病性单亲二倍体的筛查装置及存储介质和处理器 Download PDF

Info

Publication number
CN110211630A
CN110211630A CN201910491767.1A CN201910491767A CN110211630A CN 110211630 A CN110211630 A CN 110211630A CN 201910491767 A CN201910491767 A CN 201910491767A CN 110211630 A CN110211630 A CN 110211630A
Authority
CN
China
Prior art keywords
loh
screening
mutation
site
upd
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910491767.1A
Other languages
English (en)
Other versions
CN110211630B (zh
Inventor
刘晶星
赵薇薇
陈白雪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GUANGZHOU KINGMED DIAGNOSTICS GROUP Co.,Ltd.
Original Assignee
Guangzhou Jinyu Medical Inspection Group Ltd By Share Ltd
Guangzhou Kingmed Diagnostics Central Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangzhou Jinyu Medical Inspection Group Ltd By Share Ltd, Guangzhou Kingmed Diagnostics Central Co Ltd filed Critical Guangzhou Jinyu Medical Inspection Group Ltd By Share Ltd
Priority to CN201910491767.1A priority Critical patent/CN110211630B/zh
Publication of CN110211630A publication Critical patent/CN110211630A/zh
Application granted granted Critical
Publication of CN110211630B publication Critical patent/CN110211630B/zh
Priority to PCT/CN2020/094125 priority patent/WO2020244538A1/zh
Priority to US17/616,714 priority patent/US20220328131A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/10Ploidy or copy number detection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/20Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B20/00ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
    • G16B20/50Mutagenesis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H10/00ICT specially adapted for the handling or processing of patient-related medical or healthcare data
    • G16H10/40ICT specially adapted for the handling or processing of patient-related medical or healthcare data for data related to laboratory analysis, e.g. patient specimen analysis
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/70ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Genetics & Genomics (AREA)
  • Theoretical Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Evolutionary Biology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Data Mining & Analysis (AREA)
  • Biomedical Technology (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

本发明涉及一种致病性单亲二倍体的筛查装置及存储介质和处理器,属于基因检测技术领域。该筛查装置包括:数据获取模块:用于获取全外显子组测序数据;位点筛选模块:用于筛选得到预定条件的突变;LOH判断模块:用于根据上述得到的突变情况进行LOH判断;UPD判定模块:用于根据LOH判断UPD,如发生LOH的染色体数超过2条,判定为近亲关系;如发生LOH的区段为单拷贝,判定为片段缺失;其余发生LOH的区段判定为UPD。该筛查装置通过筛选出特定的突变位点,进行LOH判断并最终得到UPD的判定结果。可基于全外显子组测序的数据,在检查常规致病突变的同时提示存在致病性UPD的风险,无需额外实验和人力成本。

Description

致病性单亲二倍体的筛查装置及存储介质和处理器
技术领域
本发明涉及基因检测技术领域,特别是涉及一种致病性单亲二倍体的筛查装置及存储介质和处理器。
背景技术
基因组印记(Genomic imprinting),又称遗传印记,是通过生化途径,在一个基因或基因组域上标记其双亲来源信息的遗传学过程。这类基因称作印记基因,这类基因表达与否取决于它们所在染色体的来源(父系或母系)以及在其来源的染色体上该基因是否发生沉默(沉默机制主要为甲基化)。有些印记基因只从母源染色体上表达,而有些则只从父源染色体上表达。
正常二倍体中,一对同源染色体分别来源于父本和母本,单亲二倍体(UniParental Disomy简称UPD)是指一对同源染色体(或染色体的部分区段)来源于同一个亲本,如果这些区段包含印记基因,则会导致基因表达紊乱。
目前诊断UPD的方法主要有检测甲基化水平和SNP芯片法。其中,检测甲基化水平方法即检测一对同源染色体的同一段区间的甲基化水平是否一致,但是,甲基化的方法只能处理染色体局部的小片段,而且针对不同区域需要设计不同的实验,效率低速度慢,并不适用于全基因组范围的筛查。而使用SNP芯片检测是否存在连续的大片段纯合位点,又存在SNP芯片的方法成本较高的缺陷,且其靶标探针为多态性位点,无法同时检测其他的致病微小突变(点突变、微小插入缺失等)。
发明内容
基于此,有必要针对上述问题,提供一种致病性单亲二倍体的筛查装置,采用该筛查装置进行筛查,可基于全外显子组测序的数据,在检查常规致病突变的同时提示存在致病性UPD的风险,无需额外实验和人力成本。
一种致病性单亲二倍体的筛查装置,包括:
数据获取模块:用于获取全外显子组测序数据;
位点筛选模块:用于筛选得到预定条件的突变;
LOH判断模块:用于根据上述得到的突变情况进行LOH判断,如连续纯合位点数与其覆盖范围乘积大于预设值,则判定该区间为LOH;
UPD判定模块:用于根据LOH判断UPD,如发生LOH的染色体数超过2条,判定为近亲关系;如发生LOH的区段为单拷贝,判定为片段缺失;其余发生LOH的区段判定为UPD。
全外显子测序是目前检测基因缺陷疾病最普遍的方法,可以检测致病性点突变、微小插入缺失、拷贝数变异等,是大多数此类患者的首选项目。本发明人在长期实践经验的基础上,通过创造性的实验设计和反复摸索,在全外显子测序的基础上增加一项致病性UPD的筛查,在不增加任何成本的前提下提高了诊断阳性率。
考虑到UPD来源于一方亲本的同一条染色体的两个拷贝,因此表现为该区域所有碱基均为纯合,即杂合性丢失(loss of heterozygosity,LOH),而造成LOH的原因主要有三种:片段缺失、UPD、近亲结婚。这三种情况造成的LOH在片段大小、分布、以及临床表现上各有不同,因此可以通过检测LOH的方法来推断UPD的存在,在此理论基础之上,本发明人通过筛选出特定的突变位点,进行LOH判断并最终得到UPD的判定结果。
对于近亲关系的判定,由于UPD发生的偶然性,同时在多条染色体上发生的概率非常小,据此可用于区分近亲结婚的情况,即发生LOH的染色体数超过2条(即大于2条)的判定为近亲结婚;对于片段缺失的判定,按照常规方式判定即可,如可结合全外显子测序拷贝数变异(CNV)分析结果进行,即参考LOH区段的测序数据覆盖深度与同批次其他样本的对比,若CNV分析提示该LOH区段为单拷贝,则判定为片段缺失;特别的,很大段的缺失一般是致死性的,如果LOH区段达到整条染色体的一半以上甚至为整条染色体,若样本来源非胚胎,基本可以排除片段缺失。
在其中一个实施例中,所述预定条件的突变通过以下方法筛选得到:
筛选高质量位点:在全外显子组测序数据中筛选高质量突变位点;
除Y染色体突变:去除上述突变位点中位于Y染色体上的突变;
筛选点突变:筛选上述除Y染色体突变步骤得到突变中的点突变;
等位基因频率筛选:筛选上述点突变在群体数据库中各人种中人群等位基因频率均低于0.7的点突变位点;
突变频率筛选:去除上述点突变位点中杂合位点突变频率高于70%的位点,并且去除纯合位点突变频率低于85%的位点,即得预定条件的突变。
以上述突变进行分析,可以排除假阳性突变、体细胞突变、以及人群高频突变对LOH判定的影响,具有判定准确的优点。例如一段较大的LOH,中间掺杂有几个假阳性或体细胞的杂合突变,将会被拆成几个小的LOH,如果每个小的LOH达不到预定长度阈值(如3M),则无法识别,导致判定不准确。
上述群体数据库包括千人基因组、ESP6500、ExAC、gnomAD等,对于各人种的分类包括东亚、南亚、非洲/非裔美洲、美洲、芬兰、非芬兰欧洲等。
在其中一个实施例中,所述筛选高质量位点步骤中,所述高质量突变位点指:通过GATK-VQSR质控、总覆盖>40X、且突变频率>30%。
上述GATK-VQSR质控指在GATK软件中variant quality score recalibration得到的结果为PASS;“总覆盖>40X”是指该位点上覆盖的有效reads数超过40条。上述“突变频率>30%”是指该位点上含突变碱基的reads占所有reads的比例。
在其中一个实施例中,所述等位基因频率筛选步骤和突变频率筛选步骤之间,还包括排除假阳性位点步骤,所述排除假阳性位点步骤为:根据Hardy-Weinberg平衡在待评估区域人群频率库中排除存在的假阳性位点。所述待评估区域人群频率库指待评估个体所在地域频率库,即根据地域性特点,排除假阳性位点。
在其中一个实施例中,所述位点筛选模块中还包括质控单元,所述质控单元用于检测筛选得到的突变数量,如所述突变数量≥1万,则所述质控单元提示通过;如所述突变数量<1万,则所述质控单元提示不通过。如筛选得到的位点数过少,太少的话连续纯合位点数不够导致统计上没有显著性。
在其中一个实施例中,所述LOH判断模块中,所述连续纯合位点数≥20,所述覆盖范围≥3Mbp。
在其中一个实施例中,所述LOH判断模块中,如连续纯合位点数与其覆盖范围乘积大于200Mbp,则判定该区间为LOH。例如:连续5Mbp范围都是Hom(纯合)位点,其中Hom位点数为60,60×5>200,判断该区间是LOH。
上述200Mbp的预设值,是本发明人通过反复试验、不断测试后得到的阈值,具有判断准确,误判率低的优点。
在其中一个实施例中,所述UPD判定模块中,还包括致病风险判断单元,所述致病风险判断单元中,将判定为UPD的LOH区段进行印记基因比对,如所述LOH区段未覆盖印记基因或对应条带,提示良性UPD;如所述LOH区段覆盖印记基因或对应条带,提示致病UPD风险。
本发明还公开了一种存储介质,所述存储介质包括存储的程序,所述程序实现上述模块的功能。
本发明还公开了一种处理器,所述处理器用于运行程序,所述程序实现上述模块的功能。
与现有技术相比,本发明具有以下有益效果:
本发明的一种致病性单亲二倍体的筛查装置,通过数据获取模块、位点筛选模块、LOH判断模块和UPD判定模块的依次分析判断,通过筛选出特定的突变位点,进行LOH判断并最终得到UPD的判定结果。可基于全外显子组测序的数据,在检查常规致病突变的同时提示存在致病性UPD的风险,无需额外实验和人力成本。
附图说明
图1为实施例1中LOH在染色体上的分布示意图;
图2为图1中5号和7号染色体上LOH分布放大示意图;
图3为图1中14号、16号和19号染色体上LOH分布放大示意图;
图4为实施例2中LOH在染色体上的分布示意图;
图5为图4中15号染色体上LOH分布放大示意图;
图6为实施例3中LOH在染色体上的分布示意图;
图7为图6中5号染色体上12.57(M)LOH分布放大示意图;
其中:图1,4,6中,横坐标为各染色体号,图下半部分为连续的纯合片段占整个染色体长度的比例,图上半部分为各染色体上突变位点的分布情况;
图2,3,5,7的LOH放大示意图中,中间黑色线段为全外显子测序覆盖范围(exomebed),左边菱形点为检测到的杂合(Het)突变,右边五角星点为检测到的纯合(Hom)突变,右侧点状虚线段为印记区段(imprint location),其上四角点为印记基因范围(imprintgene)。
具体实施方式
为了便于理解本发明,下面将参照相关附图对本发明进行更全面的描述。附图中给出了本发明的较佳实施例。但是,本发明可以以许多不同的形式来实现,并不限于本文所描述的实施例。相反地,提供这些实施例的目的是使对本发明的公开内容的理解更加透彻全面。
除非另有定义,本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的,不是旨在于限制本发明。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
实施例1
一种致病性单亲二倍体的筛查,采用以下装置进行,该筛查装置包括:
数据获取模块:用于获取全外显子组测序数据;
位点筛选模块:用于筛选得到预定条件的突变;
LOH判断模块:用于根据上述得到的突变情况进行LOH判断,如连续纯合位点数与其覆盖范围乘积大于预设值,则判定该区间为LOH;
UPD判定模块:用于根据LOH判断UPD,如发生LOH的染色体数超过2条,判定为近亲关系;如发生LOH的区段为单拷贝,判定为片段缺失;其余发生LOH的区段判定为UPD。
以上述筛查装置按照下述方法运行程序,包括以下步骤:
1、数据获取
获取某例样本的全外显子组测序数据,其中突变为59312个。
2、位点筛选
2.1筛选高质量位点:
在全外显子组测序数据中筛选高质量突变位点,具体为筛选通过GATK-VQSR质控、总覆盖>40X、且突变频率>30%的突变位点,在本实施例中为45260个。
2.2除Y染色体突变:
去除上述突变位点中位于Y染色体上的突变,得到45256个突变。
2.3筛选点突变:
筛选上述除Y染色体突变步骤得到突变中的点突变,得到41273个突变。
2.4等位基因频率筛选:
筛选上述点突变在群体数据库(千人基因组、ESP6500、ExAC、gnomAD)中各人种(东亚、南亚、非洲/非裔美洲、美洲、芬兰、非芬兰欧洲)中人群等位基因频率均低于0.7的点突变位点,得到22231个突变。
2.5排除假阳性位点:
根据Hardy-Weinberg平衡在待评估区域人群频率库中排除存在的假阳性位点,得到21705个突变。
2.6突变频率筛选:
去除上述点突变位点中杂合位点突变频率高于70%的位点,并且去除纯合位点突变频率低于85%的位点,即得预定条件的突变,21644个。
3、LOH判断
上述得到的位点中,如连续纯合位点数与其覆盖范围乘积大于200Mbp,则判定该区间为LOH,其中连续纯合位点数≥20,覆盖范围≥3Mbp。
按照该标准判断,本实施例的样本中检测到5段LOH,如下表所示。
表1.LOH区间
从上述结果可以看出,上述5段LOH分别位于5条染色体上。图1为这5段LOH在染色体上的分布情况,图中椭圆圈部分表示LOH区间,图2和3分别为图1中5,7,14,16,19号染色体上的LOH放大示意图。
4、UPD判定
由于上述5段LOH分别位于5条染色体上,据此判断为近亲结婚,排除UPD致病的可能。
该病例经后期证实为近亲结婚后代。
实施例2
一种致病性单亲二倍体的筛查,以某例样本进行,采用实施例1的装置,其中:
1、数据获取
参照实施例1。
2、位点筛选
参照实施例1,得到符合预定条件的突变22210个。
3、LOH判断
上述得到的位点中,如连续纯合位点数与其覆盖范围乘积大于200,则判定该区间为LOH,其中连续纯合位点数≥20,覆盖范围≥3Mbp。
按照该标准判断,本实施例的样本中检测到1段LOH,如下表所示。
表2.LOH区间
从上述结果可以看出,上述LOH位于15号染色体上,长度为12.28M。图4为这段12.28M的LOH在染色体上的分布情况,图中椭圆圈部分表示LOH区间,图5为图4中15号染色体上的LOH放大示意图。
4、UPD判定
4.1基础判定
由于该段LOH不符合近亲关系判定规则,不符合片段缺失规则,判定为UPD。
4.2致病风险判断
上述12.28M的LOH覆盖的印记基因与Prader-Willi综合征相关。
该病例经后期证实具有Prader-Willi综合征病征。
实施例3
一种致病性单亲二倍体的筛查,以某例样本进行,采用实施例1的装置,其中:
1、数据获取
参照实施例1。
2、位点筛选
参照实施例1,得到符合预定条件的突变22947个。
3、LOH判断
上述得到的位点中,如连续纯合位点数与其覆盖范围乘积大于200,则判定该区间为LOH,其中连续纯合位点数≥20,覆盖范围≥3Mbp。
按照该标准判断,本实施例的样本中检测到两段LOH,如下表所示。
表3.LOH区间
从上述结果可以看出,上述LOH位于5号染色体上,长度分别为93.6M和12.36M。图6为其在染色体上的分布情况,图中椭圆圈部分表示LOH区间,图7分别为图6中12.57M的LOH放大示意图。
注:该例样本同时做了CMA基因芯片检测(芯片型号为CytoScan HD),检测结果为chr5:2667631-99572420和chr5:166974594-180520810两段LOH,与本方法得到结果基本一致。
4、UPD判定
4.1基础判定
由于该段LOH不符合近亲关系判定规则,不符合片段缺失规则,判定为UPD。
4.2致病风险判断
上述长度为93.6M的LOH区间覆盖了印记基因ERAP2,RNU5D-1,但目前与之相关的研究较少,不能明确为致病原因,但可提示相关风险。
实施例4
一种致病性单亲二倍体的筛查,以送检本单位的12444例全外显子测序数据为基础,中筛查致病性UPD,按照实施例1的方法进行筛查,检测到LOH为1018例,排除近亲结婚之后余800例,进行分析后发现覆盖了印记基因的为142例,其中部分案例经回访后证实与筛查结果相符度大于95%。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种致病性单亲二倍体的筛查装置,其特征在于,包括:
数据获取模块:用于获取全外显子组测序数据;
位点筛选模块:用于筛选得到预定条件的突变;
LOH判断模块:用于根据上述得到的突变情况进行LOH判断,如连续纯合位点数与其覆盖范围乘积大于预设值,则判定该区间为LOH;
UPD判定模块:用于根据LOH判断UPD,如发生LOH的染色体数超过2条,判定为近亲关系;如发生LOH的区段为单拷贝,判定为片段缺失;其余发生LOH的区段判定为UPD。
2.根据权利要求1所述的致病性单亲二倍体的筛查装置,其特征在于,所述预定条件的突变通过以下方法筛选得到:
筛选高质量位点:在全外显子组测序数据中筛选高质量突变位点;
除Y染色体突变:去除上述突变位点中位于Y染色体上的突变;
筛选点突变:筛选上述除Y染色体突变步骤得到突变中的点突变;
等位基因频率筛选:筛选上述点突变在群体数据库中各人种中人群等位基因频率均低于0.7的点突变位点;
突变频率筛选:去除上述点突变位点中杂合位点突变频率高于70%的位点,并且去除纯合位点突变频率低于85%的位点,即得预定条件的突变。
3.根据权利要求2所述的致病性单亲二倍体的筛查装置,其特征在于,所述筛选高质量位点步骤中,所述高质量突变位点指:通过GATK-VQSR质控、总覆盖>40X、且突变频率>30%。
4.根据权利要求2所述的致病性单亲二倍体的筛查装置,其特征在于,所述等位基因频率筛选步骤和突变频率筛选步骤之间,还包括排除假阳性位点步骤,所述排除假阳性位点步骤为:根据Hardy-Weinberg平衡在待评估区域人群频率库中排除存在的假阳性位点。
5.根据权利要求1所述的致病性单亲二倍体的筛查装置,其特征在于,所述位点筛选模块中还包括质控单元,所述质控单元用于检测筛选得到的突变数量,如所述突变数量≥1万,则所述质控单元提示通过;如所述突变数量<1万,则所述质控单元提示不通过。
6.根据权利要求1所述的致病性单亲二倍体的筛查装置,其特征在于,所述LOH判断模块中,所述连续纯合位点数≥20,所述覆盖范围≥3Mbp。
7.根据权利要求6所述的致病性单亲二倍体的筛查装置,其特征在于,所述LOH判断模块中,如连续纯合位点数与其覆盖范围乘积大于200Mbp,则判定该区间为LOH。
8.根据权利要求1所述的致病性单亲二倍体的筛查装置,其特征在于,所述UPD判定模块中,还包括致病风险判断单元,所述致病风险判断单元中,将判定为UPD的LOH区段进行印记基因比对,如所述LOH区段未覆盖印记基因或对应条带,提示良性UPD;如所述LOH区段覆盖印记基因或对应条带,提示致病UPD风险。
9.一种存储介质,其特征在于,所述存储介质包括存储的程序,所述程序实现权利要求1-8任一项所述模块的功能。
10.一种处理器,其特征在于,所述处理器用于运行程序,所述程序实现权利要求1-8任一项所述模块的功能。
CN201910491767.1A 2019-06-06 2019-06-06 致病性单亲二倍体的筛查装置及存储介质和处理器 Active CN110211630B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201910491767.1A CN110211630B (zh) 2019-06-06 2019-06-06 致病性单亲二倍体的筛查装置及存储介质和处理器
PCT/CN2020/094125 WO2020244538A1 (zh) 2019-06-06 2020-06-03 致病性单亲二倍体的筛查方法及应用
US17/616,714 US20220328131A1 (en) 2019-06-06 2020-06-03 Method for screening pathogenic uniparental disomy and use thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910491767.1A CN110211630B (zh) 2019-06-06 2019-06-06 致病性单亲二倍体的筛查装置及存储介质和处理器

Publications (2)

Publication Number Publication Date
CN110211630A true CN110211630A (zh) 2019-09-06
CN110211630B CN110211630B (zh) 2020-03-20

Family

ID=67791367

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910491767.1A Active CN110211630B (zh) 2019-06-06 2019-06-06 致病性单亲二倍体的筛查装置及存储介质和处理器

Country Status (3)

Country Link
US (1) US20220328131A1 (zh)
CN (1) CN110211630B (zh)
WO (1) WO2020244538A1 (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111863125A (zh) * 2020-08-04 2020-10-30 广州金域医学检验中心有限公司 基于NGS-trio的单亲二倍体检测方法及应用
WO2020244538A1 (zh) * 2019-06-06 2020-12-10 广州金域医学检验中心有限公司 致病性单亲二倍体的筛查方法及应用
CN112687336A (zh) * 2021-03-11 2021-04-20 北京贝瑞和康生物技术有限公司 确定upd类型的方法、计算设备和存储介质
CN113066529A (zh) * 2021-03-26 2021-07-02 四川大学华西医院 基于全外显子数据的近亲家系鉴定方法、装置及设备
WO2022027212A1 (zh) * 2020-08-04 2022-02-10 广州金域医学检验中心有限公司 基于NGS-trio的单亲二倍体检测方法及应用
CN115394357A (zh) * 2022-09-01 2022-11-25 杭州链康医学检验实验室有限公司 用于判断样本配对或污染的位点组合及其筛选方法和应用

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011083312A1 (en) * 2010-01-08 2011-07-14 Oxford Gene Technology (Operations) Ltd Combined cgh & allele-specific hybridisation method
US20110301854A1 (en) * 2010-06-08 2011-12-08 Curry Bo U Method of Determining Allele-Specific Copy Number of a SNP
US20130230221A1 (en) * 2011-11-29 2013-09-05 Affymetrix, Inc. Analysis of Data Obtained from Microarrays
CN106021984A (zh) * 2016-05-13 2016-10-12 万康源(天津)基因科技有限公司 一种全外显子组测序数据分析系统
KR101721480B1 (ko) * 2016-06-02 2017-03-30 주식회사 랩 지노믹스 염색체 이상 검사 방법 및 시스템

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102004005497B4 (de) * 2004-01-30 2007-01-11 Eberhard-Karls-Universität Tübingen Universitätsklinikum Diagnose von uniparentaler Disomie anhand von Single-Nukleotid-Polymorphismen
CN110211630B (zh) * 2019-06-06 2020-03-20 广州金域医学检验中心有限公司 致病性单亲二倍体的筛查装置及存储介质和处理器

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011083312A1 (en) * 2010-01-08 2011-07-14 Oxford Gene Technology (Operations) Ltd Combined cgh & allele-specific hybridisation method
US20110301854A1 (en) * 2010-06-08 2011-12-08 Curry Bo U Method of Determining Allele-Specific Copy Number of a SNP
US20130230221A1 (en) * 2011-11-29 2013-09-05 Affymetrix, Inc. Analysis of Data Obtained from Microarrays
CN106021984A (zh) * 2016-05-13 2016-10-12 万康源(天津)基因科技有限公司 一种全外显子组测序数据分析系统
KR101721480B1 (ko) * 2016-06-02 2017-03-30 주식회사 랩 지노믹스 염색체 이상 검사 방법 및 시스템

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CARL FRIEDRICH CLASSEN等: ""Dissecting the genotype in syndromic intellectual disability using whole exome sequencing in addition to genome-wide copy number analysis"", 《HUM GENET》 *
MANNY D. BACOLOD等: ""Emerging Paradigms in Cancer Genetics: Some Important Findings from High-Density Single Nucleotide Polymorphism Array Studies"", 《CANCER RES》 *
YAO-SHAN FAN等: ""Frequent detection of parental consanguinity in children with developmental disorders by a combined CGH and SNP microarray"", 《MOLECULAR CYTOGENETICS》 *
付斌: "《再生障碍性贫血临床医师诊疗手册》", 30 April 2018, 上海世界图书出版公司 *
贾静等: ""单亲二倍体染色体异常的研究进展"", 《国际生殖健康/计划生育杂志》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020244538A1 (zh) * 2019-06-06 2020-12-10 广州金域医学检验中心有限公司 致病性单亲二倍体的筛查方法及应用
CN111863125A (zh) * 2020-08-04 2020-10-30 广州金域医学检验中心有限公司 基于NGS-trio的单亲二倍体检测方法及应用
WO2022027212A1 (zh) * 2020-08-04 2022-02-10 广州金域医学检验中心有限公司 基于NGS-trio的单亲二倍体检测方法及应用
CN111863125B (zh) * 2020-08-04 2024-04-12 广州金域医学检验中心有限公司 基于NGS-trio的单亲二倍体检测方法及应用
CN112687336A (zh) * 2021-03-11 2021-04-20 北京贝瑞和康生物技术有限公司 确定upd类型的方法、计算设备和存储介质
CN112687336B (zh) * 2021-03-11 2021-06-22 北京贝瑞和康生物技术有限公司 确定upd类型的方法、计算设备和存储介质
CN113066529A (zh) * 2021-03-26 2021-07-02 四川大学华西医院 基于全外显子数据的近亲家系鉴定方法、装置及设备
CN113066529B (zh) * 2021-03-26 2023-08-18 四川大学华西医院 基于全外显子数据的近亲家系鉴定方法、装置及设备
CN115394357A (zh) * 2022-09-01 2022-11-25 杭州链康医学检验实验室有限公司 用于判断样本配对或污染的位点组合及其筛选方法和应用

Also Published As

Publication number Publication date
US20220328131A1 (en) 2022-10-13
WO2020244538A1 (zh) 2020-12-10
CN110211630B (zh) 2020-03-20

Similar Documents

Publication Publication Date Title
CN110211630A (zh) 致病性单亲二倍体的筛查装置及存储介质和处理器
Davies et al. HRDetect is a predictor of BRCA1 and BRCA2 deficiency based on mutational signatures
Sun et al. Linked-read sequencing of gametes allows efficient genome-wide analysis of meiotic recombination
Sun et al. SHOREmap v3. 0: fast and accurate identification of causal mutations from forward genetic screens
CN111863125B (zh) 基于NGS-trio的单亲二倍体检测方法及应用
Andersen et al. Genetic analyses of the human eye colours using a novel objective method for eye colour classification
Ramani et al. Ancestry prediction in Singapore population samples using the Illumina ForenSeq kit
Timmermans et al. Mimicry diversification in Papilio dardanus via a genomic inversion in the regulatory region of engrailed–invected
Hussin et al. Haplotype allelic classes for detecting ongoing positive selection
CN110648722B (zh) 新生儿遗传病患病风险评估的装置
CN106599613B (zh) 一种判定遗传性肿瘤变异位点分级的方法
Logue et al. Whole‐genome sequencing reveals absence of recent gene flow and separate demographic histories for Anopheles punctulatus mosquitoes in Papua New Guinea
CN109712671B (zh) 基于ctDNA的基因检测装置、存储介质及计算机系统
CN104232649A (zh) 基因突变体及其应用
US20230282307A1 (en) Method for detecting uniparental disomy based upon ngs-trio, and use thereof
JP4468773B2 (ja) 遺伝子情報の表示方法及び表示装置
CN114921536A (zh) 一种检测单亲二倍体和杂合性缺失的方法、装置、存储介质和设备
CN114743596A (zh) 基于二代测序数据的孟德尔遗传错误分析方法
CN112837746B (zh) 用于小麦外显子测序基因定位的探针设计方法及定位方法
CN114566213A (zh) 家系高通量测序数据的单亲二倍体分析方法及其系统
CN107545152A (zh) 一种基于Illumina数据找变异的方法
CN116209777A (zh) 基于无创产前基因检测数据的亲缘关系判定方法和装置
CN115579056B (zh) 一组用于评估精神分裂症分子分型的基因群及其诊断产品和应用
CN112735519B (zh) 一种定位偏分离性状的方法、装置及存储介质
CN112687336B (zh) 确定upd类型的方法、计算设备和存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190906

Assignee: SHANGHAI KINGMED MEDICAL DIAGNOSTICS INSTITUTE Co.,Ltd.

Assignor: GUANGZHOU KINGMED CENTER FOR CLINICAL LABORATORY|GUANGZHOU KINGMED DIAGNOSTICS GROUP Co.,Ltd.

Contract record no.: X2020980004717

Denomination of invention: Screening device, storage medium and processor of pathogenic haploid

Granted publication date: 20200320

License type: Common License

Record date: 20200805

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190906

Assignee: Zhengzhou Jinyu Clinical Laboratory Center Co.,Ltd.

Assignor: GUANGZHOU KINGMED DIAGNOSTICS GROUP Co.,Ltd.

Contract record no.: X2021980010019

Denomination of invention: Screening device, storage medium and processor for pathogenic uniparental diploid

Granted publication date: 20200320

License type: Common License

Record date: 20210928

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20211014

Address after: No. 10, helix 3 Road, International Biological Island, Guangzhou, Guangdong 510000

Patentee after: GUANGZHOU KINGMED DIAGNOSTICS GROUP Co.,Ltd.

Address before: 510335 3rd floor, 2429 Xingang East Road, Haizhu District, Guangzhou City, Guangdong Province

Patentee before: GUANGZHOU KINGMED CENTER FOR CLINICAL LABORATORY

Patentee before: GUANGZHOU KINGMED DIAGNOSTICS GROUP Co.,Ltd.

EC01 Cancellation of recordation of patent licensing contract
EC01 Cancellation of recordation of patent licensing contract

Assignee: Zhengzhou Jinyu Clinical Laboratory Center Co.,Ltd.

Assignor: GUANGZHOU KINGMED DIAGNOSTICS GROUP Co.,Ltd.

Contract record no.: X2021980010019

Date of cancellation: 20220922

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20190906

Assignee: Zhengzhou Jinyu Clinical Laboratory Center Co.,Ltd.

Assignor: GUANGZHOU KINGMED DIAGNOSTICS GROUP Co.,Ltd.

Contract record no.: X2022980016522

Denomination of invention: Screening device, storage medium and processor for pathogenic parthenodiploid

Granted publication date: 20200320

License type: Common License

Record date: 20220927