全 文 ：作物学报 ACTA AGRONOMICA SINICA 2015, 41(9): 13721383 http://zwxb.chinacrops.org/
ISSN 0496-3490; CODEN TSHPA9 E-mail: xbzw@chinajournal.net.cn

本研究由河南省科技创新杰出人才计划(114200510002), 国家转基因生物新品种培育科技重大专项(2014ZX0800402B), 河南省重点科
技攻关计划项目(132102110091)和河南省农业科学院专项资金(201315603, 201315615, 201218334)资助。
第一作者联系方式: E-mail: lhzh66666@163.com, Tel: 0371-65751589
Received(收稿日期): 2015-03-03; Accepted(接受日期): 2015-06-01; Published online(网络出版日期): 2015-06-23.
URL: http://www.cnki.net/kcms/detail/11.1809.S.20150623.1359.004.html
DOI: 10.3724/SP.J.1006.2015.01372
大豆异黄酮及其组分含量的遗传分析与 QTL 检测
梁慧珍 1 余永亮 1 杨红旗 1 许兰杰 1 董 薇 1 牛永光 1 张海洋 1
刘学义 2 方宣钧 3
1河南省农业科学院芝麻研究中心, 河南郑州 450002; 2山西省农业科学院经济作物研究所, 山西汾阳 032200; 3 海南省热带农业资源
开发利用研究所, 海南三亚 572025
摘 要: 以栽培大豆晋豆 23为母本, 以山西农家品种大豆灰布支黑豆为父本杂交衍生的 447个 RIL作为供试群体构
建遗传图谱, 利用高效液相色谱法定性、定量测定样品中的异黄酮及其组分含量。采用主基因+多基因混合遗传分离
分析法和 WinQTLCart 2.5复合区间作图法, 对大豆异黄酮及其组分含量进行混合遗传分析和 QTL定位。结果表明,
大豆苷、黄豆苷元、染料木素、染料木苷、大豆苷元和异黄酮总含量分别受 4、4、2、3、2 和 2 对主基因控制, 并
有多基因修饰。检测到 44个与大豆异黄酮及其组分含量相关的 QTL, 与大豆苷、染料木素、黄豆苷元、大豆苷元、
染料木苷和异黄酮总含量相关的 QTL 分别有 10、9、4、7、8 和 6 个。连续 2 年分别检测到与大豆苷、染料木苷、
黄豆苷元和异黄酮关联, 分别位于标记区间 satt430~satt359、satt038~satt570、satt197~sat_128和 satt249~satt285的稳
定表达 QTL, 可尝试用于分子标记辅助育种。
关键词: 大豆; 异黄酮含量; 遗传分析; QTL
Genetic Analysis and QTL Mapping of Isoflavone Contents and Its Com-
ponents in Soybean
LIANG Hui-Zhen1, YU Yong-Liang1, YANG Hong-Qi1, XU Lan-Jie1, DONG Wei1, NIU Yong-Guang1,
ZHANG Hai-Yang1, LIU Xue-Yi2, and FANG Xuan-Jun3
1 Sesame Research Center, Henan Academy of Agricultural Sciences, Zhengzhou 450002, China; 2 Industrial Crop Institute, Shanxi Academy of
Agricultural Sciences, Fenyang 032200, China; 3 Hainan Provincial Institute of Tropical Agriculture Resources, Sanya 572025, China
Abstract: A set of 447 recombinant inbred lines (RILs) derived from the cross between cultivars Jingdou 23 (female parent) and
Huibuzhi (semi-wild, male parent) was used to construct a new map. Isoflavone content and its components were quantitatively
and qualitatively evaluated by using high performance liquid chromatography (HPLC). We analyzed inheritance and detected
QTLs for isoflavone content and its components in soybean seeds using major gene plus polygene mixed inheritance analysis and
WinQTLCart 2.5 composite interval mapping. The results showed that daidzin, daidzein, genistein, genistin, glycitin, and total
isoflavone contents were controlled by four, four, two, three, two and two main-genes, respectively. However, polygene effects
were not detected in the study. Forty-four quantitative trait loci (QTLs) for isoflavone contents and its components were mapped,
including ten for daidzin, nine for genistein, four for daidzein, seven for glycitin, eight for genistin, and six for total isoflavone
content. The stable QTLs related to daidzin, genistin, glycitin, isoflavone content were respectively detected to be located in the
intervals of satt430–satt359, satt038–satt570, satt197–sat_128, and satt249–satt285 during two years, which could be used in
marker-assisted selection (MAS) for soybean breeding.
Keywords: Soybean; Isoflavone content; Genetic analysis; Quantitative trait loci
大豆异黄酮是大豆种子中积累的一类次生代谢
产物。它的组成和存在形式主要包括染料木素、黄
豆苷元、大豆苷元、染料木苷和大豆苷等。大豆异
黄酮具有特殊的生物效能, 具有弱雌激素活性、抗
第 9期 梁慧珍等: 大豆异黄酮及其组分含量的遗传分析与 QTL检测 1373


