CN109192243B - 染色体比例的修正方法、装置、介质 - Google Patents

染色体比例的修正方法、装置、介质 Download PDF

Info

Publication number
CN109192243B
CN109192243B CN201810915839.6A CN201810915839A CN109192243B CN 109192243 B CN109192243 B CN 109192243B CN 201810915839 A CN201810915839 A CN 201810915839A CN 109192243 B CN109192243 B CN 109192243B
Authority
CN
China
Prior art keywords
group
proportion
target
chromosome
ratio
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.)
Active
Application number
CN201810915839.6A
Other languages
English (en)
Other versions
CN109192243A (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.)
Chengdu Fandi Medical Laboratory Co ltd
Original Assignee
Chengdu Fandi Medical Laboratory 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 Chengdu Fandi Medical Laboratory Co ltd filed Critical Chengdu Fandi Medical Laboratory Co ltd
Priority to CN201810915839.6A priority Critical patent/CN109192243B/zh
Publication of CN109192243A publication Critical patent/CN109192243A/zh
Application granted granted Critical
Publication of CN109192243B publication Critical patent/CN109192243B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Biological Materials (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

一种染色体比例的修正方法、装置、介质,所述方法包括:一种染色体比例的修正方法,包括:根据预设的参照组内的目标染色体,判断样本组内的目标染色体是否正常;将所述样本组内正常的目标染色体形成对比组;获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例之间的差异值;根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正;当确定对所述样本组内的目标染色体的比例进行修正时,根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正。采用上述方案,可以提高无创产前检验的结果准确度。

Description

染色体比例的修正方法、装置、介质
技术领域
本发明涉及基因检测领域,尤其涉及一种染色体比例的修正方法、装置、介质。
背景技术
如今,人们通常会在胎儿出生前进行产前诊断,以判断胎儿是否存在患有21-三体综合征、18-三体综合征以及13-三体综合征的风险。
现有技术中,进行产前检验时,人们为了避免过大的创伤性,同时保证检测结果较高的准确度,人们通常选择使用无创产前胎儿DNA检测技术(Non-Invasive PrenatalTesting,NIPT)进行产前检验。NIPT技术采用Z检验进行结果筛选,通过孕妇的目标染色体,即孕妇的13号染色体、18号染色体和21号染色体比例与参照库的平均13号染色体、18号染色体和21号染色体比例进行比较得出检测结果。
然而,使用无创产前胎儿DNA检测技术时,某些情况下,由于DNA提取、建库、上机测序、操作人员、仪器和实验室环境等因素存在诸多变量,导致批次内正常样本的13号染色体、18号染色体和21号染色体的比例会与参考库存在差异。这种差异可以表示为一种随机波动,但当此种波动有显著性差异时,其会影响到无创产前检测结果的准确性。
发明内容
本发明解决的技术问题是批次内样本染色体比例与参考库差异过大时,提升该批次样本检测结果的准确性。
为解决上述技术问题,本发明实施例提供一种染色体比例的修正方法,包括:根据预设的参照组内的目标染色体数据,判断样本组内的目标染色体是否正常;将所述样本组内正常的目标染色体数据形成对比组;获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例之间的差异值;根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正;当确定对所述样本组内的目标染色体的比例进行修正时,根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正。
可选的,获取所述参照组内正常的目标染色体的平均比例;获取所述样本组内各个目标染色体的比例;根据所述参照组内目标染色体的平均比例判断所述样本组内目标染色体的比例是否正常。
可选的,所述参照组内目标染色体的平均比例meanK使用以下公式计算:
RatioKm=Ri/Rc
meanK=∑RatioKm/M;
其中Ri为参照组内第m条正常的目标染色体上片段的数目,Rc为所述第m条正常的目标染色体对应的所有常染色体上片段的数目,RatioKm为所述第m条目标染色体的比例,M为所述参照组内正常的目标染色体的条数,∑RatioKm为参照组内M条正常的目标染色体的比例的和,1≤m≤M。
可选的,当(RatioLn-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常;其中,RatioLn为所述样本组内第n条目标染色体的比例,1≤n≤N,N为所述样本组内的目标染色体的数目,SDK为所述参照组内M条正常的目标染色体的比例的标准差。
可选的,将所述参照组中正常的M条目标染色体的比例作为第一组;将所述对比组中的J条目标染色体的比例作为第二组,J为所述对比组中目标染色体的数目;1≤J≤N;对所述第一组与第二组进行假设检验,计算所述第一组与第二组之间没有差异值的概率。
可选的,若所述第一组与第二组之间没有差异值的概率小于等于预设的阈值,则判定对样本组内目标染色体的比例进行修正。
可选的,所述对比组内目标染色体的平均比例meanJ由以下公式计算:
meanJ=∑RatioLj/J,
其中,RatioLj为所述对比组中第j条目标染色体的比例,∑RatioLj为对比组中的J条目标染色体的比例的和值;
所述修正量Δr由以下公式计算:
Δr=meanJ-meank
由以下公式对样本组内目标染色体的比例进行修正:
Ratio’Ln=RatioLn-Δr,
其中,Ratio’Ln为修正后的样本组内目标染色体的比例。
可选的,使用以下公式判断所述样本组内目标染色体是否正常:
当(Ratio’Ln-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常。
本发明还提供一种染色体比例的修正装置,包括:第一判断单元,用于根据预设的参照组内的目标染色体数据,判断样本组内的目标染色体是否正常;对比单元,用于将所述样本组内正常的目标染色体数据形成对比组;差异单元,用于获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例之间的差异值;第二判断单元,用于根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正;修正单元,用于当确定对所述样本组内的目标染色体的比例进行修正时,根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正。
可选的,所述第一判断单元,用于获取所述参照组内正常的目标染色体的平均比例;获取所述样本组内各个目标染色体的比例;根据所述参照组内目标染色体的平均比例判断所述样本组内目标染色体的比例是否正常。
可选的,所述第一判断单元,用于使用以下公式计算所述参照组内目标染色体的平均比例meanK
RatioKm=Ri/Rc
meanK=∑RatioKm/M;
其中,Ri为参照组内第m条正常的目标染色体上片段的数目,Rc为所述第m条正常的目标染色体对应的所有常染色体上片段的数目,RatioKm为所述第m条目标染色体的比例,M为所述参照组内正常的目标染色体的条数,∑RatioKm为参照组内M条正常的目标染色体的比例的和,1≤m≤M。
可选的,所述第一判断单元,用于当(RatioLn-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常;其中,RatioLn为所述样本组内第n条目标染色体的比例,1≤n≤N,N为所述样本组内的目标染色体的数目,SDK为所述参照组内M条正常的目标染色体的比例的标准差。
可选的,所述差异单元,用于将所述参照组中正常的M条目标染色体的比例作为第一组;将所述对比组中的J条目标染色体的比例作为第二组,J为所述对比组中目标染色体的数目;1≤J≤N;对所述第一组与第二组进行假设检验,计算所述第一组与第二组之间没有差异值的概率。
可选的,所述第二判断单元,用于若所述第一组与第二组之间没有差异值的概率小于等于预设的阈值,则判定对样本组内目标染色体的比例进行修正。
可选的,所述修正单元,用于由以下公式计算所述对比组内目标染色体的平均比例meanJ
meanJ=∑RatioLj/J,
其中,RatioLj为所述对比组中第j条目标染色体的比例,∑RatioLj为对比组中的J条目标染色体的比例的和值;
由以下公式计算所述修正量Δr:
Δr=meanJ-meank
由以下公式对样本组内目标染色体的比例进行修正:
Ratio’Ln=RatioLn-Δr,
其中,Ratio’Ln为修正后的样本组内目标染色体的比例。
可选的,所述修正单元,还用于使用以下公式判断所述样本组内目标染色体是否正常:
当(Ratio’Ln-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常。
本发明还提供一种计算机可读存储介质,其上存储有计算机指令,所述计算机指令运行时执行上述任一种的染色体比例的修正方法的步骤。
本发明还提供一种染色体比例的修正装置,包括存储器和处理器,所述存储器上存储有计算机指令,所述计算机指令运行时所述处理器执行上述任一种的染色体比例的修正方法的步骤。
与现有技术相比,本发明实施例的技术方案具有以下有益效果:
通过将已完成测序的孕有正常胎儿的孕妇的目标染色体的比例作为参照组,使用参照组判断样本组中孕妇的目标染色体是否正常,将样本组中正常的目标染色体作为对比组,根据所述参照组内目标染色体的平均比例与所述对比组内目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内目标染色体的比例进行修正。在对样本组的目标染色体的比例修正后,使得样本组中原本正常的胎儿不会被检验得出患有21-三体综合征、18-三体综合征以及13-三体综合征的风险。
进一步,在对样本组的目标染色体的比例进行修正后,根据参照组对样本组进行检验,得到修正后降低的检验错误的概率。
附图说明
图1是本发明实施例提供的一种染色体比例的修正方法的示意图;
图2是本发明实施例提供的一种染色体比例的修正装置的结构示意图。
具体实施方式
现有技术中,进行产前检验时,人们为了避免过大的创伤性,同时保证检测结果较高的准确度,人们通常选择使用NIPT技术进行产前检验。NIPT技术采用Z检验进行结果筛选,通过孕妇的目标染色体,即孕妇的13号染色体、18号染色体和21号染色体比例与参照库的平均13号染色体、18号染色体和21号染色体比例进行比较得出检测结果。
然而,使用无创产前胎儿DNA检测技术时,某些情况下,由于DNA提取、建库、上机测序、操作人员、仪器和实验室环境等因素存在诸多变量,导致批次内正常样本的13号染色体、18号染色体和21号染色体的比例会与参考库存在差异。这种差异可以表示为一种随机波动,但当此种波动有显著性差异时,其会影响到无创产前检测结果的准确性。
本发明实施例中,通过将已完成测序的孕有正常胎儿的孕妇的目标染色体的比例作为参照组,使用参照组判断样本组中孕妇的目标染色体是否正常,将样本组中正常的目标染色体作为对比组,根据所述参照组内目标染色体的平均比例与所述对比组内目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内目标染色体的比例进行修正。在对样本组的目标染色体的比例修正后,使得样本组中原本正常的胎儿不会有被检验得出患有21-三体综合征、18-三体综合征以及13-三体综合征的风险。
为使本发明的上述目的、特征和有益效果能够更为明显易懂,下面结合附图对本发明的具体实施例做详细的说明。
参阅图1,本发明实施例提供了一种染色体比例的修正方法,其中具体步骤包括:
步骤S101,根据预设的参照组内的目标染色体数据,判断样本组内的目标染色体是否正常。
在具体实施中,参照组可以包括多名孕妇的染色体数据,这些孕妇的染色体数据为在实验室经过检验的数据,且其中包括:已确定其中孕有正常胎儿的孕妇对应的正常的目标染色体数据,以及孕有患三体综合征的胎儿的孕妇对应的目标染色体数据。
在实际应用中,正常的目标染色体即可表征孕妇的染色体呈阴性,孕妇当孕妇的染色体呈阳性时表征孕妇孕有患三体综合征的胎儿。
在具体实施中,参照组中孕妇染色体数据越多,相对应得到的无创检测结果越准确。
在实际应用中,目标染色体可以是13号染色体、18号染色体或21号染色体中的任一个。
在具体实施中,样本组可以由孕妇的染色体数据组成,其中孕妇的目标染色体正常与否未知。
在具体实施中,将样本组中的孕妇的目标染色体数据与参照组中的目标染色体数据相对比,确定样本组中正常的目标染色体。
步骤S102,将所述样本组内正常的目标染色体数据形成对比组。
在具体实施中,将样本组与参照组比对后,选出样本组中所有的正常的目标染色体,并将这些正常的目标染色体形成对比组。
步骤S103,获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例之间的差异值。
在具体实施中,染色体比例为某一染色体上的片段(Read)数目与该染色体对应的所有常染色体上的片段数目之和的比值。
例如,1号染色体的比例为1号染色体上的片段数目与1-22号所有常染色体上的片段数目之和的比值。
在具体实施中,可以将参照组内的各个目标染色体的比例数据作为一组,将对比组内的各个目标染色体的比例数据作为另一组,获取两组数据之间的差异值。
步骤S104,根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正。
在具体实施中,阈值可以由临床实验数据确定。
在具体实施中,可以根据获取的差异值与阈值之间的大小关系确定是否对样本组内的目标染色体的比例进行修正。
步骤S105,当确定对所述样本组内的目标染色体的比例进行修正时,根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正。
在具体实施中,若确定对样本组内的目标染色体比例进行修正时,说明可能由于地区差异等原因,样本组内正常目标染色体的比例与参照组内正常目标染色体的平均比例之间的偏差值过大,过大的偏差值进而导致样本组内某些正常的目标染色体被判定为阳性。
在具体实施中,为了修正样本组内正常的目标染色体的偏差值,根据参照组内的目标染色体的平均比例与对比组内的目标染色体的平均比例得到修正量。在得到修正量之后,可以将样本组内的各个目标染色体的比例减去修正量,完成修正过程。
本发明实施例中,在根据参照组内的目标染色体判断样本组内的目标染色体是否正常的过程中,可以先在参照组中选取正常的目标染色体,获取这些正常的目标染色体的平均比例,并获取样本组内各个目标染色体的比例。之后,根据参照组内正常的目标染色体的平均比例以及样本组内的目标染色体的比例,判断样本组内的目标染色体是否正常。
本发明实施例中,获取参照组内正常目标染色体的平均比例meanK的过程为,确定参照组内的M条正常的目标染色体,其中第m条正常的目标染色体的比例RatioKm的计算公式如下:
RatioKm=Ri/Rc, (1)
其中,Ri为第m条正常的目标染色体上片段的数目,Rc为所述第m条正常的目标染色体对应的所有常染色体上片段的数目。
在具体实施中,在获取参照组内M条正常目标染色体的比例后,参照组内正常目标染色体的平均比例meanK的计算公式如下:
meanK=∑RatioKm/M, (2)
其中,∑RatioKm为M条正常的目标染色体的比例之和,1≤m≤M。
在具体实施中,样本组内的N条目标染色体的比例RatioLn的获取方法可以参照上述参照组内正常的目标染色体比例的获取方法;1≤n≤N。
本发明实施例中,可以使用Z检验算法,根据参照组内的正常目标染色体的平均比例meanK来判断样本组的目标染色体是否正常。
在具体实施中,Z检验为无创产前检测方法中常用的一种判断目标染色体为阴性或阳性的方法,用户也可以根据实际情况使用其他检验方法判断目标染色体为阴性或阳性。
在具体实施中,使用Z检验前获取参照组内M条正常的目标染色体比例的标准差SDK。在本发明实施例中,可以采用如下公式计算Z检验:
(RatioLn-meanK)/SDK, (3)
当Z检验得到的值小于3时,则判定第n条目标染色体为正常。
本发明实施例中,挑选出样本组内J条正常的目标染色体作为对比组。
本发明实施例中,将参照组中M条正常的目标染色体比例作为第一组;将对比组中的J条目标染色体的比例作为第二组。
例如,参照组内有3条正常的21号染色体,对应的比例分别为0.1、0.2、0.3,则第一组内的数据为0.1、0.2和0.3。对比组内有3条21号染色体,对应的比例分别为0.1、0.3、0.4,则第二组内的数据为0.1、0.3和0.4。
本发明实施例中,计算所述第一组的数据与第二组的数据之间的差异值。
在具体实施中,可以使用秩和检验计算第一组的数据与第二组的数据之间的差异值。在实际应用中可知,秩和检验是一种假设检验,通过假设两组数据之间没有差异值,可以计算得到两组数据之间没有差异值的概率。
可以理解的是,也可以使用其他算法计算第一组的数据与第二组的数据之间的差异值。
本发明实施例中,在通过假设检验计算得到第一组的数据与第二组的数据之间没有差异值的概率后,若概率小于预设的阈值,则说明两组数据之间的差异值较大,可以对样本组内的目标染色体比例进行修正。
在具体实施中,阈值可以由用户根据实际情况设定。
本发明实施例中,确定对样本组内的目标染色体修正后,可以计算修正量。
本发明实施例中,修正量Δr可以由以下公式计算:
Δr=meanJ-meank, (4)
其中,meanJ为对比组内目标染色体的平均比例,meanJ=∑RatioLj/J,RatioLj为所述对比组中第j条目标染色体的比例,∑RatioLj为对比组中的J条目标染色体的比例的和值。
本发明实施例中,获取修正量Δr后,可以使用修正量对样本组内的目标染色体的比例进行修正。
在具体实施中,可以由以下公式对样本组内目标染色体的比例RatioLn进行修正:
Ratio’Ln=RatioLn-Δr, (5)
其中,Ratio’Ln为修正后的样本组内目标染色体的比例。
在具体实施中,在对样本组的目标染色体进行修正后,可以修正样本组内正常孕妇的目标染色体比例与参照组内正常孕妇的目标染色体比例之间的偏差值,使得偏差值在正常范围内,进而使得样本组内的正常孕妇的胎儿不会被检测出患有三体综合征。
本发明实施例中,在获取修正后的样本组内目标染色体的比例Ratio’Ln后,可以使用以下Z检验公式根据Ratio’Ln重新检测样本组内目标染色体是否正常:
(Ratio’Ln-meanK)/SDK, (6)
当(Ratio’Ln-meanK)/SDK<3时,则所述样本组内第n条目标染色体为正常。
参阅图2,本发明还提供一种染色体比例的修正装置20,包括:
第一判断单元201,用于根据预设的参照组内的目标染色体数据,判断样本组内的目标染色体是否正常;
对比单元202,用于将所述样本组内正常的目标染色体数据形成对比组;
差异单元203,用于获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例之间的差异值;
第二判断单元204,用于根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正;
修正单元205,用于当确定对所述样本组内的目标染色体的比例进行修正时,根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正。
本发明实施例中,所述第一判断单元201,用于获取所述参照组内正常的目标染色体的平均比例;获取所述样本组内各个目标染色体的比例;根据所述参照组内目标染色体的平均比例判断所述样本组内目标染色体的比例是否正常。
本发明实施例中,所述第一判断单元201,用于使用以下公式计算所述参照组内目标染色体的平均比例meanK
RatioKm=Ri/Rc
meanK=∑RatioKm/M;
其中,Ri为参照组内第m条正常的目标染色体上片段的数目,Rc为所述第m条正常的目标染色体对应的所有常染色体上片段的数目,RatioKm为所述第m条目标染色体的比例,M为所述参照组内正常的目标染色体的条数,∑RatioKm为参照组内M条正常的目标染色体的比例的和,1≤m≤M。
本发明实施例中,所述第一判断单元201,用于当(RatioLn-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常;其中,RatioLn为所述样本组内第n条目标染色体的比例,1≤n≤N,N为所述样本组内的目标染色体的数目,SDK为所述参照组内M条正常的目标染色体的比例的标准差。
本发明实施例中,所述差异单元203,用于将所述参照组中正常的M条目标染色体的比例作为第一组;将所述对比组中的J条目标染色体的比例作为第二组,J为所述对比组中目标染色体的数目;1≤J≤N;对所述第一组与第二组进行假设检验,计算所述第一组与第二组之间没有差异值的概率。
本发明实施例中,所述第二判断单元204,用于若所述第一组与第二组之间没有差异值的概率小于等于预设的阈值,则判定对样本组内目标染色体的比例进行修正。
本发明实施例中,所述修正单元205,用于由以下公式计算所述对比组内目标染色体的平均比例meanJ
meanJ=∑RatioLj/J,
其中,RatioLj为所述对比组中第j条目标染色体的比例,∑RatioLj为对比组中的J条目标染色体的比例的和值;
由以下公式计算所述修正量Δr:
Δr=meanJ-meank
由以下公式对样本组内目标染色体的比例进行修正:
Ratio’Ln=RatioLn-Δr,
其中,Ratio’Ln为修正后的样本组内目标染色体的比例。
本发明实施例中,所述修正单元205,还用于使用以下公式判断所述样本组内目标染色体是否正常:当(Ratio’Ln-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常。
本发明还提供一种计算机可读存储介质,其上存储有计算机指令,所述计算机指令运行时执行上述任一种的染色体比例的修正方法的步骤。
本发明还提供一种染色体比例的修正装置,包括存储器和处理器,所述存储器上存储有计算机指令,所述计算机指令运行时所述处理器执行上述任一种的染色体比例的修正方法的步骤。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指示相关的硬件来完成,该程序可以存储于一计算机可读存储介质中,存储介质可以包括:ROM、RAM、磁盘或光盘等。
虽然本发明披露如上,但本发明并非限定于此。任何本领域技术人员,在不脱离本发明的精神和范围内,均可作各种更动与修改,因此本发明的保护范围应当以权利要求所限定的范围为准。

Claims (14)

1.一种染色体比例的修正方法,其特征在于,包括:
根据预设的参照组内的目标染色体数据,判断样本组内的目标染色体是否正常,包括:获取所述参照组内正常的目标染色体的平均比例;获取所述样本组内各个目标染色体的比例;根据所述参照组内目标染色体的平均比例判断所述样本组内目标染色体的比例是否正常;所述获取所述参照组内正常的目标染色体的平均比例,包括:所述参照组内目标染色体的平均比例meanK使用以下公式计算:RatioKm=Ri/Rc,meanK=∑RatioKm/M;其中,Ri为参照组内第m条正常的目标染色体上片段的数目,Rc为所述第m条正常的目标染色体对应的所有常染色体上片段的数目,RatioKm为所述第m条目标染色体的比例,M为所述参照组内正常的目标染色体的条数,∑RatioKm为参照组内M条正常的目标染色体的比例的和,1≤m≤M;
将所述样本组内正常的目标染色体数据形成对比组;
获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例之间的差异值;
根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正;
当确定对所述样本组内的目标染色体的比例进行修正时,根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正。
2.根据权利要求1所述的染色体比例的修正方法,其特征在于,所述根据所述参照组内目标染色体的平均比例判断所述样本组内目标染色体的比例是否正常,包括:
当(RatioLn-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常;其中,RatioLn为所述样本组内第n条目标染色体的比例,1≤n≤N,N为所述样本组内的目标染色体的数目,SDK为所述参照组内M条正常的目标染色体的比例的标准差。
3.根据权利要求2所述的染色体比例的修正方法,其特征在于,所述获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例的差异值,包括:
将所述参照组中正常的M条目标染色体的比例作为第一组;
将所述对比组中的J条目标染色体的比例作为第二组,J为所述对比组中目标染色体的数目;1≤J≤N;
对所述第一组与第二组进行假设检验,计算所述第一组与第二组之间没有差异值的概率。
4.根据权利要求3所述的染色体比例的修正方法,其特征在于,所述根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正,包括:
若所述第一组与第二组之间没有差异值的概率小于等于预设的阈值,则判定对样本组内目标染色体的比例进行修正。
5.根据权利要求4所述的染色体比例的修正方法,其特征在于,所述根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正,包括:
所述对比组内目标染色体的平均比例meanJ由以下公式计算:
meanJ=∑RatioLj/J,
其中,RatioLj为所述对比组中第j条目标染色体的比例,∑RatioLj为对比组中的J条目标染色体的比例的和值;
所述修正量Δr由以下公式计算:
Δr=meanJ-meank
由以下公式对样本组内目标染色体的比例进行修正:
Ratio’Ln=RatioLn-Δr,
其中,Ratio’Ln为修正后的样本组内目标染色体的比例。
6.根据权利要求5所述的染色体比例的修正方法,其特征在于,在使用所述修正量对所述样本组内的目标染色体的比例进行修正之后,还包括:
使用以下公式判断所述样本组内目标染色体是否正常:
当(Ratio’Ln-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常。
7.一种染色体比例的修正装置,其特征在于,包括:
第一判断单元,用于根据预设的参照组内的目标染色体数据,判断样本组内的目标染色体是否正常,包括:获取所述参照组内正常的目标染色体的平均比例;获取所述样本组内各个目标染色体的比例;根据所述参照组内目标染色体的平均比例判断所述样本组内目标染色体的比例是否正常;所述获取所述参照组内正常的目标染色体的平均比例,包括:所述参照组内目标染色体的平均比例meanK使用以下公式计算:RatioKm=Ri/Rc,meanK=∑RatioKm/M;其中,Ri为参照组内第m条正常的目标染色体上片段的数目,Rc为所述第m条正常的目标染色体对应的所有常染色体上片段的数目,RatioKm为所述第m条目标染色体的比例,M为所述参照组内正常的目标染色体的条数,∑RatioKm为参照组内M条正常的目标染色体的比例的和,1≤m≤M;;
对比单元,用于将所述样本组内正常的目标染色体数据形成对比组;
差异单元,用于获取所述参照组内目标染色体的比例与所述对比组内的目标染色体的比例之间的差异值;
第二判断单元,用于根据预设的阈值以及所述差异值判断是否对样本组内的目标染色体的比例进行修正;
修正单元,用于当确定对所述样本组内的目标染色体的比例进行修正时,根据所述参照组内的目标染色体的平均比例与所述对比组内的目标染色体的平均比例得到修正量,使用所述修正量对所述样本组内的目标染色体的比例进行修正。
8.根据权利要求7所述的染色体比例的修正装置,其特征在于,所述第一判断单元,用于当(RatioLn-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常;
其中,RatioLn为所述样本组内第n条目标染色体的比例,1≤n≤N,N为所述样本组内的目标染色体的数目,SDK为所述参照组内M条正常的目标染色体的比例的标准差。
9.根据权利要求8所述的染色体比例的修正装置,其特征在于,所述差异单元,用于将所述参照组中正常的M条目标染色体的比例作为第一组;将所述对比组中的J条目标染色体的比例作为第二组,J为所述对比组中目标染色体的数目;1≤J≤N;对所述第一组与第二组进行假设检验,计算所述第一组与第二组之间没有差异值的概率。
10.根据权利要求9所述的染色体比例的修正装置,其特征在于,所述第二判断单元,用于若所述第一组与第二组之间没有差异值的概率小于等于预设的阈值,则判定对样本组内目标染色体的比例进行修正。
11.根据权利要求10所述的染色体比例的修正装置,其特征在于,所述修正单元,用于由以下公式计算所述对比组内目标染色体的平均比例meanJ
meanJ=∑RatioLj/J,
其中,RatioLj为所述对比组中第j条目标染色体的比例,∑RatioLj为对比组中的J条目标染色体的比例的和值;
由以下公式计算所述修正量Δr:
Δr=meanJ-meank
由以下公式对样本组内目标染色体的比例进行修正:
Ratio’Ln=RatioLn-Δr,
其中,Ratio’Ln为修正后的样本组内目标染色体的比例。
12.根据权利要求11所述的染色体比例的修正装置,其特征在于,所述修正单元,还用于使用以下公式判断所述样本组内目标染色体是否正常:
当(Ratio’Ln-meanK)/SDK<3时,判定所述样本组内第n条目标染色体为正常。
13.一种计算机可读存储介质,其上存储有计算机指令,其特征在于,所述计算机指令运行时执行权利要求1~6任一项所述的染色体比例的修正方法的步骤。
14.一种染色体比例的修正装置,包括存储器和处理器,所述存储器上存储有计算机指令,其特征在于,所述计算机指令运行时所述处理器执行权利要求1~6任一项所述的染色体比例的修正方法的步骤。
CN201810915839.6A 2018-08-13 2018-08-13 染色体比例的修正方法、装置、介质 Active CN109192243B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810915839.6A CN109192243B (zh) 2018-08-13 2018-08-13 染色体比例的修正方法、装置、介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810915839.6A CN109192243B (zh) 2018-08-13 2018-08-13 染色体比例的修正方法、装置、介质

Publications (2)

Publication Number Publication Date
CN109192243A CN109192243A (zh) 2019-01-11
CN109192243B true CN109192243B (zh) 2021-03-12

Family

ID=64921619

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810915839.6A Active CN109192243B (zh) 2018-08-13 2018-08-13 染色体比例的修正方法、装置、介质

Country Status (1)

Country Link
CN (1) CN109192243B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112712853B (zh) * 2020-12-31 2023-11-21 北京优迅医学检验实验室有限公司 一种无创产前检测装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SG191757A1 (en) * 2011-06-29 2013-08-30 Bgi Health Service Co Ltd Noninvasive detection of fetal genetic abnormality
CN104120181B (zh) * 2011-06-29 2017-06-09 深圳华大基因股份有限公司 对染色体测序结果进行gc校正的方法及装置
US10612080B2 (en) * 2014-09-22 2020-04-07 Roche Molecular Systems, Inc. Digital PCR for non-invasive prenatal testing
US20180211002A1 (en) * 2015-07-13 2018-07-26 Agilent Technologies Belgium Nv System and methodology for the analysis of genomic data obtained from a subject

Also Published As

Publication number Publication date
CN109192243A (zh) 2019-01-11

Similar Documents

Publication Publication Date Title
CN108573125B (zh) 一种基因组拷贝数变异的检测方法及包含该方法的装置
JP2016539630A5 (zh)
CN111304303B (zh) 微卫星不稳定的预测方法及其应用
EP3293270B1 (en) Reagent kit, apparatus, and method for detecting chromosome aneuploidy
CN108256292B (zh) 一种拷贝数变异检测装置
CN104156631B (zh) 染色体三倍体检验方法
EP3518974A1 (en) Noninvasive prenatal screening using dynamic iterative depth optimization
CN110268072B (zh) 确定旁系同源基因的方法和***
KR20010042824A (ko) 화학적 및 생물학적 분석의 평가방법
JP2018514234A5 (zh)
KR101678962B1 (ko) 대규모 병렬형 게놈서열분석 방법을 이용한 비침습적 산전검사 장치 및 방법
CN114292912A (zh) 一种变体核酸的检测方法
CN105825076A (zh) 消除常染色体内和染色体间gc偏好的方法及检测***
EP2924126B1 (en) Method for using microRNA (miRNA) for detection of endometriosis
CN109979529A (zh) Cnv检测装置
CN109192243B (zh) 染色体比例的修正方法、装置、介质
CN107463797B (zh) 高通量测序的生物信息分析方法及装置、设备及存储介质
CN104520437A (zh) 一种染色体非整倍性检测方法及装置
CN114592074A (zh) 一种与胎龄相关的靶标基因组合及其应用
CN113789371A (zh) 一种基于批次矫正的拷贝数变异的检测方法
CN111696622B (zh) 一种校正和评估变异检测软件检测结果的方法
RU2543155C1 (ru) Способ неинвазивной диагностики анеуплоидий плода методом секвенирования
CN110970089B (zh) 胎儿浓度计算的预处理方法、预处理装置及其应用
CN117079723A (zh) 一种与肌萎缩侧索硬化症相关的生物标志物、诊断模型及其应用
CN117153258A (zh) 校正测序数据、检测染色体非整倍体的方法和装置

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