雷雅坤^1,2,**, 刘兵强^1,*, 邸锐¹, 闫龙¹, 杨春燕¹, 郝东旭³, 张孟臣^1,*
1 河北省农林科学院粮油作物研究所 / 国家大豆改良中心石家庄分中心 / 河北省遗传育种重点实验室, 河北石家庄 050031;

² 河北省农林科学院农业信息与经济研究所, 河北石家庄 050051;

³ 河北经济管理学校, 河北石家庄 050071

^*通讯作者(Corresponding author): 张孟臣, E-mail: zhangmengchen@hotmail.com, Tel: 0311-87670653 第一作者联系方式: E-mail: yakun0628@163.com ^**同等贡献(Contributed equally to this work)
收稿日期:2015-07-09 接受日期:2015-11-20网络出版日期:2015-12-08基金:本研究由国家高技术研究发展计划(863计划)项目(2012AA101106), 国家自然科学基金项目(31000719), 国家转基因生物新品种培育重大专项(2014ZX0800402B), 国家现代农业产业技术体系建设专项(CARS-004-PS06), 河北省博士基金项目(F13E006)和河北省青年拔尖人才支持计划资助

摘要大豆油的品质取决于脂肪酸各组分在大豆中的比例, 为发掘控制大豆5种脂肪酸含量的数量性状位点(QTL), 利用冀豆12和黑豆重组自交系群体构建遗传图谱, 采用Windows QTL Cartographer 2.5和QTL Network-2.0软件的CIM和MCIM法对大豆5种脂肪酸组分进行数量性状定位。结果表明, 在石家庄和三亚各环境下共检测到16个QTL, 位于连锁群A2、B2、C2、F、G、I、L上。对2个环境联合分析, 检测到13个QTL, 其中9个用2种方法被检测到, 但这13个位点与环境互作的贡献率明显小于加性效应。其中在B2连锁群Satt168~Satt556控制硬脂酸的QTL Ste-1在河北石家庄和海南三亚均能被检测到, 贡献率均为12%, 在双尾群体和间隔挑选群体中也能检测到控制硬脂酸的QTL Ste-1, 说明这一QTL稳定存在于本组合群体中, 为今后大豆硬脂酸的QTL精细定位奠定了基础。

关键词:大豆; 脂肪酸; 数量性状位点; 环境
QTL Analysis of Fatty Acid Contents under Different Environments in Soybean
LEI Ya-Kun^1,2,**, LIU Bing-Qiang^1,**, DI Rui¹, YAN Long¹, YANG Chun-Yan¹, HAO Dong-Xu³, ZHANG Meng-Chen^1,*
1Institute of Cereal and Oil Crops, Hebei Academy of Agriculture and Forestry Sciences / Shijiazhuang Branch Center of National Center for Soybean Improvement / Hebei Genetic Breeding Laboratory, Shijiazhuang 050031, China

²Institute of Agricultural Information and Economy, Hebei Academy of Agriculture and Forestry Sciences, Shijiazhuang 050051, China

³Hebei Econmy Management School, Shijiazhuang 050071, China

Fund:This study was support by the National High Technology Research and Development Program of China (863 Program) (2012AA101106) and the National Natural Science Foundation of China (31000719), the Major Project of China on New Varieties of GMO Cultivation (2014ZX0800402B) and the China Agriculture Research Syscem (CARS-004-PS06), the Doctoral Scientific Research Foundation of Hebei Province (F13E006) and the Program for Excellent Young Talents in Hebei Province
AbstractSoybean is an important crop which contributes about 60% of the world’s oilseed production. The quality of soybean oil depends on the relative composition of fatty acid in seeds. A population of recombinant inbred lines derived from a cross of Jidou 12 × Heidou was grown in two environments and the seed samples from the environments were evaluated for fatty acid contents. Simple sequence repeat (SSR) markers were used to construct genetic linkage map. A total of 16 QTLs underlying fatty acid contents were identified by CIM and MCIM methods using Windows QTL Cartographer 2.5 and QTL Network-2.0 software at Sanya and Shijiazhuang locations, respectively. These QTLs were scattered on linkage groups A2, B2, C2, F, G, I, L. According to the two environment combined data, 13 QTLs were detected using two mapping methods. Nine of them were common in the results from two methods. Another QTL associated with stearic acid content, named as Ste-1, located on LG B2 and flanked by Satt168 and Satt556, was stable across two locations, also. QTL Ste-1 could explain 12% of the phenotypic variation at both locations. This study is helpful to improve fatty acid composition in soybean.

Keyword:Soybean; Fatty acids; Quantitative trait loci; Environment
Show Figures
Show Figures






