CN113584204A - 与油茶种子出仁率相关的dna片段、其紧密连锁的snp分子标记及其应用 - Google Patents

Abstract

本发明涉及油茶分子标记及遗传育种技术领域，具体涉及与油茶种子出仁率相关的DNA片段、其紧密连锁的SNP分子标记及其应用。本发明提供与油茶种子出仁率相关的DNA片段，其为位于油茶连锁图谱6号连锁群的rkf.06‑2，置信区间为53.7cM‑64.5cM。本发明还提供与油茶种子出仁率位点紧密连锁的SNP分子标记chr06‑88890750，其为含有如SEQ ID NO.1所示序列第156位的多态性为C/T的核苷酸序列，可以解释14.9％的油茶种子出仁率的表型方差。通过检测该SNP分子标记，可在苗期进行油茶种子出仁率表型的鉴定和辅助筛选，大大节约生产成本，显著提高了油茶种子出仁率选择育种的效率。

Description

与油茶种子出仁率相关的DNA片段、其紧密连锁的SNP分子标 记及其应用

技术领域

本发明涉及油茶分子标记及遗传育种技术领域，具体涉及与油茶种子出仁率相关的DNA片段、其紧密连锁的SNP分子标记及其应用。

背景技术

油茶(Camellia oleifera)与油菜、花生、大豆共同被作为中国四大油料，是重要的经济作物，在中国亚热带地区广泛种植。油茶籽油是一种优质的食用油，其中不饱和脂肪酸的含量可达90％以上，且富含角鲨烯、维生素E等营养成分。选育高产(油)质优的油茶优良品种一直是油茶产业健康发展的基础和保障。长期以来，以选择和杂交育种为主要手段的油茶育种工作取得了较大的进展，但油茶的常规育种存在周期长、新品种选育缓慢等问题，良种选育速度仍然不能满足油茶产业发展的需求，已成为限制油茶产业发展的重要因素之一。

与传统育种技术相比，分子标记辅助育种可从苗期开始选择，大幅缩短育种的周期，对以果实为主要目的的经济林的育种优势尤其明显。有效的分子标记是分子标记辅助育种的关键，因此，开发与油茶果实产量、种子出仁率、种仁含油率、油脂品质表型相关的分子标记，对于油茶油脂产量和品质的分子标记辅助育种及相关性状的遗传改良具有重要意义。

果实产量、鲜果出籽率、种子出仁率、种仁含油率等指标直接决定单位面积油茶产(油)量。因此，开展油茶种子出仁率性状的研究和改良，是提高油茶产量的重要途径之一，对油茶产业的提升和健康发展具有十分重要的意义。

发明内容

本发明的目的之一在于提供与油茶种子出仁率相关的DNA片段(基因位点)及与其紧密连锁的SNP分子标记，本发明的另一目的是提供所述分子标记在油茶高出仁率表型鉴定和育种中的应用。

本发明提供的油茶种子出仁率基因位点及该基因位点紧密连锁的分子标记的开发是基于已建立的油茶F1代杂交群体，开展高密度遗传连锁图构建和出仁率性状的QTL定位实现的。油茶的简化基因组序列为本发明分子标记开发的区域。

本发明中关键基因位点及连锁的SNP分子标记的开发过程基本如下：

(1)对种子出仁率差异显著的油茶无性系长林53号(47.66％)和长林81号(58.14％)控制授粉，创制种子出仁率广泛分离的油茶F1代杂交群体。

(2)采集杂交群体180个单株的完全成熟种子，测定种子出仁率。

(3)采集杂交群体180个单株及两个亲本的幼嫩叶片，采用TaKaRa MiniBEST植物基因组DNA提取试剂盒(TaKaRa,Dalian,China)提取DNA，利用EcoRI和NlaIII(Hin1II)双酶切后每样本分别构建简化基因组(ddRAD)测序文库，利用Illumina HiSeqXten平台测序。

