陈强^*, 闫龙^*, 邓莹莹, 肖二宁, 刘兵强, 杨春燕^*, 张孟臣^*
河北省农林科学院粮油作物研究所 / 国家大豆改良中心石家庄分中心 / 农业部黄淮海大豆生物学与遗传育种重点实验室 / 河北省遗传育种重点实验室, 河北石家庄 050035
^* 通讯作者(Corresponding authors): 张孟臣, E-mail: mengchenzhang@hotmail.com; 杨春燕, E-mail: chyyang66@163.com 第一作者联系方式: 陈强, E-mail: chenqiangwsm@163.com; 闫龙, E-mail: dragonyan1979@163.com
收稿日期:2016-01-26 接受日期:2016-05-09网络出版日期:2016-06-06基金:本研究由国家自然科学基金项目(31471522), 河北省自然科学基金项目(C2015301012), 国家高技术研究发展计划(863计划)项目(2012AA101106), 国家现代农业产业技术体系建设专项(CARS-004-PS06)和国家转基因生物新品种培育重大专项(2014ZX0800402B)资助

摘要大豆籽粒大小和粒形性状不仅与产量和外观品质紧密相关, 还对机械化播种有着一定的影响。本研究采用大粒栽培品种冀豆12与小粒半野生地方品种黑豆(ZDD03651)杂交衍生的包含188个重组自交系的F_6:8和F_6:9群体为材料, 对粒长、粒宽、粒厚、长宽比、长厚比和宽厚比的遗传结构进行分析, 并分别以WinQTLCart 2.5、QTLNetwork 2.1和IciMapping 4.1 3种模型对以上性状的加性效应QTL, QE互作效应及上位性互作效应进行检测。6个性状的广义遗传率介于64.01%~79.57%, 遗传力较高, 且除粒厚外的其他性状受环境影响显著。共定位到加性效应QTL 38个, 单个QTL的贡献率介于2.21%~10.71%之间, 分布在12条染色体的17个标记区间内, 且12个染色体区段至少与2种性状相关。两种及以上模型同时检测到的QTL有24个, 3种模型均能检测到的QTL共8个, 分别为 qSL-17-1、qSL-18-1、qSW-6-1、qST-2-1、qST-6-1、qSLT-2-2、qSWT-2-1和 qSWT-20-1。检测到7对上位性互作QTL, 分别涉及粒长、粒宽、长宽比、长厚比和宽厚比, 互作效应贡献率介于0.78%~6.20%之间。QE互作效应贡献率均较低, 介于0.0005%~0.3900%之间。以多种模型同时检测结果准确性较高, 可为分子标记辅助育种工作提供可靠理论基础。

关键词:大豆; 粒形; QTL; 上位性互作
Mapping Quantitative Trait Loci for Seed Size and Shape Traits in Soybean
CHEN Qiang^*, YAN Long^*, DENG Ying-Ying, XIAO Er-Ning, LIU Bing-Qiang, YANG Chun-Yan^*, ZHANG Meng-Chen^*
Institute of Cereal and Oil Crops, Hebei Academy of Agriculture and Forestry Sciences / Shijiazhuang Branch of National Soybean Improvement Center / Huang-Huai-Hai Key Laboratory of Biology and Genetic Improvement of Soybean, Ministry of Agriculture, Shijiazhuang 050035, China
Fund:This study was supported by the National Natural Science Foundation of China (31471522), the Natural Science Foundation of Heibei Province (C2015301012), the National High-tech Research and Development Program of China (2012AA101106), the China Agriculture Research System (CARS-004-PS06), and the National Major Project for Developing New GM Crops (2014ZX0800402B)
AbstractSeed size and shape not only relate to seed yield and quality, but also affect mechanical seeding in soybean ( Glycine maxL.). In this study, F_6:8 and F_6:9populations derived from Jidou 12 × Heidou were used to analyze the genetic character and detect quantitative trait loci (QTLs) for seed length, seed width, seed thickness, seed length-to-width ratio, seed length-to-thickness ratio, and seed width-to-thickness ratio. Softwares WinQTLCart 2.5, QTLNetwork 2.1 and IciMapping 4.1 were used to identify the additive, epistatic and environmentally interacted QTLs for seed size and shape related traits. As results, the heritability of the six traits varied from 64.01% to 79.57%. A total of 38 additive QTLs were identified to be located on 12 chromosomes, with the heritability varying from 2.21% to 10.71%. Eight of them ( qSL-17-1, qSL-18-1, qSW-6-1, qST-2-1, qST-6-1, qSLT-2-2, qSWT-2-1, and qSWT-20-1) were identified using three methods, simultaneously. In the meantime, seven pairs of additive × additive epistasis were detected and the heritability of epistasis pairs ranged from 0.78% to 6.20%. Additionally, the effects of QTL by environment interaction ranged from 0.0005% to 0.3900%. The QTLs identified using different mapping softwares in this study could provide a reliable theoretical basis for marker-assisted selection breeding.

Keyword:Soybean; Seed shape; Quantitative trait locus; Epistasis
Show Figures
Show Figures






