通过检测SNP位点确定等位基因的基因型的方法与流程

本发明涉及分子生物学技术领域，特别涉及一种通过检测snp位点确定等位基因的基因型的方法。

背景技术：

等位基因是指位于一对同源染色体相同位置上控制同一性状不同形态的基因，英文原文为"allele"。不同的等位基因产生例如发色或血型等遗传特征的变化。等位基因控制相对性状的显隐性关系及遗传效应，可将等位基因区分为不同的类别。在个体中，等位基因的某个形式(显性的)可以比其他形式(隐性的)表达得多。当一个生物体带有一对完全相同的等位基因时，则该生物体就该基因而言是纯合的(homozygous)或可称纯种(true-breeding)；反之，如果一对等位基因不相同，则该生物体是杂合的(heterozygous)或可称杂种(hybrid)。等位基因各自编码蛋白质产物，决定某一性状，并可因突变而失去功能。等位基因之间存在相互作用。当一个等位基因决定生物性状的作用强于另一等位基因并使生物只表现出其自身的性状时，就出现了显隐性关系。作用强的是显性，作用被掩盖而不能表现的为隐性。

单核苷酸多态性(singlenucleotidepolymorphism,snp)是一种遗传标记，位于人类基因组上。随着人类基因组计划的完成，人们发现人类的dna序列并非完全一致的，其中包含着某些变异。这些变异导致了多种遗传性疾病的产生。作为第三代遗传标志，snp位点与遗传性疾病的产生息息相关。故对snp位点进行研究及检测，将为发现遗传疾病高危群体、检测疾病的相关基因以及对生物医学的基础研究等方面提供一条新的出路。现有的snp的检测方法大致可以分为两大类：一大类是以单链构象多态性(sscp)、等位基因特异性pcr(allele-specificpcr，as-pcr)等为代表的以凝胶电泳为基础的传统经典的检测方法。另一大类是以直接测序、质谱检测技术等为代表的高通量、自动化程度较高的检测方法。多种检测方法中直接测序法准确度最高，常作为snp位点检测的金标准。但该方法所需仪器成本高，测序耗时较长，很难在一般性医院或实验室中普及。

技术实现要素：

为了解决上述技术问题，本发明提供一种通过检测snp位点确定等位基因的基因型的方法，所述方法通过as-pcr法检测待测等位基因的一个或多个snp位点，然后根据每个snp位点的检测结果，确定每个待测等位基因的基因型。

在一种实施方式中，所述方法包括以下步骤：步骤1：针对待测等位基因的每个snp位点设计三条扩增引物，其中一条扩增引物是上游引物或下游扩增引物，相应地另外二条扩增引物是下游扩增引物或上游扩增引物；所述二条扩增引物3’末端碱基与snp位点碱基分别是相同和互补，所述二条扩增引物的其它各个位置碱基相同；步骤2：构建符合所述snp位点所有可能基因型的阳性质粒对照组，根据阳性对照组qpcr结果得到各待测snp位点基因型判读区间；和，步骤3：使用步骤1设计的特异性扩增引物对待测样本进行qpcr反应，所得结果根据步骤2确定的snp位点不同基因型判读区间判别待测样本snp位点的基因型。

在一种实施方式中，所述步骤2包括以下步骤：a.对所述snp位点所有可能出现的基因型与引物组合设计至少三个不同梯度dna阳性模板实验组，每个snp位点有三种基因型，每种基因型依照步骤1的方法设计三条扩增引物，每个snp位点共有六种不同的基因型和引物组合，然后以上述组合进行qpcr反应；b.引物与所对应的纯合基因型质粒模板为特异性组，引物与非对应的纯合基因型质粒模板为非特异性组，引物与混合基因型质粒模板为混合组；c.将特异性组、非特异性组和混合组的qpcr结果换算成基因拷贝数后，按纯合基因型/纯合基因型、杂合基因型/杂合基因型比较各基因拷贝数，分别得到各基因型各浓度梯度下特异性组合与非特异性组合、混合组组合比较的倍率比；d.将所得各种基因型的倍率比按对数形式表现出来，并绘制所述snp位点基因型定量分型区域图，得到各待测snp位点基因型判读区间。

在一种实施方式中，所述待测基因是聚集素基因。

在一种实施方式中，所述snp位点是rs11136000、rs2279590和/或rs9331888。

在一种实施方式中，所述rs11136000位点在步骤1中三条扩增引物分别是：上游引物seqno.3：5’-taaagaatccactcatccttcaaga-3’，下游引物seqno.4：5’-gcaagggcccgttagagaa-3’和seqno.5：5’-gcaagggcccgttagagag-3’。

