CN104846089B - 一种孕妇外周血中胎儿游离dna比例的定量方法 - Google Patents
一种孕妇外周血中胎儿游离dna比例的定量方法 Download PDFInfo
- Publication number
- CN104846089B CN104846089B CN201510225771.5A CN201510225771A CN104846089B CN 104846089 B CN104846089 B CN 104846089B CN 201510225771 A CN201510225771 A CN 201510225771A CN 104846089 B CN104846089 B CN 104846089B
- Authority
- CN
- China
- Prior art keywords
- snp
- dna
- snp site
- conc
- err
- 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
Links
- 230000008774 maternal effect Effects 0.000 title claims abstract description 26
- 230000001605 fetal effect Effects 0.000 title claims abstract description 19
- 238000013459 approach Methods 0.000 title claims abstract description 16
- 210000003754 fetus Anatomy 0.000 claims abstract description 65
- 238000013461 design Methods 0.000 claims abstract description 17
- 239000012634 fragment Substances 0.000 claims description 27
- 238000012163 sequencing technique Methods 0.000 claims description 21
- 238000012216 screening Methods 0.000 claims description 18
- 239000007788 liquid Substances 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 11
- 239000002131 composite material Substances 0.000 claims description 10
- 108090000790 Enzymes Proteins 0.000 claims description 9
- 102000004190 Enzymes Human genes 0.000 claims description 9
- 102000003960 Ligases Human genes 0.000 claims description 8
- 108090000364 Ligases Proteins 0.000 claims description 8
- 210000000349 chromosome Anatomy 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 8
- 238000013467 fragmentation Methods 0.000 claims description 8
- 238000006062 fragmentation reaction Methods 0.000 claims description 8
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 230000003321 amplification Effects 0.000 claims description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 6
- 239000000872 buffer Substances 0.000 claims description 4
- 230000035772 mutation Effects 0.000 claims description 4
- 238000012408 PCR amplification Methods 0.000 claims description 3
- 239000007853 buffer solution Substances 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000002864 sequence alignment Methods 0.000 claims description 3
- 238000011534 incubation Methods 0.000 claims description 2
- 210000002381 plasma Anatomy 0.000 abstract description 17
- 238000007619 statistical method Methods 0.000 abstract description 3
- 238000000034 method Methods 0.000 description 22
- 239000002585 base Substances 0.000 description 21
- 210000004369 blood Anatomy 0.000 description 8
- 239000008280 blood Substances 0.000 description 8
- 108090000623 proteins and genes Proteins 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 5
- 230000011987 methylation Effects 0.000 description 3
- 238000007069 methylation reaction Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 239000003153 chemical reaction reagent Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000002156 mixing Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 239000002773 nucleotide Substances 0.000 description 2
- 125000003729 nucleotide group Chemical group 0.000 description 2
- 210000005259 peripheral blood Anatomy 0.000 description 2
- 239000011886 peripheral blood Substances 0.000 description 2
- BASFCYQUMIYNBI-UHFFFAOYSA-N platinum Chemical compound [Pt] BASFCYQUMIYNBI-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 230000006641 stabilisation Effects 0.000 description 2
- 238000011105 stabilization Methods 0.000 description 2
- 206010000234 Abortion spontaneous Diseases 0.000 description 1
- 208000031404 Chromosome Aberrations Diseases 0.000 description 1
- 208000011359 Chromosome disease Diseases 0.000 description 1
- 206010067477 Cytogenetic abnormality Diseases 0.000 description 1
- 201000010374 Down Syndrome Diseases 0.000 description 1
- 238000012896 Statistical algorithm Methods 0.000 description 1
- 210000002593 Y chromosome Anatomy 0.000 description 1
- 239000003513 alkali Substances 0.000 description 1
- 230000003322 aneuploid effect Effects 0.000 description 1
- 208000036878 aneuploidy Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007622 bioinformatic analysis Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 208000024971 chromosomal disease Diseases 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 208000015994 miscarriage Diseases 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229910052697 platinum Inorganic materials 0.000 description 1
- 102000054765 polymorphisms of proteins Human genes 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000746 purification Methods 0.000 description 1
- 239000002096 quantum dot Substances 0.000 description 1
- 108091008146 restriction endonucleases Proteins 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 208000000995 spontaneous abortion Diseases 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开一种孕妇外周血中胎儿游离DNA比例的定量方法。本发明利用设计的引物将血浆中包含设计SNP的DNA提取出来，经过一系列处理后上机，将测序产物比对到人类基因组序列后确定每个SNP位点的碱基。通过分析SNP中各碱基的比例来定量胎儿游离DNA浓度。本发明还可通过多个SNP来校正引入的误差，用统计学的方法对多个SNP同时计算，提高计算出的胎儿游离DNA浓度准确性。
Description
技术领域
本发明涉及DNA定量领域，尤其涉及一种孕妇外周血中胎儿游离DNA比例的定量方法。
背景技术
无创产前基因检测是通过采集孕妇外周血(5ml)，提取其中的游离DNA，采用高通量测序技术，结合生物信息分析，得出胎儿患染色体非整倍体（21-三体又称唐氏综合征，18-三体，13-三体）的风险。研究发现，从孕4周开始，在孕妇外周血中即可检测到胎儿的游离DNA。随着孕周增大，胎儿游离DNA含量也随之增加。孕12周后，通过抽取孕妇外周血并从中提取出胎儿游离DNA，利用新一代基因测序技术并结合生物信息学分析手段，便可准确判断胎儿是否患有染色体病。该方法最佳检测时间为孕早、中期，具有无创取样、无流产风险、高灵敏度，准确性高的特点。
无创产前基因检测技术通过分析孕妇外周血中DNA来判断胎儿是否存在染色体非整倍体，而不是针对性的分析胎儿的DNA，因此当胎儿DNA比例过低时（<3%）,可能因为胎儿DNA量太少而不能被检测出来染色体是否有异常，所以胎儿DNA的含量对检测结果的准确性和敏感性都起到了重要的作用。
在一份检测结果中，如果知道胎儿DNA的比例将是对一份结果准确性的最好的复核。通过准确的定量胎儿DNA的比例，设定该DNA比例下三体所对应的阈值，对检测结果的准确性起到复核作用，同时能够避免孕周大但胎儿游离DNA浓度低造成的假阴性。虽然现在行业内都通过接收孕12周以后的孕妇来减少这种情况的发生，但孕妇个体之间的差异还是可能存在胎儿DNA比例较低（2%的孕妇会因为胎儿DNA含量低而被要求重新抽血）。现有技术通过对男胎中含有的Y染色体信息来判断男胎的浓度，对女胎还没有有效的方法。
现有技术中，甲基化定量胎儿游离DNA浓度的方法被很多研究证实可以用来定量胎儿游离DNA浓度。采用该方法抽取外周血后，经 DNA 分离、甲基化敏感限制性内切酶消除样本中母源性( 非甲基化) DNA ，未分解的胎儿 DNA 在PCR扩增，最后定量胎儿DNA量。通过计算胎儿DNA量占外周血DNA量的比例来计算胎儿游离DNA浓度。
但该方法还存在以下两个缺点：
1、被确定的稳定的甲基化标志基因较少，它需要具备两个条件：1）孕妇和胎儿甲基化的差异要大，孕妇要完全不甲基化，胎儿要完全甲基化。2）甲基化差异必须保持稳定,不会因为个人之间在数量或在怀孕期不同时期存在差异。
2、甲基化定量胎儿游离DNA浓度的方法是通过一个基因，因为量少的原因需要经过PCR来扩增，其中引入的误差不好被校正。
另外，还有一种方法是荧光定量PCR或者绝对定量PCR法，其也被用来定量胎儿游离DNA浓度，但其有个共同的缺点就是可用的目标基因数据少，且能获得DNA量也少，很容易有偏差。
因此，现有技术还有待于改进和发展。
发明内容
鉴于上述现有技术的不足，本发明的目的在于提供一种定量孕妇外周血中胎儿游离DNA比例的方法，旨在解决现有的孕妇外周血DNA定量方法准确性低、容易有偏差的问题。
本发明的技术方案如下：
一种孕妇外周血中胎儿游离DNA比例的定量方法，其中，包括步骤：
A、设计位点：在人类基因组中，找出突变频率在0.4-0.6的SNP位点；
B、设计引物：根据设计的SNP位点，设计出覆盖SNP位点的SNP引物；
C、提取样品DNA；
D、扩增含有SNP位点的DNA片段：将样品DNA、多重PCR反应酶和SNP 引物混匀后，进行PCR反应得到PCR产物；
E、纯化扩增到的PCR产物；
F、末端修复：将PCR产物、末端修复酶与末端修复酶缓冲液混合，室温下孵育进行反应得到混合DNA液；
G、片段筛选：对步骤F所得混合DNA液进行片段筛选及纯化，得到平末端DNA片段；
H、接头连接：将所得的平末端DNA片段与连接酶、连接酶缓冲液、及指定接头混合后，室温下孵育进行反应得到混合DNA液；
I、片段筛选：对步骤H所得混合DNA液进行片段筛选及纯化，得到加接头的DNA片段；
J、PCR扩增：对所得加接头的DNA片段进行PCR 扩增及纯化，得到小片段文库；
K、文库检测：对小片段文库检测浓度和片段大小；
L、高通量测序：对检测合格的小片段文库进行高通量测序；
M、数据预处理：将高通量测序得到的数据首先经过低质量过滤，筛选出长度大于100bp的序列；
N、序列比对：将预处理后的序列与人类基因组序列进行比对；
O、序列筛选：根据比对的结果，筛选出unique read的序列；
P、SNP数据统计：根据比对结果及SNP位点的位置统计出SNP位点的覆盖深度、碱基种类、每种碱基对应数目；
Q、根据统计结果，对每个SNP位点中四种碱基的出现次数进行排序，然后计算出SNP位点的覆盖深度Depth，并根据排序结果求出SNP位点的测序错误率Error；根据每行的覆盖深度Depth和测序错误率Error，求出SNP位点的平均覆盖度和测序错误率的较大值ErrorUp；根据平均覆盖度或者用户自定的覆盖深度，对SNP位点进行过滤；或根据每个SNP位点的测序错误率Error是否大于ErrorUp来过滤；
R、计算每个SNP距离：对每一个SNP位点，计算出可能的胎儿游离DNA浓度Conc在四种SNP组合类型下所有SNP位点到相应胎儿游离DNA浓度的距离TypeDistance；
U、计算出总距离：找出每一SNP位点在四种SNP组合类型下对应最小的TypeDistance，最后求和：AllDistance=sum（TypeDistancei），i=1…N，N为SNP位点总个数；
V、选择最优的总距离：在Conc范围内计算出的AllDistance中选出其中最小的，其所对应的Conc即为实际的胎儿游离DNA浓度。
所述的孕妇外周血中胎儿游离DNA比例的定量方法，其中，所述步骤O中，筛选出错配碱基在4个以下，且只比对到一处的序列，作为unique read。
所述的孕妇外周血中胎儿游离DNA比例的定量方法，其中，所述步骤Q中，ErrorUp=(ErrorMean+ErrorSD*2), Errormean（Errori）i=1…N，N为SNP位点个数。
所述的孕妇外周血中胎儿游离DNA比例的定量方法，其中，所述步骤U中，四种SNP组合类型为：类型1：母亲和胎儿都是纯合SNP；类型2：母亲纯合SNP，胎儿杂合SNP；类型3：母亲杂合SNP，胎儿纯合SNP；类型4：母亲和胎儿都是杂合SNP。
所述的孕妇外周血中胎儿游离DNA比例的定量方法，其中，所述步骤A中，在1-13、18及21这15条染色体上共选出1035个SNP位点。
所述的孕妇外周血中胎儿游离DNA比例的定量方法，其中，Conc的范围为0.025-0.5。
所述的孕妇外周血中胎儿游离DNA比例的定量方法，其中，所述步骤R具体包括：
按以下公式计算SNP距离Distance，Distance=abs(Major * A + Minor * B) /(A**2 + B**2)**0.5,其中A，B是参数，具体如下：
Type=1,A=Err/3,B=-(1-Err);
Type=2,A=Conc/2*(1-Err)+(1-Conc/2)*Err/3,B= -((1-Conc/2)*(1-Err)+Conc/2*Err/3)；
Type=3,A=(1/2-Conc/2)*(1-Err)+(1/2+Conc/2)*Err/3,B= -((1/2+Conc/2)*(1-Err)+(1/2-Conc/2)*Err/3)；
Type=4,A=1,B=-1，每个SNP位点在四种SNP类型下对应有四个距离值，记作TypeDistance=Distance*C，C是一个参数，其在四种SNP类型下分别是：2，1,1,2；根据计算得到的四种类型对应的距离判断该SNP位点所对应的类型；其中Typej（j=1,2,3,4）表示四种类型的SNP组合，Major、Minor分别为SNP位点中四种碱基出现次数排第一及第二的值，Err=mean（Errori）i=1…1035且Errori小于ErrorUp。
有益效果：本发明利用设计的引物将血浆中包含设计SNP位点的DNA提取出来，经过一系列处理后上机，将测序产物比对到人类基因组序列后确定每个SNP位点的碱基。通过分析SNP中各碱基的比例来定量胎儿游离DNA浓度。本发明还可通过多个SNP来校正引入的误差，用统计学的方法对多个SNP同时计算，提高计算出的胎儿游离DNA浓度准确性。
附图说明
图1为本发明一种定量孕妇外周血中胎儿游离DNA比例的方法较佳实施例的流程图。
具体实施方式
本发明提供一种定量孕妇外周血中胎儿游离DNA比例的方法，为使本发明的目的、技术方案及效果更加清楚、明确，以下对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参阅图1，图1为本发明所提供的一种孕妇外周血中胎儿游离DNA比例的定量方法较佳实施例的流程图，如图1所示，其包括步骤：
一种孕妇外周血中胎儿游离DNA比例的定量方法，其特征在于，包括步骤：
S1、设计位点：在人类基因组中，找出突变频率在0.4-0.6的SNP位点；
在人类基因组中，找出突变频率（MAF）在0.4-0.6的SNP位点，在1-13,18,21这15条染色体上共选出1035个SNP位点，每条染色体上SNP位点数目相近。设计过程中位点数目可变，染色体的分布可变。SNP（Single Nucleotide Polymorphisms）指单个碱基上的变异。
本发明的SNP位点数量多，分布广泛。据估计，人类基因组中每1000个核苷酸就有一个SNP位点，它是人类可遗传的变异中最常见的一种。本发明选择1035个SNP位点，这样避免了某一条染色体异常对定量胎儿游离DNA浓度的影响。过多个SNP位点来校正引入的误差。后续可用统计学的方法对多个SNP位点同时计算，胎儿游离DNA浓度计算的准确性更高。
S2、设计引物：根据设计的SNP位点，设计出覆盖SNP位点的SNP引物；
具体可利用ION AMPLISEQ DESIGNER在线设计网站设计出覆盖目标SNP位点的引物，设计引物也可以用其他软件代替。
S3、提取样品DNA；步骤S3也可以与S1同步执行，或者调整顺序。
S4、扩增含有SNP位点的DNA片段：将样品DNA、多重PCR反应酶和SNP 引物混匀后，进行PCR反应得到PCR产物；
S5、纯化扩增到的PCR产物；
S6、末端修复：将PCR产物、末端修复酶与5X 末端修复酶缓冲液混合，室温下孵育进行反应得到混合DNA液；
S7、片段筛选：对步骤F所得混合DNA液进行片段筛选及纯化，得到平末端DNA片段；
S8、接头连接：将所得的平末端DNA片段与ligase Enzyme（连接酶）、10X ligaseBuffer（连接酶缓冲液）、及指定接头混合后，室温下孵育进行反应得到混合DNA液；
S9、片段筛选：对步骤H所得混合DNA液进行片段筛选及纯化，得到加接头的DNA片段；
S10、PCR扩增：对所得加接头的DNA片段进行PCR 扩增及纯化，得到小片段文库，具体可将加接头的DNA片段与Platinum PCR Super Mix High Fidelity与LibraryAmplification Primer Mix混合扩增纯化得到；
S11、文库检测：对小片段文库采用Qubit和Agilent Bioanalyzer2100检测浓度和片段大小；
S12、高通量测序：对检测合格的小片段文库进行高通量测序；也可采用捕获的方法实现。
S13、数据预处理：将高通量测序得到的数据首先经过低质量过滤，筛选出长度大于100bp的序列；
进行筛选过滤是因为设计的引物加上目标序列长度都大于100bp。当然，过滤的长度可根据设计的目标区域调整。
S14、序列比对：将预处理后的序列与人类基因组序列进行比对；
具体可预处理后的序列与人类基因组序列（hg19）用bowtie2进行比对。S15、序列筛选：根据比对的结果，筛选出unique read的序列：
具体筛选出错配碱基在4个以下，且只比对到一处的序列，作为unique read。序列筛选也可采用其他标准，如允许错配的碱基数目是长度的百分之3且只比对到一处；或者允许的错配碱基数目是序列上设计的SNP数目加上1且只比对到一处。
S16、SNP数据统计：根据比对结果及SNP位点的位置统计出SNP位点的覆盖深度、碱基种类、每种碱基对应数目；
具体可通过samtools软件进行统计。
S17、根据统计结果，对每个SNP位点中四种碱基的出现次数进行排序，然后计算出SNP位点的覆盖深度Depth，并根据排序结果求出SNP位点的测序错误率Error；根据每行的覆盖深度Depth和测序错误率Error，求出SNP位点的平均覆盖度和测序错误率的较大值ErrorUp；根据平均覆盖度或者用户自定的覆盖深度，对SNP位点进行过滤；或根据每个SNP位点的测序错误率Error是否大于ErrorUp来过滤；
具体来说，对每个SNP位点中四种碱基的出现次数排序后，按照从多到少的顺序分别称作Major、Minor、Third及Fourth。计算出覆盖深度Depth= Major+Minor+Third+Fourth，并利用下述公式求出测序错误率：Error=（Third+Fourth）/Depth。
测序错误率的较大值按以下公式计算：
ErrorUp= (ErrorMean+ErrorSD*2), Errormean（Errori）i=1…N，N为SNP位点个数，即Errormean为每一SNP位点的测序错误率平均值，ErrorSD是指每一个SNP位点测序错误率的方差，其中的参数2也可采用其他值，根据置信区间需要而定。
S18、计算每个SNP距离：对每一个SNP位点，计算出可能的胎儿游离DNA浓度Conc在四种SNP组合类型下所有SNP位点到相应胎儿游离DNA浓度的距离TypeDistance；
由于母亲和胎儿在一个SNP位点上可能出现四种情况，即具有四种SNP组合类型：类型1：母亲和胎儿都是纯合SNP；类型2：母亲纯合SNP，胎儿杂合SNP；类型3：母亲杂合SNP，胎儿纯合SNP；类型4：母亲和胎儿都是杂合SNP。对于每一个SNP位点，在不确定核型之前假设都存在上面四种可能。具体用Typej（j=1,2,3,4）来表示四种类型的SNP组合。
实际临床中可能的胎儿游离DNA浓度范围是0.025-0.5之间。具体用Conc来表示可能的胎儿游离DNA浓度，Conc的范围就是0.025-0.5。
按以下公式计算SNP距离Distance，Distance=abs(Major * A + Minor * B) /(A**2 + B**2)**0.5，其中A，B是参数，A**2是指A的2次方，Major与Minor在S17中提到，具体如下：
Type=1，A=Err/3,B=-(1-Err);
Type=2,A=Conc/2*(1-Err)+(1-Conc/2)*Err/3,B= -((1-Conc/2)*(1-Err)+Conc/2*Err/3)；
Type=3,A=(1/2-Conc/2)*(1-Err)+(1/2+Conc/2)*Err/3,B= -((1/2+Conc/2)*(1-Err)+(1/2-Conc/2)*Err/3)；
Type=4,A=1,B=-1；
这里Err=mean（Errori）i=1…1035，mean为求平均值，且Errori小于ErrorUp，Major、Minor分别为SNP位点中四种碱基出现次数排第一及第二的值。
每个SNP位点在四种SNP组合类型下对应有四个距离值，记作TypeDistance=Distance*C，C是一个参数，其在四种SNP类型下分别是：2（类型1），1（类型2），1（类型3），2（类型4）；根据计算得到的四种类型对应的距离判断该SNP位点所对应的类型；其中Typej（j=1,2,3,4）表示四种类型的SNP组合。
S19、计算出总距离：找出每一SNP位点在四种SNP组合类型下对应最小的TypeDistance，最后求和：AllDistance=sum（TypeDistancei），i=1…N，N为SNP位点总个数；
首先假设Conc是0.025，Err由实际数值计算得到，分别算出所有过滤后的SNP位点的四个距离值TypeDistance。找出四种类型对应下每个SNP位点最小的那个TypeDistance，最后求和：AllDistance=sum（TypeDistancei）i=1…1035且Errori小于ErrorUp，这就是0.025这个胎儿游离DNA浓度下的最小距离值。对于其他胎儿游离DNA浓度也是如此计算。
S20、选择最优的总距离：在Conc范围内计算出的AllDistance中选出其中最小的，其所对应的Conc即为实际的胎儿游离DNA浓度。
从Conc=0.025开始，每次增加胎儿游离DNA浓度可能值Conc=Conc+0.001，然后用上步同样的方法计算出AllDistance，在这475个AllDistance选出最小的，其所对应的Conc就是实际的胎儿游离DNA浓度。每次递增的大小可变。
具体实施例
本发明将DNA高通量测序法引入到孕妇游离DNA中胎儿游离DNA比例测定技术中，具体应用过程的流程结合具体实例作进一步描述。
1、材料
供试样本为未知胎儿游离DNA浓度的孕妇外周血4例，其中2例男胎，2例女胎，为体现本实施例的准确性，两例男胎都用现有的定量技术定量胎儿游离DNA浓度进行对比。
高通量测序法试剂主要由3个试剂盒组成，文库构建试剂盒（扩增文库构建试剂盒和连接文库构建试剂盒）；测序模板制备试剂盒以及测序试剂盒。
2、方法
将四个样本按照上述实验步骤经过建库和上机处理，得到样本数据，然后进行数据分析。
1）数据处理
对上步拿到的数据进行预处理，即过滤掉长度少于100bp，或者序列中测序质量在Q20的比例小于50%的序列，二者中任一条件满足就去掉。四个样本过滤情况如下：
| 原始reads数目 | 过滤后reads数目 | 有效reads比例 | |
| 样本1 | 3732137 | 3423150 | 91.7% |
| 样本2 | 3523141 | 3261449 | 92.6% |
| 样本3 | 3869887 | 3494950 | 90.3% |
| 样本4 | 6453416 | 6019736 | 93.2% |
2）将过滤后的序列用bowtie2比对到参考序列人类基因组HG19。
根据上步比对的结果，筛选出错配数目小于等于4，且比对到人类基因组一个位置的reads，四个样本过滤的情况如下：
| 质量过滤后reads数目 | 比对过滤后的reads数目 | 有效reads比例 | |
| 样本1 | 3423150 | 2385351 | 69.7% |
| 样本2 | 3261449 | 2295840 | 70.4% |
| 样本3 | 3494950 | 2150867 | 61.5% |
| 样本4 | 6019736 | 4113254 | 68.3% |
3）利用软件samtools的命令mpileup，统计出四个样本每个SNP的碱基类型，碱基比例及覆盖深度。
根据上步的统计结果，确定每个SNP的Major和Minor，即出现次数最多碱基是Major，出现次数第二多的是Minor，并根据S17计算测序错误，四个样本的测序错误分别如下：
| 样本1 | 0.19% |
| 样本2 | 0.24% |
| 样本3 | 0.14% |
| 样本4 | 0.18% |
4）将得到的上述数据，运用S18到S20的步骤计算出四个样本的胎儿游离DNA浓度如下：
| 样本编号 | 胎儿游离DNA浓度 | 男胎Y计算的胎儿游离DNA浓度 |
| 样本1 | 9.3% | 10% |
| 样本2 | 10.3% | 10.5% |
| 样本3 | 8.7% | / |
| 样本4 | 13.2% | / |
综上所述，本发明运用高通量测序的方法，结合生物信息学，统计学等技术，计算出游离DNA中胎儿游离DNA的比例。本发明不仅可以用来计算胎儿游离DNA的比例，同样适合混合血浆样本中混合比例的计算。相对于传统方法中采用一些胎儿与母亲甲基化不同的基因，或者男性特有基因等来计算胎儿游离DNA浓度，本发明采用SNP来计算胎儿游离DNA浓度，准确率高，且本发明有系统的矫正PCR引入误差的统计学算法，如动态优化算法，遍历法等，而传统的方法只是简单的通过甲基化基因的比例来判断胎儿游离DNA浓度（甲基化定量胎儿游离DNA浓度的方法），或者是用PCR定量曲线得到胎儿游离DNA浓度，没有矫正实验或者PCR引入误差的能力，同样提高了准确率。
应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (4)
1.一种孕妇外周血中胎儿游离DNA比例的定量方法，其特征在于，包括步骤：
A、设计位点：在人类基因组中，找出突变频率在0.4-0.6的SNP位点；
B、设计引物：根据设计的SNP位点，设计出覆盖SNP位点的SNP引物；
C、提取样品DNA；
D、扩增含有SNP位点的DNA片段：将样品DNA、多重PCR反应酶和SNP 引物混匀后，进行PCR反应得到PCR产物；
E、纯化扩增到的PCR产物；
F、末端修复：将PCR产物、末端修复酶与末端修复酶缓冲液混合，室温下孵育进行反应得到混合DNA液；
G、片段筛选：对步骤F所得混合DNA液进行片段筛选及纯化，得到平末端DNA片段；
H、接头连接：将所得的平末端DNA片段与连接酶、连接酶缓冲液、及指定接头混合后，室温下孵育进行反应得到混合DNA液；
I、片段筛选：对步骤H所得混合DNA液进行片段筛选及纯化，得到加接头的DNA片段；
J、PCR扩增：对所得加接头的DNA片段进行PCR 扩增及纯化，得到小片段文库；
K、文库检测：对小片段文库检测浓度和片段大小；
L、高通量测序：对检测合格的小片段文库进行高通量测序；
M、数据预处理：将高通量测序得到的数据首先经过低质量过滤，筛选出长度大于100bp的序列；
N、序列比对：将预处理后的序列与人类基因组序列进行比对；
O、序列筛选：根据比对的结果，筛选出unique read的序列；
P、SNP数据统计：根据比对结果及SNP位点的位置统计出SNP位点的覆盖深度、碱基种类、每种碱基对应数目；
Q、根据统计结果，对每个SNP位点中四种碱基的出现次数进行排序，然后计算出SNP位点的覆盖深度Depth，并根据排序结果求出SNP位点的测序错误率Error；根据每行的覆盖深度Depth和测序错误率Error，求出SNP位点的平均覆盖度和测序错误率的较大值ErrorUp；根据平均覆盖度或者用户自定的覆盖深度，对SNP位点进行过滤；或根据每个SNP位点的测序错误率Error是否大于ErrorUp来过滤；
R、计算每个SNP距离：对每一个SNP位点，计算出可能的胎儿游离DNA浓度Conc在四种SNP组合类型下所有SNP位点到相应胎儿游离DNA浓度的距离TypeDistance；
U、计算出总距离：找出每一SNP位点在四种SNP组合类型下对应最小的TypeDistance，最后求和：AllDistance=sum（TypeDistancei），i=1…N，N为SNP位点总个数；
V、选择最优的总距离：在Conc范围内计算出的AllDistance中选出其中最小的，其所对应的Conc即为实际的胎儿游离DNA浓度；
所述步骤Q中，ErrorUp= (ErrorMean+ErrorSD*2), Errormean（Errori）i=1…N，N为SNP位点个数；
所述步骤U中，四种SNP组合类型为：类型1：母亲和胎儿都是纯合SNP；类型2：母亲纯合SNP，胎儿杂合SNP；类型3：母亲杂合SNP，胎儿纯合SNP；类型4：母亲和胎儿都是杂合SNP；
所述步骤R具体包括：
按以下公式计算SNP距离Distance，Distance=abs(Major * A + Minor * B) / (A**2+ B**2)**0.5,其中A，B是参数，具体如下：
Type=1,A=Err/3,B=-(1-Err);
Type=2,A=Conc/2*(1-Err)+(1-Conc/2)*Err/3,B= -((1-Conc/2)*(1-Err)+Conc/2*Err/3)；
Type=3,A=(1/2-Conc/2)*(1-Err)+(1/2+Conc/2)*Err/3,B= -((1/2+Conc/2)*(1-Err)+(1/2-Conc/2)*Err/3)；
Type=4,A=1,B=-1，每个SNP位点在四种SNP类型下对应有四个距离值，记作TypeDistance=Distance*C，C是一个参数，其在四种SNP类型下分别是：2，1,1,2；根据计算得到的四种类型对应的距离判断该SNP位点所对应的类型；其中Typej（j=1,2,3,4）表示四种类型的SNP组合，Major、Minor分别为SNP位点中四种碱基出现次数排第一及第二的值，Err=mean（Errori）i=1…1035且Errori小于ErrorUp。
2.根据权利要求1所述的孕妇外周血中胎儿游离DNA比例的定量方法，其特征在于，所述步骤O中，筛选出错配碱基在4个以下，且只比对到一处的序列，作为unique read。
3.根据权利要求1所述的孕妇外周血中胎儿游离DNA比例的定量方法，其特征在于，所述步骤A中，在1-13、18及21这15条染色体上共选出1035个SNP位点。
4.根据权利要求1所述的孕妇外周血中胎儿游离DNA比例的定量方法，其特征在于，Conc的范围为0.025-0.5。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510225771.5A CN104846089B (zh) | 2015-05-06 | 2015-05-06 | 一种孕妇外周血中胎儿游离dna比例的定量方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510225771.5A CN104846089B (zh) | 2015-05-06 | 2015-05-06 | 一种孕妇外周血中胎儿游离dna比例的定量方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104846089A CN104846089A (zh) | 2015-08-19 |
| CN104846089B true CN104846089B (zh) | 2017-06-16 |
Family
ID=53846072
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510225771.5A Active CN104846089B (zh) | 2015-05-06 | 2015-05-06 | 一种孕妇外周血中胎儿游离dna比例的定量方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104846089B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105926043B (zh) * | 2016-04-19 | 2018-08-28 | 苏州贝康医疗器械有限公司 | 一种提高孕妇血浆游离dna测序文库中胎儿游离dna占比的方法 |
| CN106755486A (zh) * | 2017-01-19 | 2017-05-31 | 人和未来生物科技（长沙）有限公司 | 无创产前胎儿β 型地贫基因突变检测文库构建方法、检测方法和试剂盒 |
| CN106755485A (zh) * | 2017-01-19 | 2017-05-31 | 人和未来生物科技（长沙）有限公司 | 无创产前胎儿α 型地贫基因突变检测文库构建方法、检测方法和试剂盒 |
| CN107133491B (zh) * | 2017-03-08 | 2020-05-29 | 广州市达瑞生物技术股份有限公司 | 一种获取胎儿游离dna浓度的方法 |
| CN108220451B (zh) * | 2017-12-08 | 2020-10-27 | 北京科迅生物技术有限公司 | 胎儿游离核酸浓度的检测方法及试剂盒 |
| CN109979529B (zh) * | 2017-12-28 | 2021-01-08 | 北京安诺优达医学检验实验室有限公司 | Cnv检测装置 |
| CN109461473B (zh) * | 2018-09-30 | 2019-12-17 | 北京优迅医疗器械有限公司 | 胎儿游离dna浓度获取方法和装置 |
| CN110580934B (zh) * | 2019-07-19 | 2022-05-10 | 南方医科大学 | 一种基于外周血游离dna高通量测序的妊娠期相关疾病预测方法 |
| CN110993024B (zh) * | 2019-12-20 | 2023-08-22 | 北京科迅生物技术有限公司 | 建立胎儿浓度校正模型的方法及装置与胎儿浓度定量的方法及装置 |
| CN113584178A (zh) * | 2020-04-30 | 2021-11-02 | 深圳华大法医科技有限公司 | 一种无创亲子检测分析方法和装置 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2488358A (en) * | 2011-02-25 | 2012-08-29 | Univ Plymouth | Enrichment of foetal DNA in maternal plasma |
| CN103215350B (zh) * | 2013-03-26 | 2016-11-02 | 苏州贝康医疗器械有限公司 | 一种基于单核苷酸多态性位点的孕妇血浆中胎儿dna含量的测定方法 |
| CN104073554B (zh) * | 2014-04-28 | 2016-04-13 | 王弢 | 血浆特定片段游离dna定量检测试剂盒 |
-
2015
- 2015-05-06 CN CN201510225771.5A patent/CN104846089B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104846089A (zh) | 2015-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104846089B (zh) | 一种孕妇外周血中胎儿游离dna比例的定量方法 | |
| AU2021202149B2 (en) | Detecting repeat expansions with short read sequencing data | |
| CN104232777B (zh) | 同时确定胎儿核酸含量和染色体非整倍性的方法及装置 | |
| CN103874767B (zh) | 对核酸样本中预定区域进行基因分型的方法和系统 | |
| CN105441432B (zh) | 组合物及其在序列测定和变异检测中的用途 | |
| CN109767810B (zh) | 高通量测序数据分析方法及装置 | |
| WO2016011982A1 (zh) | 确定生物样本中游离核酸比例的方法、装置及其用途 | |
| CN104946773B (zh) | 一种利用snp进行产前亲权关系判定的方法 | |
| CN107949845B (zh) | 能够在多个下一代测序平台上区分胎儿性别和胎儿性染色体异常的计算机系统 | |
| CN106834502A (zh) | 一种基于基因捕获和二代测序技术的脊髓性肌萎缩症相关基因拷贝数检测试剂盒及方法 | |
| DE202013012824U1 (de) | Systeme zum Erfassen von seltenen Mutationen und einer Kopienzahlvariation | |
| CN105392894A (zh) | 确定样本基因组中是否存在拷贝数变异的方法、系统和计算机可读介质 | |
| CN113450871B (zh) | 基于低深度测序的鉴定样本同一性的方法 | |
| CN105844116A (zh) | 测序数据的处理方法和处理装置 | |
| US10294518B2 (en) | Methods and systems for ultra-sensitive detection of genomic alterations | |
| CN103114150A (zh) | 基于酶切建库测序与贝叶斯统计的单核苷酸多态性位点鉴定的方法 | |
| CN102753703A (zh) | 胎儿染色体非整倍性的检测方法 | |
| CN111091869A (zh) | 以snp为遗传标记物的亲缘关系鉴定方法 | |
| CN108595918A (zh) | 循环肿瘤dna重复序列的处理方法及装置 | |
| WO2018149114A1 (zh) | 确定胎儿染色体中微缺失微重复的方法及设备 | |
| CN109461473B (zh) | 胎儿游离dna浓度获取方法和装置 | |
| CN117524301B (zh) | 一种拷贝数变异的检测方法、装置以及计算机可读介质 | |
| CN108319817A (zh) | 循环肿瘤dna重复序列的处理方法及装置 | |
| CN109680054A (zh) | 一种低频dna突变的检测方法 | |
| CN114517223A (zh) | 一种用于筛选snp位点的方法及其应用 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20181012 Address after: 350401 A 6, 7 building, Taiwan business park, two Jin Road, North Town, Pingtan, Fujian. Patentee after: Pingtan Wan Ji medical laboratory Co.,Ltd. Address before: 361100 605, 605, 98 Xing Xing Road, Xiangan District, Xiamen, Fujian. Patentee before: XIAMEN VANGENES BIOTECHNOLOGY Co.,Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211118 Address after: 361000 room 605, qiangye building, torch high tech Zone (Xiang'an) Industrial Zone, Xiang'an District, Xiamen City, Fujian Province Patentee after: XIAMEN VANGENES BIOTECHNOLOGY Co.,Ltd. Address before: 350401 A 6, 7 building, Taiwan business park, two Jin Road, North Town, Pingtan, Fujian. Patentee before: Pingtan Wan Ji medical laboratory Co.,Ltd. |
|
| PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
| PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A quantitative method for determining the proportion of fetal free DNA in maternal peripheral blood Effective date of registration: 20230607 Granted publication date: 20170616 Pledgee: Xiamen Xiang'an District Branch of China Postal Savings Bank Co.,Ltd. Pledgor: XIAMEN VANGENES BIOTECHNOLOGY Co.,Ltd. Registration number: Y2023980042947 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right |
Granted publication date: 20170616 Pledgee: Xiamen Xiang'an District Branch of China Postal Savings Bank Co.,Ltd. Pledgor: XIAMEN VANGENES BIOTECHNOLOGY Co.,Ltd. Registration number: Y2023980042947 |
|
| PC01 | Cancellation of the registration of the contract for pledge of patent right |