氧化活性、抗溶血活性和抗真菌活性, 能有效预防
和抑制白血病、骨质疏松、结肠癌、肺癌、骨癌、
胃癌、乳腺癌、前列腺癌、妇女更年期综合症等多
种疾病的发生, 尤其是对乳腺癌和前列腺癌有良好
的预防和治疗作用[1-4]。染料木素和黄豆苷元对抑制
癌细胞的增殖具有特异作用[5-9]。自然界中异黄酮资
源十分有限 , 大豆是唯一含有异黄酮的食物资源 ,
但含量普遍较低[10], 提高大豆异黄酮含量是解决该
问题的关键。研究大豆异黄酮及其组分含量的遗传
与QTL分析, 对于推动采用分子标记辅助育种手段
选育高异黄酮含量的大豆品种, 具有十分重要的的
现实意义。国内外学者对大豆异黄酮及其组分含量
进行了遗传分析和初步QTL定位研究。Li等[11]对大
豆异黄酮与生态因子关系进行了研究, Liang等[12-14]
对大豆异黄酮及其组分含量的配合力和杂种优势以
及大豆籽粒异黄酮遗传效应和基因定位研究发现 ,
大豆鼓粒成熟期较低的均温和较大的昼夜温差有利
于异黄酮的积累; 在出苗期较多的日照和较低的均
温有利于异黄酮的形成。异黄酮含量和黄豆苷元既
受加性效应又受非加性效应的控制, 染料木苷主要
受非加性效应控制, 染料木素、大豆苷元和大豆苷
这3个成分的遗传主要受加性效应控制。大豆异黄酮
主要受制于母体遗传效应, 其次为胚基因效应, 细
胞质效应影响较小。不同遗传体系的异黄酮基因主
效应明显大于环境互作效应。利用Essex×Forrest、
AC756×RCAT Angora、Essex×PI4375、Williams 82×
Essex、Zhongdou27×Jiunong 20、Peking×Tamahomare、
Hwangkeum×IT182932、Xiaoheidou×GR8836、LHD2
×NHZ等RIL群体, 已定位了一批QTL, 基本上覆盖了
大豆所有的基因组[15-26]。其中, Meksem等[15-16]利用
Essex×Forrest 100个株系的重组自交系(Recombinant
Inbred Lines, RIL)及结合150个多态性DNA标记, 定
位了12个异黄酮及其组分含量的QTL, 分别位于第
3、第5、第9、第11和第12染色体上; 11连锁群上与
染料木苷含量紧密连锁的主效QTL可解释的遗传变
异率高达50.2%; 3连锁群上与染料木苷和大豆苷含
量紧密连锁的主效QTL可解释的遗传变异率分别为
11.1%和10.3%; 5连锁群上与大豆苷含量紧密连锁的
主效QTL可解释的遗传变异率为9.6%。Primomo
等 [17]利用AC756×RCAT Angora 207个株系F4:6的重
组自交系群体检测到17个异黄酮及其组分含量的
QTL, 分别定位在第1、第3、第5、第6、第7、第9、
第12、第13、第16和第18染色体上; Zhang等[25]利用
184个Xiaoheidou×GR8836株系的F2:9-11重组自交系
群体(RIL), 检测到22个异黄酮及其组分含量的QTL,
分别定位在第3、第4、第7、第9、第10、第12、第
13和第17染色体上 , 遗传贡献率在4.48%~7.87%之
间。检测到3对影响异黄酮含量上位性QTL, 其中 ,
qEGp1和qEGw1影响黄豆黄苷含量, qETISOp1影响
大豆异黄酮总量。张晶莹等[26]利用LHD2 × NHZ重
组自交系群体(RIL), 检测到14个异黄酮及其组分含
量的QTL, 分别定位在第3、第4、第5、第7、第9、
第11、第12、第14、第15、第17和第20染色体上。
目前, 国内外已检测到的大豆异黄酮及其组分含量
QTL, 因遗传背景和环境条件的不同 , 表现稳定的
较少 , 因此大豆异黄酮分子标记辅助育种受到限
制。本研究以晋豆23×灰布支黑豆衍生的447个RIL
为试验材料 , 利用WinQTLCart2.5[27] (http://statgen.
ncsu.edu/qtlcart/WQTLCart-1.htm)混合线性模型分
析方法, 对大豆异黄酮及其组分含量进行QTL定位分
析; 采用主基因+多基因混合遗传分离分析方法 [28],
分析大豆异黄酮及其组分含量遗传规律, 对其遗传
基础全面剖析, 为高异黄酮大豆分子标记辅助育种
提供参考依据。
1 材料与方法
1.1 供试材料
作图亲本材料为晋豆23×灰布支黑豆。母本
晋豆23是山西省主栽大豆品种 , 丰产性好 ; 父本
灰布支黑豆是山西省农家品种 , 具有较多的原始
性。作图群体采用单粒混传法(single seed mutiple
descent, SSD) RIL。1998年夏季在山西省杂交获
得F0种子 , 种植后收获一株种子产量732粒的F1
单株 , 1998年 10月在海南岛单粒种成 732个 F2,
1999年2月单株收获每个F2形成F3株系。从F3代开
始 , 每个世代每个株系随机收获10株上的共20个
豆荚 , 每荚约2粒 , 共约40粒种子 , 随机选取其中
的20粒种成20个株行 , 形成下一个世代 , 以此类
推至高世代。
田间试验采取随机区组设计 , 2次重复 , 行长
4.0 m, 行距0.4 m。2011年和2012年分别在河南省农
业科学院现代农业科技试验示范基地(河南省原阳
县)种植。田间管理按照大田模式进行, 周边设保护
行, 一播全苗。
1.2 异黄酮及其组分含量测定及数据分析
从每个家系选取 20 g大豆籽粒 , 用旋风磨
1374 作 物 学 报 第 41卷


(Retsch ZM100, Φ = 1.0 mm, Rheinische, Germany)磨
粉; 准确称取100 mg大豆粉放入带有螺帽的有机玻
璃试管, 加入5 mL含0.1% (v/v)乙酸的70% (v/v)乙醇
水溶液; 室温下振荡提取12 h; 2795×g离心10 min,
上清液经0.2 μm 滤膜过滤, 4℃冰箱保存备用。采用
津岛高效液相色谱仪(LC-6A), 定性、定量测定样品
异黄酮及其组分含量[27-28]。手动进样, 进样量为10
μL。色谱柱为YMC-Pack, ODS-AM-303, 250.0 mm ×
4.6 mm I.D.; 柱温35℃; 流动相A为含0.1% (v/v)乙
酸的超纯水 , B为含0.1% (v/v)乙酸的乙腈水溶液 ;
梯度洗脱范围为13%~35%; 运行时间为70 min; 流
速为1.0 mL min–1; 检测波长为260 nm。
根据标准样品的保留时间和最大吸收光谱定性,
以大豆苷(daidzin, D)在260 nm波长的紫外吸收值为
基础 , 参照Sun等 [29]的方法计算样品中染料木素
(genistein, G)、黄豆苷元(glycitein, GL)、大豆苷元
(daidzein, DA)、染料木苷 (genistin, GE)、大豆苷
(daidzin, D)以及异黄酮总含量(isoflavone, IF), 利用
SPSS 19.0 (SPSS, Chicago, USA)完成异黄酮及其组
分含量的方差分析、相关分析和描述统计。
1.3 遗传连锁图谱的构建与 QTL分析
王珍[30]使用该 RIL 群体为材料, 构建了一个含
有227个 SSR 标记的大豆遗传连锁图谱。以该图谱
和大豆公共图谱中的标记为基础, 梁慧珍[31]重新整
合 , 整合后图谱全长2047.6 cM, 包括27个连锁群
232个标记位点。该图谱在9(K)、13(F)和16(J) 3条染
色体上, 均出现了2个间隙, 在12(H)染色体出现了1
个间隙, 因而形成了27个连锁群。本研究采用梁慧
珍[31]整合后的大豆分子遗传图谱对大豆异黄酮及其
组分含量进行 QTL定位分析。把不同年份作为环境
因子处理 , 根据 2个环境的表型数据 , 利用
WinQTLCart 2.5软件中复合区间作图法 (CIM)
Zmapqtl方法的 Model 6, 每2个 cM对各性状进行全
基因组扫描, 以确定各性状 QTL数目及其在染色体
上的位置, 通过逐步回归指定解释给定性状最大变
异的5个标记作为余因子(co-factor), 选取临界阈值
LOD=2.5, 检测每个环境下的 QTL效应。当 LOD≥
2.5时, 认为 QTL 存在; 如果临近位点间图距小于5
cM, 就初步认定是同一个 QTL。采取 McCouch 等[32]
方式命名。
1.4 数量性状主基因+多基因混合遗传模型分析
采用主基因+多基因混合遗传分离分析法 [28],
分析晋豆 23×灰布支黑豆及其衍生的 RIL 群体, 通
过极大似然法和 IECM 对混合分布的各世代、各成
分分布的参数作出估计, 通过 AIC 准则和一组适合
性测验, 从中选出最适遗传模型, 利用其各成分分
布参数通过最小二乘法估计相应的遗传参数。主基
因遗传率 h2mg(%) = σ2mg/σ2p, 多基因遗传率 h2pg(%) =
σ2pg/σ2p; 环境遗传率 h2e(%) = σ2e/σ2p。其中 σ2p为群
体表型方差; σ2mg为主基因遗传方差; σ2pg为多基因
遗传方差; σ2e 为多环境遗传方差; 各遗传方差间的
关系为 σ2p = σ2mg+σ2pg+σ2e。
2 结果与分析
2.1 亲本及家系根系性状在RIL群体中的分离及
表型分析
观察大豆苷、染料木素、染料木苷、黄豆苷元、
大豆苷元和异黄酮总含量6个性状在分离群体中的
分离并进行方差和遗传力分析(表1)表明, 2个亲本6
个性状在后代分离群体中表现出较大的分离程度 ,
在不同环境中均存在较大差异, RIL 群体的最小值
和最大值之间差异明显。各性状的广义遗传力(h2 =
Vg/VP)分别为35.28%、23.62%、37.15%、29.67%、
21.01%和19.26%, 为 QTL分析提供了较好的遗传背
景。表1中两个环境下异黄酮及其组分含量共6个性
状表型变异系数为0.339~0.568, 方差分析均呈现极
显著差异, 表明异黄酮及其组分含量受遗传和环境
因素的共同影响。同时, 从 RIL 群体的性状变异范
围来看, 6个性状都有超越双亲的株系出现, 说明大
豆异黄酮及其组分含量遗传基因分布在双亲中, 通
过杂交可以得到超亲分离的株系。
2.2 相关分析
表2表明 , 不同环境条件下的相关分析结果基
本一致。除黄豆苷元外, 异黄酮总量与其余4个组分
含量均表现出显著或极显著正相关 , 且这4个组分
含量间均表现出显著或极显著正相关。说明各性状
之间存在着既相互影响又相互制约的关联关系。相
关分析结果与异黄酮及其组分含量在RIL群体中的
分离及遗传分析结果基本一致。
2.3 异黄酮及其组分含量的主基因+多基因混合
遗传分析
选用 2011 年和 2012 年试验数据的平均值进行
主基因+多基因混合遗传分析, 依据 AIC 值小为最
佳的准则结合 Kolmogorov试验、均匀检验性检验和
Smirnov 检验选定最适遗传模型, 用最小二乘法估
计异黄酮及其组分含量的遗传参数(表 3)。
第 9期 梁慧珍等: 大豆异黄酮及其组分含量的遗传分析与 QTL检测 1375