在一种实施方式中，所述rs2279590位点在步骤1中三条扩增引物分别是：上游引物seqno.8：5’-ggaagtcctcctgctt-3’，上游引物seqno.9：5’-ggaagtcctcctgctc-3’和下游引物seqno.10：5’-tggaaaccacttttattc-3’。

在一种实施方式中，所述rs9331888位点在步骤1中三条扩增引物分别是：上游引物seqno.13：5’-ggaagtcctcctgctt-3’，上游引物seqno.14：5’-ggaagtcctcctgctc-3’和下游引物seqno.15：5’-tggaaaccacttttattc-3’。

单核苷酸多态性(snp)指基因组水平上由单个核苷酸的变异引起的dna序列多态性，包括转换、颠换、缺失和插入。因每个snp位点通常仅含两个等位基因，在检测时能通过一个简单的“+/-”分析进行基因型分型，无需分析片段长度。现有的snp的检测方法大致可以分为两大类：一大类是以单链构象多态性(sscp)、等位基因特异性pcr(allele-specificpcr，as-pcr)等为代表的以凝胶电泳为基础的传统经典的检测方法。另一大类是以直接测序、质谱检测技术等为代表的高通量、自动化程度较高的检测方法。多种检测方法中直接测序法准确度最高，常作为snp位点检测的金标准。但该方法所需仪器成本高，测序耗时较长，很难在一般性医院或实验室中普及。本发明所采用的as-pcr法既不需要昂贵的检测仪器，也没有冗杂的实验过程，仅需设计特异性的引物，结合qpcr技术，即可实现对待测基因的snp位点的高通量测序。

本发明方法检测结果与直接测序法结果进行对比，结果显示二者完全相符，这说明该方法对snp位点基因型的检测灵敏度和金标准相同。因其不需要昂贵的检测仪器，可以在小型实验室和医院的检验室中进行，很符合临床检验灵敏度高、方便快捷的要求。故该方法可用于临床医院检验室对clu基因与ad相关snp位点基因型的检测，为临床探寻ad的易感基因提供一种简单、便捷的检测方法。

附图说明

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来说，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。

图1是rs11136000(a)位点，rs2279590(b)位点和rs9331888(c)位点基因型判读示意图；

图2是基因测序法和本发明的qpcr定量分型法检测rs11136000基因型结果图，①，②是基因测序法检测结果，基因型分别为cc,ct；③是qpcr定量分型法结果；

图3是基因测序法和qpcr定量分型法检测rs2279590基因型结果图，①，②是基因测序法检测结果，基因型分别为cc,ct；③是qpcr定量分型法结果；和

图4是基因测序法和本发明的qpcr定量分型法检测rs9331888基因型结果图，①，②，③是基因测序法检测结果，基因型分别为cc,cg,gg；④是qpcr定量分型法结果。

具体实施方式

为了使本领域技术领域人员更好地理解本申请中的技术方案，下面将结合实施例对本发明作进一步说明，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都应当属于本申请保护的范围。

实施例

1.材料

1.1.实验样本13例血液样本均来自2019年北京中医药大学东直门医院的门诊或住院患者。其中男性8人，女性5人，所有研究对象签署书面知情同意书同意后，采集外周血2ml，并保存在-20℃冰箱。正常组3例，为健康且无ad家族遗传史者，其余10例血液样本均来自ad患者。

1.2.试剂与仪器全血基因组dna提取试剂盒(qiaampdnabloodminikit)(北京兰博利德商贸有限公司)；takarataqtmhotstardna聚合酶及其他pcr反应试剂(北京兰博利德商贸有限公司)；takaradl2000dnamarker(北京兰博利德商贸有限公司)；quantinovatmgreenpcrkit(北京兰博利德商贸有限公司)；超微量分光光度仪(美谷分子仪器(上海)有限公司)；t100^tm型pcr仪(biorad)；dyy-6c电泳仪(北京六一生物科技有限公司)；quantstudio^tm6flex荧光定量pcr仪(thermofisher)；引物由深圳华大基因公司合成；测序由上海生物工程公司完成。

2.方法

2.1.基因组dna的提取基因组的提取参考qiaampdnabloodminikit说明书进行。提取的核酸利用超微量分光光度仪测定dna的含量，利用琼脂糖凝胶电泳检测dna完整性，-20℃保存备用。

2.2.特异性引物的设计从genbank中下载clu全基因序列及待测snp位点的信息，使用primerpremier5.0软件对三个snp位点设计特异性引物。针对snp位点rs11136000设计3对特异性引物，外侧引物用于常规pcr扩增及测序，2对内侧引物用来鉴别rs11136000基因的3种基因型(cc型、ct型、tt型)。针对snp位点rs2279590设计3对特异性引物，外侧引物用于常规pcr扩增及测序，2对内侧引物用来鉴别rs2279590基因的3种基因型(cc型、ct型、tt型)。针对snp位点rs9331888设计3对特异性引物外侧引物用于常规pcr扩增及测序，2对内侧引物，用来鉴别rs9331888基因的3种基因型(cc型、cg型、gg型)。引物序列见下表1。

