CN107893125B

CN107893125B - 用于鉴定桃花铃型/蔷薇型性状的单核苷酸多态性标记位点、引物对、试剂盒及应用

Info

Publication number: CN107893125B
Application number: CN201711405914.6A
Authority: CN
Inventors: 曹珂; 王力荣; 朱更瑞
Original assignee: Zhengzhou Fruit Research Institute CAAS
Current assignee: Zhengzhou Fruit Research Institute CAAS
Priority date: 2017-12-22
Filing date: 2017-12-22
Publication date: 2021-03-09
Anticipated expiration: 2037-12-22
Also published as: CN107893125A

Abstract

本发明公开了用于鉴定桃花铃型/蔷薇型性状的单核苷酸多态性标记位点、引物对、试剂盒及应用，所述的单核苷酸多态性标记位点是桃基因组第8染色体的第14,484,624位核苷酸，该核苷酸为A或C。本发明通过对大量桃种质样品的测序比对，鉴定出4063377个SNPs，利用这些SNPs对129份种质的花型(铃型/蔷薇型)性状进行全基因组关联分析，鉴定出与桃花型(铃型/蔷薇型)有显著关联的SNP位于第8染色体的第14,484,624bp处。并设计了特定的PCR引物扩增对和单碱基延伸引物，将完成单碱基延伸反应的的产物进行处理和分析，根据不同产物的分子量大小获得其基因分型结果。本发明的检测方法具有简单、快速、成本低的优点，能够实现生产中大规模的应用，具有良好的社会和经济效益。

Description

用于鉴定桃花铃型/蔷薇型性状的单核苷酸多态性标记位点、引物对、试剂盒及应用

技术领域

本发明涉及用于鉴定桃花型(铃型/蔷薇型)性状的单核苷酸多态性标记位点、引物对、试剂盒及应用，属于生物技术领域。

背景技术

选择是育种中最重要的环节之一，它是指在一个群体中选择符合要求的基因型，来进行后续的培育。但在传统育种中，由于很难获知后代的基因型，因此选择的依据通常是表现型而非基因型，这种选择方法对质量性状而言一般是有效的，但对数量性状来说，因为其表现型与基因型之间缺乏明确的对应关系，因而效率不高。此外，对于以果实性状为目标的果树来说，这些性状都有其特定的表现时期，通常需要度过3-5年甚至更长时间的童期，因而选择的时间较晚。这对于那些植株高大、占地多、生长季长的作物，特别是果树之类的园艺作物，显然是非常不利的。

近20年来迅速发展起来的基于DNA的分子标记技术，即“分子标记辅助选择”(marker-assisted selection，缩写为MAS)给育种提供了崭新的途径。它通过分析与目的基因紧密连锁的分子标记的基因型来进行育种，从而达到提高育种效率的目的。Yamamoto(2001)利用AFLP(Amplified Fragment Length Polymorphism，扩增片段长度多态性)、RAPD(random amplified polymorphic DNA，随机扩增多态性DNA)、SSR(Simple SequenceRepeats，简单重复序列)和RFLP(Restriction Fragment Length Polymorphism，限制性内切酶片段长度多态性)标记，将果皮颜色定位第6连锁群，与之连锁距离最近(3.7cM)的SSR标记为UDP96-015。由于进行分子标记辅助选择有两个重要的前提，首先是必须得到与目标性状紧密连锁的标记，即建立目标基因与分子标记的连锁关系。其次是检测的自动化，由于分子标记辅助选择要求对育种群体进行大规模检测，因而要求检测的方法要简单、快速、成本低、比较准确，以实现检测过程(包括DNA的提取、分子标记的检测、数据分析等)的自动化。但是，可以发现在过去很长一段时间内采用的分子标记，如RFLP(RestrictionFragment Length Polymorphism，限制性内切酶片段长度多态性)、RAPD(randomamplified polymorphic DNA，随机扩增多态性DNA)、AFLP和SSR等通常需要酶切或PCR后用电泳检测分型结果，很难实现这一点。

SNPs标记(single nucleotide polymorphisms，单核苷酸多态性)主要是指在基因组水平上由单个核苷酸的变异所引起的DNA序列多态性。它在基因组上分布广泛，数量众多，因此很容易检测到相对于RFLP、RAPD、AFLP、SSR标记更连锁的标记，且在检测时具有高通量、简单、快速、高灵敏度的优点，是进行分子标记辅助育种中最具潜力的标记。

