CN110607383A - 一种与茶树表儿茶素含量连锁的分子标记位点及其应用 - Google Patents

Abstract

本发明公开了一种与茶树表儿茶素(epicatechin，EC)含量连锁的分子标记位点及其应用，本发明首次发现了位于茶树基因组Scaffold2233：468642上的与茶树表儿茶素含量相关的SNP分子标记位点，其基因型与表儿茶素含量极显著相关。通过相关性和显著性验证结果表明，AA基因型样本对应的茶汤干物质中表儿茶素含量与GG和GA基因型样本相比具有极显著差异。从统计学上判断，当基因型为双突变AA时，茶树中干物质中表儿茶素含量大概率高于正常平均水平的野生型GG或单突变基因型GA。进一步建立检测该位点的检测方法，可用于评价茶树的表儿茶素含量，以进一步用于高EC(epicatechin)茶树资源筛选及分子育种，具有很大的研究价值。

Description

一种与茶树表儿茶素含量连锁的分子标记位点及其应用

技术领域

本发明涉及分子遗传育种技术领域，更具体地，涉及一种与茶树表儿茶素(epicatechin)含量连锁的分子标记位点及其应用。

背景技术

茶(Camellia sinensis(L.)O.Kuntze)属于山茶科山茶属茶组，起源于中国的西南地区，距今有5000多年的栽培历史。茶叶与咖啡、可可并称为世界三大无酒精饮料，有着重要的经济价值，并对社会和文化有着重要影响。

茶树新稍中的特征性次级代谢产物儿茶素类化合物是茶叶滋味的主要影响因子。儿茶素化合物是2-苯基苯并吡喃的衍生物，属类黄酮化合物中的黄烷-3- 醇类，占茶叶干重的12％～24％。根据其B环的羟基数目、C环上2,3位同分异构体、C环上3位是否连接没食子基团，儿茶素类化合物可分为C((+)-儿茶素，(+) -catechin)、GC(没食子儿茶素，gallocatechin)、EGC(表没食子儿茶素，epi gallocatechin)、EC(表儿茶素，epicatechin)、EGCG(表没食子儿茶素没食子酸酯，epigallocatechin-3-gallate)、GCG(没食子儿茶素没食子酸酯，Gallocatec hin gallate)、ECG(表儿茶素没食子酸酯，epicatechin-3-gallate)和CG(儿茶素没食子酸酯，Catechin gallate)，它们与茶汤的苦涩味相关。

茶叶的次生代谢产物不仅影响茶叶品质，还具有多种生理功能。表儿茶素(epicatechin，EC)是一种天然植物黄烷醇化合物，呈白色晶体，易溶于水、甲醇。表儿茶素可以作为儿茶素配方颗粒、威麦宁胶囊、龙香膏、七厘散、蠲痹抗生丸等中药复方制剂的重要活性成分，并作为评价药品质量的重要指标。表儿茶素的分子结构有多个反应活性基团和活性部位，可使表儿茶素发生多种类酚反应，现代药理研究已表明：表儿茶素等类黄酮具有抗氧化、清除自由基、加强新陈代谢、调节免疫和抗肿瘤等功能，其中表儿茶素抗氧化能力被认为是分子中的酚羟基基团和自由基的结合达到清除自由基作用，而抗肿瘤则是影响肿瘤细胞的周期进程而抑制肿瘤细胞的周期生长，同时还有诸多光谱抑菌效果，被认为具有极大的疾病预防潜力。

目前茶树育种主要通过常规方法进行，从野生群体、杂交后代中选择优良单株进行系统选育。该方法时间长、效率低，使得新品种更新换代慢，不能快速满足大众对新产品的需求。分子标记辅助育种由于可以在苗期对育种材料进行选择，可显著提高育种效率。

发掘与茶树优良性状紧密连锁的分子标记是开展茶树分子标记辅助选择育种的基础，但是目前由于传统QTL定位研究进展的局限，一直未能找到影响表儿茶素含量的SNP分子标记位点。

发明内容

本发明的目的是为了克服现有技术的不足，提供一种与茶树表儿茶素含量连锁的分子标记位点及其应用。

本发明的第一个目的是提供一种茶树表儿茶素含量数量性状连锁的分子标记。

本发明的第二个目的是提供所述分子标记在评价茶树表儿茶素含量中的应用。