表 1 RIL 群体亲本和群体家系性状描述统计与方差分析
Table 1 Descriptive statistics and analysis of variance for isoflavone contents in the RIL populations and their parents
性状
Trait
大豆苷
Daidzin
染料木素
Genistein
染料木苷
Genistin
黄豆苷元
Glycitein
大豆苷元
Daidzein
异黄酮总含量
Isoflavone
晋豆 23 Jindou 23 7.404 23.206 15.932 2.530 16.895 3.294 亲本
Parent 灰布支黑豆 Huibuzhiheidou 7.894 23.380 16.503 9.908 18.405 3.836
平均值 Mean 6.008 13.428 7.713 5.546 9.467 2.154
标准差 SD 2.036 5.525 3.421 2.081 5.378 0.827
变异系数 CV 0.339 0.411 0.444 0.375 0.568 0.384
最小值 Minimum 2.907 4.362 2.469 2.048 1.853 0.775
最大值 Maximum 9.593 23.380 16.503 9.908 23.144 3.895
偏度 Skewness 0.154 0.432 1.241 0.050 0.989 0.637
峰度 Kurtosis –1.105 –0.614 2.034 –0.613 0.559 –0.016
RIL群体
RILs
F值 F-value 26.04** 32.01** 28.38** 41.05** 10.13** 8.52**
遗传力 h2(%) 35.28 23.62 37.15 29.67 21.01 19.26
**在 P < 0.01水平显著。**Significant at P < 0.01.

表 2 6 个性状之间的相关分析
Table 2 Phenotypic correlation coefficients among six traits in the RIL populations and their parents
性状
Trait
大豆苷
Daidzin
染料木素
Genistein
染料木苷
Genistin
黄豆苷元
Glycitein
大豆苷元
Daidzein
异黄酮
Isoflavone
大豆苷 Daidzin 0.53** 0.50* –0.02 0.49* 0.52*
染料木素 Genistein 0.56** 0.76** 0.22 0.90** 0.89**
染料木苷 Genistin 0.48* 0.61** 0.33 0.90** 0.95**
黄豆苷元 Glycitein –0.07 0.36 0.41 0.34 0.38
大豆苷元 Daidzein 0.50* 0.71** 0.75** 0.29 0.97**
异黄酮 Isoflavone 0.55** 0.69** 0.91** 0.32 0.81**
*和**按双侧检验, 分别表示 0.05 和 0.01 水平相关性显著和极显著; 右上角数据为 2011年分析结果, 左下角数据位 2012年分析
结果。
* and **: significance at the 0.05 and 0.01 levels (2-tailed), respectively. Correlation coefficient in 2011 and 2012 are listed in the top
right and lower left corners, respectively.

(1) 大豆苷最适遗传模型符合4MG-AI, 即4对
主基因加性-加性×加性上位性遗传模型。4对主基因
间加性效应分别为–1.390、–0.564、0.042和–0.044;
主基因之间相互作用的上位性效应分别为2.176、
2.207、0.403、2.210、1.272和–0.201。第3对主基因
对第4对主基因(icd)上位性效应为负值, 与其他基因
间上位性效应相反, 说明不同基因间上位性效应对
大豆苷性状的影响并不完全一致。主基因遗传率为
98.34%, 没有检测到多基因效应, 环境因素引起的
变异为1.66%, 表明大豆苷主要受4对主基因影响 ,
受环境因素影响较小。
(2) 染料木素最适遗传模型符合MX2-DominanceI-
A, 即2对显性上位主基因-加性多基因混合作用遗
传模型。2对显性上位主基因间相互作用的加性效
应分别为 2.052和 10.593。多基因间加性效应为
–7.436, 其主基因遗传率为67.15%, 多基因遗传率
为32.70%, 说明染料木素存在多基因效应, 其中主
基因效应更大。
(3) 染料木苷最适遗传模型符合 MX3-AI-AI,
即3对加性-上位性主基因×加性-上位性多基因混合
遗传模型。3对主基因间加性效应分别为–3.045、
–1.072和–2.284。3对基因中每2对基因间加性效应之
间的上位性互作效应分别为0.412、–0.206和0.719, 3
对基因间上位性互作效应(iabc)为0.059, 表明3对基
因之间存在着明显上位性互作, 并且互作后的上位
性效应不同。正向的上位性效应表示两座位间互作
基因型与具有正效应的加性基因型方向相同, 负向
的上位性效应则表示两座位间互作基因型与具有正
效应的加性基因型方向相反。3对主基因遗传率为
97.75%, 多基因遗传率为2.16%, 说明染料木苷存在
多基因效应, 其中以主基因效应为主。
(4) 黄豆苷元最适遗传模型符合4MG-AI, 即4
对主基因加性-加性×加性上位性遗传模型。4对主基
因间加性效应分别为1.687、1.364、0.821和0.487。4
1376 作 物 学 报 第 41卷


表 3 6 个性状最适模型及遗传参数估计结果
Table 3 Analysis of the best models and genetic parameters for six traits
性状和遗传模型 Trait and genetic model
大豆苷
Daidzin
染料木素
Genistein
染料木苷
Genistin
黄豆苷元
Glycitein
大豆苷元
Daidzein
异黄酮
Isoflavone
参数
Parameter
4MG-AI MX2-DominanceI-A MX3-AI-AI 4MG-AI MX2-Inhibiting-A MX2-AI-A
一阶参数 1st order parameter
M 3.004 17.992 21.467 6.429 13.303 2.678
d(da) –1.390 2.052 –3.045 1.687 — –2.531
db –0.564 10.593 –1.072 1.364 — –9.147
dc 0.042 — –2.284 0.821 — —
dd –0.044 — — 0.487 — —
iab(i*) 2.176 — 0.412 1.410 4.346 3.210
iac 2.207 — –0.206 1.106 — —
iad 0.403 — — 0.734 — —
ibc 2.210 — 0.719 1.496 — —
ibd 1.272 — — 0.670 — —
icd –0.201 — — 0.287 — —
iabc — — 0.059 — — —
[d] — –7.436 — — –0.755 6.883
二阶参数 2nd order parameter
σ2p 15.143 25.502 10.706 6.547 24.838 16.522
σ2mg 14.879 17.124 10.465 6.516 13.522 8.377
σ2pg — 8.339 0.231 — 11.295 8.117
h2mg (%) 98.26 67.15 97.75 99.53 54.44 50.24
h2pg (%) — 32.70 2.16 — 45.47 49.13

