本发明采用一种引物选择合成测序分析PCR产物单体型的方法对两个或者两个以上位点变化（如SNP，突变等）的PCR产物（多个DNA模板）进行单体型分析，PCR产物被特定测序引物杂交后分成两份，各自加入一种3’端封闭的单体试剂，测序引物在DNA模板指导下延伸位点发生变化的一个碱基后，该测序引物3’端被封闭而不能继续延伸，然后对未封闭的引物继续合成测序，得到特定DNA模板的序列信息,从而确定PCR产物的单体型。
1.一种引物选择合成测序分析PCR产物单体型的方法，其特征在于：一段核酸序列的待检样本中含有两个以上的序列变化位点组成的单体型DNA模板，沿着第一个序列变化位点有选择地限制特定的测序引物进行合成测序，从而测定相应DNA模板的序列，确定其单体型，具体步骤为：步骤一扩增：利用包含一个5’端修饰活性基团的PCR引物，对待检样本进行扩增；
步骤二DNA单链制备：将PCR产物利用5’端修饰活性基团与其它经过修饰的载体结合，并通过变性将未固定的PCR另一条DNA链清除，获取固定的单链DNA模板；
步骤三测序引物的杂交：根据检测位点所包含的可能信息，在PCR单链DNA模板的第一个分析位点前，杂交完全互补的测序引物；
步骤四测序引物的封闭：根据检测位点所包含的可能信息，加入一种3’端封闭的单体试剂，即ddNTPs或3’-O-烯丙基、3’-O-氰乙基、3’-O-叠氮甲基修饰的核苷酸之一，测序引物在DNA模板指导下延伸位点发生变化的一个碱基后，该测序引物3’端被封闭而不能继续延伸；
步骤五测序引物的延伸测序：当测序引物被部分封闭后，未封闭的测序引物采用现有技术延伸，实现对DNA模板序列的测定；
步骤六DNA模板单体型的确定：根据步骤五测序信息，辨识出DNA模板中序列变化位点得到碱基信息，构建DNA单体型。
2.根据权利要求1所述的一种引物选择合成测序分析PCR产物单体型的方法，其特征在于序列变化位点包括SNP、突变、插入或缺失造成的碱基序列发生变化。
3.根据权利要求1所述的一种引物选择合成测序分析PCR产物单体型的方法，其特征在于所述PCR引物5’端修饰的活性基团为在热条件下稳定的丙烯酰胺、生物素或氨基之一，所述载体为丙烯酰胺、表面修饰亲和素或醛基或羧基基团的载体之一，使得PCR引物5’端修饰的活性基团能够与载体、或者载体表面基团发生化学反应而固定到载体上。
4.根据权利要求1所述的一种引物选择合成测序分析PCR产物单体型的方法，其特征在于所述合成测序包括单核苷酸加入焦磷酸测序、两核苷酸加入焦磷酸测序，桑格测序技术，在焦磷酸测序方法测定分析序列时，测序引物经封闭反应后，其测序反应体系可以直接用于测序反应；而在使用桑格技术测序时，应先将测序引物经过封闭反应后的反应液清除，配制桑格测序反应液、使测序引物延伸，延伸单链产物用于电泳测序分析。
5.一种对MMP7基因PCR产物进行分析的方法，该方法是根据权利要求1所述方法进行的，该PCR产物的核苷酸序列包含rs11568818 SNP位点和rs11568819 SNP位点相关联所组成的单体型，具体步骤为：a.将所有需检测的样本处理后提取DNA用作扩增模板，PCR扩增引物为SEQ ID NO:1，即
5’-AGTCTACAGAACTTTGAAAGTATGTG-3’和SEQ ID NO :2，即5’-biotin-CTATGAGAGCAGTCATTTGACTTTG-3’；
b.将PCR扩增产物与亲和素修饰的磁珠反应，使修饰生物素的DNA链固定到磁珠上，在
0.2M NaOH溶液下变性，将未固定的另一条DNA链清除，然后用洗液洗涤，得到固定的单链DNA模板，分成两份用于下列实验；
c.将测序引物SEQ ID NO:3，即：5’-CTAAAACGAGGAAGTATTACATC-3’与磁珠固定的其中一份单链DNA模板进行杂交；
d.在焦磷酸测序反应体系中加入ddATP使测序引物完成一个碱基延伸；
e.将上述d的反应体系置于焦磷酸测序仪中，按照指定的核苷酸加入顺序对未封闭测序引物杂交的DNA模板进行序列测定;f. 按照c、d、e的流程，将测序引物SEQ ID NO: 3，即：5’-CTAAAACGAGGAAGTATTACATC-
3’与磁珠固定的第二份单链DNA模板杂交，在焦磷酸测序反应体系中加入ddGTP使测序引物完成一个碱基延伸, 最后反应体系置于焦磷酸测序仪中进行序列测定;g根据e、f得到的测序信息，辨识出DNA模板中序列变化位点得碱基信息，构建DNA单体型。
一种引物选择合成测序分析PCR产物单体型的方法
技术领域
[0001] 本发明涉及生物技术领域，尤其涉及一种采用引物选择合成测序分析PCR产物单体型的方法。
背景技术
[0002] 在分子生物学研究中，DNA的序列分析是进一步研究和改造目的基因的基础，对生命科学至关重要。随着人类基因组计划和各种模式生物基因组计划的开展和完成，使人类步入了后基因时代，分子生物学相关学科得到了迅猛的发展。从基因水平上认识生命的差异，疾病发生、发展的规律，以及药物与生命体的相互作用将成为可能。在基础研究方面，研究疾病基因的遗传规律，克隆致病基因；在应用方面，直接寻找疾病的易感基因突变位点，可以获得有关与该疾病相关基因型的信息。通过DNA序列信息的分析可以为理解生命的起源、疾病发生、以及个体化医疗等提供有用的相关核酸序列信息。人的基因组中存在诸多的基因变异引起的疾病，如目前已经确定4000多种遗传疾病是由单碱基突变引起的，而一些重大疾病如癌症、糖尿病、心血管疾病、抑郁症、哮喘等，是受众多基因相关。然而基因型的信息本身是一种低外显度遗传变异，加上病例-对照研究方法的局限性，可能会出现假阴性和假阳性结果；另外,在预测功能性SNP 时,主要集中在可能引起蛋白质中氨基酸改变和影响转录调节区活性的位点上,这样势必漏掉一些其他重要的功能位点,。已有的研究发现,位于同一条染色体上的不同位点基因型之间不是孤立的，相邻的等位基因往往倾向于同时出现。如果相邻的等位基因同时出现的频率超过了随机发生的概率时,称为连锁不平衡，处于连锁不平衡的、相邻的SNP 等位位点也倾向于以一个整体遗传给后代。位于一条染色体特定区域的一组相互关联并倾向于以整体遗传给后代的基因型 的组合叫作单体型，单体型的存在可以极大地减少疾病关联研究的工作量。
[0003] 目前普遍采用间接的方法则是利用多个基因型的分型结果，以统计方法推断单体型，这种方法得到广泛了的应用，且用于构建单体型的相关软件也非常之多。然而，软件分析的结果对于临床样本而言，如果得不到直接的实验证实终究是不踏实的，临床检测也很难接受这样的分析结果。通过实验方法得到自然状态下真实的单体型情况需要进行单分子测序(如克隆测序、高通量测序等)，这些方法由于成本高，不适于大规模的个体样本分析。对于常规PCR产物中包含多个SNP位点的实验分析单体型方法，均采用一种等位基因特异引物-PCR(AS-PCR)的方法，以第一个SNP位点为扩增起始点，用两种特异引物将起始为两种的不同DNA模板分别进行扩增，然后分别进行测序，以确定其它SNP位点的碱基信息。然而，这种方法由于需要等位基因的特异扩增引物，因而扩增引物的优化设计受到限制，其扩增效率会因不同的扩增片段而不相同，甚至无法扩增出目的片段，操作繁琐、效率不高；同时需要标记两个特异扩增引物，费用也不低廉。
[0004] 现有AS-PCR方法分析单体型是将不同的DNA模板分别扩增，然后进行序列分析来实现的。如果能将不同DNA模板同时进行PCR扩增，但分别进行测序，同样可以达到单体型分析的目的，且扩增引物可以优化设计，适合所有两相邻基因型位点间的序列分析。基于上述原理，本发明提出一种PCR产物的单体型分析方法，该方法可以设计优化PCR引物扩增需要分析的核酸序列片段，然后进行序列分析，适合任意两（多）位点间的单体型分析，操作相对简便，价格低廉，可以用于临床检测。
发明内容
[0005] 本发明所解决的技术问题是：
[0006] 提供一种引物选择合成测序分析PCR产物单体型的方法，当待检样本中含有两个以上的序列变化位点组成的单体型DNA模板时，沿着第一个序列变化位点有选择地限制特定的测序引物进行合成测序，从而测定相应DNA模板的序列，确定其单体型。
[0007] 为了解决上述技术问题，本发明所采用的技术方案是：
[0008] 提供一种引物选择合成测序分析PCR产物单体型方法，针对待检样本中一段核酸序列中含有两个以上的序列变化位点组成的单体型DNA模板，可沿着第一个序列变化位点有选择地限制特定的测序引物进行合成测序，从而测定相应DNA模板的序列，确定其单体型，具体步骤为：
[0009] 步骤1）扩增：利用包含一个5’端修饰活性基团的PCR引物，对待检样本进行扩增；
[0010] 步骤2）DNA单链制备：将PCR产物利用5’端修饰活性基团与其它经过修饰的载体结合，并通过变性将未固定的PCR另一条DNA链清除，获取固定的单链DNA模板；
[0011] 步骤3）测序引物的杂交：根据检测位点所包含的可能信息，在PCR单链DNA模板的第一个分析位点前，杂交完全互补的测序引物；
[0012] 步骤4）测序引物的封闭：根据检测位点所包含的可能信息，加入一种3’端封闭的单体试剂，测序引物在DNA模板指导下延伸位点发生变化的一个碱基后，该测序引物3’端被封闭而不能继续延伸；
[0013] 步骤5）测序引物的延伸测序：当测序引物被部分封闭后，未封闭的测序引物采用现有技术延伸，实现对DNA模板序列的测定；
[0014] 步骤6）DNA模板单体型的确定：根据步骤5）测序信息，辨识出DNA模板中序列变化位点得到碱基信息，构建DNA单体型。
[0015] 上述的一种PCR产物的单体型分析方法中，序列变化位点包括SNP、突变、插入或缺失等造成的碱基序列发生变化。
[0016] 上述的一种PCR产物的单体型分析方法中， PCR引物5’端修饰的活性基团是在热条件下稳定的丙烯酰胺、生物素或氨基，所述载体为丙烯酰胺、表面修饰亲和素或醛基或羧基基团的载体，使得PCR引物5’端修饰的活性基团能够与载体、或者载体表面基团发生化学反应而固定到载体上。
[0017] 上述的一种PCR产物的单体型分析方法中，3’端羟基封闭的单体试剂包括ddNTPs，3’-O-烯丙基、3’-O-氰乙基、3’-O-叠氮甲基等修饰的核苷酸；合成测序包括包括单核苷酸加入焦磷酸测序、两核苷酸加入焦磷酸测序，桑格测序技术，在焦磷酸测序方法测定分析序列时，测序引物经封闭反应后，其测序反应体系可以直接用于测序反应；而在使用桑格技术测序时，应先将测序引物经过封闭反应后的反应液清除，配制桑格测序反应液、使测序引物延伸，延伸单链产物用于电泳测序分析。
[0018] 采用一种引物选择合成测序分析PCR产物单体型的方法来分析MMP7 (U25346)基因PCR产物，该PCR产物包含A/G (rs11568818)和C/T ( rs11568819)核酸片段的单体型，具体步骤为：
[0019] (1)将所有需检测的样本处理后提取DNA用作扩增模板，PCR扩增引物为SEQ ID NO:1，即5’-AGTCTACAGAACTTTGAAAGTATGTG-3’和SEQ ID NO :2，即5’-biotin-CTATGAGAGC AGTCATTTGACTTTG-3’；
[0020] (2) 将PCR扩增产物与亲和素修饰的磁珠反应，使修饰生物素的DNA链固定到磁珠上，在0.2M NaOH溶液下变性，将未固定的另一条DNA链清除，然后用洗液洗涤，得到固定的单链DNA模板，分成两份用于下列实验；
[0021] (3)将测序引物SEQ ID NO:3，即：5’-CTAAAACGAGGAAGTATTACATC-3’与磁珠固定的其中一份单链DNA模板进行杂交；
[0022] (4)在焦磷酸测序反应体系中加入ddATP使测序引物完成一个碱基延伸；
[0023] (5)将上述（4）的反应体系置于焦磷酸测序仪（PSQ 96MA system (Biotage AB, Uppsala, Sweden)）中，按照指定的单个核苷酸加入顺序对未封闭测序引物杂交的DNA模板进行序列测定。
[0024] (6) 按照(3),(4),(5)的流程，将测序引物SEQ ID NO: 3，即：5’-CTAAAACGAGGAAGTATTACATC-3’与磁珠固定的第二份单链DNA模板杂交，在焦磷酸测序反应体系中加入ddGTP使测序引物完成一个碱基延伸, 最后反应体系置于焦磷酸测序仪中进行序列测定。
[0025] (7)根据(5),(6)得到的测序信息，辨识出DNA模板中序列变化位点得碱基信息，构建DNA单体型。
[0026] 本发明的有益效果是：
[0027] 提供了一种引物选择合成测序分析PCR产物的单体型的方法对两个或者两个以上位点变化的PCR产物进行单体型分析，PCR产物被特定测序引物杂交后分成两份，各自加入一种3’端封闭的单体试剂，测序引物在DNA模板指导下延伸位点发生变化的一个碱基后，该测序引物3’端被封闭而不能继续延伸，然后对未封闭的引物继续合成测序，得到特定DNA模板的序列信息,从而确定PCR产物的单体型。
[0028] 与现有等位基因特异引物-PCR(AS-PCR)技术分析单体型相比，本发明对PCR引物没有限制，可以随意在包含分析片段内的区域设计、优化PCR引物，对分析的片段区域没有限制，并扩增出足够可用于序列测定的DNA模板量，适合所有片段的单体型分析。
附图说明
[0029] 图1.为一种测序引物选择合成分析PCR产物的单体型分析方法的一般流程。图中1为固相载体，2为PCR产物（DNA模板），3为互补杂交测序引物。PCR扩增产物的一条单链通过与载体1结合得到固定、并变性后制备出单链DNA模板2，与测序引物3杂交（a）； 加入ddATP，测序引物3在DNA模板第一个SNP位点碱基T的指导下发生合成反应（b），此后，该测序引物被封闭，不能继续合成；加入dNTP进行焦磷酸测序反应（c），测序信号均来自含有“C”的DNA模板，当测序测定第二个SNP位点的碱基信息时，便能知道第二个SNP位点碱基与第一个位点碱基C的关联，从而确定单体型。
[0030] 图2为MMP7 (U25346)基因PCR产物包含A/G (rs11568818)和C/T ( rs11568819)核酸片段的测序结果。（1）为单链PCR产物与测序引物杂交、且当测序引物用ddATP封闭后，按照仪器上编辑程序依次加入单个核苷酸的焦磷酸测序结果，其序列特征表明：第一个SNP位点碱基G相关联的第二个SNP位点的碱基信息为T和C; （2）为单链PCR产物与测序引物杂交、且当测序引物用ddGTP封闭后，按照仪器上编辑程序依次加入单个核苷酸的焦磷酸测序结果，其序列特征表明：第一个SNP位点碱基A相关联的第二个SNP位点的碱基信息为C，与T无关联。
具体实施方式
[0031] 下面将结合说明书附图，对本发明作进一步的说明。
[0032] 实例1： MMP7 (U25346)基因PCR产物包含A/G (rs11568818)和C/T ( rs11568819)核酸片段的单体型分析
[0033] (1)将血样本采用传统的蛋白激酶K与苯酚/氯仿抽提法提取外周血中的基因组DNA；
[0034] (2)将PCR引物（5’-AGTCTACAGAACTTTGAAAGTATGTG-3’, 5’-biotin-CTATGAGAGC AGTCATTTGACTTTG-3’， 200 ng基因组DNA，0.2mM dNTP，1 U Taq DNA聚合酶, 1×扩增缓冲液, 1.8 mM MgCl2的50μL PCR扩增体系进行扩增。扩增条件为：94℃起始变性5分钟，35个热循环：94℃变性30秒～61℃退火45秒～72℃延伸45秒，最后72℃延伸7分钟。
[0035] (3) 将PCR扩增与亲和素修饰的磁珠反应，使修饰生物素的DNA链固定到磁珠上，在0.2M NaOH溶液下变性，将未固定的另一条DNA链清除，然后用洗液(10 mM Tris–Acetate, pH 7.6)洗涤，得到固定的单链DNA模板，分成两份用于下列实验；
[0036] (4) )将测序引物（5’-CTAAAACGAGGAAGTATTACATC-3’）与磁珠固定的其中一份单链DNA模板模板在反应体系（(10 mM Tris–HCl, 2 M NaCl, 1 mM EDTA(乙二胺四乙酸钠), 0.1% Tween 20, pH 7.6) 80℃下放置5分钟，然后自然冷却至室温，完成杂交；
[0037] (5) 在焦磷酸测序反应体系(0.1 M Tris-Ac (pH 7.7), 2 mM EDTA(乙二胺四乙酸钠), 10 mM Mg(Ac)2, 0.2% BSA(小牛血清蛋白), 10 mM DTT(二巯苏糖醇), 10 mM APS(磷酰硫酸腺), 0.4 mg/mL PVP(聚乙烯吡咯烷酮), 4 mM D-luciferin(虫荧光素), 2 U/mL ATP sulfurylase(三磷酸腺苷硫酸化酶), 0.4 mM luciferase (荧光素酶), 2 U/mL apyrase VII(三磷酸腺苷双磷酸酶VII), 2 U/mL DNA聚合酶 I (Klenow fragment, exo–)中加入ddATP使测序引物在室温下完成一个碱基延伸；
[0038] (6) 将上述（5）的反应体系置于焦磷酸测序仪（PSQ 96MA system (Biotage AB, Uppsala, Sweden)）中，按照仪器上编辑程序依次加入单个核苷酸对未封闭测序引物在DNA模板指导下进行合成测序。得到图2（1）中的图谱，根据系列特征可以得出，第一个SNP位点碱基G相关联与第二个SNP位点的碱基T和C均由关联。
[0039] (7) 按照(4),(5),(6)的流程，将测序引物（5’-CTAAAACGAGGAAGTATTACATC-3’）与磁珠固定的第二份单链DNA模板模板杂交，在焦磷酸测序反应体系中加入ddGTP使测序引物完成一个碱基延伸, 最后反应体系置于焦磷酸测序仪中，按照仪器上编辑程序依次加入单个核苷酸，对未封闭测序引物在DNA模板指导下进行合成测序。得到图2（2）中的图谱，根据序列特征可以得出，第一个SNP位点碱基A相关联的第二个SNP位点的碱基信息为C，与T无关联。
[0040] (8)根据图2中的图谱，构建出的DNA单体型为下列三种：
[0041] 单体型1（SEQ ID NO:4)： 5’GTTATTGGCAGGAAGCACACAATGTATTTGTCTTTCAAA-3’[0042] 单体型 2（SEQ ID NO:5):5’GTTATTGGCAGGAAGCACACAATGTATTCGTCTTTCAAA-3’[0043] 单体型3（SEQ ID NO:6) : 5’-ATTATTGGCAGGAAGCACACAATGTATTCGTCTTTCAAA-3’[0044] 以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
