一种两核苷酸不同步合成测序分析pcr产物单体型方法
【技术领域】
[0001] 本发明属于生物技术领域，是一种PCR产物单体型分析方法，具体涉及一种两核 苷酸不同步合成测序分析PCR产物单体型方法。
【背景技术】
[0002] 人类基因组单体型图计划本身说明了单体型研宄的重要性，也为个体样本的单体 型分析提供了重要的数据库。位于一条染色体特定区域的一组相互关联并倾向于以整体遗 传给后代的SNP的组合叫作单体型（haplotype)。早期的研宄工作就表明，单体型的分析 可以极大地减少疾病关联研宄的工作量，且提供比SNP基因型分析更大的信息量，有利于 发现更多的信息。理论上说，一个包含N个独立双等位基因的SNP区域中可能存在2n种单 体型，但实际上，由于缺乏重组和（或）反复突变，SNP之间常常存在连锁不平衡的关系，因 此，单体型的种类不会超过（N+1)。对于二倍体个体（如单个人）的DNA样本分析，一般情 况下就是要从大样本所包含的可能单体型中找出具体属于哪一种或者二种单体型。如果 两个相邻SNP位点间有一个SNP位点为纯合型，那么其单体型就很容易判断。因此，单体型 分析就聚焦在两个相邻杂合型SNP间二种单体型的确定。目前对于单体型的分析，普遍采 用间接方法是利用多个SNP分型结果，以统计方法推断，这种方法得到广泛的应用，且用于 构建单体型的相关软件也非常之多。然而，软件分析的结果对于临床样本而言，如果得不 到直接的实验证实终宄是不踏实的，临床检测也很难接受这样的分析结果。通过实验方法 得到自然状态下真实的单体型情况需要进行单分子测序（如克隆测序、高通量测序等），这 些方法由于成本高，不适于大规模的个体样本分析。对于常规PCR产物中包含多个SNP位 点的实验分析单体型方法，均采用一种等位基因特异引物-PCR (AS-PCR)的方法，以第一个 SNP位点为扩增起始点，用两种特异引物将起始为两种的不同DNA模板分别进行扩增，然后 分别进行测序，以确定其它SNP位点的碱基信息。然而，这种方法既无法知道单体型的含量 (需要另外设计实验来确定），也不一定适合所有SNP位点的扩增，且操作繁琐、效率不高。 我们知道利用测序引物在不同模板上发生不同步合成反应，焦测序技术不仅能够分型SNP 位点，同时也能够对其模板含量进行定量分析。沿着这一思路，如果在一个包含N个SNP位 点的PCR产物中，这种不同步合成测序继续对各个SNP位点的碱基信息进行读取，则除了在 SNP位置提供不同DNA模板的信息外，其它位置由于SNP碱基的错位测序，也可以提供不同 DNA模板的加和信息，并从这些加和信息中构建方程组实现对不同DNA模板的含量分析。
[0003] 最近，我们实验室提出一种基于两核苷酸实时合成测序的方法，该方法通过每次 反应同时加入未标记的两种dNTPs实时测序，得到两组编码，解码这两组编码便能确定测 序片段的具体碱基信息（肖鹏峰等，中国发明专利：ZL 201210128597. 6)。该方法具有大幅 提高测序长度，且测序得到的测序信号强度比传单核苷酸合成测序的信号强等特点。然而， 该发明专利对一个DNA模板样本需要两次测序既增加了操作的琐碎程度，也增加了分析的 成本，且同步测序反应的信息不能满足对单体型的分析。
【发明内容】
[0004] 发明目的：为了克服现有技术中存在的不足，本发明提供一种两核苷酸不同步合 成测序分析PCR产物单体型的方法，为PCR产物单体型方法提供一种快速、高效、灵敏的分 析方法。
[0005] 技术方案：为实现上述目的，本发明采用的技术方案为：
[0006] -种两核苷酸不同步合成测序分析PCR产物单体型方法，步骤包括：1)DNA提取； 2) PCR扩增；3)测序；4)关联分析；
[0007] 所述步骤3)的测序方法为：同时加入两个不同的核苷酸，PCR扩增产物中不同DNA 模板的合成测序不同步，得到相邻两个SNP位点上不同的测序信息。
[0008] 进一步的，在本发明中，所述步骤3)的具体方法包括以下步骤：
[0009] 3-1)制备单链DNA模板：将PCR扩增产物与亲和素修饰的磁珠反应，修饰生物素 的DNA链固定到所述磁珠上，在0. 2MNaOH溶液下变性；将未固定的另一条DNA链清除；然 后用洗液洗涤，得到固定的单链DNA模板；
[0010] 3-2)测序引物杂交：将测序引物与所述单链DNA模板，在杂交反应体系中，80°C下 放置5min，自然冷却至室温，完成杂交；
[0011] 3-3)焦测序：配备焦测序反应体系，与所述步骤3-2)得到的杂交产物混合，按照 循环加入两核苷酸方式对DNA模板进行循环焦测序反应，得到由按照先后顺序排列的单个 测序反应的碱基片段编码信息。
[0012] 进一步的，在本发明中，步骤3-2)中，所述测序引物是一段与单链DNA模板完全互 补的一段序列；所述测序引物为5' -SEQ ID N0:l-3'所示：5' -ctaaaacgag gaagtattac atc-3'；
[0013]所述杂交反应体系的条件为：10mM Tris-HC1，2M NaCl，lmM EDTA，0. 1% Tween20? pH7. 6〇
[0014] 进一步的，在本发明中，步骤3-3)中，焦测序反应体系包括pH 7. 7的0.1M Tris-Ac,2mM EDTA,10mM Mg(Ac)2,0.2 %BSA,10mM DTT,10mM APS,0. 4mg/mL PVP,4mM D-luciferin,2U/mL ATP sulfurylase,0. 4mM luciferase，2U/mL apyrase VII，2U/mL DNA 聚合酶I。
[0015] 进一步的，在本发明中，所述循环加入两核苷酸方式包括（A+GV(C+T)，（A+C/ (G+T)，（A+TV(C+G)三组中的一组或一组以上的循环测序方式；两核苷酸加入的顺序按照 循环的方式；或按照以四种不同的两核苷酸为原料，每次加入一种两核苷酸来实施。
[0016] 进一步的，在本发明中，所述步骤3)中，每个测序反应获得的测序信息包括两核 苷酸种类信息和测序信号强度信息，所述测序信号强度与合成核苷酸数目成正比；
[0017] 所述步骤4)中，根据不同步合成测序反应所得到的两个不同SNP位点上的测序信 息进行关联分析，构建关联方程式，确定PCR产物单体型及含量；
[0018] 所述关联方程式为：对于包含任意相邻两个SNP位点的PCR产物，不超过四种不 同的DNA模板；设定这四种单体型1、…、i的含量分别为Xl、…、Xi，其中2彡i彡4,且 x:+……+Xi= 1 ;则PCR产物中混合DNA模板单个测序反应测序信号强度fn、fn'与Xl、…、 Xi存在下列关系：（f n_fn'）= …+a:n xi;其中n表示第n次测序反应；fn、fn'分别表 示第n次测序反应的测序信号强度与背景值；、…、at分别表示100%单体型DNA序列 1、…、i在第n次测序反应中合成的核苷酸数目。
[0019] 进一步的，在本发明中，所述SNP位点上的测序信息进行关联分析是指两个相邻 SNP位点均是杂合型的单个样本的关联分析；
[0020] 若PCR产物中包含若干个SNP位点的分析，将其分解成相邻两个SNP位点的单体 型关联分析。
[0021] 进一步的，在本发明中，所述步骤2)中，
[0022] 扩增条件为：94°C起始变性5min ;35个热循环为：94°C变性30s、61°C退火45s、 72°C延伸45s ;最后72°C延伸7min ;
[0023] PCR扩增所用的PCR-条引物为5'端被生物素、氨基，或者丙烯酰胺等基团修饰， 能够与表面修饰亲合素、羧基的固相载体反应，或者与丙烯酰胺单体聚合反应制备单链DNA 模板，包括：
[0024] 5' -SEQ ID NO :2_3' 所不：5' -agtctacaga actttgaaag tatgtg-3' ；
[0025] 5' -biotin-SEQ ID NO :3_3' 所不：5' -biotin-ctatgagagc agtcatttga ctttg_3,〇
[0026] 进一步的，在本发明中，所述核苷酸dNTPs合成实时合成释放的检测分子是相同 的，所述检测分子为化学发光检测的焦磷酸盐、电化学检测的氢离子或者光学检测的荧光 信号。
[0027] 进一步的，在本发明中，所述PCR扩增产物的分析指包括单个样本PCR产物和混合 样本PCR产物的分析。
[0028] 有益效果：本发明提供的两核苷酸不同步合成测序分析PCR产物单体型的方法， 应用两核苷酸同时加入到反应中的一次循环测序信息，根据PCR产物中不同DNA模板的合 成测序不同步，使得相邻两个SNP位点上不同碱基的信息是由不同测序反应所提供。根据 这些不同步合成测序反应所得到的不同SNP位点上的测序信息进行关联分析，确定PCR产 物单体型。该发明与传统AS-PCR相比，具有如下有益效果：
[0029] 1)本发明在测序长度能够保障的前提下，适合所有包括两个关联位点的PCR产物 的单体型分析。相对于对传统AS-PCR对PCR引物设计的限制，本发明对PCR引物没有限制， 任何包括相邻两位点序列的PCR产物均适合分析，拓广了分析范围。
[0030] 2)本发明只需对PCR产物进行一次循环测序，相对于对传统AS-PCR方法需要两次 扩增、两次测序，本发明操作简单、可以降低分析成本。
[0031] 3)本发明能够对分析的各种单体型进行定量分析，相对于对传统AS-PCR对PCR产 物进行的定性分析，本发明不但可以进行定性、同时可以进行定量分析，为样本的分析提供 更多的参数。
[0032] 4)本发明适合多个混合样本的单体型分析，可以直接测定混合样本中各个单体型 的比例，可以用于从大规模样本中寻找、筛选核酸标志物。
[0033] 5)本发明直接采用商品化、非标记的天然核苷酸进行合成测序，它可以在任何基 于实时合成测序的测序平台进行。
【附图说明】
[0034]图1为本发明方法分析一个样本中包含MMP7(U25346)基因A/G(rsll568818)和 C/T(rsll568819)核酸片段的PCR产物的单体型分析测序结果；图中纵坐标为信号强度 (au)，横坐标为加入两核苷酸试剂的顺序，其中A= (A+T)，T= (C+G)。
【具体实施方式】
[0035] 本发明利用两核苷酸不同步合成测序其测序长度的特点，"跟踪"不同步合成在不 同DNA模板上的测序位置，并获取这些不同SNP位点的不同步测序信息就变得可能。同时 利用这些不同测序信息提供不同DNA模板的加和信息，并从这些加和信息中构建方程组、 并求解出不同DNA模板的含量，即实现不同单体型的实验分析。从而为单体型与疾病的相 关性研宄，以及实际样本的单体型临床检测提供一种新的分析方法。
[0036] 下