对主基因加性效应值均为正值, 表明基因来源于母
本晋豆23; 4对主基因间加性效应主基因之间相互作
用的上位性效应分别为1.410、1.106、0.734、1.496、
0.670和0.287, 说明不同基因间上位性效应对黄豆
苷元性状的影响方向一致。主基因遗传率为99.53%。
没有检测到多基因效应。环境因素引起的变异为
0.47%。表明黄豆苷元主要受4对主基因影响, 受环
境因素影响较小。
(5) 大豆苷元最适遗传模型符合MX2-Inhibiting-
A, 即2对抑制作用主基因-加性多基因混合遗传模
型。因为2对主基因间的抑制作用, 只检测到两对主
基因间上位性效应(iab)为4.346, 多基因的加性效应
为–0.755, 主基因遗传率为54.44%, 多基因遗传率
为45.47%, 说明大豆苷元存在多基因效应, 其中主
基因效应和多基因效应基本相当。
(6) 异黄酮总含量最适遗传模型符合MX2-AI-A,
即2对加性-上位性主基因×加性多基因混合遗传模
型。2对主基因间加性效应分别为–2.531和–9.147。2
对基因加性效应值均为负值, 表明基因来源于父本
灰布支黑豆。2对基因间加性效应之间的上位性互作
效应为3.210, 存在着明显上位性互作。主基因遗传
率为50.24%, 多基因遗传率为49.14%, 说明异黄酮
总含量存在多基因效应, 其中主基因效应和多基因
效应基本相等。
2.4 异黄酮及其组分含量的 QTL定位
利用Windows QTL Cart 2.5软件, 共检测到42个
控制大豆异黄酮及其组分含量的QTL, 分别位于
5(A1)、8(A2)、11(B1)、14(B2)、17(D2)、13(F_1)、
18(G)、20(I)、16(J_1)、9(K_2)、19(L)、3(N)和10(O)
染色体上 (表 4和图 1), 解释的表型变异范围为
3.74%~55.20%。加性效应值正负表现不一, 加性效
应值为正, 说明该QTL大豆异黄酮及其组分含量起
正向作用, 等位基因来自母本晋豆23, 加性效应值
为负, 说明该QTL大豆异黄酮及其组分含量起负向
作用, 定位在该染色体上的等位基因来自父本灰布
支黑豆。
共检测到9个大豆苷 QTL, 分别位于5(A1)、
8(A2)、11(B1)、14(B2)、16(J_1)、9(K_2)、19(L)、
第 9期 梁慧珍等: 大豆异黄酮及其组分含量的遗传分析与 QTL检测 1377


表 4 6 个性状 QTLs 位置及其参数
Table 4 QTL positions and its parameters for six traits
性状
Trait
年份
Year
QTL
染色体
Chr.
标记区间
Marker Interval
位置
Position
LOD
加性效应
Additive
R2
(%)
2011 qD-A1-1 A1-5 satt182–satt593 51.32 2.63 0.688 5.82
qD-A1-2 A1-5 satt593–satt454 56.77 3.15 0.700 3.74
qD-A2-1 A2-8 satt315–satt612 6.01 2.54 0.758 5.67
qD-B1-1 B1-11 satt430–satt359 124.87 4.45 –1.006 14.54
qD-B2-1 B2-14 satt416–satt474 14.54 3.22 1.810 10.19
2012 qD-A2-2 A2-8 satt538–satt429 153.20 4.93 –0.947 12.31
qD-B1-1 B1-11 satt430–satt359 124.87 4.36 –1.047 15.39
qD-J_1-1 J-16 satt249–satt285 0.05 5.49 –1.182 17.05
qD-L-1 L-19 satt388–satt182 14.02 3.03 0.674 5.60
大豆苷 Daidzin
qD-O-1 O-10 satt487–satt500 0.01 4.39 –1.066 11.65
2011 qG-A2-1 A2-8 satt187–sat_129 39.24 7.03 3.875 37.43
qG-B1-1 B-11 sat_129–satt519 75.90 3.02 1.712 7.26
qG-F_1-1 F-13 satt193–gmrubp 6.01 3.27 3.895 12.02
qG-G-1 G-18 satt275–satt038 14.01 6.56 3.721 35.16
2012 qG-A2-2 A2-8 satt329–satt333 66.49 2.97 –2.422 17.06
qG-G-2 G-18 satt038–satt570 17.82 6.74 3.721 35.16
qG-O-1 O-10 satt487–satt500 0.01 2.53 –2.528 16.95
qG-O-2 O-10 satt500–satt478 13.74 2.86 –3.435 26.17
染料木素 Genistein
qG-O-3 O-10 satt500–satt478 25.74 4.40 –4.052 29.01
2011 qGE-A2-1 A2-8 satt187–sat_129 39.24 4.94 2.356 35.61
qGE-A2-2 A2-8 satt377–satt329 64.49 5.35 –1.863 25.42
2012 qGE-G-1 G-18 satt038–satt570 15.82 6.50 2.366 39.59
染料木苷 Genistin
qGE-G-2 G-18 satt038–satt570 31.82 4.92 2.590 38.92
2011 qGL-A1-1 A1-5 satt545–satt511 126.35 3.62 –0.525 4.06
qGL-A2-1 A2-8 satt327–satt538 147.43 3.59 0.877 8.99
qGL-B1-1 B1-11 satt197–sat_128 57.79 10.62 1.640 41.42
2012 qGL-A2-1 A2-8 satt538–satt429 153.20 4.05 0.739 7.32
qGL-B1-1 B1-11 satt197–sat_128 73.79 9.24 1.640 41.12
qGL-K_2-1 K-9 satt264–sat_044 35.36 8.21 1.378 23.04
黄豆苷元 Glycitein
qGL-K_2-1 K-9 sat_044–satt137 41.68 3.89 1.328 23.14
2011 qDA-A2-1 A2-8 satt162–bsc 14.18 6.90 3.057 23.24
qDA-D2-1 D2-17 satt154–satt082 65.80 7.43 –3.453 32.02
qDA-D2-2 D2-17 satt082–satt226 69.69 7.82 –3.560 30.41
qDA-D2-3 D2-17 satt488–satt461 72.54 8.80 –3.562 30.41
2012 qDA-D2-4 D2-17 satt461–satt301 80.15 3.77 –3.539 30.40
qDA-I-1 I-20 satt292–satt182 69.13 3.70 –1.579 7.18
qDA-O-1 O-10 satt487–satt500 0.01 10.72 –4.236 5.09
大豆苷元 Daidzein
qDA-O-2 O-10 satt500–satt478 7.74 7.52 –4.216 5.01
2011 qIF-A2-1 A2-8 satt538–satt429 153.20 2.94 –0.392 13.64
qIF-D2-1 D2-17 satt154–satt082 47.80 3.68 –0.381 12.26
qIF-J_1-1 J-16 satt249–satt285 0.01 7.71 –0.761 55.14
2012 qIF-D2-2 D2-17 satt301–satt171 109.47 2.79 –0.394 12.14
qIF-O-1 O-10 satt259–satt331 88.12 2.71 –0.377 12.44
异黄酮 Isoflavone
qIF-J_1-1 J-16 satt249–satt285 0.01 7.72 –0.721 55.20

1378 作 物 学 报 第 41卷



图 1 检测到的 QTL 及加性效应 QTL 在连锁群上的分布以及与以往研究结果的对比
Fig. 1 Distribution of main QTLs and additive QTLs on linkage groups and comparison of QTLs detected in this study with those in
previous studies