大豆油是主要的食用油之一, 大豆的蛋白、矿物质和维生素含量丰富, 其不饱和脂肪酸含量在油分中也具有很高的比例。因此, 应重视大豆脂肪酸各组分的比例, 以便开展大豆品质育种, 提高脂肪含量^[1]。大豆棕榈酸、硬脂酸、油酸、亚油酸与亚麻酸5种脂肪酸含量都属于数量性状^[2]。Nickell和Brummer^[3]都用RFLP标记分别在连锁群D和B2上定位到了控制棕榈酸和亚麻酸的位点。1992年Diers等^[4]检测到影响棕榈酸、油酸和亚油酸的3个QTL, 解释10%以上的遗传变异。2002年Li等^[5]通过Cook× N87-2122-4构建的重组自交系F₂和F_2:3群体, 在连锁群上检测到了控制棕榈酸的2个QTL, 分布在连锁群A1标记Satt684附近和M连锁群标记Satt175附近, F₂群体在这2个位置上可解释的遗传变异为38%和8%, F_2:3群体在这2个位置上可解释的遗传变异分别为33%和9%; 所检测到的2个连锁群上的区域的主效基因和微效基因互作效应显著。2003年Spencer等^[6]通过群体检测到影响硬脂酸的QTL, 位于B2连锁群上标记Satt070、Satt474和Satt556附近。2006年郑永战等^[7]用Essex× ZDD2315的回交群体F₁的114个单株, 检测到控制5种脂肪酸含量的16个位点, 分布于B2、C1、D1b、D2、E、H、I、L、N等9个连锁群上, 可解释的遗传变异为8.2%~39.3%, 其中在B2上检测到的控制硬脂酸的位点和Spencer检测到的位点一致。在不同群体中所定位到的控制脂肪酸含量的QTL有所差异, 通过对所报道的位点进行整合, 可以从整体上分析大豆脂肪酸的QTL, 选出比较可靠的位点用于育种工作。1996年Orf等^[8]对报道过的14个群体进行整理, 在5个群体里控制大豆脂肪酸的QTL有12个遗传贡献率超过了10%。2009年宋万坤等^[9]对大豆脂肪酸组分QTL进行了元分析, 影响大豆5种脂肪酸组分的QTL主要分布在A1、B2、D1b、D2、E、G、L这7条连锁群上, 且成簇分布。其中影响大豆棕榈酸、硬脂酸、油酸、亚油酸、亚麻酸含量的QTL分别为3、1、5、2和1个, 还有1个QTL同时影响大豆的棕榈酸、硬脂酸和油酸含量。于福宽^[10]利用鲁黑豆2号和南汇早黑豆重组自交系群体对棕榈酸、硬脂酸、油酸、亚油酸、亚麻酸、油分等性状检测到13个QTL。分布于A1、M、G、O、K、B2和N连锁群上, LOD值为2.62~12.70, 遗传贡献率为5.60%~ 21.39%。
目前, 研究者们比较了不同环境中以同一群体定位到的QTL。朱军^[11]提出的基于混合线性模型复合区间作图法为分析QTL与环境的互作效应提供了有力手段, 可以直接分析出QTL与环境互作的遗传贡献率。单大鹏等^[12]用Charleston× 东农594的RIL群体连续5年对油份进行QTL定位, 检测到11个位点, 其中有2个QTL没有与环境的互作。苗兴芬^[13]用Charleston× 东农594的群体通过3个年份3个地点的试验, 检测到5种脂肪酸的68个位点, 其中有19个存在与环境的互作。
本试验利用冀豆12和黑豆组合的F₆群体构建更完整的大豆遗传连锁图谱, 连续2年对5种脂肪酸组分含量进行QTL定位, 在三亚环境下对石家庄定位到的数量性状位点进行验证, 找到更加可靠的控制5种脂肪酸组分含量的QTL。对在石家庄和三亚检测到的脂肪酸含量进行了上位效应和QTL与环境互作效应的分析, 有利于得到受环境影响较小, 能稳定遗传的QTL。
1 材料与方法1.1 试验材料以5种脂肪酸含量均有差异的冀豆12为母本, 黑豆为父本进行杂交。通过单粒传法得到本群体F₆的211个家系。2010年冬将收获的F₆单株种植于海南三亚实验站, 每个家系1行, 行长2 m, 株行混收, 检测了211个家系脂肪酸主要组分含量。2011年夏, 将收获的F₇种植于河北石家庄堤上实验站, 对F₈进行脂肪酸主要组分含量的测定。根据F₇群体硬脂酸含量的表型数据, 从211个家系中选取含量最高和最低的各30株建立双尾群体; 每隔3株挑选1株建立间隔挑选群体。
1.2 PCR反应以本群体211个家系F₆单株的新鲜嫩叶为试验材料, 用SDS法提取大豆基因组DNA, PCR体系含5 μ L DNA、2.0 μ L 10× PCR buffer、2.5 mmol L^-1 dNTPs、1.5 μ mol L^-1引物、1 U Taq酶, 加ddH₂O补足至20 μ L。反应程序为94℃预变性5 min; 94℃变性30 s, 47℃退火30 s, 72℃延伸30 s, 35个循环; 72℃延伸5 min, 4℃保存。
1.3 脂肪酸测定方法采用安捷伦气相色谱6890N, 将供试材料用研磨仪充分研磨, 磨好后取0.04 g豆粉放入2.0 mL离心管中, 加1 mL脂肪提取液, 振荡混匀10 min, 10 625× g离心后取上清液, 加0.5 mol L^-1甲醇钠甲酯化1 mL, 静置10 min, 加饱和NaCl溶液至样品分层。吸取上清液进行气相色谱分析。选用FID检测器, FFAP弹性石英毛细管柱, 恒压的模式, 载气为高纯氮气, 流速为1 mL min^-1; H₂气流速为35 mL min^-1; 空气流速为350 mL min^-1。进样量为1 μ L。分流比为10∶ 1, 进样口的温度为250℃。检测器的温度为250℃。程序温为180℃ 3 min, 9℃ min^-1升到225℃ min^-1, 保持25 min^[14]。
1.4 QTL定位分析方法利用Map Manager QTXb20 (P=0.0001, 函数Kosambi)软件对冀豆12× 黑豆群体进行遗传连锁图谱的构建和分析。用Windows QTL Cartographer 2.5软件及复合区间作图法, 取LOD=2.5为阈值, Walk Speed=1.5 cM, 对由211个家系组成的大群体、双尾群体和间隔挑选群体的5种脂肪酸含量进行QTL检测, 同时分析各位点的遗传贡献率和加性效应。用QTL Network-2.0软件及基于混合线性模型的复合区间作图法, 对2个环境下的5种脂肪酸含量进行联合分析, 并分析所检测到位点的加性效应和与环境的互作效应。

