CN111534632A

CN111534632A - 与油茶种仁含油率相关的3个snp分子标记及其应用

Info

Publication number: CN111534632A
Application number: CN202010479322.4A
Authority: CN
Inventors: 林萍; 王开良; 任华东
Original assignee: Research Institute of Subtropical Forestry of Chinese Academy of Forestry
Current assignee: Research Institute of Subtropical Forestry of Chinese Academy of Forestry
Priority date: 2020-05-29
Filing date: 2020-05-29
Publication date: 2020-08-14
Anticipated expiration: 2040-05-29
Also published as: CN111534632B

Abstract

本发明涉及分子标记技术领域，具体涉及与油茶种仁含油率相关的3个SNP分子标记及其应用。本发明提供了3个与油茶种仁含油率高度关联的SNP分子标记PB.103793.1‑1118、PB.103793.1‑1172和PB.103793.1‑1175，可以分别解释29.07％、29.00％和28.75％的含油率表型方差。通过检测所述的SNP分子标记，可在苗期进行鉴定和辅助筛选，大大节约生产成本和提高选择效率。

Description

与油茶种仁含油率相关的3个SNP分子标记及其应用

技术领域

本发明涉及分子标记技术领域，具体涉及与油茶种仁含油率相关的3个SNP分子标记及其应用。

背景技术

油茶(Camellia oleifera Abel.)隶属山茶科(Theaceae)山茶属(Camellia L.)，是木本油料树种。油茶籽油含有丰富的营养物质，是一种优质的食用油，其不饱和脂肪酸含量达90％以上，以油酸和亚油酸为主。油茶籽油具有抗氧化、抗肿瘤、降血脂等功效，具有较高的营养保健价值。目前，以选择和杂交育种为主要手段、以产果量为主要育种目的的油茶育种已取得了重要进展，但以提高种仁含油率、改良油脂品质为目的的育种研究仍较少。油茶的常规育种周期长，新品种选育缓慢，良种选育速度还不能满足产业发展的需求，这已成为限制油茶产业发展的重要因素之一。

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

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

发明内容

本发明的目的之一在于提供与油茶种仁含油率相关的SNP分子标记，本发明的另一目的是提供所述SNP分子标记在油茶种仁含油率表型鉴定和育种中的应用。

本发明提供的与油茶种仁含油率关联位点的开发方法是基于油茶是典型的异交物种，连锁不平衡(LD)通常在较小范围内迅速消减，因此可以开展重要性状的LD作图。油茶的种仁全部转录本作为本发明标记开发的区域。在具备产生了大量明显的遗传变异的油茶自然群体的前提下，可有效开展与油茶种仁含油率变异显著相关的标记开发。

本发明中SNP分子标记的开发过程基本如下：

(1)在油茶全分布区内广泛收集油茶种质资源，建立种仁含油率广泛分离的油茶自然群体。

(2)采集自然群体500份油茶种质的完全成熟种子，用索氏抽提法测定种仁含油率。

(3)采集自然群体500个油茶单株的油脂高速合成期的种仁，采用RNAprep Pure多糖多酚植物总RNA提取试剂盒(离心柱型，TIANGEN试剂盒Code No.DP441)提取总RNA，每样本分别构建cDNA文库，利用Illumina HiSeq^TM 4000平台进行二代转录组测序。