03(N)和10(O)染色体上 , 解释的表型变异范围为
3.74%~17.05%。其中 2011年和 2012年均检测到
qD-B1-1定位在11(B1)连锁群标记satt430~satt359间
124.87 cM处 , 解释的表型变异分别为 14.54%和
15.39%, 说明这是一个能够稳定遗传的QTL。
qD-A2-1和qD-A2-2均被定位在8(A2)连锁群上 , 但2
个QTL位置相距较远, 说明8(A2)连锁群上可能存在
多个与大豆苷相关的基因。位于 16(J_1)染色体
satt405~satt249间的qD-J_1-1, 表型贡献率为17.05%,
解释的表型变异最大。
共检测到9个染料木素QTL, 分别位于8(A2)、
11(B1)、13(F_1)、18(G)和10(O)染色体上。解释的
表型变异范围为7.26%~37.43%。位于08(A2)染色体
satt187~sat_129间的qGE-A2-1, 表型贡献率为37.43%,
解释的表型变异最大。2011年检测到的qG-G-1和
2012年检测到的qG-G-2均被定位在18(G)染色体上,
第 9期 梁慧珍等: 大豆异黄酮及其组分含量的遗传分析与 QTL检测 1379


位置仅相距3.81 cM, 被定位在Satt038标记的两侧,
解释的表型变异均为 35.16%, 很可能是同一个
QTL。2011年检测到的qG-A2-1和2012年检测到的
qG-A2-2均被定位在8(A2)染色体上, 解释的表型变
异分别为37.43%和17.06%。2个QTL 位置相距较远,
很可能不是同一个QTL, 说明8(A2)连锁群上可能存
在多个与染料木素相关的基因。
共检测到4个染料木苷QTL, 分别位于8(A2)和
18(G)染色体上。解释的表型变异范围为25.42%~
39.59%。2011年检测到的qGE-G-1和2012年检测到
的qGE-G-2均被定位在18(G)染色体上, 且均被定位
在Satt038~satt570标记区间 , 解释的表型变异分别
均为39.59%和38.92%, 尽管位置相距16 cM, 仍怀
疑二者是同一QTL。qGE-G-1表型贡献率为39.59%,
解释的表型变异最大。
共检测到7个黄豆苷元QTL, 分别位于5(A1)、
8(A2)、11(B1)和9(K_2)染色体上, 解释的表型变异
范围为4.06%~41.42%。2011年检测到的qGL-B1-1和
2012年检测到的qGL-B1-2均被定位在11(B1)染色体
上, 且均被定位在satt197~sat_128标记区间, 解释的
表型变异分别均为41.42%和41.12%, 但位置相距
16 cM, 是否是同一QTL有待进一步验证。2011年检
测到的qGL-A2-1和2012年检测到的qGL-A2-2均被定
位在8(A2)染色体上, 位置仅相距5.77cM, 且被定位
在 Satt538标记的两边 , 解释的表型变异分别为
8.99%和7.32%, 有可能是同一个QTL。qGL-B1-1表
型贡献率为41.42%, 解释的表型变异最大。
共检测到8个大豆苷元QTL, 分别位于8(A2)、
17(D2)、20(I)和10(O)染色体上, 解释的表型变异范
围为5.09%~32.02%。2011年检测到的qDA-D2-1和
qDA-D2-2仅相距3.89 cM, 且被定位在satt082标记的
两边, 解释的表型变异分别为32.02%和30.41%, 很
可能是同一个QTL; qDA-D2-2和qDA-D2-3, 尽管位
置仅相距2.85cM, 但被定位在不相邻的2个标记区
间, 很可能不是同一个QTL。上面检测到的2011年3
个QTL和 2012年检测到的 qDA-D2-4均被定位在
17(D2)染色体上 , 但相距相对较远 , 很可能也不是
同一QTL。说明Gm17(D2)连锁群上可能存在多个与
大豆苷元相关的基因。 qDA-D2-1表型贡献率为
32.02%, 解释的表型变异最大。
共检测到5个异黄酮QTL, 分别位于8(A2)、
17(D2)、16(J_1)和10(O)染色体上, 解释的表型变异
范围为12.14%~55.20%。2011年和2012年均检测到
qIF-J_1-1且定位在16(J_1)连锁群satt249~satt285标
记区间0.01 cM处, 解释的表型变异分别为55.14%和
55.20%, 说明这是一个能够稳定遗传的QTL。2011
年检测到的qIF-D2-1和2012年检测到的qIF-D2-2均
被定位在17(D2)染色体上 , 但相距较远 , 很可能不
是同一QTL。说明17(D2)连锁群上可能存在多个与
异黄酮相关的基因。qIF-J_1-2表型贡献率为55.20%,
解释的表型变异最大。
3 讨论
3.1 遗传和环境对大豆异黄酮及其组分含量的
影响和基因定位
本研究通过对 2 个环境下大豆异黄酮及其组分
含量方差分析, 均呈现出极显著差异, 表明异黄酮
及其组分含量受遗传和环境因素的共同影响。两环
境下, 虽然同一家系大豆异黄酮及其组分含量有所
变化, 相关性分析中各环境下均表现出显著或极显
著正相关的趋势, 说明遗传是主要决定因素, 控制
这些性状的基因存在着紧密连锁, 异黄酮及其组分
含量间的相关比较稳定, 对其中一种组分选择, 即
可间接达到对其他组分的选择。该研究结果与 Chiari
等[33]和张晶莹等[26]相一致。同时, 从 RIL 群体的性
状变异范围来看, 6个性状都有超越双亲的株系出现,
说明 2个亲本都具有正效等位基因 , 研究结果与
Gutierrez-Gonzalez等[20]的报道一致。本研究中通过
对异黄酮及其组分含量广义遗传力分析, 发现各个
性状遗传力均不超过 40%, 同一家系试验材料在同
一地点, 不同年份间的表型有很大差异, 说明异黄
酮及其组分含量是易受环境影响的性状, 从而影响
到定位结果。在进行 QTL定位分析时, 可以进一步
采取 BLUE (Best Liner Unbiased Estimator, 最佳线
性无偏估计)或 BLUP (Best Liner Unbiased Predictor,
最佳线性无偏预测)等统计分析方法, 尽量减少环境
误差造成的影响。
3.2 大豆异黄酮及其组分含量 QTL 标记区间的
对比
基因多效性或者基因紧密连锁引起性状相关 ,
基因连锁或一因多效导致了不同性状间的相关性 ,
表现出同一个标记区间内检测出同时控制不同性状
的QTL[34]。本研究所定位到的42个与大豆异黄酮及
其组分含量相关的QTL中, 在8(A2)、17(D2)、18(G)、
16(J)和10(O)染色体上有7个标记在多环境中均被检
测到 , 异黄酮及其组分含量QTL存在共位性(表5),
1380 作 物 学 报 第 41卷


表 5 各年份重复检测到的与异黄酮及其组分含量相关的区间
Table 5 Stable marker interval for isoflavone contents and its components during two years
染色体
Chr.
标记区间
Marker interval
异黄酮及其组分含量
Isoflavone and its components
LOD R2 (%)
8(A2) satt187–sat_129 G(2011), GE(2011) 4.942–7.025 35.61–37.43
satt538–satt429 D(2012), GL(2012), IF(2011) 2.939–4.929 7.32–13.64
17(D2) satt154–satt082 DA(2011), IF(2011) 3.684–7.425 12.26–32.02
18(G) satt038–satt570 G(2012), GE(2011, 2012) 4.923–6.740 35.16–39.59
16(J) satt249–satt285 D(2012), IF(2011, 2012) 5.485–7.721 17.05–55.20
10(O) satt487–satt500 D(2012), G(2012), DA(2012) 2.531–10.719 5.09–16.95
satt500–satt478 G(2012), DA(2012) 2.861–7.521 5.01–29.01
D: 大豆苷; G: 染料木素; GE: 染料木苷; GL: 黄豆苷元; DA: 大豆苷元; IF: 异黄酮。
D: daidzin; G: genistein; GE: genistin; GL: glycitein; DA: daidzein; IF: isoflavone.

