全 文 :作物学报 ACTA AGRONOMICA SINICA 2015, 41(7): 10271038 http://zwxb.chinacrops.org/
ISSN 0496-3490; CODEN TSHPA9 E-mail: xbzw@chinajournal.net.cn
本研究由国家自然科学基金项目(31371655)资助。
* 通讯作者(Corres ponding author): 刘列钊, E-mail: liezhao2003@126.com. Tel: 023-68250701
第一作者联系方式: E-mail: wangjia0724@126.com **同等贡献(Contributed equally to this work)
Received(收稿日期): 2014-11-25; Accepted(接受日期): 2015-04-02; Published online(网络出版日期): 2015-04-14.
URL: http://www.cnki.net/kcms/detail/11.1809.S.20150414.1633.003.html
DOI: 10.3724/SP.J.1006.2015.01027
甘蓝型油菜株高、第一分枝高和分枝数的 QTL检测及候选基因筛选
王 嘉** 荆凌云** 荐红举 曲存民 谌 利 李加纳 刘列钊*
西南大学农学与生物科技学院, 重庆 400715
摘 要: 株高、分枝数及第 1分枝高是油菜重要的农艺性状。本研究利用甘蓝型油菜 GH06和 P174杂交, F2通过单
粒法连续自交至 F11构建重组自交系群体, 利用油菜 60K芯片对该群体进行基因分型, 构建高密度遗传连锁图谱。结
果表明, 该图谱包含 2795个 SNP多态性标记位点, 总长 1832.9 cM, 相邻标记间平均距离为 0.66 cM。在此图谱基础
上采用复合区间作图法(CIM), 检测到 3个农艺性状的 24个 QTL。其中 11个株高 QTL分别位于 A01、A06、A07、
A08、A10和 C06染色体, 单个 QTL解释 5.00%~15.26%的表型变异; 7个第 1分枝高 QTL分别位于 A06、C05和 C06
染色体, 单个 QTL解释 5.04%~12.99%的表型变异; 6个分枝数 QTL分别位于 A03、A07、C01、C04和 C06染色体, 单
个 QTL 解释 5.95%~8.14%的表型变异。将 156 个拟南芥株高相关基因、10 个拟南芥第 1 分枝高相关基因和 148 个
拟南芥分枝数相关基因与 QTL 对应置信区间序列进行同源比较分析(E<1E–20), 分别找出了 20个株高候选基因、3
个第 1 分枝高候选基因以及 12 个分枝数候选基因。2 个环境中在 A07 染色体上重复检测到的 QTL 置信区间检测到
与株高相关的候选基因 ATGID1B/GID1B 和 WRI1, A08 染色体上重复检测到的 QTL 置信区间检测到 SLR/IAA14 和
AXR2/IAA72 个与株高相关的候选基因。在具有部分置信区间重叠的 q2013FBH-C05-1 和 q2014FBH-C05-2 区间均检
测到第 1分枝高候选基因 PHT1;8, 在 A03和 C06染色体上的 QTL置信区间内, 分别检测到 4个分枝数候选基因, 匹
配 E值介于 0~3E–56之间。
关键词: 甘蓝型油菜; 单核苷酸多态性; 数量性状位点; 株高; 第 1分枝高度; 分枝数
Quantitative Trait Loci Mapping for Plant Height, the First Branch Height,
and Branch Number and Possible Candidate Genes Screening in Brassica
napus L.
WANG Jia**, JING Ling-Yun**, JIAN Hong-Ju, QU Cun-Min, CHEN Li, LI Jia-Na, and LIU Lie-Zhao*
College of Agronomy and Biotechnology, Southwest University, Chongqing 400715, China
Abstract: Plant height, the first branch height and branch number are important agronomic traits in rapeseed. In our study, QTL
mapping of plant height, the first branch height and branch number in Brassica napus was tested by using the high density SNP
genetic map constructed from the high generation RIL population with the rapeseed 60K chip array. The reference SNP genetic
map contains 2795 SNP markers, covering 1832.9 cM with an average distance of 0.66 cM in B. napus genome. Totally, 24 puta-
tive QTLs were identified for plant height, the first branch height and branch number by using the complex interval mapping. The
phenotypic variation was explained by individual QTL ranged from 5.00% to 15.26% for 11 QTLs of plant height, from 5.04% to
12.99% for seven QTLs of the first branch height, and from 5.95% to 8.14% for six QTLs of branch number. We collected 156
genes associated with plant height, 10 genes associated with the first branch height and 148 genes associated with branch number
in Arabidopsis thaliana and searched the homology region of the QTL confidence interval E-value<1E–20 to screen the possible
candidate genes. We found 20 genes associated with plant height, three genes associated with the first branch height and 12 genes
associated with branch number. The plant height candidate genes ATGID1B/GID1B and WRI1 were found in QTL confidence
interval on A07 and SLR/IAA14 and AXR2/IAA7 were found in QTL confidence interval on A08 in 2013 and 2014. The first
branch height candidate gene PHT1;8 was found in QTLs q2013FBH-C05-1 and q2014FBH-C05-2 that had an overlapping con-
1028 作 物 学 报 第 41卷
fidence interval. Moreover, we found four genes associated with branch number that E-value from 0 to 3E–56 on A03 and C06,
respectively.
Keywords: Brassica napus; SNP; QTL; Plant height; The first branch height; Branch number
油菜(Brassica napus L.)属于十字花科芸薹属,
为世界四大油料作物之一, 菜籽是仅次于大豆的第
二大植物油来源[1]。影响油菜产量的农艺性状主要
有株高、有效分枝数、分枝部位以及开花时间等。
由于油菜农艺性状大多为数量性状, 由微效多基因
控制 , 表现连续变异 , 受环境的影响较大 , 不能依
照质量性状的处理方法将单个基因的效应区分开 ,
因此单独依靠传统育种方法和技术在现有基础上研
究尚未有大的突破[2]。
DNA 分子标记和 QTL 作图等生物技术的发展
为复杂的数量性状的研究提供便利。Li等 [3]利用
SRAP、SSR、AFLP 标记以及功能标记对甘蓝型油
菜 F2群体在 2个环境下同时进行了 QTL定位, 检测
到 12个与一次有效分枝数相关的 QTL, 单个 QTL
能够解释 5.47%~27.39%的表型变异。Chen等[4]利用
甘蓝型油菜 2 个群体(DH 和 IF2)定位了 6 个分枝数
QTL, 其中位于 13连锁群上的 42.0~55.6 cM的 QTL
有较大的效应。易斌等[2]利用 RFLP、RAPD、SSR
以及 SRAP 等分子标记技术共检测到与油菜产量及
其相关性状有关的 QTL 17个。其中与株高相关的 3
个分别位于第 4、第 9 和第 10 连锁群上, 对性状的
解释率为 9.42%~17.58%; 与分枝部位有关的 4个分
别位于第 4、第 6和第 7连锁群上, 其中 Bp1和 Bp2
均位于第 4 连锁群 , 对性状的解释率为 8.13%~
15.20%; 与一次有效分枝有关的 2个分别位于第 1、
第 4连锁群上, 对性状的解释率为 15.29%~19.58%。
张倩[5]利用 SSR标记检测到 4个与株高相关的 QTL,
分别位于 LG1、LG2 连锁群上, 可解释表型变异的
8.54%~17.04%; 与分枝高度相关的 5 个, 分别位于
LG3、LG5、LG8 连锁群上 , 可解释表型变异的
3.92%~11.82%; 与分枝数相关的 1 个, 位于 LG3 连
锁群上, 可解释表型变异的 11.08%。Zhang等[6]利用
复合区间作图法, 在第 4和第 21连锁群上共检测到
3个与株高相关的QTL, 每个QTL分别解释 28.60%、
23.85%和 11.11%的表型变异; 在第 11、第 13、第
20 和第 21 连锁群上共检测到 4 个与分枝高度相
关的 QTL, 单个 QTL能解释 8.54%~54.59%的表型
变异。
SNP (single nucleotide polymorphism)是基于基
因组核苷酸水平上单碱基插入、缺失、转换和颠换
引起的多态性差异而开发的 DNA标记[7]。随着甘蓝
型油菜全基因组测序的完成[8]以及高通量 SNP 芯片
技术的完善和成本的下降, SNP 芯片技术将成为油
菜遗传研究的热点之一。本研究利用油菜 60K芯片
构建的高密度 SNP 遗传连锁图谱, 对株高、第 1 分
枝高和分枝数进行 QTL定位分析, 利用甘蓝型油菜
基因组数据库序列, 根据 QTL区间物理位置及拟南
芥基因组数据库信息筛选可能的候选基因。尽可能
将数量性状定位结果明确在特定的物理位置甚至候
选基因, 有利于不同材料之间定位结果的比较分析
及候选基因的克隆分析。
1 材料与方法
1.1 试验材料
以甘蓝型黄籽油菜 GH06 作母本, 甘蓝型黑籽
型双低油菜 P174作父本配制杂交组合, F2通过单粒
法再连续自交至 F11构建的重组自交系群体, 随机选
取其中的 172 个材料进行 SNP 标记分析, 构建高密
度 SNP遗传连锁图谱。所有材料均由重庆市油菜工
程技术研究中心提供。
1.2 田间试验
2012 年 9 月至 2013 年 5 月(性状统计记录为
2013)与 2013年 9月至 2014年 5月(性状统计记录为
2014)将亲本及重组自交系群体种植于重庆市北碚
区歇马镇油菜基地。随机区组设计 3个重复, 每个小
区 2行, 每行 10株。行距 0.40 m, 株距 0.24 m。田
间管理同常规生产, 确保所有样本的外部生长环境
一致, 待成熟后进行数据统计。
1.3 性状考察
油菜成熟后, 每个株系随机选取 5株测量株高、
第 1分枝高(有效分枝)、分枝数, 取平均值。株高为
从地面到植株最高部位的距离(cm); 第 1 分枝高为
地面到第 1 有效分枝的高度(cm); 分枝数为该植株
所有有效的一次分枝数。
1.4 遗传连锁图谱构建及 QTL分析
本高密度SNP遗传图谱包含2795个SNP多态性
标记位点, 图谱总长1832.9 cM, 相邻标记间平均距
离为0.66 cM [9]。采用QTL分析软件Windows QTL
Cartographer 2.5[10]及复合区间作图(composite inter-
val mapping, CIM)法对该群体的株高、第1分枝高、
第 7期 王 嘉等: 甘蓝型油菜株高、第一分枝高和分枝数的 QTL检测及候选基因筛选 1029
分枝数进行QTL定位及效应检测[11]。CIM分析时, 选
取1 cM的步长(walking speed), 按照假定检测10和
Zmapqtl模型6, 选取参数为1000次回归 , 显著水平
为0.01。LOD≥3.0时, 即认为该区间可能存在一个
QTL。运行软件后可同时给出各QTL的加性效应和
解释的表型变异。按照McCouch等[12]的方法对检测
到的QTL命名, 以“q”加上性状再加染色体编号表示,
字体为斜体。同一性状在染色体相同的位置检测重
复的QTL, 且加性效应方向一致, 认为是同一QTL。
1.5 候选基因筛选
为了筛选出与株高、第 1 分枝高及分枝数相关
的候选基因, 以检测到的 QTL置信区间在甘蓝型油
菜基因组[8]上查询对应的序列, 然后与 Shi 等[13]和
Ding等[14]搜索出的 156个拟南芥株高相关基因、148
个分枝数相关基因, Shi等[15]检测到的 10个第 1分枝
高相关基因进行 BlastN, 将 E 值设定为 E–20, 最后
筛选出每个 QTL置信区间内匹配 E值小于阈值的株
高、第 1分枝高及分枝数基因。
2 结果与分析
2.1 不同年份的表型分析
如表 1和图 1所示, RIL群体中株高、第 1分枝
高和分枝数在 2 年的均值差异都较小, 表明其是多
位点控制的典型数量性状, 且位点间加性效应在同
一亲本中具有不同的方向。各性状在 RIL 群体中均
呈连续分布和双向超亲分离, 表明这些性状均为多
基因控制的数量性状, 其偏斜度和峰值均小 1, 适合
于 QTL分析。
相关分析表明(表2), 同一性状在2013年和2014
年表现均呈极显著正相关, 且相关系数较大, 说明
株高、第1分枝高和分枝数性状遗传稳定。性状间的
相关分析还表明, 在2个生长周期里, 第1分枝高和
株高呈极显著正相关 , 相关系数分别为 0.468和
0.451; 分枝数与株高呈显著负相关或未达到显著水
平。同时第1分枝高和分枝数在2年里均呈极显著负
相关, 相关系数分别为–0.590和–0.475。
2.2 株高、第 1 分枝高及分枝数在 2 年的 QTL
分析
利用软件Windows QTL Cartographer 2.5对 2年
的株高、第 1分枝高、分枝数进行 QTL分析, 共检
测到 24个QTL (表 3和图 2)。其中 11个株高QTL, 阈
值 3 . 0 2 ~ 7 . 5 1 , 单个 Q T L 解释的表型变异为
5.00%~15.26%。在 2 年内检测到对株高有显著效应
的共同 QTL位点(q2013PH-A07和 q2014PH-A07-1、
表 1 亲本及重组自交系群体 3个农艺性状在 2年的表型分析
Table 1 Analysis of three agronomic traits in the two parents and RIL population in 2013 and 2014
亲本 Parent RIL群体 RIL population
年份和性状
Year and trait GH06
(cm)
P174
(cm)
范围
Range (cm)
均值
Mean (cm)
方差
Variance
标准差
SD
变异系数
CV
偏度
Skewness
峰度
Kurtosis
2013PH 233.5 218.5 178–271 226.91 412.12 20.32 0.09 –0.08 –0.56
2014PH 236.0 215.2 164–266 227.79 320.04 17.89 0.08 –0.12 –0.38
2013FBH 53.8 87.4 25–174 84.56 880.67 29.67 0.35 –0.27 0.16
2014FBH 70.5 92.1 26–152 84.03 483.61 21.99 0.26 –0.05 –0.06
2013BN 10.0 9.0 4–13 8.28 3.33 0.14 0.22 0.42 –0.07
2014BN 9.0 8.3 3–11 8.34 1.84 1.36 0.16 0.39 –0.41
PH: 株高; FBH: 第 1分枝高; BN: 分枝数。PH: plant height; FBH: the first branch height; BN: branch number.
表 2 甘蓝型油菜重组自交系群体农艺性状在 2013年和 2014年的相关性分析
Table 2 Correlation analysis for agronomic traits from B. napus RIL population in 2013 and 2014
2013PH 2014PH 2013FBH 2014FBH 2013BN
2014PH 0.878**
2013FBH 0.468** 0.447**
2014FBH 0.440** 0.451** 0.843**
2013BN –0.180* –0.107 –0.590** –0.456**
2014BN –0.066 –0.010 –0.403** –0.475** 0.810**
PH: 株高; FBH: 第 1分枝高; BN: 分枝数。**和*分别代表在 0.01和 0.05显著水平。
PH: plant height; FBH: the first branch height; BN: branch number. **and* denote significant correlation at the 0.01 and 0.05 probability
levels, respectively.
1030 作 物 学 报 第 41卷
图 1 甘蓝型油菜重组自交系群体株高、第 1分枝高和分枝数在
2013年和 2014年的频率分布
Fig. 1 Frequency distribution of plant height, the first branch
height, and branch number from B. napus RIL population in
2013 and 2014
q2013PH-A08和q2014PH-A08), 分别位于A07和A08
染色体上, 加性效应值均为负, 说明增效基因来自
父本P174, 在C06染色体上两年内也检测到一个具
有重叠置信区间的QTL (q2013PH-C06和q2014PH-
C06), 加性效应值大于0, 说明增效基因来自母本
GH06。对于第1分枝高在2年内检测到7个QTL, 阈值
为3.54~5.89, 单个QTL解释的表型变异为5.04%~
12.99%, 分别位于C05、C06和A06染色体上。对于分
枝数在2年内共检测到6个QTL, 阈值为3.20~4.54, 单
个QTL解释的表型变异为5.95%~8.14%。
2.3 候选基因筛选
本研究将 24 个 QTL 置信区间序列与 156 个株
高基因、10个第 1分枝高基因和 148个分枝数基因
分别进行比对, 其中有 11个与株高相关的QTL共检
测到 20个候选基因, 匹配 E值介于 0~9E–24之间, 2
个与第 1分枝高相关的QTL共检测到 3个候选基因,
5 个与分株数相关的 QTL 共检测到 12 个候选基因,
匹配 E值介于 0~3E–56之间。其中在 q2013PH-A06
和 q2013PH-A10置信区间内分别筛出 3个株高候选
基因, q2013PH-A07、q2013PH-A08、q2014PH-A01、
q2014PH-A07-1、q2014PH-A08 以及 q2014PH-C06
置信区间内分别筛出 2个株高基因。株高基因 IAA14/
SLR 和 IAA7/AXR2 在 2 年重复检测到的 QTL
(q2013PH-A08和 q2014PH-A08)区间内均检测到, 在
2年内具有重叠置信区间的 QTL q2014PH- A07-1和
q2013PH-A07 区间内均检测到株高基因 ATGID1B/
GID1B和 WRI1。第 1分枝高候选基因 PHT1;8在具
有部分置信区间重叠的 QTL q2013FBH-C05-1 和
q2014FBH-C05-2 中均检测到 , q2013BN-A03 和
q2014BN-C06 区间内分别检测到 4 个分枝数候选基
因。比对结果列于表 4, 候选基因的功能列于表 5。
3 讨论
油菜的农艺性状主要包括株高、有效分枝数、
第1分枝高度、千粒重、开花时间、主花序长、有效
角果总数、果长、每果粒数等。株高作为一个重要
的性状受到了育种家的重视, 这不仅因为株高在油
菜产量形成过程中具有重要的间接支配作用, 而且
株高直接影响能否有效机械收割。分枝不仅是油菜
角果的重要载体, 还影响其形态建成, 对油菜产量
具有十分重要的作用。一次有效分枝数作为角果数
的重要载体, 对油菜产量起着重要的作用, 在一定
范围内增加油菜分枝数可以显著增加全株角果数 ,
进而提高单株产量[44-45]。对于油菜农艺性状与产量的
关系分析, 已有大量报道, 虽然结果不尽一致, 但普
遍认同单株角果数与分枝数对产量的贡献最大[46-49]。
对油菜株高、第 1 分枝高以及分枝数的遗传分
析大多停留在经典数量遗传学研究方面, 对油菜分
枝数等农艺性状的遗传机制了解较少[50]。相对于其
他作物, 甘蓝型油菜农艺性状 QTL定位研究开展比
较晚, 而且到目前为止, 只报道了为数不多的农艺
性状QTL定位研究和少量的QTL [51-55], 而且定位结
第 7期 王 嘉等: 甘蓝型油菜株高、第一分枝高和分枝数的 QTL检测及候选基因筛选 1031
表 3 利用复合区间作图法检测到的株高(PH)、第 1分枝高(FBH)和分枝数(BN)在 2个环境中的 QTL
Table 3 Putative QTL detected by composite interval mapping for plant height (PH), the first branch height (FBH), and branch
number (BN) in two environments
性状
Trait
数量性状座位
QTL
染色体
Chromosome
标记区间
Position
置信区间
Confidence interval
LOD
加性效应
Additive
贡献率
R2 (%)
PH q2013PH-A06 A06 22.01 SNP14242A06–SNP14309A06 3.68 5.10 5.86
q2013PH-A07 A07 61.51 SNP15557A07–SNP15554A07 5.49 –6.77 10.47
q2013PH-A08 A08 75.21 SNP17561A08–SNP19318A08 5.61 –6.03 8.53
q2013PH-A10 A10 40.61 SNP22138A10–SNP22091A10 3.09 5.16 5.21
q2013PH-C06 C06 47.21 SNP41759C06–SNP36796 9.76 8.13 15.56
q2014PH-A01 A01 56.11 SNP2507A01–SNP2403A01 3.02 –4.27 5.18
q2014PH-A06 A06 58.71 SNP14986A06–SNP12820A06 3.14 4.30 5.27
q2014PH-A07-2 A07 71.31 SNP15747A07–SNP15778A07 3.11 –4.14 5.00
q2014PH-A07-1 A07 61.51 SNP15557A07–SNP15554A07 6.12 –6.62 12.80
q2014PH-A08 A08 75.21 SNP17561A08–SNP19318A08 4.86 –4.83 7.05
q2014PH-C06 C06 48.91 SNP36829–SNP36812 7.51 5.87 10.15
FBH q2013FBH-C05-2 C05 72.91 SNP37857–SNP42376C05 4.66 7.35 5.66
q2013FBH-C05-1 C05 62.71 SNP28959C05–SNP44638C05 3.55 8.29 6.92
q2013FBH-C06-2 C06 56.81 SNP24413C06–SNP24466C06 3.54 8.73 7.58
q2013FBH-C06-1 C06 49.21 SNP36829–SNP36827 5.89 11.41 12.99
q2014FBH-A06 A06 8.91 SNP13071A06–SNP13385A06 4.00 5.19 5.04
q2014FBH-C05-2 C05 66.11 SNP47509C05–SNP44638C05 4.57 5.84 6.31
q2014FBH-C05-1 C05 58.71 SNP39974C05–SNP28197 3.69 6.36 7.45
BN q2013BN-A03 A03 22.41 SNP6539A03–SNP6422A03 4.54 –0.54 7.67
q2013BN-C06 C06 36.91 SNP29430C06–SNP29760C06 3.59 –0.47 6.04
q2014BN-A07 A07 24.31 SNP7695A03–SNP17501A07 3.48 –0.35 6.13
q2014BN-C01 C01 19.01 SNP46967C01–SNP34335C01 4.09 0.46 8.41
q2014BN-C04 C04 59.21 SNP34444–SNP26297 3.59 0.35 6.37
q2014BN-C06 C06 28.11 SNP29459–SNP29410C06 3.20 –0.33 5.95
果不尽一致, 这也从侧面表明油菜农艺性状的复杂
性。因此, 研究株高、第1分枝高及分枝数具有重要
的现实意义。利用高密度SNP遗传图谱已经对甘蓝
型油菜的种皮纤维素[9]、种子硫苷含量[56]及盐、旱
胁迫下种子发芽率[57]进行了准确定位。本研究利用
高密度SNP遗传图谱对甘蓝型油菜株高、第1分枝高
及分枝数进行QTL分析, 均在同一标记区间或紧邻
的染色体区域内重复检测到一个或多个QTL, 但其
效应值均较小, 没有主效基因。
Wei等[58]利用258个DH株系以SSR和SRAP标记
技术共检测到20个与株高相关的QTL, 分布在N3、
N10、N13、N16和N17连锁群上, 单个QTL解释6.2%~
12.0%的表型变异; 检测到10个与第1分枝高相关的
QTL, 分布于N2、N8、N16和N17连锁群上 , 单个
QTL解释4.3%~16.2%的表型变异。本研究除在普遍检
测到株高QTL的染色体A07、A08和C06上[13-15,59-61]检测
到株高相关的QTL外, 还在A01和A10染色体上各检测
到1个株高QTL; 在A06、C05及C06染色体上检测到7
个与第1分枝高相关的QTL; 在A03和C06染色体上检
测到的分枝数QTL与前人研究结果 [3-4,13-15,59-61]一致 ,
另外在A07、C01及C04染色体上检测到3个与分枝数
相关的QTL, 单个QTL解释的表型变异为6.13%~
8.14%。
Tuberosa等[62]认为性状相关表现在QTL分析结
1032 作 物 学 报 第 41卷
果上, 可能存在控制不同性状QTL之间的紧密连锁
或同一QTL位点控制不同的性状现象即一因多效。
本研究中, 在C06染色体上检测到6个与株高、第1
分枝高及分枝数相关的QTL, 其中q2013PH-C06和
q2013FBH-C06-1的区域紧密连锁, 佐证了株高与第
1分枝高呈极显著正相关。而位于C06染色体上控制
株高和分枝数的QTL增效基因相反, 印证了株高与
分枝数呈显著负相关。前人对油菜的株高、第1分枝
高和分枝数等农艺性状的QTL定位结果不尽一致 ,
这是由于这些性状由微效多基因控制, 受环境影响
较大, 因此筛选出控制这些农艺性状的候选基因显
得极有意义。Shi等[15]在低磷和最佳施磷量条件下,
用30个GBM标记在检测到的QTL的置信区间发现
19个候选基因。其中在缺磷和最适条件下均在A09
染色体上的QTL置信区间检测到1个控制株高的候
选基因(LPR1), 分别在A01和A06染色体上的QTL置
信区域检测到与分枝数有关的候选基因PAP17、
FBHLH32和PHT1;8, 在A01、A02、A03、A04、A09
和C07染色体上的QTL区域共检测到10个与第1分枝
高相关的候选基因。Ding等[14]在70个QTL置信区域
内检测到136个与甘蓝型油菜农艺性状相关的候选
基因, 其中46个与株高相关, 30个为分枝数的候选
基因, 此结果与Shi等[15]搜索的部分结果一致。本研
究利用Shi等 [13]、Ding等 [14]和Shi等 [15]搜索或检测
图 2 甘蓝型油菜株高、第 1分枝高及分枝数 QTL在 SNP连锁群上分布图
Fig. 2 Putative QTL locations of plant height, the first branch height, and branch number on the SNP genetic map
为了简化起见, 包含 QTL的染色体上仅展示了 QTL区段的标记及染色体两端各 2个标记。
For simplicity, only show the markers in the QTL confidence intervals, along with the terminal two markers at each end of the
QTL-containing chromosomes.
第 7期 王 嘉等: 甘蓝型油菜株高、第一分枝高和分枝数的 QTL检测及候选基因筛选 1033
表 4 甘蓝型油菜基因组中 QTL置信区间候选基因与拟南芥相关基因的比对
Table 4 Alignment of candidate genes in QTL confidence interval in B. napus with the related genes in Arabidopsis thaliana
数量性状座位
QTL
拟南芥相关基因
Related genes in A. thaliana 性状
Trait 名称
Name
物理区间
Physical interval
预测基因
Gene prediction
基因名
Gene name
登录号
Accessions number
E值
E-value
PH q2013PH-A06 2371514–3103568 BnaA06g37850D PIF3 AT1G09530 6E–63
q2013PH-A06 2371514–3103568 BnaA06g05470D PHYA AT1G09570 0
q2013PH-A06 2371514–3103568 BnaA06g03910D GCR1 AT1G48270 1E–53
q2013PH-A07 13870810–15602504 BnaA07g19530D ATGID1B/ GID1B AT3G63010 0
q2013PH-A07 13870810–15602504 BnaA07g16350D WRI1 AT3G54320 4E–46
q2013PH-A08 5330869–5684994 BnaA08g05350D IAA14/SLR AT4G14550 1E–63
q2013PH-A08 5330869–5684994 BnaA08g05350D IAA7/AXR2 AT3G23050 9E–24
q2013PH-A10 12371885–12610538 BnaA10g17240D RGL3 AT5G17490 0
q2013PH-A10 12371885–12610538 BnaA10g16550D ATAPY2/APY2 AT5G18280 8E–87
q2013PH-A10 12371885–12610538 BnaA10g16110D BUD2/SAMDC4 AT5G18930 0
q2013PH-C06 21465525–24243988 BnaC06g19100D LUE1 AT1G80350 1E–100
q2014PH-A01 6898804–8402453 BnaA01g15540D GAMT1 AT4G26420 6E–80
q2014PH-A01 6898804–8402453 BnaA01g16300D WES1 AT4G27260 0
q2014PH-A06 5330869–5684994 BnaA06g00360D DEK1 AT1G55350 0
q2014PH-A07-2 16688281–16912983 BnaA07g21900D BZR1 AT1G75080 0
q2014PH-A07-1 13870810–15602504 BnaA07g19530D ATGID1B/GID1B AT3G63010 0
q2014PH-A07-1 13870810–15602504 BnaA07g16350D WRI1 AT3G54320 4E-46
q2014PH-A08 5330869–5684994 BnaA08g05350D IAA14/SLR AT4G14550 1E–63
q2014PH-A08 5330869–5684994 BnaA08g05350D IAA7/AXR2 AT3G23050 9E–24
q2014PH-C06 32013794–33834232 BnaC06g34300D T9L24.44 AT1G73340 2E–98
q2014PH-C06 32013794–33834232 BnaC06g32200D LPR2 AT1G71040 0
FBH q2013FBH-C05-1 8025904–12092352 BnaC05g16120D PHT1;8 AT1G20860 0
q2014FBH-C05-2 9283565–12092352 BnaC05g16120D PHT1;8 AT1G20860 0
q2014FBH-C05-2 9283565–12092352 BnaC05g18410D LPR1 AT1G23010 0
BN q2013BN-A03 19776618–24907552 BnaA03g48550D WES1 AT4G27260 0
q2013BN-A03 19776618–24907552 BnaA03g39630D TGH AT5G23080 1E–157
q2013BN-A03 19776618–24907552 BnaA03g39870D TCP5 AT5G60970 0
q2013BN-A03 19776618–24907552 BnaA03g43080D HAT1, HAT2 AT5G47370 3E–56
q2014BN-A07 10780228–11403972 BnaA07g11660D ARF5/IAA24 AT1G19850 5E–70
q2014BN-A07 10780228–11403972 BnaA07g12710D SHI AT5G66350 E–119
q2014BN-C01 24141843–26227700 BnaC01g26840D AXR4 AT1G54990 0
q2014BN-C04 34037875–38776266 BnaC04g32260D ERA1 AT5G40280 2E–65
q2014BN-C06 30464787–32535173 BnaC06g29550D BRC2, TCP12 AT1G68800 7E–66
q2014BN-C06 30464787–32535173 BnaC06g29980D AP1, AGL7 AT1G69120 6E–88
q2014BN-C06 30464787–32535173 BnaC06g31070D PP2A AT1G69960 1E–148
q2014BN-C06 30464787–32535173 BnaC06g32240D AT1G71090 AT1G71090 0
1034 作 物 学 报 第 41卷
表 5 拟南芥中株高(PH)、第 1分枝高(FBH)和分枝数(BN)候选基因功能
Table 5 Functions of candidate genes related to plant height (PH), the first branch height (FBH), and branch number (BN)
in A. thaliana
性状
Trait
基因
Gene
基因功能
Gene function
参考文献
Reference
PH PIF3 与启动子 LHY和 CCA1的 G-box形成三元复合物, 是植酸酶 B信号转导的负调控子。
Forms a ternary complex with G-box element of the promoters of LHY, CCA1.
Acts as a negative regulator of phyB signaling.
Sentandreu et al.[16]
PHYA 红色或远红光受体活性, 光反应的信号转导途径。
Red or far-red light photoreceptor activity, signal transducer activity response to light.
Oka et al.[17]
GCR1 类似 G蛋白偶联受体, 表达下调导致细胞分裂素敏感性降低。
Encodes a protein similar to G-coupled receptor, reduction of expression results in
decreased sensitivity to cytokinin.
Chen et al.[18]
GID1B 编码一个可溶的赤霉素受体, 介导植株中 GA 的信号传导, 影响植株的株高。
Encodes a soluble GA receptor, mediates signal transduction of GA in plants, affects
the plant height.
Iuchi et al.[19]
IAA14/SLR 参与侧根发育, 突变降低生长素诱导的基因表达。
Involved in lateral root development. Function mutation decreases auxin-inducible
gene expression.
Muto et al.[20]
IAA7/AXR2 控制幼苗向地生长和向光生长。
Control of gravitropic growth and development in light-grown seedlings.
Muto et al.[20]
RGL3 DELLA蛋白家族成员, 参与 GA信号转导。
DELLA subfamily member involved in GA signal transduction.
Feng et al.[21]
ATAPY2/APY2 利用 ATP酶和 ADP酶活性编码的一种酶。
Encodes an enzyme with ATPase and ADPase activity.
Chiu et al.[22]
BUD2/SAMDC4 类似 SAMDC的腺苷甲硫氨酸脱羧酶。
Adenosylmethionine decarboxylase; similar to SAMDC.
Cui et al.[23]
LUE1 突变体中赤霉素异常表达。
GA abnormal expression in mutant.
Stoppin et al.[24]
GAMT1 编码赤霉素甲基转移酶催化的赤霉素甲基化过程。
Encodes gibberellin methyltransferase catalyze a process of gibberellin methylation);
Varbanova et al.[25]
WES1 编码 IAA酰胺合酶, 调控生长素自由态向结合态转化。
Encodes an IAA-amido synthase, regulating from auxin free state to combine state
transformation.
Park et al.[26]
DEK1 玉米胚乳糊粉层细胞的发育基因, 编码蛋白酶基因超家族的膜蛋白。
A gene encoding a membrane protein of the calpain gene superfamily required for
aleurone cell development in the endosperm of maize grains.
Liang et al.[27]
BZR1 油菜素内酯信号通路中的关键转录因子, 调控下游响应基因的表达。
Encodes a positive regulator of the brassinosteroid (BR) signalling pathway that
mediates both downstream BR responses.
Wang et al.[28]
WRI1 编码 AP2转录因子, 在种子成熟时通过 LEC2启动子来调节脂肪酸积累。
Encoding AP2 transcription factor, regulating fatty acid accumulation by LEC2
promoter during seed maturation.
Kim et al.[29]
T9L24.44 参与植物的细胞伸长和分裂, 对植物生长发育起至关重要的作用。
Involved in plant cell elongation and cell division, play a vital role in plants growth
and development.
Theologis et al.[30]
LPR2 参与调控低磷胁迫下植物初生根的发育。
Involved in the regulation of plant primary roots development under low phosphorus stress.
Wang et al.[31]
第 7期 王 嘉等: 甘蓝型油菜株高、第一分枝高和分枝数的 QTL检测及候选基因筛选 1035
(续表 5)
性状
Trait
基因
Gene
基因功能
Gene function
参考文献
Reference
FBH PHT1;8 编码的磷酸盐转运蛋白 ORYsa, 提高了植株磷素吸收能力。
Improve the plant ability of phosphorus absorption by codes phosphate transporter ORYsa.
Remy et al.[32]
LPR1 参与调控生长素响应低磷胁迫。
Involved in the regulation of auxin response to low phosphorus stress.
Wang et al.[31]
BN TGH 甘油三酯水解酶, 突变体出现发育缺陷, 包括植株矮化, 畸形和导管减少。
Triglyceride hydrolase, Mutants display developmental defects, including reduced
plant height, polycotyly, and reduced vascularization.
Calderon et al.[33]
TCP5 调控植物特异的形态建成。
Regulation of plant specific morphogenesis.
Efroni et al.[34]
HAT2 在茎和根组织中调控生长素形态扮演一个相反的角色。
Plays opposite roles in the shoot and root tissues in regulating auxin-mediated
morphogenesis.
Roig et al.[35]
ARF5/IAA24 编码一种转录因子(IAA24)介导胚胎轴的形成和微管的发育。
Encodes a transcription factor (IAA24) mediating embryo axis formation and vascular
development.
Wu et al.[36]
SHI 参与赤霉素响应。
Involved in the response to gibberellic acid.
Kuusk et al.[37]
AXR4 响应机械刺激和生长素刺激, 生长素的极性运输。
Response to mechanical stimulus, response to auxin stimulus, auxin polar transport.
Dharmasiri et al.[38]
ERA1 基因缺失或抑制表达可以提高植株抗旱能力。
Involved in meristem organization and ABA-mediated signal transduction pathway,
and response to abscisic acid and drought.
Yalovsky et al.[39]
BRC2/TCP12 编码一个 TCP转录因子, 抑制腋芽的发育与生长。
Encodes a TCP transcription factor, arrests axillary bud development and prevents
axillary bud outgrowth.
Aguilar et al.[40]
AP1/AGL7 调控花序分生组织向花分生组织的转变, 调控花器官形态发育。
Regulation of transformation that inflorescence meristem into a floral meristem, and
floral organogenesis.
Yalovsky et al.[41]
PP2A 参与生长素极性运输。
Involved in auxin polar transport.
Gao et al.[42]
AT1G71090 参与生长素的极性运输和跨膜运输。
Involved in auxin polar transport and transmembrane transport.
Barbez et al.[43]
到的相关候选基因与本文检测到的QTL置信区间序
列进行同源比对, 检测到 20 个株高候选基因、3 个
第 1分枝高候选基因和 12个分枝数候选基因。2个
环境中在A07染色体上重复检测到的QTL置信区间
检测到与株高相关的候选基因 ATGID1B/GID1B 和
WRI1, A08染色体上重复检测到的 QTL置信区间检
测到 SLR/IAA14和 AXR2/IAA72个与株高相关的候选
基因。在具有部分置信区间重叠的 QTL q2013FBH-
C05-1和 q2014FBH-C05-2区间均检测到第 1分枝高
候选基因 PHT1;8, 与 Shi等[15]检测到该基因的位置
(A01 和 A06 染色体上)存在差异。在位于前人检测
到分枝数 QTL 的 A03 和 C06 染色体上[3-4,13-15,59-61]
的 QTL置信区间内, 我们分别检测到 4个分枝数候
选基因, 匹配 E 值介于 0~3E–56 之间, 暗示这两个
QTL 极有可能为控制分枝数的遗传位点, 其结果有
待进一步验证。检测到株高、第 1 分枝高及分枝数
的候选基因较少, 可能与 QTL 的贡献率小有关, 检
测到第 1 分枝高相关候选基因较少的原因还可能是
在本研究中用于搜索的第 1分枝高相关基因较少。
4 结论
同一环境下株高与第 1 分枝高呈极显著相关,
1036 作 物 学 报 第 41卷
与分枝数呈负相关, 第 1 分枝高与分枝数在 2 个环
境中均呈极显著负相关。在 2 个环境下重复检测到
3个与株高相关的 QTL, 分别位于 A07、A08和 C06
染色体上, 在 A07染色体上重复检测到的 QTL置信
区 间 检 测 到 2 个 与 株 高 相 关 的 候 选 基 因
ATGID1B/GID1B 和 WRI1, 在 A08 染色体上的 QTL
置信区间检测到 3 个与株高相关的候选基因(SLR/
IAA14、AXR2/IAA7和 BIN2/DWF12); 第 1分枝高、
分枝数各重复检测到 1 个 QTL, 分别位于 C05 和
C06 染色体上。位于 A03 和 C06 染色体上的 QTL
q2013BN-A03和 q2014BN-C06分别检测到 4个分枝
数候选基因, 极有可能为控制分枝数的遗传位点。
References
[1] James C. Global review of commercialized transgenic crops.
Curr Sci, 2003, 84: 303–309
[2] 易斌, 陈伟, 马朝芝, 傅廷栋, 涂金星. 甘蓝型油菜产量及相
关性状的 QTL分析. 作物学报, 2006, 32: 676–682
Yi B, Chen W, Ma C Z, Fu T D, Tu J X. Mapping of quantitative
trait loci for yield and yield components in Brassica napus L.
Acta Agron Sin, 2006, 32: 676–682 (in Chinese with English ab-
stract)
[3] Li W, Lin Z X, Zhang X L. A novel segregation distortion in in-
traspecific population of Asian cotton (Gossypium arboretum L.)
detected by molecular markers. J Genet Genomics, 2007, 34:
634–640
[4] Chen W, Zhang Y, Liu X P, Chen B Y, Tu J X, Fu T D. Detection
of QTL for six yield-related traits in oilseed rape (Brassica napus)
using DH and immortalized F2 populations. Theor Appl Genet,
2007, 115: 849–858
[5] 张倩. 甘蓝型油菜主要株型性状的遗传分析和 QTL初步定位.
西南大学硕士学位论文, 重庆, 2013
Zhang Q. Genetic Effects Analysis and QTL Mapping of Major
Plant-type Traits in Brassica napus L. PhD Dissertation of
Southwest University, Chingqing, China, 2013 (in Chinese with
English abstract)
[6] Zhang S H, Fu T D, Zhu J C, Wang J P, Wen Y C, Ma C Z, Jiang
Y Z. QTL mapping and epistasis analysis for plant height and
height to the first branch in rapeseed (Brassica napus L.).
BIOTECHNOLOGY: Gene Clon Funct Anal, 2007, 2: 232–235
[7] Wang D G, Fan J B, Siao C J, Berno A, Young P, Sapolsky R,
Ghandour G, Perkins N, Winchester E, Spencer J, Kruglyak L,
Stein L, Hsie L, Topaloglou T, Hubbell E, Robinson E, Mittmann
M, Morris M S, Shen N, Kilburn D, Rioux J, Nusbaum C, Rozen
S, Hudson T J, Lipshutz R, Chee M, Lander E S. Large-scale
identification, mapping and genotyping of single-nucleotide
polymorphisms in the human genome. Science, 1998, 280:
1077–1082
[8] Chalhoub B, Denoeud F, Liu S, Parkin I A, Tang H, Wang X,
Chiquet J, Belcram H, Tong C, Samans B, Corréa M, Da Silva C,
Just J, Falentin C, Koh C S, Le Clainche I, Bernard M, Bento P,
Noel B, Labadie K, Alberti A, Charles M, Arnaud D, Guo H, Da-
viaud C, Alamery S, Jabbari K, Zhao M, Edger P P, Chelaifa H,
Tack D, Lassalle G, Mestiri I, Schnel N, Le Paslier M C, Fan G,
Renault V, Bayer P E, Golicz A A, Manoli S, Lee T H, Thi V H,
Chalabi S, Hu Q, Fan C, Tollenaere R, Lu Y, Battail C, Shen J,
Sidebottom C H, Wang X, Canaguier A, Chauveau A, Bérard A,
Deniot G, Guan M, Liu Z, Sun F, Lim Y P, Lyons E, Town C D,
Bancroft I, Wang X, Meng J, Ma J, Pires J C, King G J, Brunel D,
Delourme R, Renard M, Aury J M, Adams K L, Batley J, Snow-
don R J, Tost J, Edwards D, Zhou Y, Hua W, Sharpe A G, Pater-
son A H, Guan C, Wincker P. Early allopolyploid evolution in the
post-neolithic Brassica napus oilseed genome. Science, 2014,
345: 950–953
[9] Liu L Z, Qu C M, Wittkop B, Yi B, Xiao Y, He H J, Snowdon R,
Li J. A high-density SNP map for accurate mapping of seed fibre
QTL in Brassica napus L. PLoS One, 2013, 8: e83052
[10] Wang S, Basten C J, Zeng Z B. Windows QTL Cartographer. Ver.
2.5 [computer program] department of statistics, North Carolina
State University, Raleigh, NC, 2006. http://statgen.ncsu.edu/
qtlcart/WQTL Cart.htm
[11] Lander E S, Botstein D. Mapping mendelian factors underlying
quantitative traits using RFLP linkage maps. Genetics, 1989, 121:
185–199
[12] Mccouch S R, Cho Y G, Yano M, Paul E, Blinstrub M, Morishima
H, Kinoshita T. Report on QTL nomenclature. Rice Genet Newsl,
1997, 14
[13] Shi J Q, Li R Y, Qiu D, Jiang C C, Long Y, Morgan C, Nancroft I,
Zhao J Y, Meng J L. Unraveling the complex trait of crop yield
with quantitative trait loci mapping in Brassica napus. Genetics,
2009, 182: 851–861
[14] Ding G D, Zhao Z K, Liao Y, Hu Y F, Shi L, Long Y, Xu F S.
Quantitative trait loci for seed yield and yield-related traits, and
their responses to reduced phosphorus supply in Brassica napus.
Ann Bot, 2012, 109: 747–759
[15] Shi T X, Li R Y, Zhao Z K, Ding G D, Long Y, Meng J L, Xu F S,
Shi L. QTL for yield traits and their association with functional
genes in response to phosphorus deficiency in Brassica napus.
PLOS One, 2013, 8: e54559
[16] Sentandreu M, Martín G, González-Schain N, Leivar P, Soy J,
Tepperman J M, Quail P H, Monte E. Functional profiling identi-
fies genes involved in organ-specific branches of the PIF3 regu-
latory network in Arabidopsis. Plant Cell, 2011, 23: 3974–3991
[17] Oka Y, Ono Y, Toledo-Ortiz G, Kokaji K, Matsui M, Mochizuki N,
Nagatani A. Arabidopsis phytochrome a is modularly structured
to integrate the multiple features that are required for a highly
sensitized phytochrome. Plant Cell, 2012, 24: 2949–2962
[18] Chen J G, Pandey S, Huang J, Alonso J M, Ecker J R, Assmann S
M, Jones A M. GCR1 can act independently of heterotrimeric
G-protein in response to brassinosteroids and gibberellins in
Arabidopsis seed germination. Plant Physiol, 2004, 135: 907–915
[19] Iuchi S, Suzuki H, Kim Y C, Iuchi A, Kuromori T, Ueguchi-
Tanaka M, Asami T, Yamaguchi I, Matsuoka M, Kobayashi M,
Nakajima M. Multiple loss-of-function of Arabidopsis gibberellin
receptor AtGID1s completely shuts down a gibberellin signal.
Plant J, 2007, 50: 958–966
[20] Muto H, Watahiki M K, Nakamoto D, Kinjo M, Yamamoto K T.
Specificity and similarity of functions of the Aux/IAA genes in
第 7期 王 嘉等: 甘蓝型油菜株高、第一分枝高和分枝数的 QTL检测及候选基因筛选 1037
auxin signaling of Arabidopsis revealed by promoter-exchange
experiments among MSG2/IAA19, AXR2/IAA7, and SLR/IAA14.
Plant Physiol, 2007, 144: 187–196
[21] Feng S, Martinez C, Gusmaroli G, Wang Y, Zhou J, Wang F, Chen
L, Yu L, Iglesias-Pedraz JM, Kircher S, Schäfer E, Fu X, Fan L
M, Deng X W. Coordinated regulation of Arabidopsis thaliana
development by light and gibberellins. Nature, 2008, 451:
475–479
[22] Chiu T Y, Christiansen K, Moreno I, Lao J, Loqué D, Orellana A,
Heazlewood J L, Clark G, Roux S J. AtAPY1 and AtAPY2 func-
tion as Golgi-localized nucleoside diphosphatases in Arabidopsis
thaliana. Plant Cell Physiol, 2012, 53: 1913–1925
[23] Cui X, Ge C, Wang R, Wang H, Chen W, Fu Z, Jiang X, Li J,
Wang Y. The BUD2 mutation affects plant architecture through
altering cytokinin and auxin responses in Arabidopsis. Cell Res,
2010, 20: 576–586
[24] Stoppin-Mellet V, Gaillard J, Vantard M. Katanin’s severing ac-
tivity favors bundling of cortical microtubules in plants. Plant J,
2006, 46: 1009–1017
[25] Varbanova M, Yamaguchi S, Yang Y, McKelvey K, Hanada A,
Borochov R, Yu F, Jikumaru Y, Ross J, Cortes D, Ma C J, Noel J
P, Mander L, Shulaev V, Kamiya Y, Rodermel S, Weiss D, Pich-
ersky E. Methylation of gibberellins by Arabidopsis GAMT1 and
GAMT2. Plant Cell, 2007, 19: 32–45
[26] Park J E, Seo P J, Lee A K, Jung J H, Kim Y S, Park C M. An
Arabidopsis GH3 gene, encoding an auxin-conjugating enzyme,
mediate sphytochrome B-regulated light signals in hypocotyl
growth. Plant Cell Physiol, 2007, 48: 1236–1241
[27] Liang Z, Demko V, Wilson R C, Johnson K A, Ahmad R, Perroud P
F, Quatrano R, Zhao S, Shalchian-Tabrizi K, Otegui M S, Olsen O
A, Johansen W. The catalytic domain Cyspc of the DEK1 calpain is
functionally conserved in land plants. Plant J, 2013, 75: 742–754
[28] Wang Z Y, Nakano T, Gendron J, He J, Chen M, Vafeados D,
Yang Y, Fujioka S, Yoshida S, Asami T, Chory J. Nuclear-
localized BZR1 mediates brassinosteroid-induced growth and
feedback suppression of brassinosteroid biosynthesis. Dev Cell,
2002, 2: 505–513
[29] Kim H U, Jung S J, Lee K R, Kim E H, Lee S M, Roh K H, Kim
J B. Ectopic overexpression of castor bean LEAFY
COTYLEDON2 (LEC2) in Arabidopsis triggers the expression of
genes that encode regulators of seed maturation and oil body
proteins in vegetative tissues. FEBS Open Biol, 2013, 23: 25–32
[30] Theologis A, Ecker J R, Palm C J, Federspiel N A, Kaul S, White
O, Alonso J, Altafi H, Araujo R, Bowman C L, Brooks S Y,
Buehler E, Chan A, Chao Q, Chen H, Cheuk R F, Chin C W,
Chung M K, Conn L, Conway A B, Conway A R, Creasy T H,
Dewar K, Dunn P, Etgu P, Feldblyum T V, Feng J, Fong B, Fujii
C Y, Gill J E, Goldsmith A D, Haas B, Hansen N F, Hughes B,
Huizar L, Hunter J L, Jenkins J, Johnson-Hopson C, Khan S,
Khaykin E, Kim C J, Koo H L, Kremenetskaia I, Kurtz D B,
Kwan A, Lam B, Langin-Hooper S, Lee A, Lee J M, Lenz C A,
Li J H, Li Y, Lin X, Liu S X, Liu Z A, Luros J S, Maiti R, Mar-
ziali A, Militscher J, Miranda M, Nguyen M, Nierman W C, Os-
borne B I, Pai G, Peterson J, Pham P K, Rizzo M, Rooney T,
Rowley D, Sakano H, Salzberg S L, Schwartz J R, Shinn P,
Southwick A M, Sun H, Tallon L J, Tambunga G, Toriumi M J,
Town C D, Utterback T, Van Aken S, Vaysberg M, Vysotskaia V
S, Walker M, Wu D, Yu G, Fraser C M, Venter J C, Davis R W.
Sequence and analysis of chromosome 1 of the plant Arabidopsis
thaliana. Nature, 2000, 408: 816–820
[31] Wang X, Du G, Wang X, Meng Y, Li Y, Wu P, Yi K. The function
of LPR1 is controlled by an element in the promoter and is inde-
pendent of SUMOE3 ligase SIZ1 in response to low Pi stress in
Arabidopsis thaliana. Plant Cell Physiol, 2010, 51: 380–394
[32] Remy E, Cabrito T R, Batista R A, Teixeira M C, Sá-Correia I,
Duque P. The Pht1;9 and Pht1;8 transporters mediate inorganic
phosphate acquisition by the Arabidopsis thaliana root during
phosphorus starvation. New Phytol, 2012, 195: 356–371
[33] Calderon-Villalobos L I, Kuhnle C, Dohmann E M, Li H, Bevan
M, Schwechheimer C. The evolutionarily conserved TOUGH
protein is required for proper development of Arabidopsis
thaliana. Plant Cell, 2005, 17: 2473–2485
[34] Efroni I, Blum E, Goldshmidt A, Eshed Y. A protracted and dy-
namic maturation schedule underlies Arabidopsis leaf develop-
ment. Plant Cell, 2008, 20: 2293–2306
[35] Roig-Villanova I, Bou J, Sorin C, Devlin P F, Martínez-García J F.
Identification of primary target genes of phytochrome signaling:
early transcriptional control during shade avoidance responses in
Arabidopsis. Plant Physiol, 2006, 141: 85–96
[36] Wu J, Wang F, Cheng L, Kong F, Peng Z, Liu S, Yu X, Lu G.
Identification, isolation and expression analysis of auxin response
factor (ARF) genes in Solanum lycopersicum. Plant Cell Rep,
2011, 30: 2059–2073
[37] Kuusk S, Sohlberg J J, Magnus Eklund D, Sundberg E. Function-
ally redundant SHI family genes regulate Arabidopsis gynoecium
development in a dose-dependent manner. Plant J, 2006, 47:
99–111
[38] Dharmasiri S, Swarup R, Mockaitis K, Dharmasiri N, Singh S K,
Kowalchyk M, Marchant A, Mills S, Sandberg G, Bennett M J,
Estelle M. AXR4 is required for localization of the auxin influx
facilitator AUX1. Science, 2006, 312: 1218–1220
[39] Yalovsky S, Kulukian A, Rodríguez-Concepción M, Young C A,
Gruissem W. Functional requirement of plant farnesyl transferase
during development in Arabidopsis. Plant Cell, 2000, 12:
1267–1278
[40] Aguilar-Martínez J A, Poza-Carrión C, Cubas P. Arabidopsis
BRANCHED1 acts as an integrator of branching signals within
axillary buds. Plant Cell, 2007, 19: 458–472
[41] Yalovsky S, Rodríguez-Concepción M, Bracha K, Toledo-Ortiz G,
Gruissem W. Prenylation of the floral transcription factor
APETALA1 modulates its function. Plant Cell, 2000, 12:
1257–1266
[42] Gao H B, Chu Y J, Xue H W. Phosphatidic acid (PA) binds
PP2AA1 to regulate PP2A activity and PIN1 polar localization.
Mol Plant, 2013, 6: 1692–1702
[43] Barbez E, Kubeš M, Rolčík J, Béziat C, Pěnčík A, Wang B,
Rosquete M R, Zhu J, Dobrev P I, Lee Y, Zažímalovà E, Petrášek
J, Geisler M, Friml J, Kleine-Vehn J. A novel putative auxin
carrier family regulates intracellular auxin homeostasis in plants.
Nature, 2012, 485: 119–122
[44] 王通强, 万育麟. 甘蓝型油菜数量性状的聚类和通径分析. 种
子, 1991, (5): 14–17
1038 作 物 学 报 第 41卷
Wang T Q, Wan Y L. The quantitative character’s aggregation
and path analysis of Brassica napus. Seed, 1991, (5): 14–17 (in
Chinese with English abstract)
[45] 王瑞, 李加纳, 谌利, 唐章林. 甘蓝型黄籽油菜主要性状的通
径分析. 中国农学通报, 2004, 20(6): 325–326
Wang R, Li J N, Chen L, Tang Z L. Path analysis for main char-
acters inyellow-seeded rapeseed lines (Brassica napus L). Chin
Agric Sci Bull, 2004, 20(6): 325–326 (in Chinese with English
abstract)
[46] Palosz T, Sienkowski A, Grala B. Statistic relations between
several factors of rape production. Arch Phytopathol Plant Prot,
1994, 29: 101–106
[47] 周小丽, 王通强. 双低油菜主要农艺性状的通径分析. 种子,
2005, 24(1): 70–72
Zhou X L, Wang T Q. Path analysis of main agronomic charac-
ters of double low rape. Seed, 2005, 24(1): 70–72 (in Chinese
with English abstract)
[48] 杨安中, 彭春华. 油菜单株产量与若干农艺性状的相关分析.
安徽农学通报, 2006, 12(2): 33–34
Yang A Z, Peng C H. Relevant analysis yield per plant and some
agronomic traits in rape. Anhui Agric Sci Bull, 2006, 12(2):
33–34 (in Chinese with English abstract)
[49] Tuncturk M, Ciftci V. Relationships between yield and some
yield components in rapeseed (Brassica napus L.) cultivars by
using correlation and path analysis. Pak J Bot, 2007, 39: 81–84
[50] Xu J S, Song X, Cheng Y, Zou X L, Zeng L, Qiao X, Lu G Y, Fu
G P, Qu Z, Zhang X K. Identification of QTLs for branch number
in oilseed rape (Brassica napus L.). J Genet Genomics, 2014, 41:
557–559
[51] Chen W, Zhang Y, Liu X, Chen B, Tu J, Tingdong F. Detection of
QTL for six yield–related traits in oil seed rape (Brassica napus)
using DH and irnrnortalized F2 populations. Theor Appl Genet,
2007, 115: 849–858
[52] Quijada P A, Maureira I J, Osbom T C. Confirmation of QTL
controlling seed yield in spring canola (Brassica napus L.) hy-
brids. Mol Breed, 2004, 13: 193–200
[53] Quijada P A, Udall J A, Lambert B, Osborn T C. Quantitative trait
analysis of seed yield and other complex traits in hybrid spring
rape seed (Brassica napus L.): 1. Identification of genomic re-
gions from winter germplasm. Theor Appl Genet, 2006, 113:
549–561
[54] Udall J A, Quijada P A, Lambert B, Osbom T C. Quantitative trait
analysis of seed yield and other complex traits in hyrid spring
rape seed (Brassica napus L.): 2. Identification of alleles from
unadapted germplasm. Theor Appl Genet, 2006, 113: 597–609
[55] Zhao J, Beeker H C, Zhang D, Zhang Y, Ecke W. Conditional
QTL mapping of oil content in rape seed with respect to protein
content and traits related to plant development and grain yield.
Theor Appl Genet, 2006, 113: 33–38
[56] Jian H J, Wei L J, Li J N, Xu X F, Chen L, Liu L Z. Mapping
quantitative traits loci for seed glucosinolate content in Brassica
napus using high-density SNP map. Acta Agron Sin, 2014, 40:
1386–1391
[57] 荐红举, 肖阳, 李加纳, 马珍珍, 魏丽娟, 刘列钊. 利用 SNP
遗传图谱定位盐、旱胁迫下甘蓝型油菜种子发芽率的 QTL. 作
物学报, 2014, 40: 629–635
Jian H J, Xiao Y, Li J N, Ma Z Z, Wei L J, Liu L Z. QTL mapping
for germination percentage under salinity drought stresses in
Brassica napus L. Acta Agron Sin, 2014, 40: 629–635 (in Chinese
with English abstract)
[58] Chen W, Zhang Y, Liu X P, Chen B Y, Tu J X, Fu T D. Detection
of QTL for six yield-related traits in oilseed rape (Brassica napus)
using DH and immortalized F2 populations. Theor Appl Genet,
2007, 115: 849–858
[59] Ramchiary N, Padmaja K L. Mapping of yield influencing QTL
in Brassica juncea: implications for breeding of a major oilseed
crop of dryland areas. Theor Appl Genet, 2007, 115: 807–817
[60] Butruille D V, Guries R P, Osbom T C. Linkage analysis of mo-
lecular markers and quantitative trait loci in populations of inbred
backcross lines of Brassica napus L. Genetics, 1999, 153:
949–964
[61] Yi B, Chen W, Ma C, Fu T D, Tu J X. Mapping of quantitative
trait loci for yield and yield components in Brassica napus L.
Acta Agron Sin, 2006, 32: 676–682
[62] Tuberosa R, Salvi S, Sanguineti M C, Landi P, Maecaferri M,
Conti S. Mapping QTLs regulating morpho-physiological
traits and yield: casestudies, short comings and perspectives
in drought-stressed maize. Ann Bot, 2002, 89: 941–96