油菜千粒重关联的位点qSW.C3-3的PARMS分子标记或标记组合的应用的制作方法

油菜千粒重关联的位点qsw.c3-3的parms分子标记或标记组合的应用
技术领域
1.本发明属于生物育种技术领域，具体涉及油菜千粒重关联的位点qsw.c3-3的parms分子标记或标记组合的应用。


背景技术：

2.在同一种植密度条件下，油菜单产取决于单株产量，而单株产量直接决定于全株角果数、每角粒数和粒重三个构成因子。虽然油菜产量的三个构成因子之间呈现一定程度的负相关，但其系数往往不大，因此可以通过提高产量构成因子(如粒重)来增加产量。从育种实践上看，油菜单产的增加确实能通过提高粒重等构成因子来实现。目前，油菜品种的千粒重一般都低于4g，而油菜种质资源中千粒重的最大值可以超过8克(陈苇等，2011)，说明油菜粒重仍有很大提升空间。
3.随着现代生物技术的飞快发展，分子标记技术也在不断更新迭代，由最初的rflp、ra pd/aflp、ssr逐渐发展为目前的snp标记。目前snp检测方法主要分为两大类：一大类是以单链构象多态性(sscp)、变性梯度凝胶电泳(ddge)、酶切扩增多态性序列(caps)、等位基因特异性pcr(as-pcr)等为代表的以凝胶电泳为基础的传统经典的检测方法；另一大类是以直接测序、dna芯片、变性高效液相色谱(dhplc)、质谱检测技术、高分辨率溶解曲线(hrm)等为代表的高通量、自动化程度较高的检测方法。这两类方法各有优缺点，目前已有将两者结合起来的检测方法如kasp(lister et.,2013)和parms(lu et al.,2020)等，不但操作简便、成本低廉而且能实现高通量的检测。
4.粒重是一个典型的由多基因控制的复杂数量性状，表型连续分布而且容易受环境条件的影响。粒重的遗传以加性效应为主，显性和上位性较弱，因而其杂种优势很弱，具体表现为杂种f1的粒重一般都介于双亲之间(李娜，2015)。数量遗传学和分子标记技术的结合，可将复杂的数量性状分解为单个的数量性状基因位点(quantitative trait loci,qtl)，然后像研究质量性状一样对控制数量性状的多个基因进行研究。qtl定位就是在遗传分离群体的基础上，借助分子标记和遗传图谱，利用qtl作图软件对分离群体的数量性状表型数据进行分析，从而确定数量性状基因在染色体上的位置和效应。目前，利用连锁作图(linkage ma pping)和关联作图(association mapping)的方法，在油菜种已定位到一百多个千粒重qtl(quijada et al.,2006；udall et al.,2006；radoev et al.,2008；shi et al.,2009；王峰等，2010；basunanda et al.,2010；fan et al.,2010；zhang et al.,2011；赵卫国等，2017)，但在不同群体和环境中都起作用的qtl数目有限，难以满足实际的育种需求。
5.本发明利用油菜自然群体10个环境的千粒重数据进行全基因组关联分析，旨在找到对油菜千粒重具有改良效应的新的稳定位点，并据此开发高通量低成本的实用性分子标记用于辅助选择，为后续全基因组选择育种奠定基础。


技术实现要素：