表明各性状QTL间存在着关联。如在大豆第16染色
体上的标记satt249与D和IF性状相关, 其中与IF相关
的QTL在2个环境中均被检测到。以往研究中, 该标
记与大豆蛋白质7S相关[35], 而大豆异黄酮与蛋白质
呈显著负相关[14], 说明该标记与大豆品质性状密切
相关 , 在大豆异黄酮与蛋白质形成与积累过程中 ,
二者存在相互制约的关系, 最后达到总体的平衡。
另外 , 在第11条染色体上 , 2个环境下均在标记
satt359定位到与大豆苷相关的QTL, 王英等 [36]检测
到该标记与大豆株高的光周期敏感度(PS)相连锁 ,
与短日处理比较, 长日照可增加株高, 降低蛋白质
含量。以往研究未见报道该标记与异黄酮主要组分
相关, 这可能是遗传背景不同所致。因此, 该标记可
能是控制异黄酮主要组分大豆苷的新QTL, 其作用
有待进一步研究证实。
本研究中检测出很多贡献率较高的QTL, 但其
定位的区间是相邻的。如2011年检测到的qG-G-1和
2012年检测到的qG-G-2均被定位在18(G)染色体上,
仅相距3.81 cM, 但被定位在相邻的标记区间, 解释
的表型变异均为35.16%, 是否是同一QTL尚不能确
定。下一步可以对两年群体数据进行联合分析, 进
一步探讨产生这种现象的可能原因。同时, 对有一
些两年试验同时检测的QTL, 特别效应值比较大的
QTL, 如qIF-J_1-1在2个环境中同时被检测到, 且贡
献率均超过55%, 希望进一步确认该基因位点，并精
细定位，为分子标记辅助育种奠定基础。
3.3 大豆异黄酮及其组分含量 QTL 定位的准确
性探讨
Janse[37]发现, 同时对多年间不同环境下的数据
进行分析, 能增大 QTL的检测强度, 提高 QTL的位
置和效应的准确程度, 更有利于挖掘稳定的 QTL。本
研究所采用的 RIL群体包含 474个家系, 群体数量相
对较大, QTL定位的精确性和稳定性相对较高。本次
试验中, 连续 2年检测到 satt430~satt359标记区间与
大豆苷相关, satt038~satt570标记区间与染料木苷相
关 , satt197~sat_128 标记区间与黄豆苷元相关 ,
satt249~satt285 标记区间与异黄酮相关 , 表明这些
QTL 有可能是稳定遗传的, 所定位的异黄酮及其组
分含量相关区间是相对准确和可信的。同时, 也可能
是因为标记区间范围相对较大, 在标记期间内存在
较多的相关基因, 各组分之间存在相对稳定的相关,
反映了各性状之间相关的遗传特性。同时, 本研究与
以往研究相比, 尽管试验地理环境条件不同, 研究群
体不同, G、GL、DA和 IF定位结果均部分吻合, 但
不完全统一(表 6和图 1)。例如, 与 Primomo 等[17]相
比, GL均被定位在 9(K)染色体、DA均被定位在 8(A2)
染色体、IF均被定位在16(J)染色体 satt249标记上; 与
Yoshikawa等[23]相比, G均被定位在 8(A2)染色体、GL
均被定位在 5(A1)、8(A2)和 11(B1)染色体; DA 均被
定位在 8(A2)和 17(D2)染色体、IF均被定位在 8(A2)
和 17(D2)染色体上; 与 Kassem等[38]相比, GL均被定
位在 11(B1)染色体 satt197 标记上等。不同年份间试
验均在同一条染色体上检测到 QTL, 如 IF 在 2011
年和 2012年均被定位在 16(J_1)染色体上, 且均被定
位在 satt249~satt285 区间的同一位点上; D 在 2011
年和 2012 年均被定位在 11(B1)连锁群标记
satt430~satt359间, 且均被定位在 124.87 cM处。与
前人研究结果相比 , IF 均被定位在 16(J)染色体
satt249标记、GL均定位在 11(B1)染色体 satt197标
记上, 说明这 2个 QTL能够稳定表达。

第 9期 梁慧珍等: 大豆异黄酮及其组分含量的遗传分析与 QTL检测 1381


表 6 QTL 定位结果的比较
Table 6 Comparison of QTLs detected in this study with in previous studies
QTL标记 Marker associated with QTL 性状
Trait
染色体
Chr. 本研究 This study 文献报道 In literatures
8(A2) satt187, satt329 sat_040[22], satt200[23], satt315[24]
11(B1) sat_129 sat_247[20]
13(F ) satt193 satt144[22], satt490[18], sat_197[18], satt395[25]
18(G) satt275, satt038 sat_372[23]
染料木素
Genistein
Gm10(O) satt487, satt500 satt479[25]
(A1) satt545 satt236[19], satt276[23]
8(A2) satt327, satt538 sat_215[23], satt315[24]
11(B1) satt197 satt197[16], sct_026[16], sat_247[20], satt509[23], satt359[23], BARC-031375-07086[24]
黄豆苷元
Glycitein
9(K ) satt264, sat_044 BARC-014813-01678[21], satt196[17], satt242[18-19], satt552[18-19], satt544[24]
8(A2) satt162 satt538[17], satt187[18], aw132402[23], sat_215[23], satt315[24]
17(D2) satt082, satt461 satt669[23]
20(I ) satt292 satt587[22]
大豆苷元
Daidzein
10(O) satt487, satt500
8(A2) satt538 BARC-022387-04319[18,20], sat_215[23], AW132402[23]
17(D2) satt154, satt301 satt208[22], satt669[23], satt386[14]
16(J) satt249 satt249[17], sat_339[19], satt249[14]
异黄酮
Isoflavone
10(O) satt259 satt241[22]

本研究采用主基因+多基因混合遗传分离分析
法对大豆异黄酮及其组分含量进行混合遗传分析表
明, 大豆苷、黄豆苷元、染料木素、染料木苷、大
豆苷元和异黄酮总含量分别受4、4、2、3、2和2对
主基因控制, 并有多基因修饰。王春娥[39]以科丰1号
与南农 1138-2及其构建的 184个重组自交系群体
(NJRIKY)籽粒为材料, 对12种大豆异黄酮组分含量
及其总含量进行了分离分析, 认为大豆籽粒异黄酮
总含量和各组分含量的遗传均属2~3对主基因加多
基因混合遗传模型 , 主基因遗传率为 16.42%~
45.08%, 多基因遗传率为54.78%~83.44%。对比研究
结果 , 多数结论吻合 , 少数结论不同 , 有待下一步
深入研究。但是, 本研究中也出现了 QTL定位数目
与遗传模型预测的数目不完全一致的现象。如大豆
苷元最适遗传模型符合 MX2-Inhibiting-A, 2次试验
均定位出4个主效 QTL, QTL定位分析结果与遗传分
离分析结果预测值不一致。导致这种情况的可能原
因有: (1)本研究选用的模型分离分析方法是4对主基
因加多基因混合遗传模型[40], 该模型建立时只考虑
了主基因间无连锁的情况, 对于主基因间有连锁的
情况没有分析; (2)该遗传模型建立过程中, 成分分
布数和待估参数过多, 分析过程受到限制, 需要进
一步改进参数估算方法 , 提高估算结果的准确性 ;
(3)与 QTL定位时选用的图谱质量有关。图谱本身存
在不完善的地方。尽管研究结果还存在着不足, 仍
然为大豆异黄酮及其组分含量的分子研究和遗传分
析提供了有益参考。
4 结论
大豆苷和黄豆苷元主要受4对主基因控制 ; 染
料木素受2对显性上位主基因-加性多基因混合作用
控制; 染料木苷主要受到3对主基因控制; 大豆苷元
受2对抑制作用主基因-加性多基因控制; 异黄酮总
含量受2对加性-上位性主基因×加性多基因控制。检
测到44个与大豆异黄酮及其组分含量相关的 QTL,
与大豆苷、染料木素、黄豆苷元、大豆苷元、染料
木苷和异黄酮总含量相关的 QTL分别有10、9、4、
7、8和6个。连续2年检测到 satt430~satt359标记区间
与大豆苷相关, satt038~satt570标记区间与染料木苷
相关 , satt197~sat_128标记区间与黄豆苷元相关 ,
satt249~satt285标记区间与异黄酮相关。经比对文
献，认为位于 Gm16(J)上的 satt249和位于 Gm11(B1)
上的 satt197附近存在稳定遗传的 QTL，分别控制异
黄酮和黄豆苷元含量。
1382 作 物 学 报 第 41卷