(4)采集油茶“长林4号”的根、嫩叶、成熟叶片、花瓣和未成熟种子，采用RNAprepPure多糖多酚植物总RNA提取试剂盒(离心柱型，TIANGEN试剂盒Code No.DP441)分别提取RNA，各组织RNA等比例混合，构建PacBio SMRTbell文库，在PacBio Sequel平台进行三代转录组测序。测序结果过滤掉低质量数据和冗余序列后，对所有转录本进行注释分析。过程中用到软件LoRDEC(http://www.atgc-montpellier.fr/lordec/)、CD-HIT v4.6(Fu L,NiuB,Zhu Z,Wu S,Li W,2012.CD-HIT:accelerated for clustering the next-generationsequencing data.Bioinformatics 28,3150-2.)、Coding Potential Calculator(CPC)(Kong L,Zhang Y,Ye Z-Q,et al.,2007.CPC:assess the protein-coding potential oftranscripts using sequence features and support vector machine.Nucleic AcidsResearch 35,W345.)和Coding-Non-Coding Index(CNCI)(https://github.com/www-bioinfo-org/CNCI)等是免费公开的。

(5)以(4)中获得的全长转录组序列为参考序列，采用多序列比对法，分析(3)中获得的500个样本转录组序列的SNP位点。SNP数据根据以下原则严格过滤：每个位点只有2个等位基因；基因型缺失率≤20％；最小等位基因频率≥5％；SNP质量值≥100；纯合基因型样本数超过10个；杂合基因型率≤70％。过程中用到软件bcftools v1.9software(http://www.htslib.org/doc/bcftools.html)是公开免费的。

(6)将群体的基因型数据输入GCTA v1.25.2(Jian Y,S Hong L,Goddard ME,Visscher PM,2011.GCTA:a tool for genome-wide complex trait analysis.AmericanJournal of Human Genetics 88,76-82.)软件，进行主成分分析(PCA)。

(7)将群体的基因型数据、前10个主成分(PC)数据、种仁含油率的表型数据以及Kinship矩阵数据输入TASSEL5.0(http://www.maizegenetics.net/tassel)软件中，采用统一混合线性模型方法(MLM)分析SNPs标记和油茶种仁含油率性状的连锁不平衡性。

利用上述技术措施，本发明最终获得了与油茶种仁含油率极显著关联(P<10^-5)的位于同一个转录本内的3个SNP标记PB.103793.1-1118、PB.103793.1-1172和PB.103793.1-1175(具体见表1)，它们对表型变异的贡献率分别为29.07％、29.00％和28.75％(表1)。

表1 SNP分子标记信息

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

第一方面，本发明提供与油茶种仁含油率相关的SNP分子标记，所述SNP分子标记选自PB.103793.1-1118、PB.103793.1-1172和PB.103793.1-1175中的任一个；

其中，SNP分子标记PB.103793.1-1118含有如SEQ ID NO.3所示序列第1118位的多态性为A/G的核苷酸序列；SNP分子标记PB.103793.1-1172含有如SEQ ID NO.3所示序列第1172位的多态性为G/A的核苷酸序列；SNP分子标记PB.103793.1-1175含有如SEQ ID NO.3所示序列第1175位的多态性为A/C的核苷酸序列。

进一步地，本发明的与油茶种仁含油率相关的SNP分子标记可以由核苷酸序列如SEQ ID NO.1-2所示的引物对以油茶cDNA为模板经PCR扩增获得。

所述SNP分子标记PB.103793.1-1118中，具有所述多态性的位点的基因型为A/A，对应于高含油率；基因型为A/G，对应于候选高含油率；基因型为G/G，对应于低含油率；

所述SNP分子标记PB.103793.1-1172中，具有所述多态性的位点的基因型为G/G，对应于高含油率；基因型为A/G，对应于候选高含油率；基因型为A/A，对应于低含油率；

所述SNP分子标记PB.103793.1-1175中，具有所述多态性的位点的基因型为A/A，对应于为高含油率；基因型为A/C，对应于候选高含油率，基因型为C/C，对应于低含油率。

本发明所述的含油率为种仁的含油率。

各SNP分子标记可分别单独用于或联合用于油茶种仁含油率表型的鉴定，联合使用时鉴定的准确率更高。

本发明还提供与油茶种仁含油率相关的SNP分子标记组合，其包括PB.103793.1-1118、PB.103793.1-1172和PB.103793.1-1175中的至少两个；

其中，SNP分子标记PB.103793.1-1118、PB.103793.1-1172和PB.103793.1-1175同上所述。

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

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

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

第三方面，本发明提供所述的SNP分子标记或所述的SNP分子标记组合或所述的引物或所述的试剂或试剂盒的以下任一应用：

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

(2)在油茶种质资源鉴定、改良或分子标记辅助育种中的应用，所述油茶种质资源鉴定、改良或分子标记辅助育种的性状为油茶种仁含油率；

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

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

在利用本发明提供的3个SNP分子标记进行油茶种仁中含油率表型鉴定或分子标记辅助育种时，本领域技术人员根据需要可以选择其中任意一个分子标记或多个分子标记的组合。多个联合使用，鉴定的准确率更高。

第四方面，本发明提供鉴定油茶种仁含油率表型的方法，包括：

(1)提取待鉴定油茶的总RNA，反转录合成cDNA；

(2)以cDNA为模板，利用上述的引物进行PCR扩增；

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

上述方法的步骤(1)中，所述待鉴定油茶可以为任何育种材料，包括自然群体个体和有性群体个体。

提取油茶总RNA采用RNAprep Pure多糖多酚植物总RNA提取试剂盒(离心柱型，TIANGEN试剂盒Code No.DP441)。反转录合成单链cDNA采用PrimeScript RT Master Mix试剂盒(TaKaRa,大连，中国)。

步骤(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产物。

作为一种实施方案，所述琼脂糖凝胶电泳中，琼脂糖凝胶的浓度为1.2％。胶回收使用AxyPrep DNA凝胶回收试剂盒(AxyGEN,Code No.AP-GX-50)。

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

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

若SNP分子标记PB.103793.1-1118具有所述多态性的位点的基因型为A/A，则待鉴定油茶为高含油率；若基因型为A/G，则待鉴定油茶为候选高含油率；若基因型为G/G，则待鉴定油茶为低含油率；和/或，

若SNP分子标记PB.103793.1-1172具有所述多态性的位点的基因型为G/G，则待鉴定油茶为高含油率；若基因型为A/G，则待鉴定油茶为候选高含油率；若基因型为A/A，则待鉴定油茶为低含油率；和/或，

若所述SNP分子标记PB.103793.1-1175具有所述多态性的位点的基因型为A/A，则待鉴定油茶为高含油率；若基因型为A/C，则待鉴定油茶为候选高含油率，若基因型为C/C，则待鉴定油茶为低含油率。

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

(1)提取待鉴定油茶的总RNA，反转录合成cDNA；

(2)以cDNA为模板，利用上述的引物进行PCR扩增；

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

若SNP分子标记PB.103793.1-1118具有所述多态性的位点的基因型为A/A，则待鉴定油茶为高含油率；若基因型为A/G，则待鉴定油茶为候选高含油率；和/或，

若SNP分子标记PB.103793.1-1172具有所述多态性的位点的基因型为G/G，则待鉴定油茶为高含油率；若基因型为A/G，则待鉴定油茶为候选高含油率；和/或，

若所述SNP分子标记PB.103793.1-1175具有所述多态性的位点的基因型为A/A，则待鉴定油茶为高含油率；若基因型为A/C，则待鉴定油茶为候选高含油率。

本发明的有益效果在于：本发明开发了3个与油茶种仁含油率高度关联的SNP位点，可以分别解释29.07％、29.00％和28.75％的含油率表型方差。利用这3个标记对有性油茶群体进行了辅助选择，结果表明，PB.103793.1-1118位点的基因型为A/A的单株中，78.35％的个体其种仁含油率高于群体种仁含油率平均值；PB.103793.1-1172位点的基因型为G/G的单株中，79.17％的个体其种仁含油率高于群体平均值；PB.103793.1-1175位点的基因型为A/A的单株中，79.17％的个体其种仁含油率高于群体平均值。这表明该标记用于辅助选择是切实有效的。

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

具体实施方式

以下实施例用于说明本发明，但不用来限制本发明的范围。若未特别指明，实施例中所用的技术手段为本领域技术人员所熟知的常规手段。

以下实施例中所用的自然群体材料500份单株，均由中国林业科学研究院亚热带林业研究所木本油料育种与培育研究组收集、评价，并保存于浙江金华婺城区东方红林场种质资源圃。

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

本实施例中使用普通油茶资源收集圃内500份种质资源的自然群体，其起源地涵盖我国油茶主产区的大部分，包括浙江省、湖南省、江西省、广西区、福建省、广东省等。500个体待果实完全成熟后(5％果实开裂)，采集种子，采用索氏抽提法测定种仁含油率。其操作步骤如下：

(1)制备中速滤纸包，并放入铝盒内，105℃烘至恒质量，记录铝盒和滤纸包质量(W₁)。

(2)适量油茶种子剥去硬种皮，105℃烘至恒质量，用粉碎机粉碎后，装入滤纸包内包扎好，记下铝盒、滤纸包和试样的总质量(W₂)。

(3)采用瑞士Buchi索氏提取仪B-811LSV，将称好质量的试样滤纸包放入浸提瓶中，加入约100ml石油醚，提取6h，回收石油醚，将滤纸包(内有残渣)放入铝盒内于105℃烘至恒质量，记下铝盒、滤纸包和残渣的质量(W₃)。

种仁含油率＝[(W₂-W₃)/(W₂-W₁)]×100％

油茶种仁含油率测定结果表明：自然群体种仁含油率呈正态分布，说明该性状具有数量性状特点。

实施例2油茶三代转录组测序及注释分析

1、三代测序样本RNA的提取：

采集油茶“长林4号”的根、嫩叶、成熟叶片、花瓣和未成熟种子，采用RNAprep Pure多糖多酚植物总RNA提取试剂盒(离心柱型，TIANGEN试剂盒Code No.DP441)分别提取RNA，具体步骤如下：

(1)首先在1.5ml离心管中加入500μl的裂解液SL(使用前检查是否已加入β-巯基乙醇)。取0.1g样本材料加入液氮充分研磨，迅速将研磨好的样本粉末加入到离心管中，立即旋涡剧烈震荡混匀。

(2)12000rpm离心2分钟。

(3)将上清液转移至过滤柱CS上(过滤柱CS放在收集管中)，12000rpm离心2分钟，小心吸取收集管中的上清至新的RNase-Free的离心管中，吸头尽量避免接触收集管中的细胞碎片沉淀。

(4)缓慢加入0.4倍上清体积的无水乙醇，混匀(此时可能会出现沉淀)，将得到的溶液和沉淀一起转入吸附柱CR3中，12000rpm离心15秒，倒掉收集管中的废液，将吸附柱CR3放回收集管中。

注意，若果上清液体积有损失，请相应调整乙醇的加量。

(5)向吸附柱CR3中加入350μl去蛋白液RW1，12000rpm离心15秒，倒掉收集管中的废液，将吸附柱CR3放回收集管中。

(6)DNase I工作液的配制：取10μl DNase I储存液放入新的RNase-Free的离心管中，加入70μl RDD缓冲液，轻柔混匀。

(7)向吸附柱CR3中央加入80μl的DNase I工作液，室温放置15分钟。

(8)向吸附柱CR3中加入350μl去蛋白液RW1，12000rpm离心15秒，倒掉收集管中的废液，将吸附柱CR3放回收集管中。

(9)向吸附柱CR3中加入500μl漂洗液RW(使用前请先检查是否已加入乙醇)，12000rpm离心15秒，倒掉收集管中的废液，将吸附柱CR3放回收集管中。

(10)重复步骤9。

(11)12000rpm离心2分钟，将吸附柱CR3放入一个新的RNase-Free的离心管中，向吸附膜的中间部位悬空滴加30-50μl RNase-Free ddH₂O，室温放置2分钟，12000rpm离心1分钟，得到RNA溶液。

2、三代转录组测序与注释分析：

通过纯度和浓度检测的五个组织样品总RNA，等比例混合，利用Clontech

PCR cDNA合成试剂盒进行反转录合成单链cDNA。利用KAPA HiFi PCR试剂盒以单链cDNA为模板进行第一轮PCR扩增，生成双链cDNA。生成的双链cDNA利用Blue Pippin分成0.5-2kb，2-3kb，3-6kb三个长度片段库。然后进行第二轮PCR扩增产生充足的cDNA，构建PacBio SMRTbell文库，在PacBio Sequel平台进行三代转录组测序。利用SMRTlink 5.0软件处理测序数据。测序结果过滤掉低质量数据和冗余序列后，生成CCS。根据序列是否含有5’primer，3’primer和polyA尾巴，将所有的CCS分为全长和非全长序列两大类。全长的CCS采用ICE算法在缺省参数条件下进行聚类分析产生CS。利用Arrow和LoRDEC(http://www.atgc-montpellier.fr/lordec/)软件进一步过滤CS，利用CD-HIT v4.6(Fu L,Niu B,Zhu Z,Wu S,Li W,2012.CD-HIT:accelerated for clustering the next-generationsequencing data.Bioinformatics 28,3150-2.)软件去除冗余序列。

利用Coding Potential Calculator(CPC)(Kong L,Zhang Y,Ye Z-Q,et al.,2007.CPC:assess the protein-coding potential of transcripts using sequencefeatures and support vector machine.Nucleic Acids Research 35,W345.)和Coding-Non-Coding Index(CNCI)(https://github.com/www-bioinfo-org/CNCI)软件在参数缺省条件下进行转录本的蛋白编码潜能预测。未通过蛋白编码潜能检测的转录本进一步在Swiss-Prot数据库中比对，若在Swiss-Prot数据库中仍无注释，则认为该转录本为长链非编码RNA。其它的转录本进一步在NR、Swiss-Prot、COG、KEGG和GO等数据库中比对，注释转录本。

实施例3油脂高速合成期种仁转录组测序及多态位点识别

1、500个油茶无性系油脂高速合成期种仁总RNA提取：

利用RNAprep Pure多糖多酚植物总RNA提取试剂盒(离心柱型，TIANGEN试剂盒Code No.DP441)分别提取各无性系未成熟种仁的总RNA(见实施例2)。

2、二代转录组测序：

通过纯度和浓度检测的各样品总RNA，去除其中的核糖体RNA，以最大限度地保留所有coding RNA和ncRNA。得到的RNA随机打断成短片段，再以片断化后的RNA为模板，用六碱基随机引物(random hexamers)合成cDNA第一链；接着加入缓冲液、dNTPs(dUTP代替dTTP)、RNase H和DNA polymerase I合成cDNA第二链，经过QiaQuick PCR试剂盒纯化并加EB缓冲液洗脱，经末端修复、加碱基A，加测序接头，然后通过UNG(Uracil-N-Glycosylase)酶降解第二条链。用琼脂糖凝胶电泳进行片段大小选择，进行PCR扩增。最后建好的测序文库用Illumina HiSeq^TM 4000平台进行二代转录组测序。

3、多态位点识别：

为了保证数据质量，对下机后经过初步过滤得到的clean reads进行进一步更严格的过滤，得到高质量的clean reads，用于后续的信息分析。过滤的步骤如下：

(1)去除含有接头的reads；

(2)去除全部都是A碱基的reads；

(3)去除含N比例大于10％的reads；

(4)去除低质量的reads(质量值Q≤20的碱基数占整条reads的50％以上)。

采用Tophat v2.1.1(Trapnell C,Roberts A,Goff L,et al.,2012.Differential gene and transcript expression analysis of RNA-seqexperiments with TopHat and Cufflinks.Nature protocols 7,562-78.)软件将每个样本的高质量的reads比对到参考转录组序列上(见实施例2)。剔除没有比对上的序列，其余序列利用bcftools v1.9软件(http://www.htslib.org/doc/bcftools.html)识别SNP位点。识别的SNP位点经过严格过滤，获得高质量的SNPs数据。过滤标准如下：

(1)位点上只有2个等位基因；

(2)基因型缺失率≤20％；

(3)最小等位基因频率(MAF)≥5％；

(4)SNP质量值≥100；

(5)纯合子基因型的样本数大于10个；

(6)杂合基因型样本率≤70％。

实施例4与油茶种仁含油率相关的SNP位点的筛选

1、群体结构分析：

利用GCTA v1.25.2(Jian Y,S Hong L,Goddard ME,Visscher PM,2011.GCTA:atool for genome-wide complex trait analysis.American Journal of HumanGenetics 88,76-82.)软件对油茶自然群体进行主成分分析(PCA)，利用前10个主成(PC)作为固定效应用于后续的关联分析(表2)。

表2自然群体部分个体的前10个PC值

2、关联分析：

将所有样本的SNPs位点数据、前10个PC值数据、表型数据(见实施例1)及Kinship矩阵数据导入TASSEL5.0软件中，采用MLM法分析SNPs与种仁含油率性状的连锁不平衡性，筛选与种仁含油率显著关联的分子标记，经多重检验校正，检测到了3个跟含油率存在极显著关联的位点(P<10^-5，表1)，这3位点位于PB.103793.1转录本内，均引起非同义突变，对含油率差异的贡献率分别为29.07％、29.00％和28.75％(表1)。

实施例5本发明中的3个分子标记在油茶高油育种中的应用

(1)选择一个油茶杂交F1代家系群体为材料(母本为‘长林53号’，父本为‘长林40号’，均为国家审定良种，良种号分别为“国S-SC-CO-012-2008”和“国S-SC-CO-011-2008”)，采集嫩叶提取总RNA(见实施例2)。以RNA为模板，采用Clontech cDNA合成试剂盒反转录生成单链cDNA，并稀释100倍，作为工作液。

(2)利用SEQ ID NO.1-2所示的引物对对单链cDNA工作液分别进行PCR扩增，反应体系如表3所示：

表3

PCR扩增程序为：

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

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

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

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

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

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

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

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

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

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

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

(4)分别鉴定所有个体的3个位点的基因型。对照各位点的基因型与含油率高低的关系，若PB.103793.1-1118位点的基因型为A/A，则该油茶个体为高含油油茶；若PB.103793.1-1172位点的基因型为G/G，则该油茶个体为高含油油茶；若PB.103793.1-1175位点的基因型为A/A，则该油茶个体亦为高含油率的油茶。

(5)采集所有F1代个体完全成熟种子，测定其种仁的含油率(见实施例1)。结果表明(表4)，PB.103793.1-1118位点的基因型为A/A的单株中，78.35％的个体其种仁含油率高于群体种仁含油率平均值(37.79％)；PB.103793.1-1172位点的基因型为G/G的单株中，79.17％的个体其种仁含油率高于群体平均值(37.79％)；PB.103793.1-1175位点的基因型为A/A的单株中，79.17％的个体其种仁含油率高于群体平均值(37.79％)。这表明该标记用于辅助选择是切实有效的，可用于早期鉴别或辅助鉴别，可大大节约生产成本，提高选择效率，加快油茶高油育种进程。

表4 F1单株的种仁含油率及基因型数据

注：表中“..”表示基因型缺失。

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

序列表

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

<120> 与油茶种仁含油率相关的3个SNP分子标记及其应用

<130> KHP201112304.4

<160> 3

<170> SIPOSequenceListing 1.0

<210> 1

<211> 20

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 1

aaacaccaca cacacacagc 20

<210> 2

<211> 20

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 2

aaagtttttg taaacataaa 20

<210> 3

<211> 1859

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 3

aaacaccaca cacacacagc aaggaaaaat gaaaagcatt ccagagatgt tactgtgttg 60

tagttctgat cataagccaa ttcctcttgg ggggtttgga acagctgttt atcctcttgc 120

atcctctgaa accatgaaac aatccttcct ccatgcaatc aaacttggtt acagacactt 180

cgactctgca gctttatacc agtcagagca gcctcttgga gaatcaattg ttgatgccat 240

acgcctaggc ttcattcaat ctcgccaaga cctcttcatc acctctaagc tttggtgttc 300

tgatgctcac cctcatcatg tcctccctgc tcttcaaagt tcactcaagt aattcataac 360

tgctgctctt ttaatgttat tatgctccgt aaatctttct ctttcttttt tatcagcgaa 420

agtcgcttgt tatgctccat aaatcattat ttatattctt ttgttatata ctagtttttt 480

tataaattat ttttttaatc cctctaattt gatgattttg atatacttat caatctcttg 540

gcttctatac gtaaaaaaaa taaacgaaaa tttaaacata atgatcaaat tgaaacggac 600

gtaaagtgca aggatttaag gacaaaaaat aaaaggaaga attacactga catgaatatc 660

aataaaaaaa tattatgaat gttaataatt aatgttataa actttttttt tattaattta 720

taattttttc tctctaacac cactatcagg tagttttccc aaaatataat tttgcttagt 780

tggttttttt tttttaatta ttatttatta tagaaatctt ggattggaat accttgatct 840

gtatctcatt cactggccag tgagttcaaa gccagggaaa tttgagtatc cagtgaacaa 900

gcaagagctt cttcccatgg atttcaagtc tgtttgggaa gccatggagg agtgtcagaa 960

tcttggcctc acaaaattca ttggagtcag taacttctca tgcaagaagc tccaactatt 1020

actagcaacc gcaaagatcc ctccagctgt caaccaggtc gagatgaacc cactttggca 1080

acagaagaag ctaagagagt tttgtgagaa aaaaggtatt catatcacag cttactctcc 1140

tttgggcgcc aaaggaacaa tttggggaag tgacaaagtc atggaatgtg aggtgctcaa 1200

acagattgcc aaagctagag gaaaatctgt tgcccagtca ccggttaatc tagagggaaa 1260

gaaaaaggga attaccccca attaatgccc catattagct tctctaaaat taatggatgc 1320

ttcctaagga tcaatcacgg atgagcacat tgaactattc atgtgattga gagcccacac 1380

agcccaagca caacgacgca attatcacgg cgcgctctac cactgaacta ataactcatc 1440

gtgcgagcct cctagtgggg gctcactatg ccaaaagcaa gagaaacccc catccctctc 1500

tttccttttt ttggccctca aattctatcc ttggtttgtc tcagatgggt ttatgagcaa 1560

ggggtgagtg ttttggtgaa gagcttcagt gaggagagga tgaaagagaa ccttcaaata 1620

gttgattggg agctaggcgc acaagactcc gagatgataa atcaaatttc acaatataaa 1680

ggatgtgctg gacttgattt catatcagat gaaggccctt acaaatctct ccaggattta 1740

tgggatggtg aaattgtttg atcctgtaaa cgtgtagcca aaaaccactt ggataccgtt 1800

tgataacatt ttatgcttac aacacaaatg aatgtgtgtt ttatgtttac aaaaacttt 1859

Claims

1.与油茶种仁含油率相关的SNP分子标记，其特征在于，所述SNP分子标记选自PB.103793.1-1118、PB.103793.1-1172和PB.103793.1-1175中的任一个；

2.根据权利要求1所述的SNP分子标记，其特征在于，所述SNP分子标记PB.103793.1-1118中，具有所述多态性的位点的基因型为A/A，对应于高含油率；基因型为A/G，对应于候选高含油率；基因型为G/G，对应于低含油率；

3.与油茶种仁含油率相关的SNP分子标记组合，其特征在于，所述SNP分子标记组合包括PB.103793.1-1118、PB.103793.1-1172和PB.103793.1-1175中的至少两个；

其中，SNP分子标记PB.103793.1-1118、PB.103793.1-1172和PB.103793.1-1175同权利要求1或2中所述。

4.用于扩增权利要求1或2所述的SNP分子标记或权利要求3所述的SNP分子标记组合的引物。

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

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

7.权利要求1或2所述的SNP分子标记或权利要求3所述的SNP分子标记组合或权利要求4或5所述的引物或权利要求6所述的试剂或试剂盒的以下任一应用：

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

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

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

8.鉴定油茶种仁含油率表型的方法，其特征在于，包括：

(1)提取待鉴定油茶的总RNA，反转录合成cDNA；

(2)以cDNA为模板，利用权利要求4或5所述的引物进行PCR扩增；

(3)分析PCR扩增产物中权利要求1或2所述的SNP分子标记或权利要求3所述的SNP分子标记组合的基因型，根据所述基因型判断待鉴定油茶的种仁含油率表型。

9.根据权利要求8所述的方法，其特征在于，步骤(2)中，所述PCR扩增的反应程序为：94～95℃，3～5min；94～95℃，15～30s，65～69℃，40～60s，38～45个循环；67～70℃，3～6min。

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

若SNP分子标记PB.103793.1-1175具有所述多态性的位点的基因型为A/A，则待鉴定油茶为高含油率；若基因型为A/C，则待鉴定油茶为候选高含油率，若基因型为C/C，则待鉴定油茶为低含油率。