用于单细胞全基因组DNA扩增产物质量评价的核酸组合物及方法与流程

用于单细胞全基因组dna扩增产物质量评价的核酸组合物及方法
技术领域
1.本发明涉及分子生物学领域，特别是涉及一种用于单细胞全基因组dna扩增产物质量评价的核酸组合物及方法。


背景技术：

2.胚胎植入前遗传学诊断(preimplantation genetic diagnosis，pgd)是在体外对胚胎进行遗传学检查，确定胚胎是否携带有遗传缺陷的基因，排除遗传异常胚胎，提高植入成功率，并避免基因缺陷遗传给子代。目前pgd主要集中在胚胎活检，但是由于其是侵入性诊断，可能会对胚胎造成损伤，因而对于新发患者或无先证者家系，采用精子或极体进行pgd检测前预检是一种有效的方式，既可以确定基因变异信息，又可以减少对胚胎造成的损伤。
3.为了从单精子或极体中最大限度地获取遗传信息，目前利用单细胞全基因组扩增技术(whole genome amplification，wga)对单精子或极体基因组进行扩增。但由于受到多因素干扰，无法保证细胞全基因组扩增达到100％的成功率，而质量较差的扩增产物会直接影响全面测序结果，降低pgd检测前预检的成功率。因此需要通过评价单细胞全基因组dna的扩增产物的质量，剔除不合格的样本，避免不合格样本进入下一步测序阶段，减少资源浪费。但是目前评价单细胞全基因组dna的扩增产物的质量的方法的准确度还比较低。


技术实现要素：