2 结果与分析2.1 五种脂肪酸的分布特征从表1可看出, 在不同环境中2个亲本的5种脂肪酸组分含量有一定的变化。后代群体的5种脂肪酸含量中除三亚环境的油酸最小值为14.86%, 比低亲的父本黑豆高0.81个百分点外, 其余性状都具有双向超亲的现象, 且群体后代的平均值均位于两亲本之间, 在2个环境中硬脂酸的平均值分别为3.36%和3.15%, 均倾向于低亲值, 油酸的平均值分别为17.89%和17.19%, 均倾向于高亲值。在5种脂肪酸组分中, 硬脂酸的变异幅度较小, 油酸和亚油酸的变异幅度较大, 表明在这2个性状上后代群体有一定差异, 对提高油酸和亚油酸在大豆5种脂肪酸组分中所占的比例有一定意义。5种脂肪酸组分的最大值和最小值差异较明显, 说明在后代分离群体中的表现分离程度较大, 且变异连续, 在2个环境中5种脂肪酸组分偏度的绝对值都小于1, 基本呈正态分布, 是属于多基因控制的数量性状。
表1
Table 1
表1(Table 1)

表1 亲本及后代群体脂肪酸性状统计结果 Table 1 Statistics results of fatty acids traits of parents and descendant population 性状 Trait 环境 Environment 亲本 Parent 群体 Population 冀豆12 Jidou 12 黑豆 Heidou 最大值 Max 最小值 Mix 变幅 Range 平均值 Average 标准差 SD 偏度 Skew 峰度 Kurt 棕榈酸 三亚 Sanya 13.57 11.27 15.25 10.69 4.56 12.13 0.71 0.987 2.562 Palmitic acid 石家庄 Shijiazhuang 12.56 10.51 13.02 9.67 3.35 11.21 0.72 0.172 -0.578 硬脂酸 三亚 Sanya 3.48 3.35 4.10 2.44 1.66 3.36 0.30 -0.234 0.242 Stearic acid 石家庄 Shijiazhuang 3.48 2.82 3.83 2.54 1.29 3.15 0.27 0.024 -0.569 油酸 三亚 Sanya 19.31 14.05 24.82 14.86 9.96 17.89 1.84 0.936 1.189 Oleic acid 石家庄 Shijiazhuang 18.62 15.57 21.23 10.09 11.14 17.19 1.51 -0.579 2.595 亚油酸 三亚 Sanya 53.30 57.89 59.36 51.18 8.18 55.24 1.50 -0.047 0.018 Linoleic acid 石家庄 Shijiazhuang 56.35 59.50 62.28 54.59 7.69 57.72 1.28 0.160 0.450 亚麻酸 三亚 Sanya 10.34 13.43 13.84 8.95 4.89 11.38 1.01 0.010 -0.157 Linolenic acid 石家庄Shijiazhuang 9.00 11.66 12.95 8.03 4.92 10.67 0.94 0.045 -0.217

表1 亲本及后代群体脂肪酸性状统计结果 Table 1 Statistics results of fatty acids traits of parents and descendant population

2.2 不同世代不同地点5种脂肪酸组分的相关表2表明, F₇棕榈酸和油酸、亚油酸呈极显著负相关, 硬脂酸和亚麻酸呈显著负相关, 油酸和亚油酸、亚麻酸呈极显著负相关, F₇脂肪酸间的相关除了硬脂酸和亚麻酸的相关程度外, 其余均和F₈脂肪酸间的相关性一致。棕榈酸、硬脂酸、油酸、亚油酸和亚麻酸F₇和F₈之间都呈极显著正相关。说明F₇和F₈的脂肪酸组分含量受环境因素的影响较大。
表2
Table 2
表2(Table 2)

表2 群体F7和F8两代间脂肪酸的相关 Table 2 Correlatin coefficient of five types of fatty acids between F7 and F8 populations 脂肪酸 Fatty acid F7 F8 棕榈酸Palmitic 硬脂酸Stearic 油酸 Oleic 亚油酸Linoleic 亚麻酸Linolenic 棕榈酸Palmitic 硬脂酸Stearic 油酸 Oleic 亚油酸 Linoleic 棕榈酸 Palmitic (F7) 硬脂酸 Stearic (F7) -0.096 油酸 Oleic (F7) -0.223* * -0.025 亚油酸 Linoleic (F7) -0.232* * 0.009 -0.739* * 亚麻酸 Linolenic (F7) 0.078 -0.195* -0.563* * 0.021 棕榈酸 Palmitic (F8) 0.574* * 0.079 -0.059 -0.208* * -0.009 硬脂酸 Stearic (F8) 0.017 0.549* * 0.342* * -0.343* * -0.288* * 0.297* * 油酸 Oleic (F8) -0.140* -0.121* 0.499* * -0.387* * -0.202* * -0.215* * 0.186* 亚油酸 Linoleic (F8) -0.176* -0.068 -0.426* * 0.689* * -0.104* -0.253* * -0.390* * -0.560* * 亚麻酸 Linolenic (F8) 0.051 -0.045 -0.337* * -0.009 0.607* * -0.121* -0.318* * -0.416* * -0.132* * * Significant at the 0.01 level; * significant at the 0.05 level. * * 表示P< 0.01极显著相关; * 表示P< 0.05显著相关。