(4)以二倍体油茶基因组为参考序列，分析(3)中获得的180个样本及2个亲本简化基因组的SNP位点。SNP数据根据以下原则过滤：亲本测序深度≥10X，子代测序深度≥8X；基因型缺失率≤30％；SNP质量值≥30。上述过程中使用的软件BWA是公开免费的。

(5)运用Joinmap4.0软件进行连锁图谱构建，参数设置为：Rec≤0.4,LOD≥3.0，Jump＝5，作图函数使用Kosambi；分析连锁群内标记的排列顺序，计算相邻标记间的遗传距离。

(6)运用QTL Icimapping软件进行数据分析，采用完备区间作图法(ICIM)进行QTL定位。扫描步长设置为1cM；逐步回归标记进入的概率(PIN)为0.002(POUT＝2*PIN＝0.002)；LOD值为2.5。

利用上述技术手段，本发明获得了一个位于6号连锁群(LG06)上的油茶种子出仁率基因位点rkf.06-2，该基因位点的对油茶种子出仁率的贡献率为14.9％，与该位点紧密连锁的SNP标记为chr06-88890750，其基因型为C/T(表1)。

表1基因位点及连锁SNP分子标记信息

具体地，本发明提供如下技术方案：

第一方面，本发明提供与油茶种子出仁率相关的DNA片段，其为位于油茶连锁图谱6号连锁群的rkf.06-2，置信区间为53.7cM-64.5cM。

上述DNA片段对油茶种子出仁率表型的贡献率为14.9％，可用于图位克隆和分子标记辅助选择。

与上述DNA片段紧密连锁的SNP分子标记为chr06-88890750，所述chr06-88890750位于油茶基因组6号染色体的88890750bp位，多态性为C/T。

第二方面，本发明提供与油茶种子出仁率位点紧密连锁的SNP分子标记，其包括SNP分子标记chr06-88890750，该SNP分子标记位于油茶基因组6号染色体的88890750bp位，多态性为C/T。

该SNP分子标记与油茶种子出仁率基因位点rkf.06-2紧密连锁。

具体地，SNP分子标记chr06-88890750含有如SEQ ID NO.1所示序列第156位的多态性为C/T的核苷酸序列。

优选地，SNP分子标记chr06-88890750的核苷酸序列如SEQ ID NO.1所示，SNP位点位于如SEQ ID NO.1所示第156位，多态性为C/T。

进一步地，SNP分子标记chr06-88890750可以由核苷酸序列如SEQ ID NO.2-3所示的引物对以油茶基因组DNA为模板经PCR扩增获得。

SEQ ID NO.2：5’-AAGGAAGTTCAGGGAATAAGT-3’；

SEQ ID NO.3：5’-ATTAGCTACAACTAACTAATA-3’。

SNP分子标记chr06-88890750中，具有所述多态性的位点的基因型为C/C，对应于高出仁率；基因型为C/T，对应于低出仁率。

本发明所述的出仁率为烘干至恒重的油茶种子出仁率。

第三方面，本发明提供用于扩增所述SNP分子标记的引物。

作为本发明的一种实施方式，所述引物包括如SEQ ID NO.2-3所示的引物。

本发明还提供含有所述引物的试剂或试剂盒。

第四方面，本发明提供所述与油茶种子出仁率相关的DNA片段或所述SNP分子标记或所述引物或所述试剂或试剂盒的以下任一应用：

(1)在鉴定油茶种子出仁率表型中的应用；

(2)在油茶种子出仁率改良或油茶种子出仁率分子标记辅助育种中的应用；

(3)在油茶种子出仁率的早期预测中的应用；

(4)在筛选高出仁率油茶中的应用。

第五方面，本发明提供一种鉴定油茶种子出仁率表型或筛选高种子出仁率油茶的方法，包括如下步骤：

(1)提取待鉴定油茶的基因组DNA；

(2)以基因组DNA为模板，利用序列如SEQ ID NO.2-3所示的引物进行PCR扩增；

(3)分析PCR扩增产物中所述SNP分子标记的基因型，根据所述基因型判断待鉴定油茶的种子出仁率表型。