References
[1] Patel R P, Boersma B J, Crawford J H, Hogg N, Kirk M, Kaly-
anaraman B, Paeks D A, Barnes S, Darley-Usmar V. Antioxidant
mechanisms of isoflavones in lipid systems: paradoxical effects
of peroxyl radical scavenging. Free Radical Biol Med, 2001, 31:
1570–1581
[2] Arean S, Rappa C, Del Frate E, Cenci S, Villani C. A natural al-
ternative to menopausal hormone replacement therapy. Phytoes-
trogens Minerva Ginecol, 2002, 54: 53–57
[3] Kritz-Silverstein D, Goodman-Gruen D L. Usual dietary isofla-
vone intake, bone mineral density, and bone metabolism in post-
menopausal women. J Women’s Health Gend Based Med, 2002,
11: 69–78
[4] Atkinson C, Compston J E, Day N E, Owsett M, Bingham S A.
The effects of phytoestrogen isoflavones on bone density in
women: a double-blind, randomized, lacebo-controlled trial. Am J
Clinic Nutr, 2004, 79: 326–333
[5] Lamartiniere C A, Cotroneo M S, Fritz W A, Wang J, Mentor-
Marcel R, Elgavish A. Genistein chemoprevention: timing and
mechanisms of action in murine mammary and prostate. J Nutr,
2002, 132: 552–558
[6] Liggins J, Mulligan A, Runswick S, Bingham S A. Daidzein and
genistein content of cereals. Clinic Nutr, 2002, 56: 961–966
[7] Munro I C, Harwood M, Hlywka J J, Stephen A M, Doull J,
Flamm W G, Adlercreutz H. Soy isoflavones: a safety review.
Nutr Rev, 2003, 61: 1–33
[8] Ma D F, Qin L Q, Wang P Y, Katoh R. Soy isoflavone intake in-
creases bone mineral density in the spine of menopausal women:
Metaanalysis of randomized controlled trials. Clin Nutr, 2008, 27:
57–64
[9] Nagata C. Factors to consider in the association between soy
isoflavone intake and breat cancer risk. J Epidemiol, 2010, 20:
83-89
[10] Kudou S, Fleury Y, Welti D, Magnolato D, Uchida T, Kitamura K,
Okubo K. Malonyl isoflavone glycosides in soybean seeds (Gly-
cine max L. Merr.). Agric Biol Chem, 1991, 55: 2227–2233
[11] Li W D, Liang H Z, Lu W G, Wang S F, Yang Q C, Yang C Y, Liu
Y F. Studies on the correlations between isoflavone contents in
soybean seed and the eco-physiological factors. Agric Sci China,
2004, 3: 340–348
[12] Liang H Z, Li W D, Fang X J, Cao Y N, Wang H. Genetic analy-
sis of embryo, cytoplasm and maternal effects and their environ-
ment interactions for isoflavone content in soybean [Glycine max
(L.) Merr.]. Agric Sci China, 2007, 6: 1051–1059
[13] Liang H Z, Li W D, Cao Y N, Wang H. Genetic analysis of com-
bining abilities and heterosis for the contents of soybean isofla-
vone and its components among the soybean varieties [Glycine
max (L.) Merr.]. Agric Sci China, 2005, 4: 101–105
[14] Liang H Z, Yu Y L, Wang S F, Lian Y, Wang T F, Wei Y L, Gong
P T, Liu X Y, Fang X J, Zhang M C. QTL mapping of isoflavone,
oil and protein contents in soybean (Glycine max L. Merr.). Agric
Sci China, 2010, 9: 1108–1116
[15] Meksem K, Njiti V N, Banz W J, Iqbal M J, Kassem M M, Hyten
D L, Yuang J, Winters T A, Lightfoot D A. Genomic regions that
underlie soybean seed isoflavone content. J Biomed Biotechnol,
2001, 1: 38–44
[16] Kassem M A, Shultz J, Meksem K, Cho Y, Wood A J, Iqbal M J,
Lightfoot D A. An updated ‘Essex’ by ‘Forrest’ linkage map and
first composite interval map of QTL underlying six soybean traits.
Theor Appl Genet, 2006, 113: 1015–1026
[17] Primomo V S, Poysa V, Ablett G R, Jackson C J, Gijzen M, Ra-
jcan I. Mapping QTL for individual and total isoflavone content
in soybean seeds. Crop Sci, 2005, 45: 2454–2464
[18] Gutierrez-Gonzale J J, Wu X L, Zhang J, Lee J D, Ellersieck M,
Shannon J G, Yu O, Nguyen T H, Sleper A D. Genetic control of
soybean seed isoflavone content: Importance of statistical model
and epistasis in complex traits. Theor Appl Genet, 2009, 119:
1069–1083
[19] Gutierrez-Gonzale J J, Wu X, Gillman J, Lee J, Zhong R, Yu O,
Shannon G, Ellersieck M, Nguyen H, Sleper D. Intricate envi-
ronment-modulated genetic networks control isoflavone accumu-
lation in soybean seeds. Plant Biol, 2010, 10: 105–120
[20] Gutierrez-Gonzale J J, Vuong D T, Zhong R, Yu O, Lee J, Shan-
non G, Ellersieck M, Nguyen T H, Sleper A D. Major locus and
other novel additive and epistatic loci involved in modulation of
isoflavone concentration in soybean seeds. Theor Appl Genet,
2011, 123: 1375–1385
[21] Smallwood C. Detection of quantitative trait loci for marker- as-
sisted selection of soybean isoflavone genistein. TN Res. Creat.
Exch. 2012
[22] Zeng G, Li D, Han Y, Teng W, Wang J, Qiu L, Li W. Identifica-
tion of QTL underlying isoflavone contents in soybean seeds
among multiple environments. Theor Appl Genet, 2009, 118:
1455–1463
[23] Yoshikawa T, Okumoto Y, Ogata D, Sayama T, Teraishi M, Terai
M, Toda T, Yamada K, Yagasaki K, Yamada N, Tsukiyama T,
Yamada T, Tanisaka T. Transgressive segregation of isoflavone
contents under the control of four QTLs in a cross between dis-
tantly related soybean varieties. Breed Sci, 2010, 60: 243–254
[24] Yang K, Moon J, Jeong N, Chun H, Kang S, Back K, Jeong S.
Novel major quantitative trait loci regulating the content of
isoflavone in soybean seeds. Genes Genom, 2011, 33: 685–692
[25] Zhang H J, Li J W, Liu Y J, Jiang W Z, Du X L, Li L, Li X W, Su
L T, Wang Q Y, Wang Y. Quantitative trait loci analysis of indi-
vidual and total isoflavone contents in soybean seeds. J Genet,
2014, 93: 331–338
[26] 张晶莹, 葛一楠, 孙君明, 韩粉霞, 于福宽, 闫淑荣, 杨华. 多
环境条件下大豆异黄酮主要组分的 QTL定位. 中国农业科学,
2012, 45: 3909–3920
Zhang J Y, Ge Y N, Sun J M, Han F X, Yu F K, Yan S R, Yang H.
Identification of QTLs for major isoflavone components among
multiple environments in soybean seeds. Sci Agric Sin, 2012, 45:
3909–3920 (in Chinese with English abstract)
[27] Wang S C, Basten C J, Zeng Z B. Windows QTL Cartographer
2.5 User Manual. Department of Statistics, North Carolina State
University, Raleigh, NC, 2005
[28] 曹锡文, 刘兵, 章元明. 植物数量性状分离分析 Windows 软
件包 SEA的研制. 南京农业大学学报, 2013, 36(6): 1–6
Cao X W, Liu B, Zhang Y M. SEA: a software package of segre-
gation analysis of quantitative traits in plants. J Nanjing Agric
Univ, 2013, 36(6): 1–6 (in Chinese with English abstract)
[29] Sun J M, Sun B L, Han F X, Yan S R, Yang H, Akio K. Rapid
第 9期 梁慧珍等: 大豆异黄酮及其组分含量的遗传分析与 QTL检测 1383


