与大豆蛋白质含量相关的qtl及分子标记的获得方法、分子标记和应用
【专利摘要】本发明提供两个与大豆蛋白质含量相关的QTL、分子标记、分子标记的获得方法及其应用。QTL分别位于连锁群A2、N上，分别由分子标记Satt409、Satt584定位，分子标记Satt409的5’?3’引物如SEQ ID NO.3所示，3’?5’引物如SEQ ID NO.4所示；分子标记Satt584的5’?3’引物如SEQ ID NO.5所示，3’?5’引物如SEQ ID NO.6所示。分别以相应大豆个体基因组DNA为模板，用上述5’?3’引物和3’?5’引物对上述模板进行PCR扩增，将PCR扩增得到的产物经聚丙烯酰氨凝胶电泳分离后，得到目标分子标记。得到的分子标记可用于大豆分子育种。
【专利说明】
与大豆蛋白质含量相关的QTL及分子标巧的获得方法、分子标 巧和应用
技术领域
[0001] 本发明属于分子生物学领域，具体设及与大豆蛋白质含量相关的WL及分子标记 的获得方法、分子标记和应用。
【背景技术】
[0002] 大豆是在中国乃至世界范围当中都是重要的经济作物，其原因就在于大豆是人们 日常生活中需要的植物蛋白的主要来源之一，同时大豆又是畜牧养殖业的重要的饲料作 物。因此，大豆子粒当中的蛋白质含量成为衡量大豆品种的重要指标之一。然而，目前我国 主要生产的高蛋白品种大豆无法满足人们日常生产生活需要。因此，利用现代化分子育种 技术提高现有大豆品种蛋白质含量成为当代育种科研人员的主要任务。
[0003] 大豆蛋白含量属于多基因控制的数量性状遗传，对运一性状进行WL分析和定位， 可W发掘有利用价值的SSR位点，W培育优质高产的新品种。国内外对于大豆蛋白WL做过 大量分析，但是没有对高效位点进行有效确定，也没有明确位点的真实作用，从而无法直接 应用到实际育种当中，并且可能出现重要位点的丢失。
[0004] 近年来分子标记辅助育种已经得到了一定的发展，但是却有一些不足之处。首先， 大豆蛋白质是由多基因控制的数量性状遗传，运一遗传方式是由多基因控制的，因为单一 位点对性状的影响很小，必须由多个基因对其操控才行，再加上环境因子对于表型的影响 也很大，所W单一位点的作用很难被直接评估出来，导致无法应用在实际育种当中。对于数 量性状遗传的最初步认识，是将微效多基因看成一个整体用生物统计法来分析。后来又出 现了微效多基因假说，认为数量遗传是有多个微效基因共同控制。最新的研究成果已经证 实运一假设是有缺陷的。数量性状遗传是有一个主效基因和多个微效基因共同调控，而且 基因与环境互作对表型的影响很大。可W说数量性状的表现就是有主效位点、微效位点和 环境因子共同影响下完成的。因此从基因的层面上来讲，判断出各个位点对蛋白质含量的 效应值就显得极为重要。
[0005] 另外分层评价是W卡方检验和T检验为手段的一种统计方法，用于分析一组影响 因子对某一效应产生的效应值。是一种W综合评判标准分析某一性状的统计方法。目前为 止分层评价在生物领域当中的应用并不多见。目前，国内外众多大豆育种学者都在研究大 豆的数量遗传性状。目前已发表的有关MAS和OTL定位的文章多为研究的标记和位点的潜在 应用价值，很少将其真正应用到实际育种中，更少有人对数量遗传基因的效应值给出结论， 在大豆子粒蛋白质含量运一方面也是一个空白。众所周知，大豆子粒蛋白质含量是数量性 状遗传，由多基因控制。在运些基因中必然有的基因对该性状的影响较大，有些则影响较 小。在得知各个基因对蛋白质含量运一性状的具体影响之前，很难将运些直接应用于实际 育种当中。采取分层评价是一个比较有效的解决手段。
【发明内容】
[0006] 针对上述问题，本发明通过分层评价等统计分析方法，提供5种与大豆中大豆蛋白 质含量密切相关的QTL，并准确的验证出各个大豆蛋白质含量相关WL的效应值W及作用方 式，通过变化率的计算发现3个关键WL能够增加蛋白质含量，2个关键WL能够降低大豆蛋 白质含量。从而将其直接应用于实际分子育种当中。
[0007] 本发明提供的五个9化中，Satt683属于大豆子粒蛋白质含量的主效QTLsSattAOg、 8曰1:1:180、5日1:1:584、5日1:1:181属于高效0化。
[000引其中包含的等位基因 Satt683-3、Satt409-2和Satt584-5对于大豆蛋白含量运个 性状表现的是正相关，即其对应的等位基因在大豆子粒中的作用是促进大豆蛋白质含量增 加，而Sattl81-2和Sattl80-5对于大豆蛋白含量运个性状表现为负相关，即其对应的等位 基因在大豆子粒中的作用是抑制大蛋白质含量增加。
[0009] 3个对大豆蛋白质性状相关的WL的信息见表1，其中列出的等位基因为与性状正 相关的等位基因。
[0010] 表1 3个对大豆蛋白质性状相关的QTL
[0011]
[0012] 2个对大豆蛋白质性状相关的WL的信息件表2,其中列出的等位基因为与性状负 相关的等位基因。
[0013]表2 2个对大丑蛋白质性状相关的QTL 「mi4l
[0015] 可W直接用引物扩增，鉴定待测品种是否含有该等位基因，用于分子标记辅助选 择。
[0016] ( -）、本发明提供一个与大豆蛋白质含量相关的QTL:位于连锁群N上，在公共图谱 的基因组位置为27.787cM，由分子标记Satt683定位;分子标记Satt683的5 ' -3 '引物如沈Q ID NO. 1所示，3'-5'引物如沈Q ID NO.2所示。
[0017] 该OTL包含效应值为12%的等位基因Satt683-3,该等位基因与大豆蛋白质含量正 相关。
[0018] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记的获得方法:该分子 标记的5'-3'引物如SEQ ID NO. 1所示，3'-5'引物如SEQ ID NO.2所示；W大豆品种北豆5号 为轮回亲本，与天碟蛋Z孤02114杂交，共经过两次杂交和两次回交之后获得的大豆个体基 因组DNA为模板;用上述5 '-3 '引物和3 '-5 '引物对上述模板进行PCR扩增，其PCR扩增体系为 10化，包括 1. SngAiL的DNA模板0.祉L，0.5皿01/L的 5 ' -3 ' 引物、3 ' -5 ' 各0.扣L，200mmol/L的 dNTPO. 6化，5U/iiL的化9酶0.05化，IOXBuff erl化，阳离子采用MgCl2终浓度1.5-2. Ommol/L， 其余用超纯水补足;PCR扩增程序为预变性94 °C 5分钟;变性94 °C 40秒、退火47 °C 40秒、延伸 72 °C 40秒，循环35次;延伸72 °C 5分钟；4 °C结束反应。将上述PCR扩增得到的产物经PCR产物 变形程序:94°C5分钟，电泳分离后，得到的大小为45化p、43化P或41化P的扩增片段，即为目 标分子标记。
[0019] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记:该分子标记的5'-3' 引物如沈Q ID NO. 1所示，3'-5'引物如沈Q ID NO.2所示。
[0020] 本发明还提供上述的分子标记用于大豆分子育种的方法:该分子标记的5'-3'引 物如沈Q ID NO. 1所示，3'-5'引物如沈Q ID NO. 2所示；用其5'-3'引物和3'-5'引对待测大 豆品种的基因组DNA进行PCR扩增，其PCR扩增体系为10化，包括1. SngAiL的DNA模板0.祉L， 0.5WI101/L的 5 ' -3 ' 引物、3 ' -5 ' 各0 .化L，200mmo 1 /L的dNTPO.化L，511/化的Taq酶0.05化， IOXBuffer化L，阳离子采用MgCl2终浓度1.5-2.Ommol/L，其余用超纯水补足;PCR扩增程序 为预变性94 °C 5分钟;变性94 °C 40秒、退火47 °C 40秒、延伸72 °C 40秒，循环35次;延伸72 °C 5分 钟;4°C结束反应。将上述扩增产物经PCR产物变形程序:94°C5分钟，电泳分离后，如果扩增 产物中存在大小为45化p、43化P或41化P的扩增片段，则该待测大豆品种具有与大豆蛋白质 含量相关的QTLSatt683,若扩增产物中存在大小为41化P的扩增片段，则该待测大豆品种为 具有与大豆蛋白质含量相关的QTLSatt683中的增加大豆蛋白质含量的等位基因Satt683-3 的品种，否则，待测品种不具有该OTL或该等位基因。
[0021] (二）、本发明提供的另一个与大豆蛋白质含量相关的QTL:位于连锁群A2上，在公 共图谱的基因组位置为145.565。1，由分子标记5曰《409定位，分子标记5曰《409的5'-3'弓| 物如沈Q ID NO.3所示，3'-5'引物如沈Q ID NO.4所示。
[0022] 该OTL包含效应值为10%的等位基因Satt409-2,该等位基因与大豆蛋白质含量正 相关。
[0023] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记的获得方法:该分子 标记的5'-3'引物如SEQ ID NO.3所示，3'-5'引物如SEQ ID NO.4所示；W大豆品种北豆5号 为轮回亲本，与天碟蛋Z孤02114或中黄杂交，共经过两次杂交和两次回交之后获得的大豆 个体基因组DNA为模板；用上述5 '-3 '引物和3 '-5 '引物对上述模板进行PCR扩增，其PCR扩增 体系和程序同（一）；将上述PCR扩增得到的产物经电泳分离后，得到的大小为360bp、330bp 或30化P的扩增片段，即为目标分子标记。
[0024] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记:该分子标记的5'-3' 引物如沈Q ID NO.3所示，3'-5'引物如沈Q ID NO.4所示。
[0025] 本发明还提供上述的分子标记用于大豆分子育种的方法:该分子标记的5'-3'引 物如沈Q ID NO. 3所示，3'-5'引物如沈Q ID NO.4所示；用其5'-3'引物和3'-5'引对待测大 豆品种的基因组DNA进行PCR扩增，其PCR扩增体系和程序同（一）；将上述扩增产物经电泳分 离后，如果扩增产物中存在大小为360bp、330bp或300bp的扩增片段，则该待测大豆品种具 有与大豆蛋白质含量相关的QTLSatt409,若扩增产物中存在大小为33化P的扩增片段，则该 待测大豆品种为具有与大豆蛋白质含量相关的QTLSatt409中的增加大豆蛋白质含量的等 位基因Satt409-2的品种，否则，待测品种不具有该OTL或该等位基因。
[0026] (S)、本发明提供的另一个与大豆蛋白质含量相关的Q化:位于连锁群N上，在公共 图谱的基因组位置为29.37cM，由分子标记Satt584定位，分子标记Satt584的5 ' -3 '引物如 沈Q ID NO.5所示，3'-5'引物如沈Q ID NO.6所示。
[0027] 该OTL包含效应值为10%的等位基因Satt584-5,该等位基因与大豆蛋白质含量正 相关。
[0028] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记的获得方法:该分子 标记的5'-3'引物如SEQ ID NO.5所示，3'-5'引物如SEQ ID NO.6所示；W大豆品种北豆5号 为轮回亲本，与天碟蛋Z孤02114杂交，共经过两次杂交和两次回交之后获得的大豆个体基 因组DNA为模板;用上述5 ' -3 '引物和3 ' -5 '引物对上述模板进行PCR扩增，其PCR扩增体系和 程序同（一）；将上述PCR扩增得到的产物经电泳分离后，得到的大小为40化P、35化P或30化P 的扩增片段，即为目标分子标记。
[0029] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记:该分子标记的5'-3' 引物如沈Q ID NO.5所示，3'-5'引物如沈Q ID NO.6所示。
[0030] 本发明还提供上述的分子标记用于大豆分子育种的方法:该分子标记的5'-3'引 物如沈Q ID NO. 5所示，3'-5'引物如沈Q ID NO.6所示；用其5'-3'引物和3'-5'引对待测大 豆品种的基因组DNA进行PCR扩增，其PCR扩增体系和程序同（一）；将上述扩增产物经电泳分 离后，如果扩增产物中存在大小为400bp、350bp或300bp的扩增片段，则该待测大豆品种具 有与大豆蛋白质含量相关的QTLSatt584，若扩增产物中存在大小为30化P的扩增片段，则该 待测大豆品种为具有与大豆蛋白质含量相关的QTLSatt584中的增加大豆蛋白质含量的等 位基因Satt584-5的品种，否则，待测品种不具有该OTL或该等位基因。
[0031] (四）、本发明提供的另一个与大豆蛋白质含量相关的QTL:位于连锁群Cl上，在公 共图谱的基因组位置为105.092CM，由分子标记Satt 181定位，分子标记Satt 181的5 ' -3 '引 物如沈Q ID NO.7所示，3'-5'引物如沈Q ID NO.8所示。
[0032] 该WL包含效应值为8%的等位基因Sattl81-2,该等位基因与大豆蛋白质含量负 相关。
[0033] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记的获得方法:该分子 标记的5'-3'引物如SEQ ID NO.7所示，3'-5'引物如SEQ ID NO.8所示；W大豆品种北豆5号 为轮回亲本，与大明白豆子或艾卡166杂交，共经过两次杂交和两次回交之后获得的大豆个 体基因组DNA为模板;用上述5 '-3 '引物和3 '-5 '引物对上述模板进行PCR扩增，其PCR扩增体 系和程序同（一）；将上述PCR扩增得到的产物经电泳分离后，得到的大小为3(K)bp、280bp或 ISObp的扩增片段，即为目标分子标记。
[0034] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记:该分子标记的5'-3' 引物如沈Q ID NO.7所示，3'-5'引物如沈Q ID NO.8所示。
[0035] 本发明还提供上述的分子标记用于大豆分子育种的方法:该分子标记的5'-3'引 物如沈Q ID NO. 7所示，3'-5'引物如沈Q ID NO.8所示；用其5'-3'引物和3'-5'引对待测大 豆品种的基因组DNA进行PCR扩增，其PCR扩增体系和程序同（一）；将上述扩增产物经电泳分 离后，如果扩增产物中存在大小为300bp、280bp或ISObp的扩增片段，则该待测大豆品种具 有与大豆蛋白质含量相关的QTLSattlSl,若扩增产物中存在大小为28化P的扩增片段，则该 待测大豆品种为具有与大豆蛋白质含量相关的QTLSattlSl中的减少大豆蛋白质含量的等 位基因Sattl81-2的品种，否则，待测品种不具有该OTL或该等位基因。
[0036] (五）、本发明提供的另一个与大豆蛋白质含量相关的Q化:位于连锁群H上，在公共 图谱的基因组位置为87.597cm，由分子标记SattlSO定位，分子标记SattlSO的5'-3'引物如 沈Q ID NO.9所示，3'-5'引物如沈Q ID NO. 10所示。
[0037] 该WL包含效应值为8%的等位基因Sattl80-5,该等位基因与大豆蛋白质含量负 相关。
[0038] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记的获得方法:该分子 标记的5'-3'引物如SEQ ID NO.9所示，3'-5'引物如SEQ ID NO. 10所示；W大豆品种北豆5 号为轮回亲本，与Biogedulote或艾卡166或瑞丰2号杂交，共经过两次杂交和两次回交之后 获得的大豆个体基因组DNA为模板；用上述5'-3'引物和3'-5'引物对上述模板进行PCR扩 增，其PCR扩增体系和程序同（一）；将上述PCR扩增得到的产物经电泳分离后，得到的大小为 350bp、340bp或3(K)bp的扩增片段，即为目标分子标记。
[0039] 本发明还提供上述与大豆蛋白质含量相关的WL的分子标记:该分子标记的5'-3' 引物如沈Q ID NO.9所示，3'-5'引物如沈Q ID NO. 10所示。
[0040] 本发明还提供上述的分子标记用于大豆分子育种的方法:该分子标记的5'-3'引 物如SEQ ID NO. 9所示，3'-5'引物如SEQ ID NO. 10所示；用其5'-3'引物和3'-5'引对待测 大豆品种的基因组DNA进行PCR扩增，其PCR扩增体系和程序同（一）；将上述扩增产物经电泳 分离后，如果扩增产物中存在大小为350bp、34化P或30化P的扩增片段，则该待测大豆品种 具有与大豆蛋白质含量相关的QTLSattlSO,若扩增产物中存在大小为3(K)bp的扩增片段，贝U 该待测大豆品种为具有与大豆蛋白质含量相关的QTLSattlSO中的减少大豆蛋白质含量的 等位基因Sattl80-5的品种，否则，待测品种不具有该OTL或该等位基因。
[0041] 上述电泳分离，为聚丙締酷胺凝胶电泳或毛细管电泳。了解了蛋白相关基因的效 应值，我们就可W通过检验某大豆品种当中是否有高效位点来判断该品种的大豆是否是蛋 白含量高的优质大豆，运样就实现该基因在分子辅助育种当中的实际应用。
[0042] 本发明在已知的大量大豆蛋白质相关WL中确定重要的、效应大的QTL，并确定各 个基因对大豆蛋白质形状的具体影响：是正相关还是负相关，并利用运些WL实现W分子标 记为手段进行育种。
【附图说明】
[0043] 图1大豆Satt683在表1所示的250份资源大豆中的编号1-61中的带型；
[0044] 图2卡方分析资源等位基因的分层分析；
[0045] 图3 T-测验的资源等位基因的分层分析；
[0046] 图4等位基因在群体的SSR分析；
[0047] 图5卡方分析群体等位基因的分层分析；
[0048] 图6 T-测验的群体等位基因的分层分析；
[0049] 图7各个等位基因在资源和北豆群体中变化率；
[0化0] 图8 Satt683的引物对母本为北豆5号父本为天碟蛋ZDD02114的群体的基因组DNA 扩展产物毛细管电泳图；
[0化1]图9 Satt409的引物对母本为北豆5号父本为天碟蛋ZDD02114或中黄的群体的基 因组DNA扩展产物毛细管电泳图；
[0化2] 图10 Sat巧84的引物对母本为北豆5号父本为天碟蛋ZDD02114的群体的基因组 DNA扩展产物毛细管电泳图；
[0053]图11 SattlSl的引物对母本为北豆5号父本为大明白豆子或艾卡166的群体的基 因组DNA扩展产物毛细管电泳图；
[0化4] 图12 SattlSO的引物对母本为北豆5号父本为Biogedulote或艾卡166或瑞丰2号 的群体的基因组DNA扩展产物毛细管电泳图。【具体实施方式】
[0055] W下实施例中：大豆蛋白质含量测定全部使用FOSS Infratec 1241谷物分析仪； 毛细管电泳分析全部使用全自动毛细管电泳仪，Advanced Anal}ftical Technologies(美 国AATI公司），型号Rragment Analyzer?。
[0056] 实施例1实验材料的种植和蛋白质含量检测：
[0057] (1)实验材料的种植和采样
[0058] 资源材料:选取的是250份中国地区的主栽大豆品种。表1是250份资源的品种名称 及子粒蛋白质含量。
[0059] 种植场所:吉林省农科院。
[0060] 种植方法:Im行长，每行21粒种子。
[0061 ] 管理方法：同大田管理。
[0062] 采样方法：营养生长阶段，取植株顶端最幼嫩处叶片用于提取DNA、PCR反应。收获 时按单株收取大豆子粒用于测定蛋白质含量。
[0063] 验证材料：是W黑龙江省主栽大豆品种北豆5号作为轮回亲本，与多马卡?托利 萨、坡黄、化rosoy、绿75、中特1号、NOVAXenUi巧、东农163、杜纳吉卡、Amsoy、中豆27、95- 5383、艾卡166、大明白豆子、东山69、瑞丰2号、中黄、Biogedulote和天碟蛋ZDD02114运19个 品种的大豆分别杂交，总共经过两次杂交和两次回交之后收获BC2F2。群体总共330个个体。
[0064] 种植场所:黑龙江省农星科研育种中屯、。
[0(?日]种植方法:5m行长，每行100粒种子。
[0066] 管理方法：同大田管理。
[0067] 采样方法：营养生长阶段，取植株顶端最幼嫩处叶片用于提取DNA、PCR反应。收获 时按单株收取大豆子粒用于测定蛋白质含量。
[0068] (2)蛋白质含量的测定:分别测定上述资源材料和验证材料中每种大豆的子粒蛋 白质含量。
[0069] 使用步骤
[0070] 1)打开In打atec 1241谷物分析仪的电源。
[0071] 2)当仪器进行完自检预热好之后，从应用模型列表栏中通过上下键找到所需的模 型。应用模型名称中有附加的"STM"的。
[0072] 3)将大豆子粒倒入样品槽，之后支架装入对接单元，注意有齿轮条的一侧在右边。 向下轻轻按支架，W确保齿轮条和齿轮正确晒合。将样品槽支架放入对接单元，运时候要小 屯、，W免损坏内部的小齿轮。
[0073] 4)在键盘区按下分析键开始分析。
[0074] FOSS Infratec 1241谷物分析仪测出的数据是".crcT格式。需要先转换成EXCEL 格式才能打开。导出数据步骤：
[00巧]IFT文件夹中包含Diskl ,Disk2两个目录和一个说明文件，运行Diskl目录下的 setup.exe。
[0076] 开始一 FOSS^IFT;
[0077] Result Data/CRD一TXT-fOrmat(Columns)。
[0078] 在Source folder中选择20071030.crd的路径，并在文件前选中V。再选择输出路 径，如Desktop，点击Convert selected。
[0079] 为文件起名，保存。退出IFT。得到一个TXT文件。
[0080] 新建一个Excel,打开TXT文件，在文本导入向导中选择"分隔符号"点击下一步，选 择叮ab键"和"逗号"，点击下一步，点击完成。就可W将FOSS Infratec 1241谷物分析仪测 得的大豆蛋白质含量的数据用Excel格式打开。
[0081] 表1 250份大豆资源名称及其子粒蛋白质含量
[0083]
[0084] 实施例2 DNA提取和PCR扩增：
[0085] (1)大豆叶片DNA的提取（250份资源大豆与验证材料330个个体均提取，每种大豆 分别提取）
[0086] 1)将实施例1中取得的用于提取DNA、PCR反应的新鲜叶片置于1.5ml离屯、管中，将 离屯、管放在浸于液氮中的离屯、管架上，用塑料钻头低溫钻磨，直至样品变为白色细粉末为 止。
[0087] 2)打开水浴锅预热至65°C，同时将CTAB提取液(按照40ml: 1化1加入琉基乙醇）加 热至65°C。溫度达到之后将每个离屯、管中加入0.7mlCTAB提取液，充分混匀之后，水浴40分 钟，并且间歇轻摇。
[0088] 3)将水浴之后的离屯、管中加入等体积氯仿，然后8000(转/分)离屯、20分钟。取上清 液，再次加入等体积氯仿。离屯、8000(转/分)20分钟。
[0089] 4)取上清液加至预先预冷的0.7ml异丙醇中，直至DNA析出。
[0090] 5)用无水乙醇冲洗DNA，之后惊干，WSOiil超纯水溶解。
[0091] 6) 1 %琼脂糖凝胶电泳检测。电压1OOV电流1OOA运行20分钟。
[0092] 注意事项：1)移液器黄色枪头的尖端部分减去一小部分，防止剪切力破坏DNA结 构）
[0093] 2)移液W及冲洗DNA的时候一定要缓慢轻柔，W免损坏DNA结构
[0094] 表2用于提取大豆叶片DNA的CTAB提取液配方
[00951
[0
[0097] (2)PCR 反应
[0098] 选取分布在大豆DNA13个连锁群上的29对引物（见表3)，利用选取的29对SSR引物 对250份资源大豆的基因组DNA进行PCR扩增，并通过电泳分析:在表4所示体系中，按W下程 序分别进行扩增，引物由上海生工合成。
[0099] PCR程序:预变性94 °C 5分钟
[0100] 变性 94 °C 40 秒
[0101] 退火 47 °C 40 秒
[0102] 延伸 72°C40 秒
[0103] 循环35次
[0104] 延伸72 °C 5分钟
[0105] 4°C结束反应。
[0106] 表3在大豆DNA13个连锁群上选取的29对引物位置和序列
[0109] 表4 PCR反应（1化1体系）
[0110]
[Cn11」 选取原则：1)巧重巧13个连锁稱巧杀上郡选取SSK引物
[0112] 2)-条连锁群上选取的几个引物遗传距离尽量远
[0113] (3)聚丙締酷胺凝胶电泳:上述步骤(2)得到的扩增产物分别进行聚丙締酷胺凝胶 电泳检测，得到表1所示的250个资源品种大豆和验证资源330个个体的DNA的29对引物的 PCR扩增产物带型，部分带型见图1。
[0114] 表5聚丙締酷胺凝胶电泳药品
[0115]
[0116]
[0117] 聚丙締酷胺凝胶电泳步骤
[0118] 将大板用酒精擦干净，再均匀涂抹亲和硅烷(酒精:亲和硅烷:冰乙酸=5ml: 50山： 5化1)
[0119] 将耳板放入抽风楓内，用酒精擦干净，再均匀涂抹剥离硅烷。
[0120] 将干净的边条放在大板两侧，再盖上耳板。两块板涂抹硅烷的一面要朝里。用夹子 固定。
[0121] 将配置好的6%PA胶内加入20化IAPS和50iilTEMED。均匀灌入大板和耳板之间，注 意不要出现气泡。
[0122] 将齿梳平滑一面朝里插入两板之间5mm左右。调水平。静置一小时左右即可使用
[0123] 上板，夹紧，将电泳槽内倒入缓冲液TBE(上槽为1倍TBE，下槽为1/3倍TBE)。
[0124] 点样槽里滴入少量染色剂，1800V电压预电泳10分钟左右。期间对PCR产物内加入 染色剂并作变性处理。
[0125] 将齿梳插入点样槽（齿朝下），点样，1800V电压运行90分钟。下板。取下耳板和边 条。
[01%] 银染：1500ml蒸馈水中加入ISOOiil硝酸银溶液。将大板浸泡在银染溶液里20~30 分钟。
[0127] 蒸馈水冲洗半分钟。
[0128] 10)显影：1500ml蒸馈水，加入1500山甲醒、3 %化OH和6.25g化2C03。将大板浸泡在 银染溶液里20~30分钟。
[0129] 11)统计带型，记录数据，扫描。
[0130] 表6分层中的部分引物扩增出的片段大小
[0131]
[0132] 实施例3大豆蛋白质含量相关等位基因分层评价：
[0133] (1)大豆蛋白质含量相关等位基因分层评价
[0134] 将实施例1得到的表1所示的250个资源品种中每个品种的蛋白质含量的数据与实 施例2(3)中得到的其带型数据一一对应，选取蛋白质含量前40%和后40%利用MATLAB R2009a卡方检验分析等位基因的效应值。Wal地a = 0.Ol和0.1作为分层的阔值，将结果分 为3个层次。其中当检验P值小于0.01时，我们认为该位点为主效位点，而P值大于0.1时，认 为是微效位点，参考图2，P值小于0.01水平(极显著）的等位基因有Satt683-3、Satt578-4、 Satt409-2 和 Satt079-1。显著水平在 0.1 和 0.01 之间的有 Sattl81-2、Sattl80-5、Sattl26- 1、Satt578-2、Satt584-5、Satt409-10、Sat_311-2、Sat_242-2、Satt578-6、Satt584-6、 8过的713-1、53的012-5、53古_216-2、53古_216-3、53的584-2、53的700-2和53^369-2，其余的 等位基因其显著水平均大于0.1。
[0135] 经过Minitabl4.12T检验之后参照图3,显著水平小于0.01的共有十个等位基因， 与之前的卡方检验到的4个极显著等位基因对比之下可W发现共有=个等位基因在卡方检 验和T检验检测下，均表现为极显著，分别是Satt683-3、Satt578-4和Satt409-2。除了运S 个极显著水平的等位基因，其余的等位基因要么表现为一般显著或是不显著，要么就是在 两种检验方法中的显著水平不一致。可W初步区分出主效基因和微效基因。
[0136] (2)关键等位基因在北豆5号群体中的验证
[0137] 根据资源分层数据，选出极显著水平的全部引物W及其他显著水平较高的引物。 进行卡方检验和T检验。同时为了验证分层效果的成立，再选择一个显著性水平极低的引 物。同样进行两种分析。
[0138] 选取极显著水平的等位基因Satt683-3、Satt578-4和Satt409-2,选取卡方检测一 般显著而T检测极显著的等位基因Sattl80-5、Sattl81-2和Satt584-5,选取卡方检测和T检 测都不显著的等位基因Sat_174-2共7个等位基因的引物(表3)，W北豆5号为轮回亲本的群 体的330个个体的基因组DNA为模板，在表4所示体系中，按W下程序分别进行PCR扩增，并通 过电泳分析。引物由上海生工合成。部分电泳图见图4。统计带型。
[0139] PCR程序:预变性94°C 5分钟；
[0140] 变性 94 °C 40 秒，
[0141] 退火 47 °C 40 秒，
[0142] 延伸 72°C40 秒，
[0143] 循环35次；
[0144] 延伸72 °C 5分钟，
[0145] 4。(：结束反应。
[0146] 与本实施例（1)对资源的带型数据处理方法相同，对7对引物在北豆5号群体中的 带型数据进行卡方检验和T检验：
[0147] 将北豆5号导入系群体个体(验证资源330个个体）的蛋白质含量的数据与带型数 据一一对应
[0148] W蛋白质含量为关键字全排序。选出前40%和后40%
[0149] 利用MATLAB R2009a进行卡方检验
[0150] 打开Min;Uabl4.12进行T检验。验证卡方检验结果。
[0151] 结果如图5和图6。
[0152] 从W上两张图可W看出，Satt683-3在两种检测中依然都是极显著，说明该等位基 因确实属于大豆子粒蛋白质含量的主效基因。Satt409-2在T检验中仍是极显著，尽管卡方 检验结果有所下降，但是仍然保持在较高的水平。Sattl80-5和Satt584-5却表现出在两种 检验中都为极显著，再加上同上一轮分层评价表现一致的Sattl81-2,可W判断出运四种等 位基因应该属于高效基因，只是其作用效果不如Satt683-3明显。而Sat_174-2则与上一轮 的结果完全一致，都为不显著，说明该基因确实在增加蛋白质含量运方面效应值很低。意外 的是在资源中显著水平极高的引物Satt578的各等位基因在北豆群体中显著水平极低，甚 至未发现等位基因Satt578-4。运个极个别现象产生原因可能是北豆5号导入系群体的20个 亲本的基因组当中没有该基因，另外测量仪器不够准确等原因产生的系统误差也可能造成 运种极个别现象发生。
[0153] 因此，虽然两个群体中的分层评价略有出入，但差异并不大，总体上的层次还是比 较明显，那就是在资源中显著水平较高的，在北豆5号群体中依然较高，不同引物的等位基 因在两个群体中的作用是一致的。运样就可W认为，对于运29个引物的等位基因分层评价 其结论是准确的。
[0154] 实施例4利用表型变化率验证分析等位基因加性效应：
[0155] 本发明通过变化率来表示是否含有该等位基因的变化，运种变化率是基于资源或 群体中个体成熟期后表型值差异，若某个等位基因对大豆子粒蛋白质含量有正方向的加性 效应，那么含有该等位基因的资源或群体个体的表型数据要比不含有该等位基因的资源或 个体的要高，但是不排除标记之间的互作和等位基因间的互作的影响。上述变化率的计算 方法为：
[0156]
[0157] 其中Rate Of change是变化率，化Iue A代表含有A等位基因的个体或资源的蛋白 质含量表型值的平均值，Value B代表含有B等位基因或不含有A等位基因的个体或资源的 蛋白质含量表型值的平均值。
[0158] 利用表型变化率进一步验证关键等位基因加性效应
[0159] 从理论上讲，存在某等位基因的大豆品种其蛋白质含量如果高于不存在该等位基 因的大豆品种，那么则说明该等位基因对于蛋白质含量运一数量性状表现为正向调控的正 相关基因。反之，该等位基因则是负相关的基因。本发明选择了经过两次分层评价筛选得到 的5条极显著等位基因，选出含有运5条等位基因的北豆5号小群体(表7)。通过重要引物等 位基因在资源和北豆5号群体中对蛋白质含量产生的变化率来进行验证。结果如图7。筛选 之后得到几个引物及等位基因相应的验证群体。群体的母本是北斗5号，父本见表7。
[0162]从表中可W明确看出运5个等位基因对资源和群体中大豆子粒蛋白质含量的影响 效果是一致的。Satt683-3、Satt409-2和Satt584-5运S条等位基因在两个群体中的变化率
[01601 車7絡4去妈至Il掛"C田革全島#竺:仿五:巧 [0161 都是正的，说明他们对于大豆蛋白含量运个性状表现的是正相关，即该等位基因在大豆子 粒中的作用是促进大豆蛋白质含量增加。而Sattl81-2和Sattl80-5在两个群体中的变化率 都是负的，证明运两个等位基因对于大豆蛋白含量运个性状表现为负相关，即该等位基因 在大豆子粒中的作用是抑制大蛋白质含量增加。
[0163] 实施例5等位基因位点在分子标记育种中的应用
[0164] 9化位点Satt683的引物（SEQ ID N0.1、SEQ ID 腸.2)、9化位点Satt409的引物 (沈9 10側.3、沈9 10側.4)、9化位点8曰的584的引物（569 10側.5、沈9 10側.6)、9化位 点SattlSl的引物（SEQ ID N0.7、SEQ ID 側.8)、9化位点SattlSO的引物（SEQ ID N0.9、SEQ ID NO. 10)，W北豆5号为轮回亲本的群体的330个个体的基因组DNA为模板，在表4所示体系 中，按W下程序分别进行PCR扩增，并通过毛细管电泳分析(步骤和参数参照仪器说明）。引 物由上海生工合成。
[01化]PCR程序:预变性94°C 5分钟；
[0166] 变性 94°C40 秒，
[0167] 退火 47 °C 40 秒，
[016引延伸72 °C 40秒，
[0169] 循环35次；
[0170] 延伸72 °C 5分钟；
[0171] 4°C结束反应。
[0172] 毛细管电泳分析结果表明，Satt683的引物对母本为北豆5号父本为天碟蛋 Z孤02114的群体的基因组DNA扩展出了大小为41化P的分子标记片段（图8)，说明该群体基 因组DNA中含有提高大豆蛋白含量的Satt683-3等位基因；
[0173] Satt409的引物对母本为北豆5号父本为天碟蛋ZDD02114或中黄的群体的基因组 DNA扩展出了大小为33化P的分子标记片段（图9)，说明该群体基因组DNA中含有提高大豆蛋 白含量的Satt409-2等位基因；
[0174] Satt584的引物对母本为北豆5号父本为天碟蛋Z孤02114的群体的基因组DNA扩展 出了大小为3(K)bp的分子标记片段(图10)，说明该群体基因组DNA中含有提高大豆蛋白含量 的Satt584-5等位基因；
[0175] SattlSl的引物对母本为北豆5号父本为大明白豆子或艾卡166的群体的基因组 DNA扩展出了大小为28化P的分子标记片段（图11)，说明该群体基因组DNA中含有降低大豆 蛋白含量的Sattl81-2等位基因；
[0176] SattlSO的引物对母本为北豆5号父本为Biogedulote或艾卡166或瑞丰2号的群体 的基因组DNA扩展出了大小为3(K)bp的分子标记片段（图12)，说明该群体基因组DNA中含有 降低大豆蛋白含量的Sattl80-5等位基因。
【主权项】
1. 一个与大豆蛋白质含量相关的QTL，其特征在于:位于连锁群A2上，在公共图谱的基 因组位置为145.565cM，由分子标记Satt409定位，分子标记Satt409的5 ' -3 '引物如SEQ ID NO.3所示，3'-5'引物如SEQ ID NO.4所示。2. 根据权利要求1所述的QTL，其特征在于：该QTL包含效应值为10 %的等位基因 Satt409-2,该等位基因与大豆蛋白质含量正相关。3. -种权利要求1所述与大豆蛋白质含量相关的QTL的分子标记的获得方法，其特征在 于:该分子标记的5'_3'引物如SEQ ID NO.3所示，3'-5'引物如SEQ ID NO.4所示；以大豆品 种北豆5号为轮回亲本，与天鹅蛋ZDD02114或中黄杂交，共经过两次杂交和两次回交之后获 得的大豆个体基因组DNA为模板；用上述5 '-3 '引物和3 '-5 '引物对上述模板进行PCR扩增， 将PCR扩增得到的产物经电泳分离后，得到的大小360bp、330bp或300bp的扩增片段，即为目 标分子标记。4. 一种权利要求1所述与大豆蛋白质含量相关的QTL的分子标记，其特征在于：该分子 标记的5'-3'引物如SEQ ID NO.3所示，3'-5'引物如SEQ ID NO.4所示。5. -种权利要求4所述的分子标记用于大豆分子育种的方法，该分子标记的5'-3'引物 如SEQ ID NO. 3所示，3'-5'引物如SEQ ID NO.4所示；用其5'-3'引物和3'-5'引对待测大豆 品种的基因组DNA进行PCR扩增，将扩增产物经电泳分离后，如果扩增产物中存在大小为 360bp、330bp或300bp的扩增片段，则该待测大豆品种具有与大豆蛋白质含量相关的 QTLSatt409,若扩增产物中存在大小为330bp的扩增片段，则该待测大豆品种为具有与大豆 蛋白质含量相关的QTLSatt409中的增加大豆蛋白质含量的等位基因 QTLSatt409-2的品种， 否则，待测品种不具有该QTL或该等位基因。6. -个与大豆蛋白质含量相关的QTL，其特征在于:位于连锁群N上，在公共图谱的基因 组位置为29.37cM，由分子标记Satt584定位，分子标记Satt584的5'-3'引物如SEQ ID N0.5 所示，3'-5'引物如SEQ ID NO.6所示。7. 根据权利要求6所述的QTL，其特征在于：该QTL包含效应值为10 %的等位基因 Satt584-5,该等位基因与大豆蛋白质含量正相关。8. -种权利要求6所述与大豆蛋白质含量相关的QTL的分子标记的获得方法，其特征在 于:该分子标记的5'_3'引物如SEQ ID NO.5所示，3'-5'引物如SEQ ID NO.6所示；以大豆品 种北豆5号为轮回亲本，与天鹅蛋ZDD02114杂交，共经过两次杂交和两次回交之后获得的大 豆个体基因组DNA为模板;用上述5 ' -3 '引物和3 ' -5 '引物对上述模板进行PCR扩增，将PCR扩 增得到的产物经电泳分离后，得到的大小400bp、350bp或300bp的扩增片段，即为目标分子 记 D9. 一种权利要求6所述与大豆蛋白质含量相关的QTL的分子标记，其特征在于：该分子 标记的5'-3'引物如SEQ ID NO.5所示，3'-5'引物如SEQ ID NO.6所示。10. -种权利要求9所述的分子标记用于大豆分子育种的方法，该分子标记的5 ' -3 '弓丨 物如SEQ ID NO. 5所示，3'-5'引物如SEQ ID NO.6所示；用其5'-3'引物和3'-5'引对待测大 豆品种的基因组DNA进行PCR扩增，将扩增产物经电泳分离后，如果扩增产物中存在大小为 400bp、350bp或300bp的扩增片段，则该待测大豆品种具有与大豆蛋白质含量相关的 QTLSatt584,若扩增产物中存在大小为300bp的扩增片段，则该待测大豆品种为具有与大豆 蛋白质含量相关的QTLSatt584中的增加大豆蛋白质含量的等位基因 Satt584-5的品种，否 贝1J，待测品种不具有该QTL或该等位基因 D
【文档编号】C12Q1/68GK105925722SQ201610541719
【公开日】2016年9月7日
【申请日】2016年7月11日
【发明人】齐照明, 陈庆山, 辛大伟, 蒋洪蔚, 武小霞, 胡振帮, 杜翔宇, 齐慧冬, 李伟
【申请人】东北农业大学