上述方法的步骤(1)中，所述待鉴定油茶具体选自长林53号×长林81号的杂种F1代个体。

步骤(2)中，所述PCR扩增的反应程序为：94～95℃，3～5min；94～95℃，15～30s，65～69℃，40～60s，38～45个循环；67～70℃，3～6min。优选为，95℃，3min，1个循环预变性；95℃，15s变性，68℃，45s延伸，40个循环；68℃，5min，1个循环彻底延伸。

步骤(2)中，在所述扩增后，通过琼脂糖凝胶电泳检测并回收所得到的PCR产物。

步骤(3)中，分析SNP分子标记的基因型可采用本领域常规技术手段，例如测序等，可以SEQ ID NO.2-3为测序引物进行测序。

步骤(3)中所述判断待鉴定油茶的种子出仁率表型的方法为：

若SNP分子标记chr06-88890750具有所述多态性的位点的基因型为C/C，则待鉴定油茶为高出仁率；若基因型为C/T，则待鉴定油茶为低出仁率。

本发明提供鉴定高出仁率油茶的方法，包括：

(1)提取待鉴定油茶的基因组DNA；

(2)以DNA为模板，利用SEQ ID NO.2-3所示的引物进行PCR扩增；

(3)分析PCR扩增产物中所述SNP分子标记的基因型，根据所述基因型判断待鉴定油茶是否为高出仁率油茶。

若SNP分子标记chr06-88890750具有所述多态性的位点的基因型为C/C，则待鉴定油茶为高出仁率；若基因型为C/T，则待鉴定油茶为低出仁率。

本发明的有益效果在于：本发明提供了一个油茶种子出仁率基因位点，对种子出仁率表型的贡献率为14.9％，并开发了与基因位点紧密连锁的SNP分子标记。利用该SNP分子标记对长林53号×长林81号的杂种F1代个体进行了辅助选择，结果表明，该位点为高出仁率基因型的单株中，70.42％的个体的种子出仁率高于群体种子出仁率平均值(59.16％)；低出仁率基因型的个体中，69.10％的个体的种子出仁率低于群体平均值(59.16％)，表明该标记用于种子出仁率表型的辅助选择是切实有效的。

在油茶常规选择育种中，种子出仁率性状的鉴定需要幼苗造林5-6年才能鉴定，费时费力。本发明中的SNP位点位置明确，检测方法方便快速，不受环境影响，目的性更强，工作量小，效率更高，成本低。因此，通过检测所述SNP位点，可在苗期进行鉴定和辅助筛选，大大节约生产成本和提高选择效率。在油茶育种中，可选择本发明的分子标记及其检测方法鉴定出高出仁率油茶进行育种，可提高油茶育种的选择效率，加快育种进程。

附图说明

图1为本发明实施例3中油茶种子出仁率基因位点rkf.06-2的位置示意图。

具体实施方式

以下实施例用于说明本发明，但不用来限制本发明的范围。

以下实施例中所用的长林53号×长林81号的杂种F1代单株，均由中国林业科学研究院亚热带林业研究所木本油料育种与培育研究组收集、评价，并保存于浙江金华婺城区东方红林场种质资源圃。

实施例1油茶种子出仁率分离群体的构建及性状测定

本实施例中使用长林53号和长林81号分别作为母本和父本，采用控制授粉技术，创制了经济性状广泛分离的杂种F1代群体。F1群体采用随机区组3次重复设计保存在浙江省婺城区东方红林场。180个子代个体待果实完全成熟后(5％果实开裂)，采集种子，测定种子出仁率。其操作步骤如下：

(1)每样本取10-15粒种子，105℃烘至恒质量，记录种子质量(W₁)。

(2)油茶种子剥去硬种皮，记下种仁的总质量(W₂)。

(3)计算每个样本的种子出仁率：

种子出仁率＝W₂/W₁×100％。

油茶种子出仁率测定结果表明：杂交群体种子出仁率显著分离，说明该性状具有数量性状特点。

实施例2油茶连锁图谱构建

1、基因组DNA提取