表2 群体F7和F8两代间脂肪酸的相关 Table 2 Correlatin coefficient of five types of fatty acids between F7 and F8 populations

2.3 遗传图谱构建参照Cregan等^[15]公布的大豆遗传连锁图谱, 在20条连锁群上, 约每10 cM挑选2~3对引物, 共筛选了495对引物, 其中有多态的为221对, 从中选取130对分布均匀的引物, 利用构建遗传连锁图谱, 能得到条带清晰的为118对, 与冀豆12带型一致的记为A, 与黑豆带型一致的记为B, 杂合的带型记为H, 缺失的记为-。
利用SSR标记对冀豆12× 黑豆群体构建大豆遗传连锁图谱, 共包括21条连锁群。与大豆公共图谱相比较, 连锁群与公共图谱有很好的对应, 标记顺序与公共图谱完全一致。遗传连锁图谱的标记为117个SSR标记, 每个连锁群上平均包括6个标记, 图谱总长度为1501.0 cM, 标记间平均距离为15.6 cM^[16]。
2.4 不同环境5种脂肪酸组分的QTL定位采用复合区间作图法分析石家庄和三亚2个环境大豆5种脂肪酸含量的数据, 共检测到16个QTL (图1和表3)。
在石家庄环境可以检测到4个控制硬脂酸的QTL, 分别位于连锁群B2、C2和I, 在三亚环境还可以检测到2个QTL, 分别位于连锁群B2和C2 (表3), LOD值为2.50~ 7.03, 遗传贡献率为11%~14%。其中QTL Ste-1在2个环境中都被检测到, 增效基因均来自母本冀豆12, 且可解释的遗传变异都为12%。
图1
Fig. 1

	Figure OptionViewDownloadNew Window
	图1 冀豆12× 黑豆中检测出的5种脂肪酸含量QTL在连锁群的分布Fig. 1 Distrbution of QTL for five types of fatty acids detected from Jidou 12× Heidou on linkage map

表3
Table 3
表3(Table 3)

表3 2年中用CIM法定位的5种脂肪酸含量QTL Table 3 QTLs mapping of five fatty acid contents with CIM method in two years 性状 Trait 位点 QTL 连锁群 Chr. 区间 Marker interval 加性效应 Additive effect LOD值 LOD score 贡献率 PEV (%) 位置 Position (cM) 棕榈酸 Palmitic acid (F8) Pal-1 F Satt663-Satt335 -0.17 2.69 6 40.26 硬脂酸 Stearic acid (F7) Ste-1 B2 Satt168-Satt556 -0.13 6.80 12 54.40 Ste-2 C2 Satt643-Satt286 0.16 6.49 11 66.74 硬脂酸 Stearic acid (F8) Ste-1 B2 Satt168-Satt556 -0.09 4.13 12 48.40 Ste-3 B2 Satt556-Satt534 -0.09 5.12 12 59.80 Ste-4 C2 Satt286-Satt460 0.10 7.03 14 77.00 Ste-5 I Satt496-Satt049 -0.07 2.50 6 37.03 油酸 Oleic acid (F7) Ole-1 L Satt278-Satt462 -0.52 2.53 5 38.43 油酸 Oleic acid (F8) Ole-2 I Satt162-Satt148 -0.47 4.23 9 75.44 Ole-3 I Sat_418-Satt162 -0.49 3.85 10 68.20 亚油酸 Linoleic acid (F8) Lei-1 C2 Satt457-Satt322 -0.45 4.85 12 42.94 Lei-2 I Sat_418 0.32 3.09 6 60.20 亚麻酸 Linolenic acid (F7) Len-1 I Satt162-Satt148 0.37 3.92 8 81.44 Len-2 C2 Sat_062-Satt281 0.28 2.44 5 14.01 亚麻酸 Linolenic acid (F8) Len-3 A2 Satt409-Satt429 0.29 3.57 10 36.94 Len-4 C2 Satt457-Satt322 0.28 3.07 7 131.50 Len-5 G Satt612-Sat_372 0.33 5.07 15 78.93 Negative values represent gene additive effect from the female parent Jidou 12, positive values from the male parent Heidou Heidou parent. 加性效应A中的负值代表增效基因来自母本冀豆12, 正值代表增效基因来自父本黑豆。

表3 2年中用CIM法定位的5种脂肪酸含量QTL Table 3 QTLs mapping of five fatty acid contents with CIM method in two years