HPLC method for determination of 12 isoflavone components in
soybean seeds. Agric Sci China, 2011, 10: 70–77
[30] 王珍. 大豆 SSR 遗传图谱构建及重要农艺性状 QTL分析. 广
西大学硕士学位论文, 广西南宁, 2004
Wang Z. Construction of Soybean SSR Based Map and QTL
Analysis of Important Agronomic Traits. MS Thesis of Guangxi
University, Nanning, China, 2004 (in Chinese with English ab-
stract)
[31] 梁慧珍. 大豆子粒性状的遗传及 QTL 分析. 西北农林科技大
学博士学位论文, 陕西杨凌, 2006
Liang H Z. Genetic Analysis and QTL Mapping of Seed Traits in
Soybean [Glycine max (L.) Merr]. PhD Dissertation of Northwest
A＆F University. Yangling, China, 2006 (in Chinese with English
abstract)
[32] McCouch S R, Cho Y G, Yano M, Paul E, Blinstrub M, Mor-
ishima H, Kinoshita T. Report on QTL nomenclature. Rice Genet
Newsl, 1997, 14: 11–14
[33] Chiari L, Naoe L K, Piovesan N D, José I C, Cruz C D, Moreira
M A, Barros E G. Genetic parameters relating isoflavone and
protein content in soybean seeds. Euphytica, 2004, 138: 55–60.
[34] Hagiwara W E, Onish K, Takamure I, Sano Y. Transgressive seg-
regation due to linked QTLs for grain characteristics of rice.
Euphytica, 2006, 150: 27–35
[35] Panthee D R, Kwanyuen P, Sams C E, West D R, Saxton A M,
Pantalone V R. Quantitative trait loci for β-conglycinin (7S) and
glycinin (11S) fractions of soybean storage protein. J Am Oil
Chem Soc, 2004, 81: 1558–9331
[36] 王英, 程立锐, 冷建田, 吴存祥, 侯文胜, 韩天富. 开花后不
同光周期条件下大豆农艺性状和品质性状的 QTL 分析. 作物
学报, 2010, 36: 1092–1099
Wang Y, Cheng L R, Leng J T, Wu C X, Hou W S, Han T F. QTL
mapping of agronomic and quality traits in soybean under differ-
ent post-flowering photoperiods. Acta Agron Sin, 2010, 36:
1092–1099 (in Chinese with English abstract)
[37] Jansen R C, Van Ooijien J M, Stam P, Lister C, Dean C. Geno-
type-by-environment interaction in genetic mapping of multiple
quantitative trait loci. Theor Appl Genet, 1995, 91: 33–37
[38] Kassem A, Meksem K, Iqbal M, Njiti V, Banz W, Winters T,
Wood A, Lightfoot D. Definition of soybean genomic regions that
control seed phytoestrogen amounts. J Biomed Biotechnol, 2004,
1: 52–60
[39] 王春娥. 我国大豆资源豆腐豆乳得率和异黄酮含量的遗传变
异及两类性状的 QTL 分析. 南京农业大学博士学位论文, 江
苏南京, 2008
Wang C E. Genetic Variability of Toftu and Soymilk Output and
Isoflavone Content in Soybean Germplasm from China and Qtl
Mapping of the Two Kinds of Traits. PhD Dissertation of Nanjing
Agricultural University, Nanjing, China, 2008 (in Chinese with
English abstract)
[40] 王金社, 李海旺, 赵团结, 盖钧镒. 重组自交家系群体 4 对主
基因加多基因混合遗传模型分离分析方法的建立. 作物学报,
2010, 36: 191–201
Wang J S, Li H W, Zhao T J, Gai J Y. Establishment of segrega-
tion analysis of mixed inheritance model with four major genes
plus polygenes in recombinant inbred lines population. Acta
Agron Sin, 2010, 36: 191–201 (in Chinese with English abstract)



欢迎订阅 2016 年《麦类作物学报》
《麦类作物学报》是由西北农林科技大学、中国作物学会、国家小麦工程技术研究中心联合主办(教育部主
管)的专业性学术期刊, 也是全国唯一的一份麦类作物专刊。主要刊载麦类作物(小麦、大麦、燕麦、黑麦等)遗
传育种、生理生化、栽培管理、食品加工、产品贸易等方面有创见性的学术论文、领先水平的科研成果、学术
报告、有新意的文献综述以及学术动态等。读者对象为国内外农业科技人员、农业院校师生及高级农业技术推
广和管理人员。
《麦类作物学报》为“农业科学中文核心期刊”、“中国科技核心期刊”、“中国科技精品期刊”, 现已被英国
《国际农业与生物技术文摘》数据库(CABI)、美国《化学文摘》数据库(CA)、美国《剑桥科学文摘》数据库(CSA)、
俄罗斯《文摘杂志》(AJ)数据库、日本《科学技术》数据库(JST)、波兰《哥白尼索引》数据库(IC)、《中国科学
引文数据库》(核心库)等国内外多家权威性检索系统收录。影响因子排名已连续 3年居全国农业期刊前 10位。
《麦类作物学报》为月刊, 每月中旬出版, 定价 30.00元/册, 全年 360元, 国内刊号: CN 61-1359/S, 国际刊号:
1009-1041。全国各地邮局均可订阅, 邮发代号: 52-66。漏订者可直接汇款至编辑部补订。国外总发行: 北京 中
国国际图书贸易总公司, 代号: 1479B。
热忱欢迎国内外专家随时指导和赐稿, 亦欢迎各有关课题组、单位和个人出版专辑、刊登广告。
联系人: 华千勇 电话: 029-87082642 (兼传真)
通讯地址: 陕西杨凌邰城路 3号《麦类作物学报》编辑部, 邮编: 712100
网址: http://mlzwxb.alljournals.ac.cn/; http://www.tcrop.net/(两个网址是同一个网站。有在线投稿、在线
查询稿件等在线办公系统)