表1snp位点常规pcr和qpcr引物序列

2.3.基因测序分型使用设计的常规pcr引物及taqtmhotstardna聚合酶进行pcr扩增。pcr反应体系组成如下：10xpcrbuffer2.5ul，dntpmixture(2.5mm)2ul，takarataq酶(5u/ul)0.125ul，上下游引物(10um)各1ul，dna模板1ul，灭菌水17.375ul，反应总体系25ul。扩增反应条件：95℃3min，95℃30s，58℃30s，72℃30s，共35个循环；终延伸72℃5min，4℃保存。pcr产物通过琼脂糖凝胶电泳进行检测分析，dl2000dnamarker作为参考，回收并纯化所需dna条带，送由上海生物工程公司完成测序。

2.4.阳性对照组质粒模板的构建使用北京市东直门医院采集的3例健康者的血液作为样本，从血液中提取基因组dna，按上述方法进行pcr扩增，回收所需dna条带并送测，pcr产物作为模板构建质粒载体。因s11136000，rs2279590和rs9331888位点各对应2种碱基，三种基因型。故根据测序结果，以已经构建好的质粒载体为模板进行点突变实验，获得每个snp位点相对应的另一种碱基，并构建新的质粒载体。将已经构建好的两种纯合子质粒载体等比例混合，获得杂合质粒载体，从而构建符合各snp位点基因型的阳性对照组质粒模板。载体的构建使用pcr产物作为模板，pmdtm18-tvector(takara)作为载体，步骤参考takarapmdtm18-tvectorcloningkit说明书。点突变实验采用已经测序的质粒载体作为模板，点突变引物自行设计，选择quikchangelightningmultisite-directedmutagenesiskit中的pcr反应试剂，pcr反应体系组成如下：10xreactionbuffer5ul，dntpmixture1ul，上下游突变引物(10um)各1ul，quicksolutionreagent1.5ul，质粒dna模板1ul，灭菌水14.5ul，随后加入quickchangelighteningenzyme1ul，反应总体系26ul。扩增反应条件：95℃2min，95℃20s，60℃10s，68℃30s，共18个循环；终延伸68℃5min，4℃保存。pcr产物根据说明书构建质粒载体。

2.5.qpcr定量分型采用设计的qpcr引物，选择sybrgreenqpcrmaster进行qpcr检测。qpcr反应体系组成如下：2×sybrgreenmastermix10ul，上下游引物(10um)1ul，dna模板1ul，灭菌水7ul，反应总体系25ul。rs11136000、rs9331888扩增反应条件:95℃2min，95℃5s，64℃20s，共40个循环。rs2279590扩增反应条件:95℃2min，95℃5s，55℃20s，共40个循环。将实验组和阳性对照组所得结果进行比较，对各组临床样本进行snp分型。

2.6.基因型判读qpcr实验前对每个snp位点设计4个不同梯度的dna模板(10⁹拷贝数、10⁸拷贝数、10⁷拷贝数、10⁶拷贝数)阳性对照组。因一个snp位点有三种基因型，两对特异性引物，故每个snp位点共有6种不同的基因型和引物组合。其中引物与所对应的纯合基因型质粒模板为阳性对照组，引物与非对应的纯合基因型质粒模板为非特异性组，引物与混合基因型质粒模板为混合组。反应结束后，以ct值为横坐标，特异性引物与相对应的纯合子基因型组起始拷贝数浓度的对数为纵坐标，绘制定量标准曲线(相关系数r²>0.99)。通过定量标准曲线将混合组和非特异性组ct值等比例换算成拷贝数，换算结果按纯合基因型/纯合基因型，杂合基因型/杂合基因型比较，所得结果按倍率形式体现，再将倍率按logistic对数形式转换，并对结果进行分析。根据阳性对照组qpcr结果得到所需的snp位点基因型判读区间为：对于rs11136000,当结果落在[1,+∞]时基因型为cc；[-1,1]基因型为ct；[-1,-∞]基因型为tt。对于rs2279590,当结果落在[1,+∞]时基因型为cc；[-1,1]基因型为ct；[-1,-∞]基因型为tt。对于rs9331888,当结果落在[1,+∞]时基因型为cc；[-1,1]基因型为cg；[-1,-∞]基因型为gg。各snp位点基因型判读区间如下图1所示。以clur1rs11136000作为具体例子对qpcr定量分型的判读区间构建过程进行说明。