本发明的第三个目的是提供所述分子标记的引物。

本发明的第四个目的是提供所述引物在评价茶树表儿茶素含量中的应用。

本发明的第五个目的是提供一种评价茶树表儿茶素含量的试剂盒。

本发明的第六个目的是提供一种评价茶树表儿茶素含量的方法。

本发明的第七个目的是提供所述分子标记SNP位点、所述引物、或所述试剂盒任一在分子辅助育种中的应用。

为了实现上述目的，本发明是通过以下技术方案予以实现的：

发明人经过长期探索性的研究发现了一个与表儿茶素(epicatechin)连锁的 SNP位点分子标记。进一步利用其建立检测该位点的检测方法，可用于评价茶树的表儿茶素含量，以进一步用于资源筛选和分子育种。

因此本发明要求保护一种茶树表儿茶素含量数量性状连锁的分子标记，所述分子标记为位于茶树基因组Scaffold2233：468642的SNP位点，即SEQ ID NO:1 所示核苷酸序列的第501个碱基。

茶树基因组Scaffold2233：468642，该位点为G或者为A，其基因型与茶树干物质中表儿茶素含量极显著相关，通过相关性分析和显著性验证表明，AA基因型样本对应的茶汤干物质中表儿茶素含量与GG和GA基因型样本相比具有极显著差异。从统计学上判断，当基因型为双突变AA时，茶树中干物质中表儿茶素含量大概率高于正常平均水平的野生型GG或单突变基因型GA。本发明所述茶树表儿茶素含量具体为茶鲜叶干物质表儿茶素的比例。

所述分子标记SNP位点在评价茶树表儿茶素含量中的应用，也属于本发明的保护范围。

本发明还要求检测所述分子标记的引物，所述引物，其核苷酸序列如SEQ ID NO：2～3所示。

引物F：CAATGCATAATGCTCTACCC(SEQ ID NO：2)；

引物R：TGTGTGCGAATCGTTGAAGC(SEQ ID NO：3)。

所述引物在评价茶树表儿茶素含量中的应用，也属于本发明的保护范围。

进一步，本发明要求保护一种评价茶树表儿茶素含量的试剂盒，包括检测所述分子标记SNP位点的试剂。

优选地，所述试剂为所述引物，其核苷酸序列如SEQ ID NO：2～3所示。

最优选地，所述试剂盒含有核苷酸序列如SEQ ID NO：2～3所示的引物、 2×TaqPCR Master Mix、ddH₂O。

其使用方法为：

(1)采用CTAB法提取茶树嫩芽总DNA，并确保每个DNA样品的A260/A280在1.8～2.0之间，浓度大于100μg/μl；

(2)PCR扩增

PCR体系(10μl)如下：

2×Taq PCR Master Mix	5μl
		引物	各0.5μl
DNA template	1μl
		ddH<sub>2</sub>O	3μl

PCR扩增程序如下：

(3)产物纯化

将PCR扩增产物进行凝胶电泳，之后使用市售凝胶电泳DNA回收试剂盒进行回收纯化。

(4)测序及结果判读

将回收纯化的产物送测序公司进行Sanger法测序，在Scaffold2233：468642 位点，从统计学上判断，当基因型为双突变AA时，茶树中干物质中表儿茶素含量大概率高于正常平均水平的野生型GG或单突变基因型GA。

同时，本发明要求保护一种评价茶树表儿茶素含量的方法，检测所述的分子标记SNP位点的基因型。

优选地，利用所述引物检测所述的分子标记SNP位点的基因型。

所述分子标记、所述引物、所述试剂盒、或所述试剂盒在分子辅助育种中的应用，也属于本发明的保护范围。

与现有技术相比，本发明具有如下有益效果：

本发明首次发现了与茶树表儿茶素含量相关的SNP分子标记位点，其位于茶树基因组Scaffold2233：468642上，其基因型与表儿茶素含量极显著相关，AA 基因型样本对应的茶汤干物质中表儿茶素含量与GG和GA基因型样本相比具有极显著差异。从统计学上判断，当基因型为双突变AA时，茶树中干物质中表儿茶素含量大概率高于正常平均水平的野生型GG或单突变基因型GA。进一步建立检测该位点的检测方法，可用于评价茶树的表儿茶素含量，以进一步用于茶树资源筛选及分子育种。这是开展茶树分子标记辅助选择育种的基础，有很大的研究价值。