在双尾群体和间隔挑选群体中, 均可定位到控制硬脂酸的QTL Ste-1, 遗传贡献率为15%~27%, 在双尾群体中定位到位于B2和C2上的QTL与石家庄环境中所检测到的位点一致(图2), 其中位于B2上的QTL Ste-3在间隔挑选群体中也被检测到, 说明分布在连锁群B2和C2上的控制硬脂酸的QTL比较稳定。在双尾群体I连锁群标记Satt148附近和在间隔挑选群体H连锁群标记Satt541- Sat-118检测到了控制硬脂酸的位点, 这2个QTL与冀豆12× 黑豆群体定位到的位点不同, 有待进一步验证。
图2
Fig. 2

	Figure OptionViewDownloadNew Window
	图2 复合区间作图法对硬脂酸的QTL定位似然比图谱Fig. 2 Likelihood ratio maps of QTL tagging for stearic acid using CIM

在石家庄环境检测到1个棕榈酸含量QTL (QTL Pal-1), 位于F连锁群, LOD值为2.69, 可解释6%的遗传变异。在石家庄环境检测到2个油酸含量QTL, 都分布在I连锁群的Satt162附近, 遗传贡献率分别为9%和10%, 其增效基因均来自母本冀豆12, 2个QTL距离较近, 在连锁群I上可能有控制油酸的主效QTL存在, 在海南环境检测到1个QTL, 位于连锁群L。在石家庄环境下只检测到2个亚油酸含量QTL, 分别位于连锁群C2和I, 可解释的遗传变异分别为12%和6%。在石家庄环境检测到3个控制亚麻酸的QTL, 分别位于连锁群A2、C2和G上, LOD值为2.44~5.07, 遗传贡献率为7%~15%。其中位于连锁群G标记Satt612-Sat_372之间的QTL Len-5, 在石家庄环境本群体的F₄中也曾被检测到, 其增效基因来自父本黑豆。在三亚环境检测到2个控制亚麻酸的QTL, 分别位于连锁群C2和I, 可解释的遗传变异为5%和8%。
2.5 联合分析2个环境下5种脂肪酸含量QTL与环境的互作利用F₆群体的分子数据及石家庄和三亚2个环境下5种脂肪酸含量的数据, 用QTLNetwork-2.1软件基于混合线性模型的复合区间作图法对2个环境的大豆脂肪酸组分含量进行联合分析(表4), 共检测到13个QTL, 可解释的遗传变异范围为3.8%~15.4%, 和CIM法相比, 2种方法都能检测到的QTL有9个, 但贡献率较CIM法有所下降。
13个QTL中有12个与环境具互作效应, 其效应值为0.0001~0.4458, 2个环境对单个位点的效应值大致相等, 但作用相反。位点与环境的互作均没有达到显著水平。其中分布在I连锁群标记Satt049~Sat_418控制亚油酸的QTL与环境互作效应的贡献率为0, 说明该位点与环境没有发生互作效应。控制亚油酸的QTLN Lei-2与环境互作的遗传贡献率最大, 为0.89%。13个QTL与环境互作效应总的遗传贡献率为3.04%, 明显小于自身加性效应的遗传贡献率。说明本群体所定位的QTL受环境影响较小, 比较稳定、可靠。从表4可知, 除了3个控制油酸位点的加性效应都为正值, 增效基因均来自父本外, 控制其余4种脂肪酸位点的加性效应同时存在正值和负值, 这也说明了5种脂肪酸组分除F₇三亚环境的油酸没有负向超亲外其余均有双向超亲分布, 且最大值和最小值之间差异显著。
表4
Table 4
表4(Table 4)

表4 2年中用MCIM法定位的5种脂肪酸含量QTL及其与环境的互作 Table 4 QTLs mapping of five fatty acid contents QE interaction with MCIM in two years 性状 Trait 位点 QTL 标记区间 Marker interval F值 F-value 加性效应 Additive effect H2A H2AEi AEi1 AEi2 棕榈酸 Palmitic acid NPal-1 Satt409-Satt429 10.62 -0.29 0.0554 0.0006 -0.0002 0.0002 NPal-2 Satt335-Satt144 20.21 0.27 0.1088 0.0008 0.0003 -0.0003 NPal-3 Satt607-Satt361 8.25 -0.13 0.0338 0.0005 0.0002 -0.0003 硬脂酸Steatic acid NSte-1 Satt286-Satt460 39.92 -0.19 0.1538 0.0006 0.0163 -0.0168 NSte-2 Satt556-Satt534 30.27 0.11 0.1149 0.0056 -0.0001 0.0001 油酸 Oleic acid NOle-1 Satt462-Satt481 9.33 0.44 0.0407 0.0019 0.0002 -0.0002 NOle-2 Satt317-Sct_193 9.80 1.81 0.0392 0.0060 0.4442 -0.4458 NOle-3 Sat_418-Satt162 11.92 0.74 0.0573 0.0007 0.0001 -0.0001 亚油酸 Linoleic acid NLei-1 Satt457-Satt322 11.06 0.29 0.0380 0.0007 -0.0006 0.0006 NLei-2 Satt317-Sct_193 13.03 -0.47 0.0529 0.0089 -0.1395 0.1406 NLei-3 Satt049-Sat_418 10.97 -0.47 0.0584 0.0000 0.0000 -0.0000 亚麻酸 Linolenic acid NLen-1 Satt409-Satt429 12.64 0.44 0.0663 0.0036 -0.0409 0.0410 NLen-2 Satt457-Satt322 17.95 -0.36 0.0771 0.0005 0.0001 -0.0001 AEi1, AEi2 indicate the additive effect of interaction of QTL with environments, in Sanya of 2010 and in Shijiazhuang of 2011, respectively. AEi1、AEi2分别表示2010年三亚和2011年石家庄2个环境下QTL与环境互作的加性效应。