下表为clur1(rs11136000)阳性质粒六种不同基因组和引物组合的qpcr结果。其中clu-p1-for/p3c-rev引物与clur1的cc纯合质粒模板组合、clu-p1-for/p3t-rev引物与clur1的tt纯合质粒模板组合为特异性组；clu-p1-for/p3c-rev引物与clur1的tt纯合质粒模板组合、clu-p1-for/p3t-rev引物与clur1的cc纯合质粒模板组合为非特异性组；剩余两组为混合组。根据上述实验方法，以特异性组作为阳性标准设置四个梯度(10⁹拷贝数、10⁸拷贝数、10⁷拷贝数、10⁶拷贝数)，qpcr结果构建定量标准曲线(相关系数r²>0.99)，再通过定量标准曲线将混合组和非特异性组ct值等比例换算成拷贝数，以得到下表。

表2clur1(rs11136000)特异性引物与质粒模板qpcr结果

下图为六种组合(特异性组、非特异性组、混合组)qpcr结果换算拷贝数后，按纯合基因型/纯合基因型、杂合基因型/杂合基因型比较拷贝数，分别得到cc、tt、ct基因型各梯度下特异性组合与非特异性组合、混合组组合比较的倍率比。

表3clur1(rs11136000)各基因型qpcr定量分型倍率比

最后将所得三种基因型的倍率比按对数形式表现出来，并绘制基因型定量分型区域图。如图1a所示。

3.结果

3.1.基因组dna的提取10例ad患者的外周血采用dna提取试剂盒提取基因组dna，经超微量分光光度仪检测吸光度(a)，a260nm/a280nm比值均在1.8-2.0之间，表明提取dna纯度符合后续实验要求。

3.2.基因测序法和qpcr定量分型法结果比较对10组样本的rs11136000，rs2279590和rs93318883个snp位点进行基因测序和qpcr定量分型结果如图2，图3，图4所示。两种方法检测10例样本rs11136000，rs2279590和rs9331888位点基因型的符合率为100％，结果见表4-6。

表4两种方法检测rs11136000基因型分布(n＝10)

表5两种方法检测rs2279590基因型分布(n＝10)

表6两种方法检测rs9331888基因型分布(n＝10)

应该理解到披露的本发明不仅仅限于描述的特定的方法、方案和物质，因为这些均可变化。还应理解这里所用的术语仅仅是为了描述特定的实施方式方案的目的，而不是意欲限制本发明的范围，本发明的范围仅受限于所附的权利要求。

本领域的技术人员还将认识到，或者能够确认使用不超过常规实验，在本文中所述的本发明的具体的实施方案的许多等价物。这些等价物也包含在所附的权利要求中。

序列表

<110>北京中医药大学

<120>通过检测snp位点确定等位基因的基因型的方法

<160>15

<170>siposequencelisting1.0

<210>1

<211>22

<212>dna

<213>人工序列(artificialsequence)

<400>1

agaatccactcatccttcaaga22

<210>2

<211>25

<212>dna

<213>人工序列(artificialsequence)

<400>2

atttttgtatttttagtagagatgg25

<210>3

<211>25

<212>dna

<213>人工序列(artificialsequence)

<400>3

taaagaatccactcatccttcaaga25

<210>4

<211>19

<212>dna

<213>人工序列(artificialsequence)

<400>4

gcaagggcccgttagagaa19

<210>5

<211>19

<212>dna

<213>人工序列(artificialsequence)

<400>5

gcaagggcccgttagagag19

<210>6

<211>24

<212>dna

<213>人工序列(artificialsequence)

<400>6

cctggagacttctacagggaccgt24

<210>7

<211>26

<212>dna

<213>人工序列(artificialsequence)

<400>7

aatgtttggaaaccacttttattctt26

<210>8

<211>16

<212>dna

<213>人工序列(artificialsequence)

<400>8

ggaagtcctcctgctt16

<210>9

<211>16

<212>dna

<213>人工序列(artificialsequence)

<400>9

ggaagtcctcctgctc16

<210>10

<211>18

<212>dna

<213>人工序列(artificialsequence)

<400>10

tggaaaccacttttattc18

<210>11

<211>27

<212>dna

<213>人工序列(artificialsequence)

<400>11

ggagtctcagcaccacctgaacttgac27

<210>12

<211>25

<212>dna

<213>人工序列(artificialsequence)

<400>12

tcgtcctgttggttctgtgatggca25

<210>13

<211>25

<212>dna

<213>人工序列(artificialsequence)

<400>13

gagcaggcttcccagagaaagtccg25

<210>14

<211>25

<212>dna

<213>人工序列(artificialsequence)

<400>14

gagcaggcttcccagagaaagtccc25

<210>15

<211>25

<212>dna

<213>人工序列(artificialsequence)

<400>15

tcgtcctgttggttctgtgatggca25