CN111676307B

CN111676307B - 与油茶种子油脂中软脂酸含量相关的snp分子标记及其应用

Info

Publication number: CN111676307B
Application number: CN202010477521.1A
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: 2023-03-21
Anticipated expiration: 2040-05-29
Also published as: CN111676307A

Abstract

本发明涉及分子标记技术领域，具体涉及与油茶种子油脂中软脂酸含量相关的SNP分子标记及其应用。本发明提供了1个与油茶软脂酸含量高度关联的SNP分子标记PB.34365.1‑621，可以解释9.66％的软脂酸含量表型方差。通过检测该SNP分子标记，可在苗期进行油茶油脂中软脂酸含量的鉴定和辅助筛选，大大节约生产成本和提高选择效率。

Description

与油茶种子油脂中软脂酸含量相关的SNP分子标记及其应用

技术领域

本发明涉及分子标记技术领域，具体涉及与油茶种子油脂中软脂酸含量相关的SNP分子标记及其应用。

背景技术

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

相比于传统育种手段，分子标记辅助育种可从苗期开始选择，大幅缩短育种的进程，对以果实为主要目的的经济林育种优势尤其明显。分子标记辅助育种离不开有效的分子标记，因此，开发与油茶籽油品质表型相关的分子标记，对于油茶籽油品质的分子标记辅助育种以及油茶籽油品质的改良具有重要意义。

软脂酸(palmitic acid，PA)，学名为“十六烷酸”，又称棕榈酸，分子式CH₃(CH₂)₁₄COOH，是一种广泛存在于植物油脂中的饱和中长链脂肪酸。目前普遍认为，饱和脂肪酸的过量摄入是造成高血压、高血脂、高胆固醇、非酒精性脂肪肝等高发疾病的重要原因之一。以降低植物油脂中的饱和脂肪酸含量，提高不饱和脂肪酸含量为主要目标，开展食用植物油品质的遗传改良一直是科学研究的重要方向。研究油茶油脂中软脂酸含量的遗传调控规律，开展低软脂酸含量油茶良种选育研究，对于提高茶油的营养价值，促进人类饮食健康具有重要意义。

发明内容

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

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

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

(1)在油茶全分布区内广泛收集油茶种质资源，建立油脂中软脂酸含量发生广泛分离的油茶自然群体。

(2)采集自然群体500份油茶种质的完全成熟种子，用气象色谱法测定成熟种子油脂的6种脂肪酸成分含量，包括硬脂酸、软脂酸、油酸、亚油酸、亚麻酸、二十碳烯酸，具体方法按照GB/T 17376《动植物油脂脂肪酸甲酯制备》和GB/T17377《动植物油脂脂肪酸甲脂的气相色谱分析》执行。