于3月份采集长林53号×长林81号家系180个个体及其双亲本的春稍嫩叶，采用KAC法(TaKaRa试剂盒Code No.9768)提取基因组总DNA。具体步骤如下：

(1)首先在1.5ml离心管中加入500μl的Buffer HS II；准确称取100mg油茶嫩叶进行液氮研磨；迅速将研磨好的粉末加入到离心管中混匀，然后加入10μl的RNase A(10mg/ml)，充分振荡混匀，于56℃水浴温育10分钟。

(2)加入62.5μl的Buffer KAC，充分混匀。冰上放置5分钟，12,000rpm离心5分钟。取上清，加入与上清液等体积的Buffer GB，充分混匀。

(3)将Spin Column安置于Collection Tube，溶液移至Spin Column中(由于溶液较多，一般需要分两次过柱，每次过柱的体积量不要超过700μl)，12,000rpm离心1分钟，弃滤液。

(4)将500μl的Buffer WA加入至Spin Column中，12,000rpm离心1分钟，弃滤液。

(5)将700μl的Buffer WB加入至Spin Column中，12,000rpm离心1分钟，弃滤液。

(6)重复操作步骤(5)。

(7)将Spin Column安置于Collection Tube上，12,000rpm离心2分钟。

(8)将Spin Column安置于新的1.5ml的离心管上，在Spin Column膜的中央处加入30～50μl的Elution Buffer或灭菌蒸馏水，室温静置5分钟。

(9)12,000rpm离心2分钟洗脱DNA。

2、dd-RAD简化基因组测序

各样本基因组DNA经过最佳酶切组合EcoRI和NlaIII双酶切后连接接头，接头包括3部分，分别为测序引物、分子识别序列(barcode)和与内切酶酶切基因组后产生的粘性末端互补的序列。然后对每个样品进行进行PCR扩增。扩增程序：98℃、2min；98℃、30s，60℃、30s，72℃、15s，13个循环；72℃、5min。PCR产物用2％琼脂糖凝胶电泳，300～500bp长度的片段用AxyPrep DNA凝胶回收试剂盒从凝胶中回收、纯化。将带有不同barcode的纯化后的DNA样本，每12个个体混池为一个样品，构建dd-RAD测序文库，利用Illumina HiSeqXten平台测序。

3、SNP位点识别和基因分型

(1)测序数据过滤，测序获得的原始序列数据首先按照如下步骤过滤：

1)根据序列上的barcode，用一个自定义的Perl脚本快速的将12个样本的混合数据按个体分开；

2)具有barcode，且barcode后面紧跟内切酶的识别位点的序列保留，其余序列舍弃；

3)缺失的核苷酸数＞3个的序列，舍弃；

4)其他低质量的、污染的序列使用NGS QC Toolkit软件包(Patel R K,JainM.NGS QC Toolkit:A Platform for Quality Control of Next-Generation SequencingData[M].Springer US,2015.)进一步过滤。

(2)SNP识别与过滤：将每个样本的高质量的reads比对到参考基因组序列上。剔除没有比对上的序列，其余序列识别SNP位点。识别的SNP位点经过严格过滤，获得高质量的SNPs数据。过程中用到软件Tophat v2.1.1、bcftools v1.9和BWA是公开免费的。SNPs过滤标准如下：

1)亲本测序深度≥10X，子代测序深度≥8X；

2)基因型缺失率≤30％；

3)SNP质量值≥30。

4、遗传图谱构建

(1)标记分离模式检测：运用JoinMap4.0软件中的CP功能，对检测到的所有SNP进行分析，标记分离比率用卡方检验计算，确定各个标记的分离模式，如ab×cd,ef×eg,hk×hk,lm×ll,nn×np，cc×ab，ab×cc等，过滤、剔除显著不正常分离(P<0.05)或含有非正常碱基的标记。ef×eg、hk×hk、lm×ll和nn×np四种类型的标记用于后续连锁图构建。