在大豆中, 粒长、粒宽、粒厚等籽粒大小性状是衡量百粒重的重要指标, 与产量密切相关^[1]。而大豆粒形性状直接关系大豆的品质和商品性, 既是机械化播种的重要限制因素, 又在一定程度上决定大豆在国际贸易中的商业价值^[2]。籽粒大小和粒形性状为多基因控制的数量性状^[3], 籽粒大小性状遗传力一般为59%~79%, 粒形性状的遗传力一般在19%~56%间^[4]。通过QTL定位分析, 检测籽粒大小和粒形相关QTL位点, 对分子标记辅助选择的高产、高品质大豆育种具有重要的指导意义。
在水稻中籽粒大小和粒形的研究较为深入^{[5, 6, 7, 8, 9]}, 定位了较多QTL, 并克隆了控制粒长和粒重的基因GS3^[5], 控制粒宽和粒重的基因GW2^[6], 控制粒宽、长宽比和粒重的基因qSW5^[7]等多个基因。大豆籽粒大小和粒形相关性状QTL定位的研究相对较少^{[10, 11, 12, 13]}。大豆数据网站SoyBase (http://www.soybase.org/)中公布了粒长QTL 18个、粒宽QTL 16个、粒厚QTL14个。近期, 牛远等^{[14, 15]}在初定位基础上, 采用Lasso法和复合区间作图法, 同时将O连锁群中Satt331- Satt592区段内多效性QTL分解后精细定位, 并对目标区段内候选基因进行预测。Xie等^[16]对位于C2连锁群Satt640-Satt422区间内的粒形QTL精细定位分析, 分别检测到粒长、粒宽、粒厚及百粒重相关联的SNP标记, 并预测了结果区段内可能的效应基因。Hu等^[17]在分离群体中检测到12个籽粒大小和粒形相关QTL, 在关联分析中检测到41个SNP标记与6个性状相关联。而随着研究的深入, QTL间的上位性互作及QE互作逐渐受到研究人员的重视^{[18, 19, 20]}。Niu等^[21]以257份大豆栽培品种为材料, 采用关联分析方法检测到59个主效QTL和31个QTL与环境互作效应。梁慧珍等^[22]利用QTLNetwork 2.0软件检测到3对粒形性状上位性互作QTL。闫宁等^[23]采用Bayes分层广义线性模型对粒形性状的主效QTL、QTL与环境互作效应、QQ上位性互作及QTL与细胞质互作效应分别进行检测。苏成付等^[24]和李杰勤等^[25]研究结果显示不同的定位程序适用的遗传模型范围不同, 建议使用多种作图法共同作图, 并优先标定共同发现的QTL。
目前, 大豆粒形相关性状上位性互作及QE互作研究相对较少, 且前人研究多采用单一模型检测QTL。本研究利用重组自交系群体为材料, 分别采用WinQTLCart 2.5^[26]、QTLNetwork 2.1^[27]和IciMapping 4.1^{[28, 29]}软件对籽粒大小和粒形各性状的加性QTL, QTL间上位性互作及QE互作进行检测, 以得到可靠的QTL信息, 为进一步精细定位及辅助育种提供依据。
1 材料与方法1.1 材料及田间试验以冀豆12× 黑豆杂交构建的包含188个家系的F_6:8和F_6:9重组自交系为材料。冀豆12为河北省农林科学院粮油作物研究所育成的百粒重为22 g的大粒品种, 黑豆(ZDD03651)为陕西省半野生地方品种, 百粒重为7 g。两亲本粒长、粒宽等籽粒大小相关性状间存在较大差异。将该群体F_6:8(2011年)和F_6:9(2013年)重组自交系播种于石家庄田间, 采用3 m行长, 3行区, 随机区组试验设计, 行距50 cm, 株距10 cm。材料自然成熟后从每小区收取中间行中间10株。
1.2 性状调查自每家系籽粒中随机选取10粒, 采用数显游标卡尺测量粒长(seed length, SL)、粒宽(seed width, SW)、粒厚(seed thickness, ST), 长宽比(seed length- to-width, SLW) = 粒长/粒宽、长厚比(seed length-to- thickness, SLT) = 粒长/粒厚、宽厚比(seed width-to- thickness, SWT) = 粒宽/粒厚^[21], 重复2次, 以平均数为观测值。
1.3 统计分析采用SPSS 17.0软件中相应程序对各性状进行方差、相关性和偏相关等分析。
广义遗传力h²= σ ²_g/(σ ²_g+σ ²_ge/e+σ ²_e/re) × 100%, 其中σ ²_g表示基因型方差, σ ²_ge表示基因型与环境互作的方差, σ ²_e为误差项, e为环境数, re为重复次数^{[30, 31]}。
1.4 QTL定位在陈强^[32]和雷雅坤等^[33]分子遗传连锁图谱的基础上, 采用WinQTLCart 2.5^[34]中的复合区间作图法(CIM)检测每个性状主效QTL位点, 选择模型6, 向前回归方法, 选择被检区间每侧10 cM的控制窗口, 进行1000次排列获得QTL显著的LOD阈值(P< 0.05)。利用QTLNetwork 2.1软件中基于混合线性模型的复合区间作图法(MCIM)分析QTL 之间的互作效应^{[35, 36]}, 将Permutation次数设置为1000次, 用于计算F值。以1 cM的步长对全基因组进行QTL扫描, 选择结果中P值显著的互作确定上位性互作及环境互作QTL。以IciMapping 4.1软件利用MET模块中ICIM-ADD和ICIM-EPI进行多环境联合分析, 以ICIM法进行上位性QTL定位, 常用LOD值为2.0, 由于本研究包含2个环境, 再结合定位结果, 故选取LOD值为4.0, 超过阈值就认为存在QTL^[37]。同一性状的LOD峰值相距在5 cM以内者视为同一QTL。采取McCouch等^[38]的方式命名QTL。

2 结果与分析2.1 表型变异及相关分析两年的表型数据分析结果(表1)显示, 冀豆12的粒长、粒宽和粒厚高于黑豆, 而黑豆的长宽比、长厚比和宽厚比高于冀豆12。群体中粒形相关6个性状均存在较大程度分离, 变异系数介于5.30%~ 8.71%之间, 且各性状均存在超双亲分离现象。各性状遗传力介于64.01%~79.57%之间。基因型与环境间的双因素方差分析结果显示, 粒长、粒宽和宽厚比受基因型、环境和基因型× 环境的影响均极显著(P< 0.01); 粒宽受基因型的影响极显著(P< 0.01); 长宽比受基因型和环境的影响均极显著(P< 0.01); 长厚比受基因型和环境的影响极显著(P< 0.01), 而基因与环境互作的影响显著(P< 0.05)。各性状的表型值均连续变异, 且为正态分布(图1)。
表1
Table 1
表1(Table 1)

表1 亲本及群体籽粒大小与粒形性状在2011和2013年中的表现 Table 1 Characteristics of seed size and shape traits in RILs derived from the cross between Jidou 12 and Heidou in 2011 and 2013 性状 Trait 年份 Year 亲本 Parent 重组自交系群体 RIL 双因素方差分析P值 P-value of double factor variance analysis 冀豆12 Jidou12 黑豆 Heidou 平均值 Mean 最大值 Max. 最小值 Min. 标准差 SD 变异系数 CV(%) 基因型 Genotype (G) 环境 Environment (E) 基因型× 环境 G × E 粒长SL (mm) 2011 8.68± 0.15 7.33± 0.19 7.92 9.63 6.56 0.55 6.98 0.000* * 0.000* * 0.000* * 2013 7.60± 0.12 7.05± 0.22 7.49 8.62 6.17 0.47 6.27 粒宽SW (mm) 2011 7.08± 0.18 4.83± 0.09 5.80 7.10 4.57 0.39 6.67 0.000* * 0.000* * 0.000* * 2013 6.70± 0.16 4.77± 0.09 5.75 6.70 4.65 0.35 6.08 粒厚ST (mm) 2011 5.52± 0.10 3.19± 0.12 4.35 5.58 3.17 0.36 8.17 0.000* * 0.467 0.974 2013 5.56± 0.11 3.27± 0.17 4.38 5.58 3.26 0.37 8.44 长宽比SLW 2011 1.23± 0.01 1.52± 0.02 1.37 1.69 1.13 0.08 6.09 0.000* * 0.000* * 0.189 2013 1.13± 0.01 1.48± 0.01 1.30 1.54 1.13 0.07 5.38 长厚比SLT 2011 1.57± 0.02 2.30± 0.02 1.83 2.34 1.44 0.16 8.71 0.000* * 0.000* * 0.040* 2013 1.37± 0.02 2.16± 0.02 1.71 2.15 1.36 0.14 8.18 宽厚比SWT 2011 1.28± 0.01 1.52± 0.02 1.34 1.65 1.11 0.08 6.26 0.000* * 0.000* * 0.000* * 2013 1.21± 0.01 1.46± 0.03 1.32 1.61 1.07 0.07 5.30 * Significant at the 0.05 probability level. * * Significant at the 0.01 probability level. SL: seed length; SW: seed width; ST: seed thickness; SLW: seed length-to-width ratio; SLT: seed length-to-thickness ratio; SWT: seed width-to-thickness ratio. * 表示0.05的显著水平; * * 表示0.01的显著水平。

表1 亲本及群体籽粒大小与粒形性状在2011和2013年中的表现 Table 1 Characteristics of seed size and shape traits in RILs derived from the cross between Jidou 12 and Heidou in 2011 and 2013

图1
Fig. 1

Figure OptionViewDownloadNew Window

图1 群体粒形相关性状频率分布图 A: 群体粒长表型分布; B: 群体粒宽表型分布; C: 群体粒厚表型分布; D: 群体长宽比表型分布; E: 群体长厚比表型分布; F: 群体宽厚比表型分布Fig. 1 Frequency distribution of seed size traits of soybean populations in 2011 and 2013 A: frequency distribution of seed length; B: frequency distribution of seed width; C: frequency distribution of seed thickness; D: frequency distribution of seed length-to-width ratio; E: frequency distribution of seed length-to-thickness ratio; F: frequency distribution of seed width-to-thickness ratio.

偏相关分析结果显示(表2), 粒长与粒宽、长宽比间均显著正相关, 与粒厚极显著负相关, 与长厚比显著负相关, 与宽厚比在2011年相关性不显著, 2013年极显著负相关。粒宽与粒厚、宽厚比极显著正相关, 与长宽比极显著负相关; 与长厚比在2011年相关性不显著, 2013年极显著负相关。粒厚与长宽比极显著正相关, 与宽厚比极显著负相关, 与长宽比相关性不显著。长宽比与长厚比极显著正相关, 与宽厚比极显著负相关。长宽比与长厚比间为极显著正相关。粒长与宽厚比2011年简单相关分析为极显著正相关, 偏相关分析无显著相关性, 2013年简单相关无显著相关性, 偏相关分析为极显著负相关。性状间相关性不稳定, 表明该相关受其他性状的影响较大。
表2
Table 2
表2(Table 2)

表2 两年中粒形性状间的简单相关和偏相关分析 Table 2 Simple and partial correlation coefficients for seed traits in soybean in two years 年份 Year 性状 Trait 粒长 SL 粒宽 SW 粒厚 ST 长宽比 SLW 长厚比 SLT 宽厚比 SWT 2011 粒长SL 0.876* * -0.388* * 0.697* * -0.176* 0.030 粒宽SW 0.574* * 0.782* * -0.441* * -0.080 0.393* * 粒厚ST 0.357* * 0.746* * -0.062 0.388* * -0.722* * 长宽比SLW 0.443* * -0.468* * -0.415* * 0.827* * -0.650* * 长厚比SLT 0.397* * -0.296* * -0.698* * 0.747* * 0.917* * 宽厚比SWT 0.209* * 0.219* * -0.629* * 0.086 0.706* * 2013 粒长SL 0.863* * -0.415* * 0.623* * 0.138 -0.292* * 粒宽SW 0.577* * 0.815* * -0.366* * -0.341* * 0.631* * 粒厚ST 0.384* * 0.707* * -0.063 0.468* * -0.814* * 长宽比SLW 0.471* * -0.380* * -0.303* * 0.687* * -0.488* * 长厚比SLT 0.344* * -0.284* * -0.685* * 0.715* * 0.893* * 宽厚比SWT -0.001 -0.043 -0.687* * 0.067 0.710* * In the panel of 2011 or 2013, the data at lower left corner are coefficients of correlation between traits, and those at upper right corner are coefficients of partial correlation between traits. * Significant at the 0.05 probability level. * * Significant at the 0.01 probability level. SL: seed length; SW: seed width; ST: seed thickness; SLW: seed length-to-width ratio; SLT: seed length-to-thickness ratio; SWT: seed width-to-thickness ratio. 表中左下角为各性状间单相关系数, 右上角为各性状间偏相关系数; * 表示0.05的显著水平; * * 表示0.01的显著水平。

表2 两年中粒形性状间的简单相关和偏相关分析 Table 2 Simple and partial correlation coefficients for seed traits in soybean in two years

2.2 主效QTL定位及QE互作效应分析WinQTLCart 2.5和QTLNetwork 2.1进行1000次排列获得QTL显著的LOD阈值均为2.5。两年中6性状共检测到38个主效应QTL, 分布在12条染色体中(图2), 单个QTL的贡献率介于2.21%~10.71%之间。24个QTL能在2种及以上的方法中检测到, 其中8个QTL可被3种方法同时检测到(表3)。
图2
Fig. 2

	Figure OptionViewDownloadNew Window
	图2 检测到的主效QTL在连锁群中的分布Fig. 2 Distribution of main-effect QTLs for seed shape traits on linkage groups

表3
Table 3
表3(Table 3)

表3 不同方法同时检测到的大豆粒形性状加OTL及QE互作效应 Table 3 QTL for seed shape tratis in soybean

表3 不同方法同时检测到的大豆粒形性状加OTL及QE互作效应 Table 3 QTL for seed shape tratis in soybean

检测到4个控制粒长的加性QTL, 贡献率范围为4.04%~8.76%。其中qSL-17-1和qSL-18-1可被3种方法同时检测到, 贡献率最大分别为4.67%和8.72%, 加性效应最大分别为-0.10 mm和0.15 mm, 增效基因分别来自于冀豆12和黑豆, 与环境互作的贡献率最大分别为0.10%和0.29%。检测到3个控制粒宽的加性QTL, 其中qSW-6-1在3种方法中均能被检测到, 贡献率最大为9.59%, 加性效应最大值为-0.12 mm, 与环境互作的贡献率最大值为0.02%; qSW-17-1和qSW-2-1分别被两种方法检测到, 贡献率最大分别为6.37%和7.44%, 加性效应值为-0.12 mm和-0.15 mm, 环境互作贡献率分别为0.03%和0.04%。检测到7个控制粒厚的加性QTL, 其中qST-6-1和qST-2-1在3种方法中可同时被检测到, 贡献率最大分别为6.80%和5.72%, 加性效应值最大分别为-0.10 mm和-0.09 mm, 环境互作贡献率最大值分别为0.20%和0.09%; qST-8-1、qST-12-1、qST-13-1和qST-17-1同时被2种方法检测到, 贡献率最大分别为5.61%、8.01%、6.73%和6.25%, 加性效应最大分别为-0.07、-0.09、-0.10和-0.10 mm, 与环境互作贡献率范围为0.01%~0.27%。
检测到3个控制长宽比的加性QTL, 其中qSLW-2-1在两种方法中同时被检测到, 贡献率最大为8.67%, 加性效应最大分别为0.07, 与环境互作贡献率为0.39%。共检测到9个控制长厚比的加性QTL, 其中qSLT-2-2在3种方法中同时被检测到, 贡献率最大为8.69%, 加性效应最大为0.07, 与环境互作贡献率为0.15%。qSLT-8-1和qSLT-13-1分别在两种方法中检测到, 贡献率最大分别为5.07%和8.23%。共检测到8个控制宽厚比的加性QTL, qSWT-2-1和qSWT-20-1可同时在3种方法中被检测到, 贡献率最大分别为7.87%和7.17%, 加性效应均为0.02, 与环境互作贡献率分别为0.29%和0.23%; qSWT-6-1、qSWT- 8-1、qSWT-13-1和qSWT-17-1可分别被两种方法检测到, 贡献率分别为7.05%、10.71%、8.51%和6.46%。
2.3 粒形性状上位性互作效应结果共检测到7对上位性互作QTL, 分别涉及粒长、粒宽、长宽比、长厚比和宽厚比(表4)。粒长上位性互作QTL发生在非主效QTL间, 贡献率为6.20%; 粒宽和长宽比上位性互作均发生在主效QTL与非主效QTL间, 贡献率分别为4.58%和4.06%, 效应值分别为-0.10 mm和-0.02, 表明该互作为亲本型大于重组型; 长厚比中检测到2对发生在主效QTL间的上位性互作, 为qSLT-13-1分别与qSLT-6-1和qSLT-2-2发生的互作, 贡献率分别为1.41%和0.94%, 效应值均为0.02。宽厚比上位性QTL为qSWT-20-1分别与qSWT-8-1和qSWT-13-3互作, 效应值均为0.01, 但效应相反。
表4
Table 4
表4(Table 4)

表4 粒形性状上位互作效应QTL Table 4 Epistatie effects QTL for seed shape traits

表4 粒形性状上位互作效应QTL Table 4 Epistatie effects QTL for seed shape traits

3 讨论3.1 与前人研究结果的比较在QTL定位研究中, 前人建议使用多种作图法, 并优先标定共同发现的QTL ^{[24, 25]}。本研究采用的3种作图方法, 共检测到籽粒大小和粒形相关性状的加性效应QTL 38个, 其中3种方法同时检测到的QTL有8个, 分别是粒长QTLqSL-17-1和 qSL-18-1, 粒宽QTLqSW-6-1, 粒厚QTLqST-2-1和qST-6-1, 长厚比QTLqSLT-2-2, 宽厚比QTLqSWT-2-1和qSWT-20-1。参照Song等^[39]的大豆公共图谱标记信息与前人研究结果比较发现, 以上8个位点中有5个与前人报道基本一致。本研究检测到的位于第18染色体的粒长QTL, 标记区间为Satt570-Sat_290, 与Moongkanna等^[10]检测到粒长QTL区间Satt303-Sat_ 315部分重叠; 位于第17染色体的粒长QTL标记区间为Satt301-Satt310, 与Yu等^[18]检测到粒长QTL区间Satt266-Sat_354部分重叠。本研究检测到的位于第6染色体的粒宽、粒厚QTL标记区间为Sat_062-Satt28, 与Yu等^[18]检测到粒宽QTL区间Satt640-Satt422重叠, 且Xie等^[16]在该区段内对粒长、粒宽和粒厚精细定位分析, 并分别检测到关联标记, Salas等^[11]在该区段内同样检测到粒宽和粒厚QTL。本研究检测到的位于第20染色体的宽厚比QTL标记区间为Satt496-Satt049, 与Yu等^[18]检测到粒宽QTL区间Satt419-Satt354部分重叠。该结果表明多种方法同时检测到的位点有着较高的准确性。而对于粒厚QTLqST-2-1、长厚比QTLqSLW-2-2、宽厚比QTLqSWT-2-1则未见报道。
本研究用2种方法检测到的QTL有14个, 仅用1种方法检测到的QTL有14个。其中在第4染色体中的粒长QTLqSL-4-1标记区间为Satt607-Satt361, Niu等^[21]检测到的粒长相关联的标记AW277661位于该区段内; 在第8染色体中粒厚QTLqST-8-2标记区间Satt424-Satt341内同样含有Niu等^[21]检测到的粒厚相关联的标记AW132402; 在第5染色体长厚比和宽厚比的QTLqSLT-5-1和qSWT-5-1标记区间Sat_267-Satt200, 与Salas等^[11]检测到的粒形QTL标记区间A975_1~K636_2部分重叠。
本研究中检测到的位于在第2染色体Satt266- satt611区间、第2染色体Sat_069-Sat_183区间、第12染色体Satt422-Satt541区间、第13染色体Satt144-Satt656区间和第17染色体Satt669-Sat_292区间, 共检测到13个QTL, 均未见前人报道。
3.2 染色体区段多效性分析本研究检测到的38个QTL集中分布在17个标记区间内, 其中13个标记区段存在同时与多个性状相关的现象。前人研究表明, 作物相关性状QTL的集中分布是一种普遍现象, 有不少QTL可能具有多效性。Mansur等^[40]认为与性状间的高度遗传力和遗传相关一致, 可能是一个上位基因作用的结果, 也可能是控制这些性状的基因紧密排列所致。在水稻中, 杜景红等^[41]应用剩余杂合体衍生群体, 在初定位基础上, 把水稻第6染色体短臂上控制水稻谷壳硅含量的qHUS-6分解为qHUS6-1和qHUS6-2, 并进一步将qHUS6-2分解为一个效应较大的qHUS6-2a和效应较小的qHUS6-2b。在大豆中, 牛远等^{[14, 15]}通过构建剩余杂合系和重组自交系, 将O染色体Satt331-Satt592 区间剖分为与标记S21-S22和O23~O19/O19~O21关联的2个多效性QTL。Xie等^[16]对第6染色体Satt640-Satt422间的粒形性状效应分析, 分别检测到粒长、粒宽和粒厚相关联的SNP标记。梁慧珍等^[22]对粒形性状进行QTL定位分析, 同样发现同一区段同时与多个性状相关的现象, 通过对不同性状间增效基因的来源对区段的多效性进行了初步推测。
本研究第2染色体Sat_069-Sat_183区段和第6染色体Sat_062-Satt281区段均同时与粒宽、粒厚、长宽比和长厚比相关, 第4染色体Satt607-Satt361区段同时与粒长、长宽比和长厚比相关, 第13染色体Satt144-Satt656区段同时与粒厚、长厚比和宽厚比相关。以上4个区段内的QTL, 对于籽粒大小性状(粒长、粒宽或粒厚)增效基因来自冀豆12, 而对于粒形性状(长宽比、长厚比或宽厚比)则相反, 推测其可能为多个QTL的聚合。位于第17染色体Satt301-Satt310区段同时与粒长和粒宽相关, 且增效基因均来自冀豆12, 推测其可能为多效性QTL。具体是一个主效基因还是由多个QTL 紧密连锁的聚合体还需进一步研究。

4 结论共检测到籽粒大小和粒形相关性状中加性效应QTL 38个, 以3种方法同时检测到的位点8个, 分别为qSL-17-1、qSL-18-1、qSW-6-1、qST-2-1、qST-6-1、qSLT-2-2、qSWT-2-1和qSWT-20-1, 可靠性较高, 其中qST-2-1、qSLT-2-2和qSWT-2-1为首次报道。环境互作效应贡献率均较低, 而上位性互作位点7对, 效应相对较高, 其中对粒形性状影响显著, 不能忽略。
The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.


参考文献View Option
原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

[1]	Liang H Z, Li W D, Wang H, Fang X J. Genetic effects on seed traits in soybean. Acta Genet Sin, 2005, 32: 1199-1204[本文引用:1]
[2]	Wilson D O Jr. Storage of orthodox seeds. In: Basra A S ed. Seed Quality: Basic Mechanisms Agricultural Implications. New York: Food Products Press, 1995. pp 173-208[本文引用:1]
[3]	Nelson R L, Wang P. Variation and evaluation of seed shape in soybean. Crop Sci, 1989, 29: 147-150[本文引用:1]
[4]	Cober E R, Voldeng H D, Fregeau-Reid J A. Heritability of seed shape and seed size in soybean. Crop Sci, 1997, 37: 1767-1769[本文引用:1]
[5]	Fan C C, Xing Y Z, Mao H L, Lu T T, Han B, Xu C G, Li X H, Zhang Q F. GS3, a major QTL for grain length and weight and minor QTL for grain width and thickness in rice, encodes a putative transmembrane protein. Theor Appl Genet, 2006, 112: 1164-1171[本文引用:2]
[6]	Song X J, Huang W, Shi M, Zhu M Z, Lin H X. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat Genet, 2007, 39: 623-630[本文引用:2]
[7]	Wan X Y, Weng J F, Zhai H Q, Wang J K, Lei C L, Liu X L, Guo T, Jiang L, Su N, Wan J M. QTL analysis for rice grain width and fine mapping of an identified QTL allele gw-5 in a recombination hotspot region on chromosome 5. Genetics, 2008, 179: 2239-2252[本文引用:2]
[8]	Zhang X J, Wang J F, Huang J, Lan H X, Wang C L, Yin C F, Wu Y Y, Tang H J, Qian Q, Li J Y, Zhang H S. Rare allele of OsPPKL1 associated with grain length causes extra-large grain and a significant yield increase in rice. Proc Natl Acad Sci USA, 2012, 109: 21534-21539[本文引用:1]
[9]	Sun L J, Li X J, Fu Y C, Zhu Z F, Tan L B, Liu F X, Sun X Y, Sun X W, Sun C Q. GS6, a member of the GRAS gene family, negatively regulates grain size in rice. J Integr Plant Biol, 2013, 55: 938-949[本文引用:1]
[10]	Moongkanna J, Nakasathien S, Novitzky W P, Kwanyuen P, Sinchaisri P, Srinives P. SSR markers linking to seed traits and total oil content in soybean. Thai J Agric Sci, 2011, 44: 233-241[本文引用:2]
[11]	Salas P, Oyarzo-Llaipen J C, Wang D, Chase K, Mansur L. Genetic mapping of seed shape in three populations of recombinant inbred lines of soybean (Glycine max L. Merr. ). Theor Appl Genet, 2006, 113: 1459-1466[本文引用:3]
[12]	刘晓芬. 大豆栽培品种群体粒形性状及百粒重的关联分析. 南京农业大学博士学位论文, 江苏南京, 2012 Liu X F. Association Analysis for Seed Shape Traits and 100-seed Weight in Soybean (Glycine max L. Merr. ). PhD Dissertation of Nanjing Agricultural University, Nanjing, China, 2012 (in Chinese with English abstract)[本文引用:1]
[13]	Xu Y, Li H N, Li G J, Wang X, Cheng L G, Zhang Y M. Mapping quantitative trait loci for seed size traits in soybean (Glycine max L. Merr. ). Theor Appl Genet, 2011, 122: 581-594[本文引用:1]
[14]	牛远, 徐宇, 李广军, 王云清, 刘晓芬, 李河南, 魏世平, 章元明. 大豆籽粒大小和粒形的驯化研究. 大豆科学, 2012, 31: 522-528 Niu Y, Xu Y, Li G J, Wang Y Q, Liu X F, Li H N, Wei S P, Zhang Y M. Domestication of seed size and shape traits in soybean. Soybean Sci, 2012, 31: 522-528 (in Chinese with English abstract)[本文引用:2]
[15]	牛远, 谢芳腾, 布素红, 谢尚潜, 韩世凤, 耿青春, 刘兵, 章元明. 大豆粒形性状QTL的精细定位. 作物学报, 2013, 39: 609-616 Niu Y, Xie F T, Bu S H, Xie S Q, Han S F, Geng Q C, Liu B, Zhang Y M. Fine mapping of quantitative trait loci for seed shape traits in soybean. Acta Agron Sin, 2013, 39: 609-616 (in Chinese with English abstract)[本文引用:2]
[16]	Xie F T, Niu Y, Zhang J, Bu S H, Zhang H Z, Geng Q C, Feng J Y, Zhang Y M. Fine mapping of quantitative trait loci for seed size traits in soybean. Mol Breed, 2014, 34: 2165-2178[本文引用:3]
[17]	Hu Z B, Zhang H R, Kan G Z, Ma D Y, Zhang D, Shi G X, Hong D L, Zhang G Z, Yu D Y. Determination of the genetic architecture of seed size and shape via linkage and association analysis in soybean (Glycine max L. Merr. ). Genetica, 2013, 141: 247-254[本文引用:1]
[18]	Yu S B, Li J X, Xu C G. Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci USA, 1997, 94: 922-931[本文引用:4]
[19]	孙亚男, 仕相林, 蒋洪蔚, 孙殿军, 辛大伟, 刘春燕, 胡国华, 陈庆山. 大豆百粒重QTL的上位效应和基因型×环境互作效应分析. 中国油料作物学报, 2013, 34: 598-603 Sun Y N, Shi X L, Jiang H W, Sun D J, Xin D W, Liu C Y, Hu G H, Chen Q S. Epistatic effects and QE interaction effects of QTLs for 100-seed weight in soybean. Chin J Oil Crop Sci, 2013, 34: 598-603 (in Chinese with English abstract)[本文引用:1]
[20]	Xing Y Z, Tan Y F, Hua J P, Sun X L, Xu C G, Zhang Q F. Characterization of the main effects, epistatic effects and their environmental interactions of QTLs on the genetic basis of yield traits in rice. Theor Appl Genet, 2002, 105: 248-257[本文引用:1]
[21]	Niu Y, Xu Y, Liu X F, Yang S X, Wei S P, Xie F T, Zhang Y M. Association mapping for seed size and shape traits in soybean cultivars. Mol Breed, 2013, 31: 785-794[本文引用:4]
[22]	梁慧珍, 余永亮, 杨红旗, 张海洋, 董薇, 杜华, 崔暐文, 刘学义, 方宣钧. 大豆粒形性状主效QTL环境互作和上位性检测. 中国农业科学, 2013, 46: 5081-5088 Liang H Z, Yu Y L, Yang H Q, Zhang H Y, Dong W, Du H, Cui W W, Liu X Y, Fang X J. Main, environmentally interacted and epistatic QTL for seed shape traits in soybean. Sci Agric Sin, 2013, 46: 5081-5088 (in Chinese with English abstract)[本文引用:2]
[23]	闫宁, 谢尚潜, 耿青春, 徐宇, 李广军, 刘兵, 汪霞, 李其刚, 章元明. 利用Bayes分层广义线性模型剖析大豆籽粒性状的遗传基础. 作物学报, 2013, 39: 258-268 Yan N, Xie S Q, Geng Q C, Xu Y, Li G J, Liu B, Wang X, Li Q G, Zhang Y M. Genetic basics of seed traits in soybean with bayes hierarchical gengralized linear model method. Acta Agron Sin, 2013, 39: 258-268 (in Chinese with English abstract)[本文引用:1]
[24]	苏成付, 赵团结, 盖钧镒. 不同统计遗传模型QTL定位方法应用效果的模拟比较. 作物学报, 2010, 36: 1100-1107 Su C F, Zhao T J, Gai J Y. Simulation comparisons of effectiveness among QTL mapping procedures of different statistical genetic models. Acta Agron Sin, 2010, 36: 1100-1107 (in Chinese with English abstract)[本文引用:2]
[25]	李杰勤, 张启军, 叶少平, 赵兵, 梁永书, 彭勇, 吴发强, 王世全, 李平. 四种不同QTL作图方法的比较研究. 作物学报, 2005, 31: 1473-1477 Li J Q, Zhang Q J, Ye S P, Zhao B, Liang Y S, Peng Y, Wu F Q, Wang S Q, Li P. Comparative research on four mapping methods of QTLs. Acta Agron Sin, 2005, 31: 1473-1477 (in Chinese with English abstract)[本文引用:2]
[26]	Wang S C, Basten C J, Zeng Z B. Windows QTL Cartographer 2. 5. North Carolina State University, Raleigh, NC, USA, 2001-2006[本文引用:1]
[27]	Yang J, Zhu J, Williams R W. Mapping the genetic architecture of complex traits in experimental populations. Bioinformatics, 2007, 23: 1527-1536[本文引用:1]
[28]	Li H H, Ye G Y, Wang J K. A modified algorithm for the improvement of composite interval mapping. Genetics, 2007, 175: 361-374[本文引用:1]
[29]	Meng L, Li H H, Zhang L Y, Wang J K. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J, 2015, 3: 269-283[本文引用:1]
[30]	Nyquist W E, Baker R J. Estimation of heritability and prediction of selection response in plant populations. Crit Rev Plant Sci, 1991, 10: 235-322[本文引用:1]
[31]	Lu W G, Wen Z X, Li H C, Yuan D H, Li J Y, Zhang H. Identification of the quantitative trait loci (QTL) underlying water soluble prote incontent in soybean. Theor Appl Genet, 2013, 126: 425-433[本文引用:1]
[32]	陈强. 大豆籽粒相关性状QTL定位分析. 河北科技师范学院硕士学位论文, 河北秦皇岛, 2014 Chen Q. QTL Mapping for seed related traits in soybean (Glycine max L. Merr. ). MS Dissertation of Hebei Normal University of Science & Technology, Qinhuangdao, China, 2014 (in Chinese with English abstract)[本文引用:1]
[33]	雷雅坤, 闫龙, 杨春燕, 宋晓坤, 张孟臣, 黄占景. 大豆公共遗传图谱C1连锁群SSR标记空白区段的填补. 华北农学报, 2012, 27(6): 5-10 Lei Y K, Yan L, Yang C Y, Song X K, Zhang M C, Huang Z J. Complete the blank section with SSR markers on linkage group C1 of public genetic map in soybean. Acta Agric Boreali-Sin, 2012, 27(6): 5-10 (in Chinese with English abstract)[本文引用:1]
[34]	Wang S C, Basten C J, Zeng Z B. Windows QTL Cartographer V2. 5. Department of Statistics. North Carolina: North Carolina State University, 2007[本文引用:1]
[35]	Yang J, Zhu J. Predicting superior genotypes in multiple environments based on QTL effects. Theor Appl Genet, 2005, 110: 1268-1274[本文引用:1]
[36]	Yang J, Hu C C, Hu H, Yu R D, Xia Z, Ye X Z, Zhu J. QTL Network: mapping and visualizing genetic architecture of complex traits in experimental populations. Bioinformatics, 2008, 24: 721-723[本文引用:1]
[37]	张红梅, 李海朝, 文自翔, 顾和平, 袁星星, 陈华涛, 崔晓艳, 陈新, 卢为国. 大豆籽粒维生素E含量的QTL分析. 作物学报, 2015, 41: 187-196 Zhang H M, Li H C, Wen Z X, Gu H P, Yuan X X, Chen H T, Cui X Y, Chen X, Lu W G. Identification of QTL associated with vitamin e content in soybean seeds. Acta Agron Sin, 2015, 41: 187-196 (in Chinese with English abstract)[本文引用:1]
[38]	McCouch S, Chen X, Panaud O, Temnykh S, Xu Y, Cho Y G, Huang N, Ishii T, Blair M. Microsatellite marker development, mapping and applications in rice genetics and breeding. Plant Mol Biol, 1997, 35: 89-99[本文引用:1]
[39]	Song Q J, Marek L F, Shoemaker R C, Lark K G, Concibido V C, Delannay X, Specht J E, Cregan P B. A new integrated genetic linkage map of the soybean. Theor Appl Genet, 2004, 109: 122-128[本文引用:1]
[40]	Mansur L M, Orf J, Chase K, Jarvik T, Cregan P B, Lark K G. Genetic mapping of agronomictraits using recombi-nant inbred lines of soybean. Crop Sci, 1996, 36: 1327-1336[本文引用:1]
[41]	杜景红, 樊叶杨, 吴季荣, 庄杰云. 水稻第6染色体短臂产量性状QTL簇的分解. 中国农业科学, 2008, 41: 939-945 Du J H, Fan Y Y, Wu J R, Zhuang J Y. Dissection of QTLs for yield traits on the short arm of rice chromosome 6. Sci Agric Sin 2008, 41: 939-945 (in Chinese with English abstract)[本文引用:1]