附图说明

图1为不同季节的表儿茶素含量。

图2为Scaffold2233：468642位点与引物示意图，N表示Scaffold2233：468642 位置上的待测碱基，加粗并下划线部分为上下游引物。

图3为样本2-73在Scaffold2233：468642位点的基因型SNaPshot测序结果。

图4为样本2-74在Scaffold2233：468642位点的基因型SNaPshot测序结果。

图5为样本2-87在Scaffold2233：468642位点的基因型SNaPshot测序结果。

图6为试剂盒检测Scaffold2233：468642位点基因型测序图。

具体实施方式

下面结合说明书附图和具体实施例对本发明作出进一步地详细阐述，所述实施例只用于解释本发明，并非用于限定本发明的范围。下述实施例中所使用的试验方法如无特殊说明，均为常规方法；所使用的材料、试剂等，如无特殊说明，为可从商业途径得到的试剂和材料。

实施例1

一、实验样本

采集位于广东省茶树种质资源库(广东，英德，113.3OE,24.3ON)的191 份茶树材料，其中广东124份、福建20份、广西15份，浙江9份，湖南6份，云南6份，江西1份，贵州1份，台湾1份。另外8份肯尼亚茶种后代，1份格鲁吉亚种后代，所选材料具有广泛代表性。

所选用的资源随机分布在资源库中。采用双行单株种植，每行4m，行距1.5 m，株距35cm。资源库进行常规水肥管理。在2016年年末对资源进行修剪并深坑施基肥，每亩4吨有机肥、0.75吨花生麸和10斤复合肥。2017年春茶及夏茶之后进行修剪并在根外追肥，每亩30斤复合肥和60斤尿素。分别在2017年的3月15日、6月25日和9月28日采摘茶树新梢一芽二叶，制作蒸青样，并按照水浸提法制备茶汤。

二、表型数据分析

1、实验步骤

利用高效液相色谱法对茶汤中与茶树滋味相关的表儿茶素，[epicatechin] (EC)进行检测，参照国标法进行检测。

利用SPSS软件对表儿茶素含量的指数大小范围、平均值、标准差和变异系数进行分析。将数量性状以0.5个标准差，将数据分为10个等级，用于计算性状的Shannon-Wiener多样性指数。使用最佳线性无偏预测(Best Liner Unbiased Prediction,BLUP)方法，采用一年多点模型估计育种值，同时估算广义遗传力。

2、实验结果

表儿茶素含量见表1.

表1不同季度不同资源EC占干物质的百分比：

群体的表儿茶素含量变异情况见表2和图1。

表2 EC性状(表儿茶素含量)表型变异：

三、基因型和性状关联分析

1、实验步骤