表4 2年中用MCIM法定位的5种脂肪酸含量QTL及其与环境的互作 Table 4 QTLs mapping of five fatty acid contents QE interaction with MCIM in two years

3 讨论3.1 两年QTL定位结果比较由于大豆脂肪酸组分含量属于数量性状, 环境和年份对其影响较大, 从2年的脂肪酸组分含量定位结果可知, 控制5种脂肪酸组分含量的QTL较集中分布在B2、C2和I连锁群。其中, 在B2连锁群标记Satt168-Satt556控制硬脂酸的QTL在2个环境中可被重复检测到, 且2年的遗传贡献率都为12%。这一位点和Spencer等^[6]利用普通硬脂酸含量品种Dare和高硬脂酸含量品种FAM94-41的重组自交系F₂、F_2:3报道过的位点一致, 郑永战等^[7]也曾用高含油量品种Essex和低含油量品种ZDD2315的回交群体报道过这个位点。在不同年份地点、环境下得到的相同的基因位点, 是硬脂酸的重要基因位点。在石家庄环境B2连锁群上还可以检测到一个控制硬脂酸的QTL位点, 和2006年郑永战等^[7]利用回交群体报道的位点一致。本试验定位的区间和曾报道的相比, 遗传距离在公共图谱上缩短了1.13 cM。三亚环境下检测到的在L连锁群上控制油酸的位点和Maria等^[17]在2008年用G99~G725和N00~3350的F_2:3定位到的位点都与标记Satt462相连锁。在C2连锁群标记Satt286-Satt460控制硬脂酸的位点和在G连锁群上控制亚麻酸的位点曾在本群体F_2:4中被检测到^[18], 且这一位点和邹筱^[19]报道的与脂肪、硬脂酸、油酸、亚油酸含量相关的Satt079标记的遗传距离较近。由于数量性状受环境和群体的影响, 所以, 在不同群体和同一群体不同环境中重复检测到的位点比较可靠, 挖掘这些有利基因有着重要的意义。
3.2 遗传距离较近的QTL本试验检测到控制5种脂肪酸组分含量的QTL位点中, 在不同年份所检测到的控制同一性状的位点分布在同一个连锁群不同标记之间且遗传距离很近, 如在C2连锁群上控制硬脂酸的QTL, 在2008年三亚环境下和2009年石家庄环境下被定位在标记Satt316-Satt307, 在2010年三亚环境下被定位在标记Satt643-Satt286, 在2011年石家庄环境下被定位在标记Satt286-Satt460, 遗传贡献率分别为7%, 11%, 14%。在G连锁群上控制亚麻酸的QTL位点, 在2008年石家庄环境下被定位在标记Satt517- Satt612, 在2011年石家庄环境下被定位在标记Satt612- Sat_372, 遗传贡献率分别为20%和15%。
以上几对分布在连锁群C2和G上的QTL在公共图谱中的遗传距离很近, 由于环境对定位造成的影响, 这几个QTL可能是一个, 但在QTL的初定位中没有明确的结论, 是否为同一个QTL只能通过这一区域引物的加密, 在不同环境反复试验继续验证。Fulton等^[20]认为, 在多种环境中能同时检测到的QTL可能比那些只能在一种环境中检测到且效应值更高的QTL更有利用价值, 这样的QTL在转移到新的遗传背景或在不同的环境条件下评价时更稳定。同时分析多个环境下的数据, 能增大QTL的检测强度, 准确估计QTL的位置和效应^[21]。对于遗传距离较近的QTL, 单大鹏等^[10]将其归为基本一致的QTL, 但对可以归为一个QTL的遗传距离没有做出限制, 对于本试验中在硬脂酸和亚麻酸上检测到的这2个标记区域可以作进一步研究, 以利于说明这个问题。
3.3 大豆脂肪酸组分含量QTL的富集区与相关性在数量性状基因定位时经常将控制不同性状的基因定位到连锁群的同一区域, 本研究发现, 控制亚麻酸和与其相关的油酸、亚油酸3个性状的QTL有部分重叠现象(表3和图1)。在I连锁群标记Satt162-Satt148同时控制油酸和亚麻酸的含量, 在I连锁群标记Sat_418附近检测到控制油酸和亚油酸这2个性状的QTL。由表2可知, 在F₈中, 油酸和亚油酸、亚麻酸都呈极显著负相关, 与于福宽^[10]和邹筱等^[19]报道的结果一致, 说明这一相关在本群体F₇和不同群体中都稳定存在^{[22, 23]}, 试验结果表明, 在连锁群I上很可能存在控制这一相关的基因位点。在F₈中, 亚油酸和亚麻酸存在极显著负相关, 这一相关徐豹等^[24]曾报道过, C2连锁群Satt457-Satt322的QTL同时控制亚油酸和亚麻酸的含量, 说明2个性状间重叠的QTL与相关性具有一致性。这些位点从分子层面上对油酸、亚油酸和亚麻酸的相关性做出了更好的解释, 也证明了基因一因多效的现象。通过脂肪酸组分间的相关, 对其中一种组分的选择可以影响到遗传距离较近的另一种脂肪酸组分。本研究中具有多效性脂肪酸QTL集中在C2和I连锁群上。同时控制多个数量性状的可能是同一个基因位点, 也可能是多个基因紧密连锁共同对一个性状起作用^[25]。目前的QTL定位中精细定位的研究还很少, 随着分子技术的发展尤其是SNP的运用, 遗传连锁图谱的不断加密和完善, 可以把QTL准确定位到大豆遗传连锁图谱中, 把在同一个区域中定位到的距离较近的QTL分解开, 可以通过精细定位和图位克隆技术, 最终将有价值的基因应用到分子育种上^[26]。