4.本发明研究发现，lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1这八个基因可以作为检测单细胞全基因组dna的扩增产物的质量的靶基因，利用扩增这八个基因的引物对来评价单细胞全基因组dna扩增产物的质量，可以提高评价质量的准确率。基于此，有必要提供一种lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1在评价单细胞全基因组dna扩增产物质量中的应用。
5.此外，还提供了一种可以提高评价单细胞全基因组dna扩增产物质量的准确率的核酸组合物，一种准确率高的单细胞全基因组dna扩增产物质量评价方法和检测试剂盒。
6.其中一个实施例中，单细胞为生殖细胞或极体。
7.一种核酸组合物，包括用于扩增lhcgr的引物对、用于扩增trex1的引物对、用于扩增evc的引物对、用于扩增musk的引物对、用于扩增abcc8的引物对、用于扩增daoa的引物对、用于扩增aldoa的引物对和用于扩增slc25a1的引物对中的至少六种。
8.其中一个实施例中，用于扩增lhcgr的引物对的核苷酸序列如seq id no.1～2所示，用于扩增trex1的引物对的核苷酸序列如seq id no.3～4所示，用于扩增evc的引物对的核苷酸序列如seq id no.5～6所示，用于扩增musk的引物对的核苷酸序列如seq id no.7～8所示，用于扩增abcc8的引物对的核苷酸序列如seq id no.9～10所示，用于扩增daoa的引物对的核苷酸序列如seq id no.11～12所示，用于扩增aldoa的引物对的核苷酸
序列如seq id no.13～14所示，用于扩增slc25a1的引物对的核苷酸序列如seq id no.15～16所示。
9.一种单细胞全基因组dna扩增产物质量的评价方法，包括以下步骤：
10.采用上述核酸组合物，利用实时荧光定量pcr扩增单细胞全基因组dna的扩增产物，并根据扩增结果判断单细胞全基因组dna的扩增产物是否合格；其中，在核酸组合物中的至少六个扩增引物对有对应的扩增曲线且扩增曲线对应的ct值小于等于30，及扩增引物对应有单一的熔解曲线峰型且熔解曲线峰型与阳性对照一致时，单细胞的全基因组dna的扩增产物合格。
11.其中一个实施例中，单细胞全基因组dna的扩增产物是单细胞基因组dna经多重置换扩增(multiple displacement amplification，mda)而得到的扩增产物。
12.其中一个实施例中，单细胞为精子或极体。
13.一种单细胞全基因组dna的扩增产物的质量的评价方法，包括以下步骤：
14.采用上述核酸组合物，扩增单细胞全基因组dna扩增产物，得到待电泳的扩增产物；及将待电泳的扩增产物进行琼脂糖凝胶电泳，根据电泳结果判断所述单细胞全基因组dna的扩增产物是否合格；其中，在核酸组合物中的至少六个扩增产物引物对有对应电泳条带时，单细胞的全基因组dna扩增产物合格。
15.一种用于评价单个细胞全基因组dna扩增产物质量的检测试剂盒，包括上述核酸组合物。
16.一种用于单细胞基因组dna测序的试剂盒，包括上述核酸组合物和辅助试剂，辅助试剂包括单细胞全基因组dna提取试剂和单细胞全基因dna扩增试剂中的至少一种。
附图说明
17.图1～8分别为极体1全基因组dna扩增产物、阳性对照和阴性对照所对应的lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1引物的qpcr扩增曲线；
18.图9～16分别为极体1全基因组dna扩增产物和阳性对照所对应的lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1引物的熔解曲线。
具体实施方式
19.为使本发明的上述目的、特征和优点能够更加明显易懂，下面对本发明的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明。但是本发明能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似改进，因此本发明不受下面公开的具体实施例的限制。
20.除非另有定义，本文所使用的所有的技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中在本发明的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本发明。本文所使用的术语“和/或”包括一个或多个相关的所列项目的任意的和所有的组合。
21.本发明研究发现，lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1这八个基因可以作为检测单细胞全基因组dna扩增产物质量的靶基因，利用扩增这八个基因的引物对来评价单细胞全基因组dna扩增产物的质量，可以提高评价质量的准确率。基于此，本申
请一实施方式提供一种lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1在评价单细胞的全基因组dna的扩增产物的质量中的应用。
22.具体地，lhcgr是位于人源2号染色体上的保守区段，trex1是位于人源3号染色体上的保守区段，evc是位于人源4号染色体上的保守区段，musk是位于人源9号染色体上的保守区段，abcc8是位于人源11号染色体上的保守区段，daoa是位于人源13号染色体上的保守区段，aldoa是位于人源16号染色体上的保守区段，slc25a1是位于人源22号染色体上的保守区段。
23.在其中一些实施例中，单细胞为生殖细胞或极体。
24.此外，本技术一实施方式还提供了一种核酸组合物，该核酸组合物包括用于扩增lhcgr的引物对、用于扩增rex1的引物对、用于扩增evc的引物对、用于扩增musk的引物对、用于扩增abcc8的引物对、用于扩增daoa的引物对、用于扩增aldoa的引物对和用于扩增slc25a1的引物对中的至少六种。
25.具体地，扩增lhcgr的引物对扩增的是lhcgr的部分区段，该区段的具体位置为2号染色体上的48956152～48956555；扩增rex1的引物对扩增的是rex1的部分区段，该区段的具体位置为3号染色体上的48508102～48508845；扩增evc的引物对扩增的是evc的部分区段，该区段的具体位置为4号染色体上的5755303～5755786；扩增musk的引物对扩增的是musk的部分区段，该区段的具体位置为9号染色体上的113449276～113449475；扩增abcc8的引物对扩增的是abcc8的部分区段，该区段的具体位置为11号染色体上的17417090～17417273；扩增daoa的引物对扩增的是daoa的部分区段，该区段的具体位置为13号染色体上的106113890～06114214；扩增aldoa的引物对扩增的是aldoa的部分区段，该区段的具体位置为16号染色体上的30060465～30060617；扩增slc25a1的引物对扩增的是slc25a1的部分区段，该区段的具体位置为22号染色体上的19164197～19164714。
26.在其中一个实施例中，用于扩增lhcgr的引物对的核苷酸序列如seq id no.1～2所示。具体地，seq id no.1的核苷酸序列为：5
’‑
tgcccatcatggctcctatt
‑3’
；seq id no.2的核苷酸序列为：5
’‑
gcctgcccagacctagtaat
‑3’
。可以理解的是，在其他实施例中，用于扩增lhcgr的引物对不限于上述，还可以是其他可以扩增lhcgr的引物对。
27.在其中一个实施例中，用于扩增rex1的引物对的核苷酸序列如seq id no.3～4所示。具体地，seq id no.3的核苷酸序列为：5
’‑
tttcgacatggaccggactg
‑3’
；seq id no.4的核苷酸序列为：5
’‑
aggctgggactagtgttcct
‑3’
。可以理解的是，在其他实施例中，用于扩增rex1的引物对不限于上述，还可以是其他可以扩增rex1的引物对。
28.在其中一个实施例中，用于扩增evc的引物对的核苷酸序列如seq id no.5～6所示。具体地，seq id no.5的核苷酸序列为：5
’‑
tctctcaagacgtggaggct
‑3’
；seq id no.6的核苷酸序列为：5
’‑
ctgacacatgggaagtaggca
‑3’
。可以理解的是，在其他实施例中，用于扩增evc的引物对不限于上述，还可以是其他可以扩增evc的引物对。
29.在其中一个实施例中，用于扩增musk的引物对的核苷酸序列如seq id no.7～8所示。具体地，seq id no.7的核苷酸序列为：5
’‑
tcaagtttcctctcacctttattc
‑3’
；seq id no.8的核苷酸序列为：5
’‑
gccatcatcactgtcttccac
‑3’
。可以理解的是，在其他实施例中，用于扩增musk的引物对不限于上述，还可以是其他可以扩增musk的引物对。
30.在其中一个实施例中，用于扩增abcc8的引物对的核苷酸序列如seq id no.9～10
所示。具体地，seq id no.9的核苷酸序列为：5
’‑
gtgctccagatctgataaggct
‑3’
；seq id no.10的核苷酸序列为：5
’‑
ccatgcagggcacatcatc
‑3’
。可以理解的是，在其他实施例中，用于扩增abcc8的引物对不限于上述，还可以是其他可以扩增abcc8的引物对。
31.在其中一个实施例中，用于扩增daoa的引物对的核苷酸序列如seq id no.11～12所示。具体地，seq id no.11的核苷酸序列为：5
’‑
tgtctcttctaaatgtcccattcct
‑3’
；seq id no.12的核苷酸序列为：5
’‑
ctgattctctctgtaggtggg
‑3’
。可以理解的是，在其他实施例中，用于扩增daoa的引物对不限于上述，还可以是其他可以扩增daoa的引物对。
32.在其中一个实施例中，用于扩增aldoa的引物对的核苷酸序列如seq id no.13～14所示。具体地，seq id no.13的核苷酸序列为：5
’‑
gcatagtgagtcctgtgtt
‑3’
；seq id no.14的核苷酸序列为：5
’‑
cctattcggacttcattcaa
‑3’
。可以理解的是，在其他实施例中，用于扩增aldoa的引物对不限于上述，还可以是其他可以扩增aldoa的引物对。
33.在其中一个实施例中，用于扩增slc25a1的引物对的核苷酸序列如seq id no.15～16所示。具体地，seq id no.15的核苷酸序列为：5
’‑
cccctgttgggtggatatagag
‑3’
；seq id no.16的核苷酸序列为：5
’‑
aagttcatccacgaccagacc
‑3’
。可以理解的是，在其他实施例中，用于扩增slc25a1的引物对不限于上述，还可以是其他可以扩增slc25a1的引物对。
34.在其中一个实施例中，核酸组合物包括用于扩增lhcgr的引物对、用于扩增rex1的引物对、用于扩增evc的引物对、用于扩增musk的引物对、用于扩增abcc8的引物对、用于扩增daoa的引物对、用于扩增aldoa的引物对和用于扩增slc25a1的引物对；用于扩增lhcgr的引物对的核苷酸序列如seq id no.1～2所示，用于扩增trex1的引物对的核苷酸序列如seq id no.3～4所示，所述用于扩增evc的引物对的核苷酸序列如seq id no.5～6所示，所述用于扩增musk的引物对的核苷酸序列如seq id no.7～8所示，所述用于扩增abcc8的引物对的核苷酸序列如seq id no.9～10所示，所述用于扩增daoa的引物对的核苷酸序列如seq id no.11～12所示，所述用于扩增aldoa的引物对的核苷酸序列如seq id no.13～14所示，所述用于扩增slc25a1的引物对的核苷酸序列如seq id no.15～16所示。通过将上述核酸组合物设置为这八对引物对可以提高评价单细胞全基因组dna扩增产物质量的准确率。
35.上述核酸组合物包括用于扩增lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1的引物对。上述核酸组合物可以提高评价单细胞全基因组dna扩增产物质量的准确率。
36.此外，本技术一实施方式还提供一种单细胞全基因组dna扩增产物质量的评价方法，该评价方法包括以下步骤：采用上述核酸组合物，利用实时荧光定量pcr扩增单细胞全基因组dna的扩增产物，并根据扩增结果判断单细胞全基因组dna的扩增产物是否合格；其中，核酸组合物中至少六对扩增引物对有对应的扩增曲线且扩增曲线对应的ct值小于等于30，及扩增引物对应有单一的熔解曲线峰型且熔解曲线峰型与阳性对照一致时，单细胞的全基因组dna的扩增产物合格。
37.在其中一个实施例中，单细胞全基因组dna扩增产物是单细胞基因组dna经多重置换扩增(multiple displacement amplification，mda)而得到的扩增产物。
38.在其中一个实施例中，单细胞为精子或极体。
39.上述单细胞全基因组dna扩增产物质量的评价方法，是通过qpcr进行，通过对qpcr结果中扩增曲线的监测，能够直接反映该引物对应的靶基因在qpcr扩增过程中是否有被扩
增；通过对ct值监测，能够判断样本中该引物对应的靶基因扩增产物的浓度是否正常；通过对熔解曲线的监测，能够判断该引物的扩增产物是否是唯一及扩增过程有无异常。这三个指标结合能够反映该引物对应的靶基因在qpcr中的扩增质量如何，从而能推断出原样本的扩增产物质量如何。质量合格的样本才能进入下一步进行文库构建和测序，进行遗传学分析。
40.此外，本技术一实施方式还提供了另一种单细胞全基因组dna扩增产物质量的评价方法，该方法包括以下步骤：采用上述核酸组合物，扩增单细胞全基因组dna的扩增产物，得到待电泳检测的扩增产物；及将待电泳的扩增产物进行琼脂糖凝胶电泳，根据电泳结果判断所述单细胞全基因组dna扩增产物是否合格；其中，在核酸组合物中的至少六个扩增产物引物对有对应电泳条带时，单细胞的全基因组dna的扩增产物合格。
41.上述单细胞全基因组dna扩增产物质量的评价方法，是通过普通pcr对单细胞全基因组dna的扩增产物进行扩增，在通过对pcr扩增产物进行电泳，观察引物对应的靶基因的电泳条带的有无可知该靶基因在原扩增产物中是否被覆盖到，反映出全基因组dna的扩增产物质量。
42.此外，本技术一实施方式还提供一种单细胞全基因组dna测序方法，该方法包括以下步骤：裂解单细胞，获得单细胞全基因组dna；以单细胞全基因组dna为模板扩增，获得单细胞全基因组dna的扩增产物；采用上述单细胞全基因组dna扩增产物质量的评价方法评价单细胞全基因组dna扩增产物的质量，获得合格的单细胞全基因组dna的扩增产物；将合格的单细胞全基因组dna的扩增产物构建文库并测序。
43.此外，本技术一实施方式还提供一种用于评价单细胞全基因组dna扩增产物质量的检测试剂盒，该检测试剂盒包括上述核酸组合物。
44.在其中一个实施例中，单细胞为生殖细胞或极体。该检测试剂盒可以应用于胚胎植入前遗传学诊断，用于减少测序成本。
45.上述用于评价单细胞全基因组dna扩增产物质量的检测试剂盒，利用了lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1在评价单细胞全基因组dna的扩增产物的质量中的应用，能够实现对单个细胞全基因组dna扩增产物的质量的评价，并且准确率高。
46.此外，本技术一实施方式还提供一种用于单细胞基因组dna测序的试剂盒，包括上述核酸组合物和辅助试剂，辅助试剂包括单细胞全基因组dna提取试剂和单细胞全基因dna扩增试剂中的至少一种。
47.上述用于单细胞基因组dna测序的试剂盒包括上述核酸组合物，具有上述核酸组合物相应的优点。
48.具体实施例
49.以下结合具体实施例进行详细说明。以下实施例如未特殊说明，则不包括除不可避免的杂质外的其他组分。实施例中采用试剂和仪器如非特别说明，均为本领域常规选择。实施例中未注明具体条件的实验方法，按照常规条件，例如文献、书本中所述的条件或者生产厂家推荐的方法实现。
50.实施例1
51.1.取样
52.采用由合作医院提供的精液和废弃极体。将采集到的来自同一人的精液，采用非
连续密度梯度离心法进行洗涤，洗涤后的精子转移至含有1ml精子冲洗培养液(coopersurgical，产品型号：art
‑
1005)的15ml离心管中。显微操作挑取单个精子、分取极体，分别命名为单精子1～12及极体1～2，每个样本分别放至2μl含有pbs缓冲液的pcr管底部，立即于
‑
80℃保存，避免反复冻融。
53.2.采用qiagen repli
‑
g single cell kit对单个精子和极体进行全基因组扩增，具体方法如下：
54.(1)裂解
55.1)准备dlb buffer：首次使用时加入500μl的无核酸酶水，充分混合并短暂离心。
56.2)配制buffer d1：33μl dlb buffer加上3μl dtt；该体积可满足12次反应，若一次实验未完全用完，可于
‑
20℃保存3个月。
57.3)配制buffer d2：用试剂盒提供的pbs将buffer d1体积补足为4μl，得到buffer d2。
58.4)向在上述取样步骤中得到的有单个精子或极体的pcr管中加入3μl buffer d2，小心轻弹管壁混匀，瞬时离心。
59.5)上述步骤4)中的pcr管置于65℃孵育10min(pcr仪热盖使用105℃)。
60.6)向上述孵育完成的样品中加入3μl stop solution，小心轻弹管壁混匀，瞬时离心，置于冰盒上。
61.(2)等温扩增
62.1)取出repli
‑
g sc reaction buffer、无核酸酶水置于室温，融化后涡旋混匀，瞬时离心，取出repli
‑
g sc dna polymerase置于冰盒上。
63.2)根据表1配制反应混合液，当加入无核酸酶水和repli
‑
g sc reaction buffer后，涡旋混匀，瞬时离心，然后再加入repli
‑
g sc dna polymerase。配制好的反应混合液放置冰上并立即使用。
64.表1
65.组分体积(μl)无核酸酶水9repli
‑
g sc reaction buffer29repli
‑
g sc dna polymerase2总体积40
66.3)将上述混合好的40μl反应混合液加入到经过裂解变性的样本管中。
67.4)设置pcr仪反应程序，按照表2条件进行孵育，得到单精子和极体全基因组dna的扩增产物：
68.表2
[0069][0070]
3.采用实时荧光定量pcr(qpcr)对上述步骤2中得到的单精子和极体全基因组dna的扩增产物分别进行扩增，扩增所需酶、buffer和染料选用kapa sybr fast qpcr kits
(roche)试剂盒，具体步骤如下：
[0071]
(1)将单精子或极体扩增产物用无核酸酶水稀释10倍，作为qpcr反应模板；选用外周血基因组dna作为阳性对照；将无核酸酶水作为阴性对照。
[0072]
(2)按照表3配制qpcr反应体系，引物组中的每对引物单独为一个扩增体系，引物对序列如表4所示：
[0073]
表3
[0074]
试剂名称使用量kapa sybr mix 2
×
10μl模板dna2μlf引物3μlr引物3μlhigh rox0.4μlddh2o1.6μl总体积20μl
[0075]
表4
[0076]
[0077][0078]
(3)按照表5参数进行qpcr扩增：
[0079]
表5
[0080][0081]
4.结果及判断：
[0082]
qpcr结果如图1～图16及表6所示，其中：
[0083]
图1～图8所示为极体1全基因组dna扩增产物的qpcr扩增曲线结果。在图1～图8中，极体1全基因组dna的扩增产物的扩增曲线对应于标号2的曲线，阳性对照的扩增曲线对应于标号为1的曲线，阴性对照的扩增曲线对应于标号为3的曲线。由图1～图8可知，极体1全基因组dna的扩增产物所对应的八对引物均有扩增曲线。
[0084]
图9～图16为极体1全基因组dna扩增产物的qpcr熔解曲线结果。在图9～图16中，极体1全基因组dna的熔解曲线对应标号2的曲线，阳性对照的熔解曲线对应于标号为1的曲线，阴性对照的熔解曲线对应于标号为3的曲线。图9～图16所示的熔解曲线表明，八对引物对所对应的极体1扩增产物的熔解曲线均是单一峰型且与阳性对照一致。
[0085]
表6为单精子和极体全基因组dna扩增产物对应的各引物对qpcr的ct值及tm值统计结果。由表6可知，极体1扩增产物所对应的八对引物的ct值均小于30。
[0086]
综上，极体1全基因组dna的扩增产物合格。
[0087]
表6
[0088][0089][0090]
根据表6对其他样本均根据此方法对其扩增产物质量进行判断，即核酸组合物中至少六对扩增引物对有对应的扩增曲线且扩增曲线对应的ct值小于等于30及扩增引物对应有单一的熔解曲线峰型且熔解曲线峰型与阳性对照一致，则判断该样本的扩增产物质量合格。结果表明：在本次实验中，用以试验的单精子1～12和极体1～2中，其中，单精子1、单精子3、单精子4、单精子5、单精子6、单精子7、单精子8、单精子9、单精子10、极体1和极体2扩
增产物质量合格。
[0091]
5.测序及单体型构建：
[0092]
(1)dna文库构建及测序：采用ion ampliseq
tm library kit 2.0对合格样本的扩增产物进行多重pcr文库构建，所构建文库采用da8600进行上机测序，测序质控数据如表7所示。其中，覆盖度(coverage)说明可捕获到的snp位点比例，frac_100x表示测序深度大于100x的snp位点比例。
[0093]
表7
[0094][0095][0096]
(2)构建单体型：对以上测序完成的单精子进行单体型构建，通过变异位点上下游1m范围内snp位点连锁分析，共检测到56个用于构建男女双方单体型的snp位点，以上扩增产物质量评价合格的单精子均成功构建单体型。所构建的单体型如表8所示：
[0097]
表8
[0098]
[0099][0100]
注：“miss”代表该位点漏检，“cc”、“gt”和“cta”是由于碱基的插入或缺失造成的，即等位基因上碱基的一对多或多对一。
[0101]
对比例1
[0102]
本对比例与实施例1的步骤相同，采用的样本有所不同，其不同在于：步骤5中对步骤4中不合格的样本(单精子2、单精子11和单精子12)进行测序和单体型构建，测序质控结果如表9所示，单体型构建结果如表10所示。
[0103]
表9
[0104][0105]
表10
[0106]
[0107]
[0108][0109]
注：“miss”代表该位点漏检，“cc”、“gt”和“cta”是由于碱基的插入或缺失造成的，即等位基因上碱基的一对多或多对一。
[0110]
对比例2
[0111]
本对比例与实施例1的步骤大致相同，其不同在于：无实施例1中的步骤3和4；对比例1中样本dna扩增产物测序质控结果如表11所示：
[0112]
表11
[0113]
样本名称raw_readsmapped_readsmean_dpcoveragecapture_ratefrac_30xfrac_100x男方3208332319514323493.5310.731410.9997女方3185364317173023128.0810.731210.9995单精子21307174330573847712.610.16890.29090.06230.0392单精子223234180321651716287.360.29940.53680.1550.1097单精子233985407396383625660.620.4270.64890.26460.2449单精子242941109291927718185.970.16330.64170.08690.0626单精子252798494278814012757.160.16520.47070.06720.0597单精子262820574281052418389.610.20350.65310.10220.0724单精子27295917929487818216.120.31950.28480.12520.0768单精子283173426315163710702.050.22670.40790.11490.083单精子292431629242726715061.160.19410.5960.06410.0383单精子302630311262249114073.670.54390.55610.26340.2102单精子31226125022505208504.280.56520.41340.09410.0728单精子32307174330573847712.610.16890.29090.06230.0392单精子342445548243747716967.671.00000.69660.96780.9419单精子352485999247710116108.690.92830.66250.85010.7721单精子362343328233679713551.460.84600.59630.61110.4650
[0114]
对比例3
[0115]
本对比例与实施例1的步骤大致相同，其不同在于：无实施例1中的步骤3和4；对比例2中样本dna扩增产物测序质控结果如表12所示：
[0116]
表12
[0117]
[0118][0119]
从表7和表8中可知，实施例1中通过本发明的方法进行了全基因组dna的扩增产物的质量评价后筛选出的合格样本，测序深度大于100x的snp位点比例平均值为0.87，coverage的平均值为0.97，一般情况下能满足构建单体型的要求。结合表8的单体型构建结果可知，筛选出的合格的单精子构建的单体型，均满足致病位点上下游附近各至少需要2个有效snp位点的要求，可用于后续的遗传连锁分析检测。此外，不同的单精子所构建的单体型一致，相互验证了构建的单体型的准确性。经过本发明的方法筛选出的合格样本单体型构建成功率为100％。
[0120]
从表9和表10中可知，对比例1中，通过本发明的方法筛选出的不合格的样本，测序深度大于100x的snp位点比例平均值为0.25，coverage的平均值为0.61，与男女双方的gdna的frac_100x值和coverage值相差太远，结合表10可知，用这三个不合格的样本构建出的单体型无法满足致病位点上下游附近各至少需要2个有效snp位点的要求，不可用于后续的遗传连锁分析检测。经过本发明的方法筛选出的不合格样本单体型构建成功率为0。
[0121]
上述结果说明lhcgr、trex1、evc、musk、abcc8、daoa、aldoa和slc25a1在评价单细胞全基因组dna的扩增产物的质量中的应用的准确度很高。
[0122]
此外，从表11中可知，对比例2中，未经本发明的方法筛选的样本中测序深度大于100x的snp位点比例平均值只有0.38，coverage平均值只有0.41，相比较于男女双方的gdna的frac_100x值和coverage值严重偏低。一般难以用于单体型构建，即使能构建出单体型也难以满足后续遗传连锁分析，各样本之间构建出的单体型也无法相互验证。
[0123]
从表12中可知，对比例3中，未经本发明的方法筛选的样本中测序深度大于100x的snp位点比例平均值只有0.54，coverage平均值只有0.85，不是所有样本都能满足构建单体型的要求，且有些样本frac_100x值和coverage值都很低，如单精子52、单精子53和单精子54，与男女双方的gdna的frac_100x值和coverage值相差太大，难以满足后续遗传连锁分析。
[0124]
显然，如果对全基因组dna的扩增产物的质量不经过筛选即进行测序及后续遗传
分析，会大量浪费测序及遗传分析的资源，浪费时间。
[0125]
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0126]
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准，说明书及附图可以用于解释权利要求的内容。