桃(学名：Amygdalus persica L.)：蔷薇科、桃属植物。花型有蔷薇型和铃型两种。然而，现有的SSR标记与桃花型(铃型/蔷薇型)性状连锁距离远，在杂交后代的早期鉴定中准确率较低；而现有的与目标性状连锁的SNPs标记来自芯片鉴定的结果，由于芯片位点较少(少于9000个)，因此不能保证鉴定出的SNPs是与目标性状最关联或连锁的位点。

发明内容

针对现有技术的不足，本发明的目的是提供用于鉴定桃花铃型/蔷薇型性状的单核苷酸多态性标记位点、引物对、试剂盒及应用，该标记位点在鉴定桃花型性状时具有较高的准确率。

为了实现上述目的，本发明所采用的技术方案是：

用于鉴定桃花铃型/蔷薇型性状的单核苷酸多态性标记位点，所述的单核苷酸多态性标记位点是桃基因组第8染色体的第14,484,624位核苷酸，该核苷酸为A或C。

用于鉴定桃花铃型/蔷薇型性状的PCR扩增引物对，所述引物对中的上游引物是根据桃基因组第8染色体的第14,484,624位核苷酸及其上游序列进行设计，所述引物对中的下游引物是根据桃基因组第8染色体的第14,484,624位核苷酸的下游序列进行设计。

所述的上游引物的核苷酸序列如SEQ ID NO.2所示，下游引物的核苷酸序列如SEQID NO.3所示。

用于鉴定桃花铃型/蔷薇型性状的单碱基延伸引物，所述单碱基延伸引物为如SEQID NO.4、SEQ ID NO.5和SEQ ID NO.6所示的核苷酸序列。

用于鉴定桃花铃型/蔷薇型性状的试剂盒，包括PCR扩增引物对和单碱基延伸引物。

一种单核苷酸多态性标记位点在鉴定或辅助鉴定桃花铃型/蔷薇型性状中的应用。

一种单核苷酸多态性标记位点在桃花铃型/蔷薇型性状分子标记辅助选择育种方面的应用。

一种试剂盒在辅助选择育种方面的应用。

一种利用单核苷酸多态性标记位点鉴定桃花铃型/蔷薇型性状的方法，包括以下步骤：

(1)PCR扩增：以待测桃的基因组DNA为模板，用PCR扩增引物对进行PCR扩增，获得PCR扩增产物；

(2)SAP反应：配制SAP反应体系，加入至步骤(1)PCR扩增反应后的反应体系中，除去PCR扩增反应中未反应的dNTP；

(3)单碱基延伸反应：配制单碱基延伸反应体系，加入至步骤(2)SAP反应后的反应体系中；

(4)基因分型：将完成单碱基延伸反应的产物进行树脂脱盐处理，点在芯片上，由芯片扫描仪扫描，检测分析，根据不同产物的分子量大小进行基因分型，预测桃花铃型/蔷薇型性状。

本发明通过对大量桃种质样品的测序比对，鉴定出4063377个SNPs，利用这些SNPs对129份种质的花型(铃型/蔷薇型)性状进行全基因组关联分析，鉴定出与桃花型(铃型/蔷薇型)有显著关联的SNP位于第8染色体的第14,484,624bp处。

本发明针对桃基因组第2版的第8染色体的第14,484,624位核苷酸多态性的特性，设计了特定的PCR引物扩增对和单碱基延伸引物，将完成单碱基延伸反应的的产物进行树脂脱盐处理，点在芯片上，由芯片扫描仪扫描，进行MALDI-TOF质谱检测，Typer4.0软件分析实验数据，根据不同产物的分子量大小获得其基因分型结果。

本发明对7个杂交群体共170份待测桃单株进行SNP准确率的验证，结果表明，本发明对显性表型(花型铃型)中具有较高的准确率，可达94.48％；对隐性表型(花型蔷薇型)中准确率为97.50％。这说明，利用本发明的单核苷酸标记位点进行检测具有简单、快速、成本低的优点，能够实现生产中大规模的应用。

具体实施方式

以下结合实施例对本发明的具体实施方式作进一步详细说明。

实施例1SNPs标记位点的获得