采用CTAB法提取191个茶树资源嫩芽总DNA，并确保每个DNA样品的 A260/A280在1.8～2.0之间，浓度大于100μg/μL。将提取好的DNA样品，检测分别位于“舒茶早”CSS栽培种茶树基因组(http：//tpia.teaplant.org/index.html)的 SNP位点(Scaffold2233：468642)的基因型，进行性状和标记的关联分析，关联的显著性水平以P值进行判断，p值小于1.25E-05为显著性水平。

2、实验结果

该SNP位点不同季节的P值见表3。

表3：不同季节Scaffold2233：468642位点的P值

实施例2分子标记在另外一个群体中的验证

一、实验方法

将位于Scaffold2233：468642的SNP位点在另一个含98个种质的群体中进行验证。

1、检测各样本的表儿茶素含量。具体检测方法同实施例1。

2、利用SnaPShot技术平台检测各样本的Scaffold2233：468642的SNP位点的基因型。

该方法针对不同突变位点设计不同长度的引物SNaPshot反应后，产物通过电泳分离、五色荧光检测、Gene mapper分析，可在一个测序反应中检测多个SNP 位点。应用SNaPshot进行定点的序列分析，其基本原理遵循了DNA直接测序中的双脱氧终止法，所不同的是PCR反应中只有不同荧光标记的ddNTP。由于每个SNP位点的引物3′端都紧靠SNP点，因此每一种引物在聚合酶作用下，根据模板的序列，只延伸一个核苷酸。然后用先进的荧光检测系统，检测延伸的那个核苷酸的种类。

(1)引物设计

根据Scaffold2233：468642在基因组的位置设计引物并进行合成。其中，Scaffold2233：468642上下游各延伸500bp。其核苷酸序列如SEQ ID NO：1所示(图2，其中N表示Scaffold2233：468642位置上的待测碱基)。

PCR引物：

F：CAATGCATAATGCTCTACCC(SEQ ID NO：2)；

R：TGTGTGCGAATCGTTGAAGC(SEQ ID NO：3)。

单碱基延伸引物：

actgactgactgactgactgactgactgCTTTCCATTCTACTCTGTTGCTGCTA。

(2)PCR扩增

PCR体系(10μl)如下：

2×Taq PCR Master Mix	5μl
		PrimerMix(按扩增情况配比)	1μl
DNA template	1μl
		ddH<sub>2</sub>O	3μl

PCR扩增程序如下：

(3)PCR产物纯化

使用虾碱酶纯化法进行纯化。虾碱酶MIX(EX-SAP)的主要功能成分为SAP 及ExoI.SAP酶，可以将残余dNTPs去磷酸化，ExoI降解游离单链引物。取4μl 的PCR产物，加入2μl的EX-SAP酶。具体反应体系如下所示：

消化体系组分	体积(μl)
		ddH<sub>2</sub>O	0.75
SAP(1U/μl)	0.5
		ExoI(5U/μl)	0.15
10*SAPbuffer	0.6
		PCR产物	4
总体积	6

之后在PCR仪上进行消化温育：37℃40min，85℃5min，4℃forever。

(4)SNaPshot反应

PCR产物当做模板进行SNaPshot反应。

SNaPshot反应体系如下所示：

试剂	用量(μl)
		SNaPshot Mix	0.5
Pooled PCR Products	3
		Pooled Primers	1
dH<sub>2</sub>O	0.5
		总体积	5

SNaPshot反应程序为：

之后，对SNaPshot产物进行纯化，直接向SNaPshot反应产物中加入2ul的 SAPmix，具体反应体系如下：

组分	体积(μl)
		水	0.9
SAP(1U/ul)	0.5
		10*SAP buffer	0.6
总计	2

在PCR仪上进行SNaPshot产物消化反应，反应程序为：37℃40min，75℃15 min，4℃forever。

(5)上机检测

取2μl消化后的SNaPshot反应产物加入到8μl含有0.4％LIZ120的去离子甲酰胺中，95℃变性5min，然后放-20℃骤冷，然后上3730XL测序。

(6)结果分析

用GeneMarker分析得到的.fsa结果，导出峰图及表格文件，并统计出每个样品的SNP突变型。

二、实验结果

各样本的表儿茶素含量及Scaffold2233：468642的SNP位点的基因型见表4，部分样本SNaPshot测序结果见图3到图5。

表4验证群体中资源EC(epicatechin)的占干物质的含量及基因型：

显著性分析结果显示，Scaffold2233：468642的基因型与表儿茶素含量极显著相关，相关系数为0.37，p-value为2.03×10^-4，F值(6.91/3.94)为14.9，为隐性突变，AA基因型样本对应的茶汤干物质中表儿茶素含量与GG和GA基因型样本相比具有极显著差异。从统计学上判断，当基因型为双突变AA时，茶树中干物质中表儿茶素含量大概率高于正常平均水平的野生型GG或单突变基因型GA。

实施例3一种评价茶树表儿茶素含量的试剂盒

一、组成

其核苷酸序列如SEQ ID NO：2～3所示的引物、2×Taq PCR Master Mix、 ddH₂O。

其中，引物F：CAATGCATAATGCTCTACCC(SEQ ID NO：2)；

引物R：TGTGTGCGAATCGTTGAAGC(SEQ ID NO：3)。

二、使用方法

(1)采用CTAB法提取茶树嫩芽总DNA，并确保每个DNA样品的 A260/A280在1.8～2.0之间，浓度大于100μg/μl；

(2)PCR扩增

PCR体系(10μl)如下：

2×Taq PCR Master Mix	5μl
		引物	各0.5μl
DNA template	1μl
		ddH<sub>2</sub>O	3μl

PCR扩增程序如下：

(3)产物纯化

将PCR扩增产物进行凝胶电泳，之后使用市售凝胶电泳DNA回收试剂盒进行回收纯化。

(4)测序及结果判读

将回收纯化的产物送测序公司进行Sanger法测序，将测序结果与SEQ ID NO：1所示核苷酸序列进行比对，根据图2所示(加粗并下划线部分为上下游引物)，Scaffold2233：468642位点位于扩增产物的第115个碱基。从统计学上判断，当基因型为双突变AA时，茶树中干物质中表儿茶素含量大概率高于正常平均水平的野生型GG或单突变基因型GA。

实施例4一种评价茶树表儿茶素含量的试剂盒的应用

一、实验方法

用实施例3的试剂盒检测实施例2中的98个茶树样本。

二、实验结果

检测结果与实施例2采用SnaPShot技术平台检测的结果一致，本试剂盒可以用于评价茶树表儿茶素含量。部分样本测序峰图如图6所示。

序列表

<110> 广东省农业科学院茶叶研究所

<120> 一种与茶树表儿茶素含量连锁的分子标记位点及其应用

<160> 3

<170> SIPOSequenceListing 1.0

<210> 1

<211> 1001

<212> DNA

<213> Camellia sinensis

<400> 1

ttttaattga gtgattaata tatatagtac agccgaccca tttaatggga tttaaggctt 60

ggtttggttt gattaatata tataggaaaa actttgctga acacattccc agaataaatt 120

aaccagacaa tctgaaatcc aagcccccca aaaatgcaag accatacaca caatttttct 180

gtattgctga gtctgctgac atcgcatttg ctcagttaca tgcctctcta ttttgttttt 240

gttttacttt tttaattttt ttcccccttt ctacttcttt tcatcttttc tgccaggaat 300

tcaatgtata aaagctatat ctagatgttt ttccaatcca atttagattg aaaaattgag 360

tcaatacacc ttttgtacca ggaattcaat gcataatgct ctacccaatt tagacaacac 420

aacatctatt taaaaggcta taaaaaatag caaacaaatt tgtctagaga atttaaagtt 480

gacatattta tggacaagca ntagcagcaa cagagtagaa tggaaagaaa tgacaatgaa 540

actcaccgtg cagtggtgag tggtgaacca gcacaccccc taggattcat gacaacagga 600

aaacaaccgc gcttcaacga ttcgcacaca aaccaccgaa tattcttatc catactccca 660

tcagcagttc cgggaacaag caaaagcgtg gtatccaatc cacgctcctc cgccaagtcc 720

aaattcgccg gccaatccaa agaaatgact cctccatcat cactactcac acacactctc 780

tgatacacaa gcttctcttc cgcctcctcc tcctcctcat cactgtccca aactcgaatc 840

ctcccactgt tcaacctcac gaaatgccta tcctccttca ctagcttctc attcacgtcc 900

tccaattcga acgatatcga aggacaacga agaatcacga aacggttaaa cggcgtggga 960

ctggtgaaca aaatccactc gccgatgtca gagagtcgag g 1001

<210> 2

<211> 20

<212> DNA

<213> Camellia sinensis

<400> 2

caatgcataa tgctctaccc 20

<210> 3

<211> 20

<212> DNA

<213> Camellia sinensis

<400> 3

tgtgtgcgaa tcgttgaagc 20

Claims (8)

1.一种茶树表儿茶素含量数量性状连锁的分子标记，其特征在于，所述分子标记为位于茶树基因组Scaffold2233：468642的SNP位点，即SEQ ID NO:1的第501个碱基。

2.权利要求1所述分子标记在评价茶树表儿茶素含量中的应用。

3.检测权利要求1所述分子标记的引物，其核苷酸序列如SEQ ID NO：2～3所示。

4.权利要求3所述引物在评价茶树表儿茶素含量中的应用。

5.一种评价茶树表儿茶素含量的试剂盒，其特征在于，包括检测权利要求1所述分子标记SNP位点的试剂。

6.根据权利要求5所述的试剂盒，其特征在于，所述试剂为权利要求3所述引物。

7.一种评价茶树表儿茶素含量的方法，其特征在于，检测权利要求1所述的分子标记SNP位点的基因型。

8.权利要求1所述分子标记、权利要求3所述引物、或权利要求5所述试剂盒任一在分子辅助育种中的应用。