(2)遗传连锁图构建：用JoinMap4.0软件构建连锁图，参数设置为：Rec≤0.4,LOD≥3.0，Jump＝5，作图函数使用Kosambi；分析连锁群内标记的排列顺序，计算相邻标记间的遗传距离。构建的连锁图共有15个连锁群，上图标记2780个，共覆盖3327cM，平均间距为1.20cM。

实施例3油茶种子出仁率基因位点及连锁SNP位点挖掘

运用QTL Icimapping软件进行数据分析，采用完备区间作图法(ICIM)进行种子出仁率QTL定位。扫描步长设置为1cM；逐步回归标记进入的概率(PIN)为0.002(POUT＝2*PIN＝0.002)；LOD值为2.5。LOD显著性阈值通过运行1000次permutation检验确定。在油茶6号染色体上定位到油茶种子出仁率基因位点rkf.06-2，贡献率为14.9％，与其紧密连锁的SNP分子标记为chr06-88890750(表1，图1)。

实施例4出仁率基因位点与连锁SNP分子标记在油茶育种中的应用

1、从长林53号×长林81号油茶杂交F1代家系群体中随机挑选124个单株为材料，采集嫩叶提取总DNA(见实施例2)。DNA溶液稀释100倍，作为工作液。

2、利用SEQ ID NO.2-3所示的引物对对DNA工作液进行PCR扩增，反应体系如表2所示：

表2 PCR扩增的反应体系

PCR扩增程序为：

3、PCR扩增产物进行凝胶检测和纯化回收并测序、基因分型。凝胶检测和纯化回收按照AxyPrep DNA凝胶回收试剂盒(AxyGEN,Code No.AP-GX-50)说明书进行，其流程如下：

(1)配制1.2％的琼脂糖凝胶，将50μl扩增产物全部上样，电泳电压为5V/cm，电泳约20分钟至上样缓冲液中二甲苯青达到距离凝胶前端1cm处时停止电泳。

(2)在紫外灯下切下含有目的DNA的琼脂糖凝胶，用纸巾吸尽凝胶表面的液体并切碎。计算凝胶重量，该重量作为一个凝胶体积(例如100mg＝100μl体积)。

(3)加入3个凝胶体积的Buffer DE-A，混合均匀后于75℃加热，每2～3分钟间断混合，直至凝胶块完全熔化。

(4)加入0.5个Buffer DE-A体积的Buffer DE-B，混合均匀。

(5)将上述溶液转移到DNA制备管中，12000rpm离心1分钟，弃滤液。

(6)加入500μl Buffer W1，12000rpm离心30秒，弃滤液。

(7)加入700μl Buffer W2，12000rpm离心30秒，弃滤液。以同样的方法再用700μlBuffer W2洗涤一次，12000rpm离心1分钟，弃滤液。

(8)将制备管放回离心管中，12000rpm离心1分钟。

(9)将制备管置于洁净的1.5ml离心管中，在制备膜中央加25～30μl去离子水，室温静置1分钟。12000rpm离心1分钟洗脱DNA。

(10)凝胶回收DNA，以对应的扩增引物为测序引物，采用一代测序测定扩增产物核苷酸序列，用Chromas软件判读测序峰图上SNP位点的基因型。

4、分别鉴定所有个体的chr06-88890750位点的基因型。对照各位点的基因型与出仁率高低的关系，若该位点的基因型为C/C，则该油茶个体为高出仁率油茶；若基因型为C/T，则该油茶个体为低出仁率油茶。

5、采集124个F1代个体及双亲本完全成熟种子，测定其种子出仁率(见实施例1)。结果表明(表3)，chr06-88890750位点为高出仁率基因型的单株中，70.42％的个体其出仁率高于群体种子出仁率平均值(59.16％)；基因型为低出仁率基因型的个体中，69.10％的个体其种子出仁率低于群体平均值(59.16％)。这表明该标记用于种子出仁率表型的辅助选择是切实有效的，可用于油茶种子出仁率表型的早期鉴别或辅助鉴别，可大大节约生产成本，提高选择效率，加快油茶育种进程。

表3母本长林53号、父本长林81号及F1单株的种子出仁率及基因型数据