本发明以从中国农业科学院郑州果树研究所桃种质资源圃随机获得的129份桃种质为样品，采用常规CTAB法提取样品DNA，并通过Illumina HiSeq 2000测序仪对129份桃种质进行重测序，得到121Gb数据，平均覆盖桃基因组89.28％，平均测序深度为4.21×左右。根据测序得到的50-150bp的reads，与桃参考基因组第2版(http://www.rosaceae.org/node/355)进行比对，鉴定出4063377个SNPs。利用这些SNPs对129份种质的表型性状进行全基因组关联分析，鉴定出与桃花型(铃型/蔷薇型)性状有显著关联的SNP位于第8染色体的第14,484,624bp处。

实施例2利用SNP标记鉴定桃花型(铃型/蔷薇型)性状的方法

1、DNA提取

采用常规CTAB法提取待测桃样品组织的DNA，去除RNA，DNA样品总体积不低于15μl。用紫外光光度计测定DNA样品在260nm、280nm处的OD值，计算DNA含量和OD_260/280的比值。DNA样品纯度OD_260/280值应在1.8-2.0之间，浓度稀释至10ng/μl。

2、设计引物

根据桃基因组第2版的第8染色体的第14,484,624位处左右各150bp序列(具体的核苷酸序列见表1)，设计引物。

表1SNP旁侧序列信息

其中，M表示A或C。

引物由生物技术公司合成后，稀释PCR扩增上下游引物至浓度均为0.5μM。稀释单碱基延伸引物至延伸引物1、2、3浓度分别为8μM、10μM、15μM，备用。引物序列见表2。

表2引物序列

3、PCR反应体系及Mass ARRAY分析

按照SEQUENOM-iPLEX标准操作流程进行PCR、SAP、单碱基延伸反应，主要步骤如下：

(1)PCR扩增反应

首先，配制PCR扩增体系，具体见表3。

表3Sequenom MassArray系统基因分型扩增PCR反应体系

其中，Primer Mix为浓度0.5μM的上游引物和下游引物各加入0.5μl。

将上述试剂转移到PCR管或板中，

PCR扩增程序如下：94℃15min；94℃20s，56℃30s，72℃1min，共45个循环；72℃3min；4℃保存。

配置1％琼脂糖凝胶，对PCR反应后的产物进行电泳，如果条带明亮且单一，则进行后续试验。

(2)SAP反应

利用该反应除去未用尽的dNTP，首先配制反应体系，见表4：

表4SAP酶促反应体系

将配好的上述溶液2μl加入到完成PCR反应的PCR管或板中。按如下程序进行SAP反应。

SAP反应程序：37℃40min；85℃5min；4℃保存。

(3)单碱基延伸反应

首先配制单碱基延伸反应体系，见表5：

表5单碱基延伸反应体系

其中，iPLEX Extend Primer Mix为单碱基延伸引物1、2、3的混合引物，引物1、2、3的浓度分别为8μM、10μM、15μM，单碱基延伸引物1、2、3的加入量相同，共0.94μl。

然后将配好的溶液2μl加入到SAP反应后的PCR管或板中，进行延伸反应，程序如下：94℃30s；94℃5s，(52℃5s，80℃5s，5个循环)，共40个循环；72℃3min；4℃保存。

将完成延伸反应的的产物进行树脂脱盐处理，点在芯片上，由芯片扫描仪扫描，进行MALDI-TOF质谱检测，Typer4.0软件分析实验数据，根据不同产物的分子量大小获得其基因分型结果。当完成延伸反应的产物分子量为7247.6左右时，Chr8：14,484,624bp处为纯合位点，分型为A/A，对应种质的花型为蔷薇型。分子量为7207.7左右时，Chr8：14,484,624bp处为纯合位点，分型为C/C；分子量为7227.6左右时，Chr8：14,484,624bp处为杂合位点，分型为A/C；C/C和A/C这两种分型对应的种质花型为铃型。

实施例3 7个杂交群体利用桃花型(铃型/蔷薇型)性状SNP标记进行表型性状的盲测验证

1、实验材料的选择

以郑州果树研究所资源圃中种植的常规桃品种为实验材料，从中选取调查过表型性状的7个杂交群体共170个单株，具体见表6。

表6供试杂交群体的名称及群体大小信息

2、利用桃花型关联SNP标记的鉴定方法

以本发明桃基因组第8染色体的第14,484,624位作为核苷酸多态性标记位点，对7个杂交群体共170个桃单株花型(铃型/蔷薇型)性状进行盲测鉴定，具体的鉴定方法参照实施例2的方法。

3、杂交群体中分型结果对表型的预测(表7)

表7显示了花型(铃型/蔷薇型)关联SNP在不同杂交群体的鉴定结果。

表7桃花型(铃型/蔷薇型)关联SNP在7个杂交群体的分型结果及表型预测准确率

综合分析，本发明在7个群体鉴定表型的准确率如表8所示，可以看出本发明针对显性表型准确率为94.48％，针对隐性表型准确率为97.50％，平均准确率为95.99％。

表8本发明涉及的SNP在对7个杂交群体鉴定后的表型预测能力

综上所述，本发明的SNP位点能够帮助实现利用桃幼苗DNA早期预测桃成龄后花型(铃型/蔷薇型)性状的目的，该方法是利用重测序得到的400万以上SNPs中全基因组关联分析得到最关联的位点从而实现的，由于原始SNPs数量庞大，从而保证了得到的关联标记关联性高，有效提高了利用SNPs预测的准确性。

以上所述仅为本发明最佳的实施例，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

序列表

<110> 中国农业科学院郑州果树研究所

<120> 用于鉴定桃花铃型/蔷薇型性状的单核苷酸多态性标记位点、引物对、试剂盒及应用

<160> 6

<170> SIPOSequenceListing 1.0

<210> 1

<211> 301

<212> DNA

<213> 桃(Amygdalus persica L.)

<220>

<223> m=a或c

<400> 1

agatgcatgt atatatataa gtaaaaaaga gagagagaga gagagagaga gagagagaga 60

ccgctcatcg tttgtcggga ggggggaccc aacagatccg tcgttattta tgaatccaac 120

tgtacatata cgtgtacaat agatagtcca matgattttg cataggtgtt acccaacaga 180

ccaaatcagt agttatttat tacataatgc ctaattaggg ttaaggagct tcggaaaata 240

attttgctgc ttgcggcagg actttaatac tgattccgag agatgataaa gaacttcaat 300

a 301

<210> 2

<211> 30

<212> DNA

<213> 人工序列()

<400> 2

acgttggatg tccaactgta catatacgtg 30

<210> 3

<211> 30

<212> DNA

<213> 人工序列()

<400> 3

acgttggatg ctactgattt ggtctgttgg 30

<210> 4

<211> 23

<212> DNA

<213> 人工序列()

<400> 4

cctcaacacc tatgcaaaat cat 23

<210> 5

<211> 24

<212> DNA

<213> 人工序列()

<400> 5

cctcaacacc tatgcaaaat catg 24

<210> 6

<211> 24

<212> DNA

<213> 人工序列()

<400> 6

cctcaacacc tatgcaaaat catt 24

Claims

1.用于鉴定桃花铃型/蔷薇型性状的SNP分子标记，其序列如SEQ ID NO.1所示，其中M为A或C。

2.用于鉴定桃花铃型/蔷薇型性状的PCR扩增引物对，其特征在于，所述引物对中上游引物的核苷酸序列如SEQ ID NO.2所示，下游引物的核苷酸序列如SEQ ID NO.3所示。

3.用于鉴定桃花铃型/蔷薇型性状的单碱基延伸引物，其特征在于，所述单碱基延伸引物为如SEQ ID NO.4、SEQ ID NO.5和SEQ ID NO.6所示的核苷酸序列。

4.用于鉴定桃花铃型/蔷薇型性状的试剂盒，其特征在于，所述的试剂盒包括权利要求2所述的PCR扩增引物对和权利要求3所述的单碱基延伸引物。

5.一种如权利要求1所述的SNP分子标记在鉴定或辅助鉴定桃花铃型/蔷薇型性状中的应用。

6.一种如权利要求1所述的SNP分子标记在桃花铃型/蔷薇型性状分子标记辅助选择育种方面的应用。

7.一种如权利要求4所述的试剂盒在辅助选择育种方面的应用。

8.一种利用权利要求1所述SNP分子标记鉴定桃花铃型/蔷薇型性状的方法，其特征在于，包括以下步骤：

(4)基因分型：将完成单碱基延伸反应的产物进行树脂脱盐处理，点在芯片上，由芯片扫描仪扫描，检测分析，根据不同产物的分子量大小进行基因分型，预测桃花铃型/蔷薇型性状；

其中，PCR扩增引物中上游引物的核苷酸序列如SEQ ID NO.2所示，下游引物的核苷酸序列如SEQ ID NO.3所示；

单碱基延伸反应体系中单碱基延伸引物为如SEQ ID NO.4、SEQ ID NO.5和SEQ IDNO.6所示的核苷酸序列。