The authors have declared that no competing interests exist.

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


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

[1]	宋晓昆, 张颖君, 闫龙, 杨春燕, 郑艳艳, 蒋春志, 荆慧贤, 张孟臣, 黄占景. 大豆脂肪酸组份相关、变异特点分析. 华北农学报, 2010, 25(增刊): 68-73 Song X K, Zhang Y J, Yan L, Yang C Y, Zheng Y Y, Jiang C Z, Jing H X, Zhang M C, Huang Z J. A study on correlation and variability of fatty acid composition contents of soybean cultivars. Acta Agric Boreali-Sin, 2010, 25(suppl): 68-73 (in Chinese with English abstract)[本文引用:1]
[2]	吕景良, 邵荣春, 吴百灵, 梁岐, 吴桂荣. 东北地区大豆品种资源脂肪酸组成的分析研究. 作物学报, 1990, 16: 349-355 Lü J L, Shao R C, Wu B L, Liang Q, Wu G R. Fatty acid compo-sition analysis of soybean germplasm resources in Northeast. Acta Agron Sin, 1990, 16: 349-355 (in Chinese with English abstract)[本文引用:1]
[3]	Brummer E C, Graef G L, Orf J, Wilcox J R, Shoemaker R C. Mapping QTL for seed protein and oil content in eight soybean populations. Crop Sci, 1997, 37: 370-378[本文引用:1]
[4]	Diers B W, Shoemaker R C. Restriction fragment length polymorphism analysis of soybean fatty acid content. J Am Oil Chem Soc, 1992, 69: 1242-1244[本文引用:1]
[5]	Li Z L, Wilson R F, Rayford W E, Boerma H R. Molecular mapping genes conditioning reduced palmitic acid content in N87-2122-4 soybean. Crop Sci, 2002, 42: 373-378[本文引用:1]
[6]	Spencer M M, Pantalone V R, Meyer E J. Mapping the Fas locus controlling stearic acid content in soybean. Theor Appl Genet, 2003, 106: 615-619[本文引用:2]
[7]	郑永战, 盖钧镒, 卢为国, 李卫东, 周瑞宝, 田少君. 大豆脂肪及脂肪酸组分含量的QTL定位. 作物学报, 2006, 32: 1823-1830 Zheng Y Z, Gai J Y, Lu W G, Li W D, Zhou R B, Tian S J. QTL mapping for fat and fatty acid composition contents in soybean. Acta Agron Sin, 2006, 32: 1823-1830 (in Chinese with English abstract)[本文引用:3]
[8]	Orf J H, Diers B W, Boerma H R. Genetic improvement: conventional and molecular based strategies. In: Boerma H R, Specht J E, eds. Soybeans: Improvement, Production, and Uses, 3rd Edn. Madison, WI: American Society of Agronomy Inc. , 2004. pp 417-450[本文引用:1]
[9]	宋万坤, 王晶, 朱命喜, 齐照明, 刘春燕, 陈庆山, 胡国华. 大豆脂肪酸组分相关QTL元分析. 大豆科学, 2009, 28: 774-778 Song W K, Wang J, Zhu M X, Qi Z M, Liu C Y, Chen Q S, Hu G H. Meta-analysis of fatty acid QTLs in soybean. Soybean Sci, 2009, 28: 774-778 (in Chinese with English abstract)[本文引用:1]
[10]	于福宽. 大豆种质脂肪酸主要组分鉴定与QTL标记定位. 中国农业科学院硕士学位论文, 北京, 2011 Yu F K. Mapping QTLs and Screening Germplasm for the Major Fatty Acid Components in Soybean Seeds. MS Thesis of Graduate School of Chinese Academy of Agricultural Sciences, Beijing, China, 2011 (in Chinese with English abstract)[本文引用:3]
[11]	朱军. 运用混合线性模型定位复杂数量性状基因的方法. 浙江大学学报(自然科学版), 1999, 33: 327-335 Zhu J. Methods for locating complex quantitative trait genes using hybrid linear model. J Zhejiang Univ (Nat Sci), 1999, 33: 327-335 (in Chinese with English abstract)[本文引用:1]
[12]	单大鹏, 齐照明, 邱红梅, 单彩云, 刘春燕, 胡国华, 陈庆山. 大豆油分含量相关的QTL间的上位效应和QE互作效应. 作物学报, 2008, 34: 952-957 Shan D P, Qi Z M, Qiu H M, Shan C Y, Liu C Y, Hu G H, Chen Q S. Epistatic effects of QTLs and QE interaction effects on oil content in soybean. Acta Agron Sin, 2008, 34: 952-957 (in Chinese with English abstract)[本文引用:1]
[13]	苗兴芬. 多年多点条件下大豆脂肪酸含量QTL分析. 东北农业大学硕士学位论文, 黑龙江哈尔滨, 2010 Miao X F. QTL Analysis of Fatty Acids Content of Soybean in Multiple Years and Sites. MS Thesis of Northeast Agricultural University, Harbin, China, 2010 (in Chinese with English abstract)[本文引用:1]
[14]	张颖君, 高慧敏, 蒋春志, 胡梦芸, 刘兵强, 李辉. 大豆种子脂肪酸含量的快速测定. 大豆科学, 2008, 27: 859-862 Zhang Y J, Gao H M, Jiang C Z, Hu M Y, Liu B Q, Li H. Fast analysis on fatty acids of soybean seed by gas chromatography. Soybean Sci, 2008, 27: 859-862 (in Chinese with English abstract)[本文引用:1]
[15]	Cregan P B, Jarvik T, Bush A L, Shoemaker R C, Lark K G, Kahler A L, Kaya N, van Toai T T, Lohnes D G, Chung J. An integrated genetic linkage map of the soybean genome. Crop Sci, 1999, 39: 1464-1490[本文引用:1]
[16]	雷雅坤, 闫龙, 杨春燕, 宋晓昆, 张孟臣, 黄占景. 大豆公共遗传图谱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]
[17]	Maria J M, Joseph W B, Roger H B. Molecular mapping and confirmation of QTLs associated with oleic acid content in N00-3350 soybean. Crop Sci, 2008, 48: 2223-2234[本文引用:1]
[18]	宋晓昆. 大豆脂肪酸组分遗传及QTL分析. 河北师范大学硕士学位论文, 河北石家庄, 2010 Song X K. Genetic Analysis and QTL Mapping of Fatty Acid Composition Contents in Soybean. MS Thesis of Hebei Normal University, Shijiazhuang, China, 2010 (in Chinese with English abstract)[本文引用:1]
[19]	邹筱. 大豆脂肪和脂肪酸主要组分含量QTL定位. 中国农业科学院研究生院硕士学位论文, 北京, 2014 Zou X. Quantitative Trait Loci of Oil and Main Fatty Acid[本文引用:2]
20	Components in Soybean. MS Thesis of Graduate School of Chinese Academy of Agricultural Sciences, Beijing, China, 2014 (in Chinese with English abstract)[本文引用:1]
[20]	Fulton T M, Beck Bunn T, Emmatty D, Eshed Y, Lopez J, Petiard V, Uhlig J, Zamir D, Tanksley S D. QTL analysis of an advanced backcross of Lycopersicon peruvianumto the cultivated tomato and comparisons with QTLs found in other wild species. Theor Appl Genet, 1997, 95: 881-894[本文引用:1]
[21]	Jansen R C, Van Ooijen J W, Stam P, Lister C, Dean C. Genotype-by-environment interaction in genetic mapping of multiple quantitative trait loci. Theor Appl Genet, 1995, 91: 33-37[本文引用:1]
[22]	陈霞. 黑龙江省主栽大豆品种脂肪、脂肪酸组份的测定及其相关性的分析. 大豆科学, 1996, 15: 92-95 Chen X. Determination of fat and fatty acid composition of soybean cultivars in Heilongjiang province and analysis of their correlation. Soybean Sci, 1996, 15: 92-95 (in Chinese with English abstract)[本文引用:1]
[23]	李文滨, 郑宇宏, 韩英鹏. 大豆种质资源脂肪酸组分含量及品质性状的相关性分析. 大豆科学, 2008, 27: 740-745 Li W B, Zheng Y H, Han Y P. Analysis of fatty acid composition and other quality traits in soybean varieties developed in Heilongjiang province. Soybean Sci, 2008, 27: 740-745 (in Chinese with English abstract)[本文引用:1]
[24]	徐豹, 庄炳昌, 路琴华, 胡传璞. 中国大豆主要生产品种蛋白质、脂肪及其组份的相关分析. 大豆科学, 1988, 7: 175-183 Xu B, Zhuang B C, Lu Q H, Hu C P. Correlation analysis between protein content and its composition, fat content and its composition in current soybean cultivar seeds in China. Soybean Sci, 1988, 7: 175-183 (in Chinese with English abstract)[本文引用:1]
[25]	张忠臣, 战秀玲, 陈庆山, 滕卫丽, 杨庆凯, 李文滨, 刘迎雪, 郭强, 张丽娜. 大豆QTL定位的研究进展. 大豆科学, 2004, 23: 222-227 Zhang Z C, Zhan X L, Chen Q S, Teng W L, Yang Q K, Li W B, Liu Y X, Guo Q, Zhang L N. A review of QTL research in soybean. Soybean Sci, 2004, 23: 222-227[本文引用:1]
[26]	单大鹏, 朱荣胜, 陈立君, 齐照明, 刘春燕, 胡国华, 陈庆山. 大豆蛋白质含量相关QTL间的上位效应和QE互作效应. 作物学报, 2009, 35: 41-47 Shan D Y, Zhu R S, Chen L J, Qi Z M, Liu C Y, Hu G H, Chen Q S. Epistatic effects and QE interaction effects of QTLs for protein content in soybean. Acta Agron Sin, 2009, 35: 41-47 (in Chinese with English abstract)[本文引用:1]