(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.9 software(http://www.htslib.org/doc/bcftools.htm1)是公开免费的。

(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标记和油茶油脂中软脂酸含量性状的连锁不平衡性，检测到1个位点与软脂酸含量极显著关联(P＝7.88×10^-6)，对表型变异的贡献率为9.66％(表1)。

表1 SNP分子标记信息

利用上述技术措施，本发明最终获得了与油茶种子油脂中软脂酸含量极显著关联的SNP标记PB.34365.1-621，这个标记位于油茶PB.34365.1转录本(核苷酸序列如SEQ IDNO.3所示)的621bp处，碱基为C/C、C/T或T/T。该位点的突变形成编码蛋白的第177个氨基酸错义突变(表1)。

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

第一方面，本发明提供与油茶种子油脂中软脂酸含量相关的SNP分子标记，包括PB.34365.1-621；

SNP分子标记PB.34365.1-621含有如SEQ ID NO.3所示序列第621位的多态性为C/T的核苷酸序列。

进一步地，本发明的与油茶油脂中软脂酸含量相关的SNP分子标记可以由引物P1：5’-GACATAGACACAAATCTAGA-3’(SEQ ID NO.1)和P2：5’-CCAAACGGGGCCTTACATCAG-3’(SEQID NO.2)以油茶cDNA为模板经PCR扩增获得。

所述SNP分子标记PB.34365.1-621中，具有所述多态性的位点的基因型为T/T，对应于低软脂酸含量，基因型为C/C或C/T，对应于候选高软脂酸含量。

本发明所述的软脂酸含量为种子油脂中的软脂酸含量。

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

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

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

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

(1)在鉴定油茶种子油脂中软脂酸含量表型中的应用；

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

(3)在油茶种子油脂中软脂酸含量的早期预测中的应用；

(4)在筛选低软脂酸含量油茶中的应用。

第四方面，本发明提供鉴定油茶种子油脂中软脂酸含量表型的方法，包括：

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

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

(3)分析PCR扩增产物中所述的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.34365.1-621具有所述多态性的位点的基因型为T/T，则待鉴定油茶为低软脂酸含量；若基因型为C/C或C/T，则待鉴定油茶为候选高软脂酸含量。

本发明提供鉴定低软脂酸含量油茶的方法，包括：

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

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

(3)分析PCR扩增产物中所述的SNP分子标记的基因型，根据所述基因型判断待鉴定油茶是否为低软脂酸含量油茶；若SNP分子标记PB.34365.1-621具有所述多态性的位点的基因型为T/T，则待鉴定油茶为低软脂酸含量油茶。

本发明的有益效果在于：本发明开发了1个与油茶种子油脂中软脂酸含量高度关联的SNP位点，可以解释9.66％的软脂酸含量表型方差。利用该标记对有性油茶群体进行了辅助选择，结果表明，该位点基因型为T/T的单株中，85.19％的个体其种子油脂中软脂酸含量低于群体软脂酸含量平均值。这表明该标记用于辅助选择是切实有效的。

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

具体实施方式

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

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

实施例1油茶种子油脂中软脂酸含量分离群体的构建及性状测定

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

(1)适量油茶种子用烤箱80℃烘烤过夜至恒重，剥去硬种皮。

(2)种仁用粉碎机粉碎后，用中速滤纸包好，加入适量石油醚浸泡抽提过夜。

(3)待石油醚完全挥发后，利用Agilent6890N气相色谱仪按照GB/T17376-2008、GB/T17377-2008方法测定脂肪酸成分及含量。

脂肪酸成分含量测定结果表明：自然群体种子油中软脂酸含量均呈正态分布，说明该性状具有数量性状特点。

实施例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:a tool 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与软脂酸含量性状的连锁不平衡性，筛选与油脂中软脂酸含量显著关联的分子标记，经多重检验校正，检测到了一个跟软脂酸含量存在极显著关联的位点PB.34365.1-621(具体见表1)，其基因型为C/C、C/T或T/T，该标记引起非同义突变，对软脂酸含量差异的贡献率为9.66％(表1)。

实施例5分子标记PB.34365.1-621在低软脂酸油茶育种中的应用

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

(2)利用P1和P2引物对单链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)分别鉴定所有个体的PB.34365.1-621位点的基因型。若该位点的基因型为T/T，则油茶个体为低软脂酸含量油茶；若该位点基因型为C/C或C/T，则油茶个体为高软脂酸含量或候选高软脂酸含量油茶。

(5)采集所有F1代个体完全成熟种子，测定其种子油脂中软脂酸含量(见实施例1)。结果表明(见表4)，该位点基因型为T/T的单株中，85.19％的个体其种子油脂软脂酸含量低于群体软脂酸含量平均值(8.3353％)。这表明该标记用于辅助选择，特别是筛选低软脂酸含量的油茶个体是切实有效的，可用于早期鉴别或辅助鉴别，可大大节约生产成本，提高选择效率，加快油茶高油育种进程。

表4 F1单株的种子油脂中软脂酸含量及PB.34365.1-621位点基因型数据

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

序列表

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

<120> 与油茶种子油脂中软脂酸含量相关的SNP分子标记及其应用

<130> KHP201112312.3

<160> 3

<170> SIPOSequenceListing 1.0

<210> 1

<211> 20

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 1

gacatagaca caaatctaga 20

<210> 2

<211> 21

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 2

ccaaacgggg ccttacatca g 21

<210> 3

<211> 1989

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 3

gacatagaca caaatctaga gtgagaaatg caaaatccaa aacagcagaa gaggagggca 60

ggcattagat tagagagaga gagagagaga gaatggggaa ggggccattc tgttttggcc 120

ggaagatgag ccgtcgccgc cctaagaagt ctagacagac ggcggcgacg gcgcctccgt 180

cgccgaataa taagtcgatg tcggcgatgg tggctgggtc gtcgagcaag acgaagaaga 240

aggcaggagg ggcgaggctg tggatgaggt tggataagtt tgggcaatcg gagctcatcg 300

agtgcgacag gagcacaatc atcaagcgcg tctctattcc tgttagggat ttgaggatcc 360

ttggcccttt cttctctcac tcctccaata tccttgctag ggagaaagct atggttgtca 420

acttagagtt tataaaggct atagttactg ctgaagaagt gctgttgctt gatcctctct 480

gccaggaggt tcttccgttt gtggatcagc tgcggcaaca acttcctcat aaaagcgtgt 540

ccaaaagtca aggagggggt caaatggatg tacaagataa tgaaatgcaa ttatcatctg 600

ggggacaatg gttacccgtt cctgaagccg ttgaaggttt gcagtgtgac ctcccatttg 660

agtttcaagt tctcgagatt gcattgaagg ttgtctgtac atatttggac tctagtgtag 720

cagaccttga gagagatgct taccctgtat tggatgaact ggcaaggaat gttagcacca 780

agaatcttga gcacgtgagg agtttgaaaa gcaatctgac ccgtttgctt gcacgggtac 840

aaaaggtgcg tgatgaaatt gaacatcttt tagatgacaa tgaagacatg gcacatttat 900

acttaacaag gaaatggatt cagagtcagc aatttgaggc ttatttagca ggtacagctt 960

ctaacagcat tgtcagtgct ggacatcaac ttcgccgtct tagttctaac aggagtggaa 1020

gtttattgtt gagcaaccat ttgaatgacc atgatgtgga agatctagaa atgttgctcg 1080

aggcttattt catgcagctg gatggtaccc gtgacaagat actatctgta agtgttatta 1140

ttattgacaa tagcatgagc ataatccaaa ataaactaat tacattaaac catcttatga 1200

cattgttcct acattgaatt ttcactgtac attgtatata ttgaaggatg gcaagctaat 1260

ttcatagaat taaggttaag cattttgtga agttggatta ttcaacccaa ccaaccacct 1320

atgtgactct aggaaacaca actttggttg attgcagaaa tattgcagct gctgtacgat 1380

cttctgagtc aatttctata gagacaaaag tcaatagctc ctcacaaatt tgtcctttct 1440

atctgtagag cactttccat ggaacggttc ttgatgctct gtggattgta gaaacaaatg 1500

ttcccacccc accaccatcc caccatcact aatacaacca cttgtacaat tgtttcaagg 1560

acgcaaccaa acataggaaa atgctatttg tgaaaaacat tttccgaaaa aaattctgcg 1620

ggaattattt tttatcaaaa caaacgtaat gctagtttca tgcaaatgag aagtttgaag 1680

ttcagtatta atctttcaag tttttgctag gcaagacaac gtggacactt gagcctctaa 1740

tattgaggtc atgaaattaa ttggatccac ctgccttcaa gaccattggt aagatagatc 1800

ttaattgggc aatgacagtt ggtagaaaga aatcaattga gtcccactta gcaacgttat 1860

gaaggctact tcattgagat tttcattctg cctctccatt cttgccccat ccagacatag 1920

tctgtgtatt tagtttcata ttttgagatg gaagaatgtc gtttaaagct gatgtaaggc 1980

cccgtttgg 1989

Claims

1.与油茶种子油脂中软脂酸含量相关的SNP分子标记，其特征在于，所述SNP分子标记为PB.34365.1-621；

SNP分子标记PB.34365.1-621如SEQ ID NO.3所示，其中，第621位的多态性为C/T。

2.根据权利要求1所述的SNP分子标记，其特征在于，所述SNP分子标记PB.34365.1-621中，具有所述多态性的位点的基因型为T/T，对应于低软脂酸含量，基因型为C/C或C/T，对应于候选高软脂酸含量。

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

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

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

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

（1）在鉴定油茶种子油脂中软脂酸含量表型中的应用；

（2）在油茶种质资源鉴定、改良或分子标记辅助育种中的应用，所述油茶种质资源鉴定、改良或分子标记辅助育种的性状为油茶种子油脂中软脂酸含量；

（3）在油茶种子油脂中软脂酸含量的早期预测中的应用；

（4）在筛选低软脂酸含量油茶中的应用。

7.鉴定油茶种子油脂中软脂酸含量表型的方法，其特征在于，包括：

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

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

（3）分析PCR扩增产物中权利要求1或2所述的SNP分子标记的基因型，根据所述基因型判断待鉴定油茶种子油脂中的软脂酸含量表型。

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

9.根据权利要求7或8所述的方法，其特征在于，步骤（3）中所述判断待鉴定油茶种子油脂中的软脂酸含量表型的方法为：

10.鉴定低软脂酸含量油茶的方法，其特征在于，包括：

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

（3）分析PCR扩增产物中权利要求1或2所述的SNP分子标记的基因型，根据所述基因型判断待鉴定油茶是否为低软脂酸含量油茶；若SNP分子标记PB.34365.1-621具有所述多态性的位点的基因型为T/T，则待鉴定油茶为低软脂酸含量油茶。