虽然，上文中已经用一般性说明及具体实施方案对本发明作了详尽的描述，但在本发明基础上，可以对之作一些修改或改进，这对本领域技术人员而言是显而易见的。因此，在不偏离本发明精神的基础上所做的这些修改或改进，均属于本发明要求保护的范围。

序列表

<110> 中国林业科学研究院亚热带林业研究所

<120> 与油茶种子出仁率相关的DNA片段、其紧密连锁的SNP分子标记及其应用

<130> KHP211117126.2

<160> 3

<170> SIPOSequenceListing 1.0

<210> 1

<211> 317

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 1

aaggaagttc agggaataag tggtcaaccg gccaccagta gtggctgatt ggccacctgt 60

attagccgac cggccaccaa tggtagccga cagccaccat tggtagccga tcggccagaa 120

atgaaatttt gcacaagtgt gttcctgagt tggagcagag gtgttttgaa agttggcagc 180

actgttggtg ttcaagcaag tgcaaagaaa ggctgcaaag gcaaggcagt ggcattttcg 240

ggaggagaat tcaagcatgc atgcatccct cgtccctcga ggccgtttaa ttgctttatt 300

agttagttgt agctaat 317

<210> 2

<211> 21

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 2

aaggaagttc agggaataag t 21

<210> 3

<211> 21

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 3

attagctaca actaactaat a 21

Claims (10)

1.与油茶种子出仁率相关的DNA片段，其特征在于，其为位于油茶连锁图谱6号连锁群的rkf.06-2，置信区间为53.7cM-64.5cM。

2.与油茶种子出仁率位点紧密连锁的SNP分子标记，其特征在于，包括SNP分子标记chr06-88890750，所述chr06-88890750含有如SEQ ID NO.1所示序列第156位的多态性为C/T的核苷酸序列。

3.根据权利要求2所述的SNP分子标记，其特征在于，由序列如SEQ ID NO.2-3所示的引物对以油茶基因组DNA为模板经PCR扩增获得。

4.根据权利要求2或3所述的SNP分子标记，其特征在于，所述SNP分子标记chr06-88890750中，具有所述多态性的位点的基因型为C/C，对应于高出仁率；基因型为C/T，对应于低出仁率。

5.用于扩增权利要求2～4任一项所述的SNP分子标记的引物。

6.根据权利要求5所述的引物，其特征在于，包括序列如SEQ ID NO.2-3所示的引物。

7.含有权利要求5或6所述的引物的试剂或试剂盒。

8.权利要求1所述的与油茶种子出仁率相关的DNA片段或权利要求2～4任一项所述的SNP分子标记或权利要求5或6所述的引物或权利要求7所述的试剂或试剂盒的以下任一种应用：

(1)在鉴定油茶种子出仁率表型中的应用；

(2)在油茶种子出仁率改良或油茶种子出仁率分子标记辅助育种中的应用；

(3)在油茶种子出仁率的早期预测中的应用；

(4)在筛选高出仁率油茶中的应用。

9.一种鉴定油茶种子出仁率表型或筛选高种子出仁率油茶的方法，其特征在于，包括如下步骤：

(1)提取待鉴定油茶的基因组DNA；

(2)以基因组DNA为模板，利用序列如SEQ ID NO.2-3所示的引物进行PCR扩增；

(3)分析PCR扩增产物中权利要求2～4任一项所述的SNP分子标记的基因型，根据所述基因型判断待鉴定油茶种子出仁率表型；

优选地，步骤(2)中，所述PCR扩增的反应程序为：94～95℃，3～5min；94～95℃，15～30s，65～69℃，40～60s，38～45个循环；67～70℃，3～6min。

10.根据权利要求9所述的方法，其特征在于，步骤(3)中，所述判断待鉴定油茶种子出仁率表型的方法为：

若SNP分子标记chr06-88890750具有所述多态性的位点的基因型为C/C，则待鉴定油茶为高出仁率；若基因型为C/T，则待鉴定油茶为低出仁率。

Images