6.本发明的目的是在于提供了油菜c03染色体上和粒重显著关联的位点qsw.c3-3，其峰值snp标记bn-scaff_22067_1-p219321位于油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基。
7.本发明的另一个目的在于提供了与油菜千粒重关联的位点组合qsw.a3-3和qsw.c3-3，其峰值snp标记bn-a03-p14176469和bn-scaff_22067_1-p219321分别位于油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基。
8.本发明还有一个目的是提供了针对上述峰值snp的parms标记引物。
9.本发明的最后一个目的在于提供了上述parms标记引物在油菜千粒重选择育种中的应用。
10.为了达到上述目的，本发明采取以下技术措施：
11.与油菜千粒重相关的snp位点的获得：
12.(1)提取油菜关联群体331份品系的总dna，利用油菜60k snp芯片对每个样本进行基因型分析(li et al.,2020)。
13.(2)利用illumina beadstudio基因分型软件(http://www.illumina.com/)计算群体材料在每个位点的缺失率(missing rate)、杂合率(heterozygous rate)、次要等位基因频率(minor allele frequency)。
14.(3)将关联群体331份品系分别在10个环境种植，成熟期挑选每个小区中长势均匀的10个单株进行收获，按单株分别脱粒后将种子晒干进行千粒重的考察，计算各品系平均值。
15.(4)利用关联群体的高密度snp基因型数据、10个环境的千粒重表型数据、群体结构和亲缘关系，利用tassel 5.0软件(bradbury et al.,2007)进行全基因组关联分析。最终获得与粒重显著关联的位点qsw.a3-3和qsw.c3-3，其峰值snp标记bn-a03-p14176469和bn-scaff_22067_1-p219321位于油菜darmorv4.1参考基因组a03染色体第13373806和c03染色体第20819936位碱基处，均可以在两个环境中被重复检测到。
16.(5)提取油菜a03染色体第13,373,806位碱基上下游各100bp的序列，按照parms引物设计原则获得引物序列：qsw.a3-3f：gaattaagtagaagctcgcaaacac、qsw.a3-3ra：gaaggtgaccaagttcatgcttgttgatggttgacttctggtga、qsw.a3-3rg：gaaggtcggagtcaacggatttgttgatggttgacttctggtgg。
17.提取油菜c03染色体第20,819,936位碱基上下游各100bp的序列，按照parms引物设计原则获得引物序列：qsw.c3-3f：atttcattcttcgttatagcaggc、qsw.c3-3ra：gaaggtgaccaagttcatgctgagcaaatctggtcatagaaagaga、qsw.c3-3rg：gaaggtcggagtcaacggattgagcaaatctggtcatagaaagagg。
18.本发明的保护内容包括：
19.检测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的系统在油菜千粒重选择育种中的应用，所述的系统为试剂、试剂盒或仪器。
20.油菜千粒重选择育种的装置，其特征在于：所述的装置包括数据接收模块和数据处理模块，所述的数据接收模块被配置为接收待测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的基因型；所述的数据处理模块用于将来自于接受模块的基因型转换
为待测油菜的判定结果。
21.油菜千粒重选择育种的计算机可读存储介质，其特征在于：所述计算机可读存储介质使计算机运行如下步骤：检测待测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的基因型，将所述基因型转换成判定待测油菜的判定结果。
22.以上所述的应用中，所述的试剂为检测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的引物。
23.以上所述应用中，所述的引物为：qsw.c3-3f：atttcattcttcgttatagcaggc、qsw.c3-3ra：gaaggtgaccaagttcatgctgagcaaatctggtcatagaaagaga、qsw.c3-3rg：gaaggtcggagtcaacggattgagcaaatctggtcatagaaagagg。
24.以上所述的应用中，所述的试剂盒为parms检测试剂盒。
25.以上所述的应用中，当待测油菜的基因型为aa时，判定为千粒重较高的植株，当待测油菜的基因型为gg时，判定为千粒重较低的植株。
26.本发明的保护范围还包括：
27.检测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20819936位碱基的系统在油菜千粒重选择育种中的应用，所述的系统为试剂、试剂盒或仪器。
28.油菜千粒重选择育种的装置，其特征在于：所述的装置包括数据接收模块和数据处理模块，所述的数据接收模块被配置为接收待测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基的基因型；所述的数据处理模块用于将来自于接受模块的基因型转换为待测油菜的判定结果。
29.油菜千粒重选择育种的计算机可读存储介质，其特征在于：所述计算机可读存储介质使计算机运行如下步骤：检测待测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基的基因型，将所述基因型转换成判定待测油菜的判定结果。
30.以上所述的应用中，所述的试剂为检测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基的引物。
31.以上所述应用中，所述的引物为qsw.a3-3f：gaattaagtagaagctcgcaaacac、qsw.a3-3ra：gaaggtgaccaagttcatgcttgttgatggttgacttctggtga、qsw.a3-3rg：gaaggtcggagtcaacggatttgttgatggttgacttctggtgg。qsw.c3-3f：atttcattcttcgttatagcaggc、qsw.c3-3ra：gaaggtgaccaagttcatgctgagcaaatctggtcatagaaagaga、qsw.c3-3rg：gaaggtcggagtcaacggattgagcaaatctggtcatagaaagagg。
32.以上所述的应用中，所述的试剂盒为parms检测试剂盒。
33.以上所述的应用中，当待测油菜a03染色体的第13,373,806位碱基的基因型为aa且c03染色体的第20,819,936位碱基的基因型为aa，判定为千粒重较高的植株；当待测油菜a03染色体的第13,373,806位碱基的基因型为gg且c03染色体的第20,819,936位碱基的基因型为gg，判定为千粒重较低的植株。
34.与现有技术相比，本发明具有以下优点：
35.(1)本发明获得了与油菜千粒重显著关联的位点qsw.a3-3，且在两个环境中均可重复检测，平均可解释6.3％的表型方差，平均效应为0.30克，可有效应用于油菜粒重的遗
传改良。
36.(2)本发明获得了与油菜千粒重显著关联的位点qsw.c3-3，且在两个环境中均可重复检测，平均可解释8.6％的表型方差，平均效应为0.22克，可有效应用于油菜粒重的遗传改良。
37.(3)本发明获得了与油菜千粒重显著关联的parms标记，其检测方法操作简便、成本低廉且通量高，可对千粒重显著关联的基因组单倍型区域进行选择，提高选择效率和准确性。
附图说明
38.图1为油菜关联群体331份材料在10个环境中的千粒重表型频率分布示意图。
具体实施方式
39.本发明所述技术方案，如未特别说明，均为本领域的常规技术；所述试剂或材料，如未特别说明，均来源于商业渠道。本发明中使用的甘蓝型油菜基因组的版本号是darmor-bzh(version 4.1)。
40.实施例1：油菜千粒重显著关联snp标记的获得
41.(1)采集油菜关联群体331份核心株系的叶片，用ctab法提取总dna，利用油菜60k snp芯片对每个样本进行基因型分析(li et al.,2020)。
42.(2)利用illumina beadstudio基因分型软件(http://www.illumina.com/)计算每个位点的标记杂合率(heterozygous rate)、缺失率(missing rate)、最小等位基因频率(minor allele frequency)。以缺失率≤0.2、杂合率≤0.2、最小等位基因频率》0.05以及snp标记在甘蓝型油菜基因组中唯一匹配为筛选标准进行snp标记的过滤，最终获得24,508个高质量snp标记用于全基因组关联分析。利用structure2.3软件对331份核心株系进行群体结构分析(pritchard et al.,2000)，利用spagedi软件计算331份核心株系间亲缘关系(hardy and vekemans,2002)。
43.(3)将关联群体的331个株系分别于2012年(w12)、2013年(z13、w13和n13)、2014年(z14、w14和n14)、2015年(w15)、2016年(w16)和2017年(w17)在河南省农业科学院原阳基地(z)、中国农业科学院油料作物研究所阳逻基地(w)和江西省红壤研究所进贤基地(n)种植，合计10个环境。成熟期挑选每个小区中长势均匀的10个单株进行收获，按单株分别脱粒后将种子晒干进行千粒重的考察，计算各株系平均值。通过对关联群体各个环境中的千粒重数据进行分析，结果表明所有10个环境中的千粒重数据均呈正态分布而且表现很大的变异幅度(图1)，可以用于后续的全基因组关联分析。
44.(4)结合关联群体的基因型数据和千粒重表型数据，利用tassel 5.0软件(bradburyet al.2007)的四个模型(glm-pca，glm-q，mim-pca+k，mlm-q+k)进行全基因组关联分析。通过对不同环境和模型中检测到的显著性关联snp标记进行整合，获得与千粒重关联的位点qsw.a3-3和qsw.c3-3，分别可以在两个环境中被重复检测到，其峰值snp标记bn-a03-p14176469和bn-scaff_22067_1-p219321分别位于a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基，平均可解释6.3％和8.6％的表型方差，平均效应为0.30和0.22克。
45.表一、粒重关联位点qsw.a3-3和qsw.c3-3的信息
[0046][0047][0048]
实施例2：一种与千粒重显著关联parms标记的开发
[0049]
实施例1中所获得的关联标记为来源于snp芯片，只有用于分子杂交的探针序列信息。油菜snp芯片一次可以检测数以万计的位点，但其操作比较繁琐而且需要用到特殊设备。另外，若用油菜snp芯片对大量的育种中间材料进行检测则成本昂贵，因此需要将其转化为操作简单而且成本低廉的基于pcr扩增的检测方法，如parms(penta-primer amplifica tion refractory mutation system，五引物扩增受阻突变系统)标记。该标记系统包括一对荧光通用引物(fam和hex作为报告荧光)、一对snp等位基因特异引物及一条反向共用引物，可快速简单地进行snp等位基因型检测。
[0050]
(1)针对与qsw.a3-3关联的峰值snp标记bn-a03-p14176469，提取油菜a03染色体第13,373,806位碱基上下游各100bp序列。按照parms引物设计原则获得引物序列为：qsw.a3-3f：gaattaagtagaagctcgcaaacac、qsw.a3-3ra：gaaggtgaccaagt tcatgcttgttgatggttgacttctggtga、qsw.a3-3rg：gaaggtcggagtcaac ggatttgttgatggttgacttctggtgg。
[0051]
针对与qsw.c3-3关联的峰值snp标记bn-scaff_22067_1-p219321，提取油菜c03染色体第20,819,936位碱基上下游各100bp序列。按照parms引物设计原则获得引物序列为：qsw.c3-3f：atttcattcttcgttatagcaggc、qsw.c3-3ra：gaaggtgaccaagttcatgctgagcaaatctggtcatagaaagaga、qsw.c3-3rg：gaaggtcggagtcaacggattgagcaaatctggtcatagaaagagg。
[0052]
(2)以油菜关联群体的基因组dna为模板，利用上述引物进行荧光定量pcr扩增，采用tecanf200扫描fam、hex和rox信号并输出结果，最后转换成基因型。
[0053]
实施例3：parms标记在油菜千粒重筛选中的应用
[0054]
在331份关联群体(li et al.,2020)中，以实施例2提供的parms标记qsw.a3-3检测出基因型为aa的材料有79份，基因型为gg的为240份(表二)。两种基因型的千粒重差异在九个环境中达到显著性水平，aa和gg两种基因型的千粒重差值从0.18到0.34，平均0.23克。
[0055]
表二、parms标记qsw.a3-3两种基因型在关联群体十环境千粒重的比较
[0056]
基因型n13n14w12w13w14w15w16w17z13z14aa(n＝79)4.914.374.483.954.094.353.783.943.874.19gg(n＝240)4.724.094.143.773.894.103.563.723.684.09aa-gg0.19*0.28**0.34**0.18*0.20*0.25*0.22*0.22*0.18*0.10
[0057]
*和**分别代表显著性水平p＝0.05和0.01。
[0058]
在331份关联群体中，以实施例2提供parms标记qsw.c3-3检测出基因型为aa的材料有66份，基因型为gg的为262份(表三)。两种基因型的千粒重差异在七个环境中达到显著性水平，aa和gg两种基因型的千粒重差值从0.19到0.41，平均0.30克。
[0059]
表三、parms标记qsw.c3-3两种基因型在关联群体十环境千粒重的比较
[0060]
基因型n13n14w12w13w14w15w16w17z13z14aa(n＝66)4.774.354.573.994.204.473.893.943.814.14gg(n＝262)4.644.114.163.763.894.103.543.753.694.09aa-gg0.130.24*0.41***0.22*0.31**0.37**0.34**0.19*0.120.05
[0061]
*，**和***分别代表显著性水平p＝0.05，0.01和0.001。
[0062]
在331份关联群体中，以实施例2提供parms标记qsw.a3-3和qsw.c3-3的引物进行联合检测，qsw.a3-3中基因型为aa且qsw.c3-3中基因型为aa的材料有23份，qsw.a3-3中基因型为gg且qsw.c3-3中基因型为gg的材料有175份，两种基因型的千粒重差异在九个环境中都达到显著性水平。其中，5个环境中的差异达到极显著水平，差值从0.41到0.56，平均0.46克。
[0063]
表四、parms标记qsw.a3-3和qsw.c3-3两种基因型在关联群体十环境千粒重的比较
[0064][0065]
*，**和***分别代表显著性水平p＝0.05，0.01和0.001。
[0066]
以上结果足以说明我们制备的parms标记qsw.a3-3和qsw.c3-3与油菜的千粒重是高度关联的，无论是单独还是联合使用均具有很好的选择效果。

技术特征：
1.检测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的系统在油菜千粒重选择育种中的应用，所述的系统为试剂、试剂盒或仪器。2.油菜千粒重选择育种的装置，其特征在于：所述的装置包括数据接收模块和数据处理模块，所述的数据接收模块被配置为接收待测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的基因型；所述的数据处理模块用于将来自于接受模块的基因型转换为待测油菜的判定结果。3.油菜千粒重选择育种的计算机可读存储介质，其特征在于：所述计算机可读存储介质使计算机运行如下步骤：检测待测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的基因型，将所述基因型转换成判定待测油菜的判定结果。4.根据权利要求1所述的应用，所述的试剂为检测油菜darmorv4.1参考基因组c03染色体第20,819,936位碱基的引物。5.根据权利要求4所述的应用，所述的引物为：qsw.c3-3f：atttcattcttcgttatagcaggc、qsw.c3-3ra：gaaggtgaccaagttcatgctgagcaaatctggtcatagaaagaga和qsw.c3-3rg：gaaggtcggagtcaacggattgagcaaatctggtcatagaaagagg。6.检测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基的系统在油菜千粒重选择育种中的应用，所述的系统为试剂、试剂盒或仪器。7.油菜千粒重选择育种的装置，其特征在于：所述的装置包括数据接收模块和数据处理模块，所述的数据接收模块被配置为接收待测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基的基因型；所述的数据处理模块用于将来自于接受模块的基因型转换为待测油菜的判定结果。8.油菜千粒重选择育种的计算机可读存储介质，其特征在于：所述计算机可读存储介质使计算机运行如下步骤：检测待测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基的基因型，将所述基因型转换成判定待测油菜的判定结果。9.根据权利要求6所述的应用，所述的试剂为检测油菜darmorv4.1参考基因组a03染色体第13,373,806位碱基和c03染色体第20,819,936位碱基的引物。10.根据权利要求9所述的应用，所述的引物为：qsw.a3-3f：gaattaagtagaagctcgcaaacac、qsw.a3-3ra：gaaggtgaccaagttcatgcttgttgatggttgacttctggtga和qsw.a3-3rg：gaaggtcggagtcaacggatttgttgatggttgacttctggtgg；qsw.c3-3f：atttcattcttcgttatagcaggc、qsw.c3-3ra：gaaggtgaccaagttcatgctgagcaaatctggtcatagaaagaga和qsw.c3-3rg：gaaggtcggagtcaacggattgagcaaatctggtcatagaaagagg。

技术总结
本发明属于生物育种技术领域，具体涉及一油菜千粒重关联的PARMS分子标记或标记组合的应用，申请人筛选到两个千粒重关联位点qSW.A3-3和qSW.C3-3的峰值SNP标记Bn-A03-p14176469和Bn-scaff_22067_1-p219321分别位于油菜A03染色体的第13,373,806和C03染色体的第20,819,936位碱基处，平均可解释6.3%和8.6%的表型方差，平均效应为0.30和0.22克。利用这两个SNP开发的PARMS标记对油菜关联群体进行检测，发现其操作简便且分型清晰，而且两种纯合基因型千粒重的差值分别为0.23和0.30克，两者联合则高达0.46克，有良好的选择效果。有良好的选择效果。


技术研发人员：师家勤 王汉中 王新发 刘贵华 梁华兵 詹杰鹏
受保护的技术使用者：权利要求书1页说明书6页序列表（电子公布）附图1页
技术研发日：2023.05.07
技术公布日：2023/8/9