秦伟伟, 李永祥, 李春辉, 陈林, 吴迅, 白娜, 石云素, 宋燕春, 张登峰, 王天宇^*, 黎裕^*
中国农业科学院作物科学研究所, 北京 100081
^* 通讯作者(Corresponding authors): 王天宇, E-mail: wangtianyu@caas.cn; 黎裕, E-mail: liyu03@caas.cn 第一作者联系方式: E-mail: 13811531756@139.com
收稿日期:2015-03-09 接受日期:2015-06-01网络出版日期:2015-07-01基金:本研究由国家重点基础研究发展计划(973计划)项目(2011CB100105, 2014CB138200), 国家自然科学基金项目(91335206)和中国农业科学院科技创新工程项目资助

摘要籽粒大小及百粒重是决定玉米产量的重要因素。为解析籽粒性状遗传基础, 本研究以玉米自交系黄早四(HZS)和Mo17为亲本, 构建包含130个重组自交系(recombination inbred line, RIL)的RIL群体。基于GBS (genotyping- by-sequencing)技术获得的高密度多态性SNP (single nucleotide polymorphism)位点, 构建了包含1262个Bin标记的高密度遗传图谱。采用完备区间作图法, 对5个环境条件下的粒长、粒宽、百粒重、粒长/粒宽4个性状分别进行QTL (quantitative trait locus)定位, 共检测到30个QTL。利用5个环境性状均值, 共检测到11个QTL。其中粒长主效QTL qklen1、粒长/粒宽主效QTL qklw1在3个单环境条件下均被检测到, 且定位在第1染色体相邻区域, 物理位置分别为210~212 Mb、207~208 Mb, 表型贡献率分别为22.60%和26.79%, 被认为是控制玉米籽粒形状的主效位点。针对第1染色体207~212 Mb区间, 采用成组法 t检验, 对黄早四(受体)和Mo17 (供体)构建的BC₃F₁回交群体进行单标记分析。结果表明, 在BC₃F₁群体中 qklen1和 qklw1同样具有显著的遗传效应。本研究结果不仅为分子标记辅助选择籽粒性状提供了实用标记, 而且为主效基因的进一步精细定位和候选基因挖掘奠定了基础。

关键词:玉米( Zea mays L.); 籽粒性状; 数量性状位点(QTL); 高密度遗传图谱
QTL Mapping for Kernel Related Traits Based on a High-Density Genetic Map
QIN Wei-Wei, LI Yong-Xiang, LI Chun-Hui, CHEN Lin, WU Xun, BAI Na, SHI Yun-Su, SONG Yan-Chun, ZHANG Deng-Feng, WANG Tian-Yu^*, LI Yu^*
Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing 100081, China

AbstractKernel size and weight are major determinants of grain yield. For understanding molecular mechanisms of kernel related traits, a recombinant inbred line (RIL) mapping population including 130 families was developed from the cross of two maize elite inbreds, Huangzaosi (HZS) and Mo17. By using the approach of GBS (genotyping-by-sequencing), the high-density genetic map with 1262 bin markers was constructed. QTLs for kernel related traits were identified by stepwise regression (RSTEP-LRT) using Windows QTL ICI-Mapping software under five environments. In total, 30 QTLs were detected under single environment and 11 QTLs were detected under the joint environment. The kernel length major QTL qklen1and the length/width major QTL qklw1 were found in the adjoining regions with a strong genetic effect under three environments. QTL qklen1was located in a region of 210-212 Mb on chromosome 1 with explained 22.60% of phenotypic variance, and qklw1 was located in a region of 207-208 Mb on chromosome 1 with explained 26.79% of phenotypic variance. In addition, for further verification, another introgression population of BC₃F₁ was developed with Mo17 as the donor parent and HZS as the recurrent parent. The result of the single marker analysis suggested that qklen1 and qklw1 also had a significantly genetic effect, which is similar to that in the RIL population. The present study provides a good basis for studying genetic mechanism and molecular marker assisted selection for the improvement of kernel related traits in maize.

Keyword:Maize; Kernel related trait; QTL; High-density genetic map
Show Figures
Show Figures








玉米作为重要粮、经、饲兼用作物, 据预测全球需求量到2050年将增加一倍^[1]。因此, 高产仍是玉米育种的主要目标。玉米产量性状是最复杂的性状之一^[2]。将产量性状分解为若干构成因子, 不仅可提高产量相关基因挖掘效率^[3], 而且有利于阐述影响产量的遗传基础^{[4, 5]}。玉米的产量构成因子可分为百粒重、穗粒数、单位面积有效穗数。其中, 百粒重具有较高的遗传力, 是玉米育种过程中的重要选择性状。百粒重主要取决于籽粒的大小和形状^[6]。在作物的选择驯化过程中, 粒形和粒重也被认为是产量的重要构成因子^[7]。此外, 玉米籽粒性状(粒重、粒长、粒宽、粒厚)与产量均表现出较高的相关性^[8], 且相对于产量而言, 具有受环境影响小的特点^[9]。因此, 将籽粒相关性状作为研究目标是解析玉米产量遗传基础的重要途径。
玉米基因组学的迅速发展^[10]和新型标记的开发, 推动了玉米籽粒相关数量性状位点(QTL)的挖掘进程。如Liu等^[11]利用F_2:3群体及256个SSR标记定位到28个籽粒相关QTL; Zhang等^[12]利用永久性F₂群体及217个SSR标记定位到54个籽粒相关QTL; Nikolic等^[13]利用F₃群体及234个分子标记定位到18个籽粒相关QTL。但是, 往往由于构建遗传连锁图谱时用到的标记数量较少, 使遗传图谱分辨率较低, 定位到的QTL置信区间大都在10~20 cM之间, 导致其结果难以用于分子标记辅助育种, 也限制了玉米产量遗传基础的深度解析。因此, 有必要利用高密度分子标记构建遗传连锁图谱, 提高籽粒相关性状QTL定位的精确度。
近年来, 随着单核苷酸多态性(SNP)标记的开发成本和使用成本大幅降低, 利用SNP标记鉴定遗传作图群体中个体的基因型并构建高密度遗传图谱得以实现。如Colasuonno等^[14]利用5019个SNP标记构建的高密度遗传连锁图谱对硬粒小麦黄色素含量进行QTL定位, 标记间平均距离为0.36 cM。Guo等^[15]利用玉米永久性F₂群体整合18 840个SNP标记, 构建了含3184个区块(Bin)的高密度遗传连锁图谱, 标记间平均距离缩小到0.84 cM, 并对杂种优势进行遗传解析。Zou等^[16]利用高粱RIL群体, 通过构建的含3418个Bin的高密度遗传连锁图谱, 标记间平均距离为0.50 cM, 对花期和株型相关性状进行了QTL定位。这些研究在很大程度上缩小了QTL置信区间, 对开发紧密连锁的分子标记和下一步的QTL精细定位提供了更多的遗传信息。本研究在不同环境条件下对玉米粒长、粒宽、百粒重和粒长/粒宽进行了表型鉴定, 并借助GBS (genotyping-by- sequencing)技术获得的SNP标记构建高密度遗传图谱, 进行玉米籽粒性状的QTL定位, 旨在利用高密度遗传图谱检测玉米籽粒性状QTL, 从而为分子标记辅助育种和玉米产量相关性状的主效QTL精细定位及候选基因挖掘提供新的信息和科学依据。
1 材料与方法1.1 试验材料以玉米自交系黄早四和Mo17为亲本配制组合, 从F₂代开始采用“ 单粒传法” 连续自交至F₈代, 构建成一套包含130个重组自交系(recombination inbred line, RIL)群体。另外, 以玉米自交系黄早四为受体亲本, 以Mo17为供体亲本, 从BC₁开始每个世代取混合花粉与黄早四连续回交3代, 获得含152个单株的BC₃F₁高代回交群体。
1.2 田间试验及表型鉴定对RIL群体2009年在北京、新疆、河南及2010年在新疆、河南共5个环境条件下进行表型鉴定。在每个环境下均采用随机区组设计, 单行区, 2次重复, 行长3.0 m, 每行11株, 行距0.6 m。植株成熟后, 从每个小区收获5株果穗用于籽粒相关性状的考察。所考察的性状包括粒长、粒宽和百粒重。粒长和粒宽为每个果穗中间部分连续籽粒的10籽长和10粒宽, 单位为cm。百粒重为每个小区果穗脱粒后, 3次随机取样100粒的平均值, 单位为g。粒长/粒宽为二级性状, 即粒长与粒宽的比值。
2012年在北京对BC₃F₁群体进行表型鉴定。对所有材料均采用单粒播种, 行长3.0 m, 每行定苗11株, 行距0.6 m, 株距25 cm。玉米拔节期单株编号挂牌, 待植株完全成熟后, 按编号单株收获。取果穗中部籽粒, 测量粒长、粒宽和百粒重, 每个果穗重复测量3次。
全部田间试验中的施肥、灌溉、防虫和除草同当地大田管理一致。
1.3 DNA的提取和基因型鉴定对RIL群体, 在植株长到六叶时混合取每个家系的幼嫩叶片, 采用CTAB法^[17]提取基因组DNA。利用简化基因组测序(genotyping-by-sequencing, GBS)技术^[18]鉴定家系基因型。
对BC₃F₁群体, 在玉米拔节期取亲本和回交群体单株新鲜叶片, 利用CTAB法提取DNA^[17]。参考中国农业大学赖锦盛教授课题组提供的黄早四和Mo17重测序结果, 针对RIL群体中定位到的QTL qklen1和qklw1的基因组区段, 利用Primer5.0设计开发InDel标记。通过PCR产物扩增, 8.0%非变性聚丙烯酰胺凝胶电泳, 0.1% AgNO₃染色, 甲醛和NaOH显色, 进行BC₃F₁群体基因型鉴定。
1.4 高密度遗传连锁图谱的构建基于GBS测序^[18], 在基因组范围内筛选多态性SNP标记。从原始基因型数据中过滤掉等位基因频率AF< 0.05的SNP位点。采用重组最大简约法(maximum parsimony of recombination, MPR), 根据RIL群体的基因型推测双亲基因型。根据双亲SNP基因型将群体基因型数值化。然后将连续多个SNP标记且具有相同的基因型概括为一个区块(Bin), 利用其构建RIL群体的重组区块图谱(recombination bin map)^[19]。以每个重组区块作为一个标记, 利用JoinMap4.0^[20]软件构建其高密度遗传图谱。
1.5 QTL检测及验证利用QTL IciMapping V3.3软件中的完备区间作图法对RIL群体的籽粒相关性状进行QTL作图及遗传效应分析^{[21, 22]}, 作图步长为1.00 cM, PIN设置为0.001, LOD阈值设置为2.50。在BC₃F₁群体中, 利用SAS软件中基于单因子的成组法t检验对精细定位候选区段的遗传效应作进一步验证。

2 结果与分析2.1 玉米籽粒相关性状的表型分析黄早四籽粒属于硬粒型, 籽粒宽、饱满; 而Mo17籽粒属于马齿型, 籽粒长、扁平, 在5种不同环境下亲本间籽粒性状的表型具有明显的差异(表1)。因此, 利用这两个自交系构建的群体后代存在较丰富的表型变异(表2)。例如, 在RIL群体中, 家系粒长均值为9.03 cm, 变异范围为7.67~10.26 cm。粒长、粒宽、粒长/粒宽、百粒重等4个籽粒相关性状的表型频率均呈正态分布或近似正态分布(图1), 可作为数量性状进行QTL作图。此外, 对这4个籽粒相关性状多环境联合方差分析(表3)显示, 粒长、粒宽、粒长/粒宽、百粒重的遗传力分别为0.861、0.860、0.869、0.855。
表1
Table 1
表1(Table 1)

表1 5个不同环境下亲本间籽粒性状的表型差异 Table 1 Difference of the kernel related traits between parents under five environments 性状 Trait 2009北京 2009 Beijing 2009河南 2009 Henan 2009新疆 2009 Xinjiang 2010河南 2010 Henan 2010新疆 2010 Xinjiang 黄早四 HZS Mo17 黄早四 HZS Mo17 黄早四 HZS Mo17 黄早四 HZS Mo17 黄早四 HZS Mo17 粒长KLEN (cm) 8.83 9.72 7.79 9.87 9.03 10.28 8.42 8.63 8.83 11.00 粒宽KWID (cm) 8.83 8.60 8.53 8.75 8.65 9.21 9.00 7.95 8.64 9.33 百粒重KWEI (g) 20.53 27.35 20.00 28.95 26.93 29.70 21.72 24.55 24.87 33.05 粒长/粒宽KLEN/KWID 1.00 1.13 0.91 1.13 1.04 1.12 0.94 1.09 1.02 1.18 KLEN: 10-kernel length; KWID: 10-kernel width; KWEI: 100-kernel weight; KLEN/KWID: kernel length/kernel width; HZS: Huangzaosi.

表1 5个不同环境下亲本间籽粒性状的表型差异 Table 1 Difference of the kernel related traits between parents under five environments

表2
Table 2
表2(Table 2)

表2 RIL群体和BC3F1群体中籽粒性状的表现 Table 2 Performance of kernel related traits in RIL and BC3F1 populations 群体 Population 性状 Trait 变异范围 Range 均值 Mean 标准差 SD 偏度 Skewness 峰度 Kurtosis RIL 粒长KLEN (cm) 7.67-10.26 9.03 0.561 -0.035 -0.473 粒宽KWID (cm) 7.28-10.11 8.63 0.502 -0.088 -0.023 百粒重KWEI (g) 14.82-31.52 23.76 3.110 -0.070 -0.184 粒长/粒宽KLEN/KWID 0.85-1.22 1.05 0.067 0.091 0.033 BC3F1 粒长KLEN (cm) 9.00-11.70 9.89 0.535 0.563 0.277 粒宽KWID (cm) 7.60-10.50 8.82 0.445 0.543 0.299 百粒重KWEI (g) 13.41-33.42 23.26 3.476 -0.102 0.387 粒长/粒宽KLEN/KWID 1.05-1.32 1.12 0.427 1.091 2.428 KLEN: 10-kernel length; KWID: 10-kernel width; KWEI: 100-kernel weight; KLEN/KWID: kernel length/kernel width.

表2 RIL群体和BC3F1群体中籽粒性状的表现 Table 2 Performance of kernel related traits in RIL and BC3F1 populations

图1
Fig. 1

	Figure OptionViewDownloadNew Window
	图1 RIL群体中粒长、粒宽、百粒重和粒长/粒宽表型频率分布Fig. 1 Frequency distributions of the KLEN, KWID, KWEI, KLEN/WEID in the RIL population KLEN: 10-kernel length; KWID: 10-kernel width; KWEI: 100-kernel weight; KLEN/KWID: kernel length/kernel width.

表3
Table 3
表3(Table 3)

表3 籽粒相关性状在5个环境下的联合方差分析 Table 3 Analysis of variance (ANOVA) for kernel-related traits under five environments 性状 Trait 方差来源 Source of variation 均方 MS F值 F-value 广义遗传力 H2 粒长KLEN 环境 Environment 75.724 243.649* * 0.861 基因型 Genotype 2.497 8.034* * 环境× 基因型 Environment× Genotype 0.726 2.335* * 粒宽KWID 环境 Environment 40.277 245.109* * 0.860 基因型 Genotype 1.775 10.802* * 环境× 基因型 Environment× Genotype 0.369 2.248* * 百粒重KWEI 环境 Environment 4573.754 683.667* * 0.855 基因型 Genotype 71.835 10.738* * 环境× 基因型 Environment× Genotype 21.766 3.254* * 粒长/粒宽KLEN/KWID 环境 Environment 0.312 83.856* * 0.869 基因型 Genotype 0.036 9.567* * 环境× 基因型 Environment× Genotype 0.008 2.218* * * * indicates significant level at P< 0.01. KLEN: 10-kernel length; KWID: 10-kernel width; KWEI: 100-kernel weight; KLEN/KWID: kernel length/kernel width. * * 表示P< 0.01, 差异极显著。

表3 籽粒相关性状在5个环境下的联合方差分析 Table 3 Analysis of variance (ANOVA) for kernel-related traits under five environments

2.2 高密度遗传图谱构建利用GBS技术, 检测RIL群体的基因型, 共得到78 507个高质量SNP。基于这些SNP, 构建RIL群体的重组区块图谱共得到1262个Bin标记。利用这些Bin构建的遗传连锁图谱总长度为1524.53 cM, 标记间最大距离为24.14 cM, 最小距离0.38 cM, 标记间的平均距离为1.21 cM。
2.3 RIL群体中的QTL定位分析在单环境作图中, 共检测到4个籽粒相关性状30个QTL, 分布在第1、第2、第4、第6、第7、第8、第10染色体上, 单个QTL的表型贡献率变化范围在6.28%~27.63%之间(表4)。利用多环境平均值进行QTL定位时, 共检测到所有性状11个QTL, 分布在第1、第2、第4、第6、第9、第10染色体上, 单个QTL解释的表型变异在6.01%~26.79%之间(表5)。单个环境和多个环境检测到的QTL在遗传图谱上的位置如图2所示。
2.3.1 粒长QTL 在单环境条件下, 共定位到7个粒长QTL, 分布在第1、第4、第7、第10染色体上。利用多环境平均值进行QTL定位时, 在第1染色体物理位置210~212 Mb处检测到一个粒长主效QTL并命名为qklen1, 其LOD值为7.15, 表型贡献率达22.60% (表5)。该位点可在2009北京、2009新疆、2010新疆各单环境下被检测到, 相应表型贡献率分别为23.35%、21.00%、10.63% (表4)。
表4
Table 4
表4(Table 4)

表4 单环境下籽粒相关性状QTL分析 Table 4 QTL analysis of kernel related traits under single environment 性状 Trait 染色体 Chr. 标记区间 Marker interval LOD值 LOD score 贡献率 PVE (%) 加性效应 Additive 环境 Environment 粒长KLEN 1 M1210096805-M1211894566 7.537 23.354 0.334 2009北京 2009 Beijing 10 M10147002195-M10147354711 2.835 8.185 -0.196 2009北京 2009 Beijing 1 M1209234039-M1210096805 7.347 21.000 0.367 2009新疆 2009 Xinjiang 7 M7165815642-M7166162680 3.810 10.220 -0.256 2009新疆 2009 Xinjiang 10 M10122808985-M10123448730 2.586 6.765 0.208 2009新疆 2009 Xinjiang 4 M4039558232-M4039794707 2.749 12.363 0.254 2010河南 2010 Henan 1 M1210096805-M1211894566 3.297 10.628 0.244 2010新疆 2010 Xinjiang 粒宽KWID 6 M6112810876-M6113109410 4.188 12.889 -0.215 2009北京 2009 Beijing 10 M10147467011-M10147845941 4.373 12.762 -0.213 2009北京 2009 Beijing 6 M6098574401-M6102160956 3.837 11.342 -0.216 2009新疆 2009 Xinjiang 10 M10147002195-M10147354711 3.212 9.675 -0.194 2009新疆 2009 Xinjiang 10 M10147002195-M10147354711 2.892 13.083 -0.230 2010河南 2010 Henan 10 M10142715899-M10143442026 3.070 10.785 -0.185 2010新疆 2010 Xinjiang 百粒重KWEI 1 M1258796378-M1259239220 4.051 10.564 1.197 2009北京2009 Beijing 4 M4239687946-M4240246524 2.559 6.574 0.954 2009北京2009 Beijing 7 M7156739212-M7158128129 2.533 6.367 0.941 2009北京 2009 Beijing 10 M10147002195-M10147354711 4.240 10.970 -1.220 2009北京 2009 Beijing 4 M4233993423-M4234603338 4.323 15.056 1.659 2009河南 2009 Henan 8 M8172252275-M8172457715 4.031 14.115 -1.637 2009河南 2009 Henan 10 M10147002195-M10147354711 2.711 9.879 -1.287 2009新疆 2009 Xinjiang 粒长/粒宽 KLEN/KWID 1 M1207528316-M1207726397 5.947 19.574 0.038 2009北京 2009 Beijing 1 M1207528316-M1207726397 5.085 14.750 0.035 2009河南 2009 Henan 2 M2191772105-M2193504116 4.529 12.521 -0.032 2009河南 2009 Henan 1 M1202169469-M1202644766 9.582 27.633 0.047 2009新疆 2009 Xinjiang 2 M2185112743-M2185399317 2.532 6.283 -0.023 2009新疆 2009 Xinjiang 6 M6102160956-M6102459999 3.973 10.408 0.029 2009新疆 2009 Xinjiang 1 M1196195341-M1198194592 3.074 12.035 0.029 2010河南 2010 Henan 2 M2154493476-M2159081032 4.647 19.355 -0.038 2010河南 2010 Henan 1 M1206305038-M1207528316 4.400 14.258 0.029 2010新疆 2010 Xinjiang 8 M8140007454-M8141932923 2.697 8.285 -0.022 2010新疆 2010 Xinjiang KLEN: 10-kernel length; KWID: 10-kernel width; KWEI: 100-kernel weight; KLEN/KWID: kernel length/kernel width; PVE: phenotypic variation explained.

表4 单环境下籽粒相关性状QTL分析 Table 4 QTL analysis of kernel related traits under single environment

表5
Table 5
表5(Table 5)

表5 多环境条件下籽粒相关性状表型平均值QTL定位分析 Table 5 QTL analysis of kernel related traits under joint-environments using phenotypic means 性状 Trait QTL 标记区间 Marker interval LOD值 LOD score 贡献率 PVE (%) 加性效应 Additive 环境数 Environment 粒长KLEN qklen1 M1210096805-M1211894566 7.1510 22.5986 0.2690 3 粒宽KWID qkwid4 M4237508696-M4237523196 2.6397 6.0111 0.1225 0 qkwid6 M6104672903-M6105080338 5.0018 12.2122 -0.1795 1 qkwid9 M9020404359-M9025815280 3.4446 8.0292 -0.1431 0 qkwid10 M10143442026-M10143738783 6.6194 15.9065 -0.1999 1 百粒重KWEI qkwei2 M2014268014-M2014614195 2.5272 6.2901 -0.7886 0 qkwei4-1 M4187601372-M4187992567 4.8716 12.7470 -1.1073 0 qkwei4-2 M4233993423-M4234603338 6.1193 16.0548 1.2394 1 qkwei10 M10147002195-M10147354711 3.6302 9.3084 -0.9431 2 粒长/粒宽 KLEN/KWID qklw1 M1207528316-M1207726397 9.1968 26.7880 0.0346 3 qklw2 M2152921409-M2154493476 5.2433 13.9095 -0.0257 1 KLEN: 10-kernel length; KWID: 10-kernel width; KWEI: 100-kernel weight; KLEN/KWID: kernel length/kernel width; PVE: phenotypic variation explained.

表5 多环境条件下籽粒相关性状表型平均值QTL定位分析 Table 5 QTL analysis of kernel related traits under joint-environments using phenotypic means

2.3.2 粒宽QTL 在单环境条件下, 共定位到6个粒宽QTL, 分布在第6、第10染色体上。其中在第10染色体物理位置147 Mb处存在一粒宽QTL, 在2009北京、2009新疆、2010河南各单环境条件下均能被检测到(表4)。以多环境平均值进行QTL定位时, 在第4、第6、第9和第10染色体上各检测到1个粒宽QTL, 分别命名为qkwid4、qkwid6、qkwid9、qkwid10, 其表型贡献率在6.01%~15.91%之间(表5)。
2.3.3 百粒重QTL 在各个单环境下共检测到7个粒重QTL, 分布在第1、第4、第7、第8、第10染色体上。以多环境平均值进行QTL定位时, 共检测到4个QTL。其中qkwei10表型贡献率为9.31%, LOD值为3.63, 物理位置位于第10染色体147 Mb左右(表5)。同时, 该位点在2009北京、2009新疆2个单环境下均能检测到(表4)。
2.3.4 粒长/粒宽QTL 在各个单环境下共检测到10个粒长/粒宽QTL, 分布在第1、第2、第6、第8染色体上。以多环境平均值进行QTL定位时, 共检测到2个QTL, 分布在第1、第2染色体上。其中, qklw1物理位置在第1染色体207~208 Mb处, 其LOD值为9.20, 表型贡献率可达26.79% (表5)。同时该主效位点在2009北京、2009河南、2010新疆3个单环境下均能检测到(表4)。
2.4 精细定位的候选区段RIL群体的籽粒相关性状QTL定位结果表明, qklen1和qklw1不仅是两个环境钝感的主效QTL, 而且其物理位置在第1染色体上相邻(图3)。因此, 这两个QTL被认为是精细定位的重要候选区段。针对该区间, 利用5个InDel标记(表6), 采用成组法t检验, 对BC₃F₁群体粒长、粒长/粒宽进行单标记分析。其中, 针对粒长的单标记分析结果表明, 主效区段内5个InDel标记对应的P值均达到极显著水平(P< 0.01)。另外, 在粒长/粒宽的单标记分析中, 标记ML162、ML178和ML93对应的P值同样达到极显著水平(表6)。以上表明, 主效位点qklen1和qklw1在BC₃F₁群体中同样具有显著的遗传效应。
图2
Fig. 2

Figure OptionViewDownloadNew Window

图2 籽粒相关性状的QTL在遗传图谱上的分布 彩色图案为单个环境监测到的QTL, 黑色图案为多环境联合定位到的QTL; 2009BJKLEN: 2009年北京粒长; 2009XJKLEN: 2009年新疆粒长; 2010HNKLEN: 2010年河南粒长; 2010XJKLEN: 2010年新疆粒长; 2009BJKWID: 2009年北京粒宽; 2009XJKWID: 2009年新疆粒宽; 2010HNKWID: 2010年河南粒宽; 2010XJKWID: 2010年新疆粒宽; 2009BJKWEI: 2009年北京粒重; 2009HNKWEI: 2009年河南粒重; 2009XJKWEI: 2009年新疆粒重; 2009BJKLEN/KWID: 2009年北京粒长/粒宽; 2009HNKLEN/KWID: 2009年河南粒长/粒宽; 2009XJKLEN/KWID: 2009年新疆粒长/粒宽; 2010HNKLEN/KWID: 2010年河南粒长/粒宽; 2010XJKLEN/KWID: 2010年新疆粒长/粒宽。Fig. 2 Distributions of identified QTL for kernel related traits on genetic linkage maps The color symbols stand for QTLs detected in a single environment, and the black symbols stand for QTLs detected across five environments; 2009BJKLEN: 10-kernel length in Beijing in 2009; 2009XJKLEN: 10-kernel length in Xinjiang in 2009; 2010HNKLEN: 10-kernel length in Henan in 2010; 2010XJKLEN: 10-kernel length in Xinjiang in 2010; 2009BJKWID: 10-kernel width in Beijing in 2009; 2009XJKWID 10-kernel width in Xinjiang in 2009; 2010HNKWID: 10-kernel width in Henan in 2010; 2010XJKWID: 10-kernel width in Xinjiang in 2010; 2009BJKWEI: 100-kernel weight in Beijing in 2009; 2009HNKWEI: 100-kernel weight in Henan in 2009; 2009XJKWEI: 100-kernel weight in Xinjiang in 2009; 2009BJKLEN/KWID: kernel length/kernel width in Beijing in 2009; 2009HNKLEN/KWID: kernel length/kernel width in Henan in 2009; 2009XJKLEN/KWID: kernel length/kernel width in Xinjiang in 2009; 2010HNKLEN/KWID: kernel length/kernel width in Henan in 2010; 2010XJKLEN/KWID: kernel length/kernel width in Xinjiang in 2010; KLEN: 10-kernel length; KWID: 10-kernel width; KWEI: 100-kernel weight; KLEN/KWID: kernel length/kernel width.

表6
Table 6
表6(Table 6)

表6 InDel标记引物及相应P值 Table 6 Primers of InDel markers and the corresponding P-value InDel标记 InDel marker 物理位置 Physical position (Mb, RefGen_v2) 正向引物 Forward primer (5′ -3′ ) 反向引物 Reverse primer (5′ -3′ ) P值 P-value 粒长 KLEN 粒长/粒宽 KLEN/KWID ML74 207.17 CTGCAAGACAGCAGCAAGACT ATTGTGCGAATTTGGGAAGAT 0.0053 0.0746 ML194 207.61 TTCATCACTGCCTCAACCAC GTTGGGAACTTCGGAACAAT 0.0053 0.0567 ML162 209.29 TTGCGTAGATCGTCGGTGTC AGTGGGATTGGATGGGTTTG 0.0077 0.0055 ML178 211.33 CAGCAGCGGAAGAGGAACCA CGGAGCAAGAAAGGCAAGCA 0.0023 0.0018 ML93 214.24 GGAGCACCAATAGGCAGTCA TCTCGCAGTAGCCGACAAAG 0.0036 0.0042 KLEN: 10-kernel length; KLEN/KWID: kernel length/kernel width.

表6 InDel标记引物及相应P值 Table 6 Primers of InDel markers and the corresponding P-value

图3
Fig. 3

	Figure OptionViewDownloadNew Window
	图3 粒长和粒长/粒宽相关主效QTL在第1染色体上的分布Fig. 3 Distributions of KLEN and KLEN/KWID related major QTL on chromosome 1 KLEN: 10-kernel length; KLEN/KWID: kernel length/kernel width.

3 讨论玉米是世界上最重要的粮食作物之一, 高产仍是其育种的首要目标^[23]。而籽粒产量性状遗传基础复杂, 是由相互关联的多基因控制的一系列生理生化过程的最终体现^[24]。产量性状的复杂性导致直接解析其遗传基础的一定难度。但是, 与产量显著相关的玉米籽粒性状, 例如, 粒长、粒宽、粒长/粒宽和百粒重等, 遗传力较高, 表型评价更加准确, 可作为产量性状的重要构成因子。在此条件下, 加强产量构成因子的遗传基础的解析和功能分子标记开发对于玉米高产育种具有重要意义。
本研究利用亲本黄早四和Mo17构建的RIL群体, 借助高密度遗传图谱, 对多环境下的籽粒相关性状进行QTL定位分析。结果显示, 一些QTL在不同环境下定位在基因组的同一区域, 如粒长主效QTL qklen1、粒长/粒宽QTL qklw1在3个不同环境条件下都被检测到。这说明这些QTL在不同的环境条件下都能稳定地表达, 是进一步精细定位和挖掘主效基因的候选项区段。另外, 不同性状QTL被定位于基因组中同一位点或相近区域, Tuberosa等^[25]认为“ 一因多效” 很可能是导致该现象的主要原因; 控制不同性状的基因存在紧密连锁关系也可引起这种富集现象的发生。
通过遗传连锁作图, 本文构建了高密度遗传连锁图谱, 图谱共包含1262个Bin标记, 长度为1524.53 cM, 标记间的平均距离仅为1.21 cM。与前人研究相比, QTL作图精度显著提高。Li等^[26]利用179个分子标记对本文中的同一RIL群体的籽粒性状进行QTL分析, 在多环境联合条件定位到的Mqkwei4、Mqklen1、Mqkwid6、Mqkwid10分别和本研究中定位到的qkwei4-2、qklen1、qkwid6、qkwid10位于染色体同一区域。而本研究中定位到的qkwid4、qkwid9、qkwei2、qkwei4-1、qkwei10在Li等^[26]研究中并未检测到, 分析其原因, 在低密度遗传图谱中, 较大的区间标记将导致部分重组无法被有效检测, 这将引起表型与基因型关联之间出现偏差, 进而导致一些主效QTL检测的缺失。另外, Li等^[26]在RIL群体中定位到的一些QTL置信区间一般为10~30 cM, 而本文检测到的QTL位点置信区间一般为0.8~3.0 cM, 作图精度比低密度图谱提高了10倍左右。例如, 本文中粒长QTL位点qklen1置信区间为1.2 cM, 对应在第1染色体物理位置210~212 Mb处。Li等^[26]也在此区域定位到一粒长主效QTL, 但置信区间达到了36.0 cM。进一步显示了高密度遗传图谱在数量性状遗传定位中的巨大优势。
本研究发现, 在第1染色体上临近的2个主效QTL位点, qklen1 (210~212 Mb)和qklw1 (207~208 Mb), 在多环境下都能稳定表达。Li等^[26]、Austin等^[27]、张向歌等^[28]、张伟强等^[29]也在相同或邻近该区域内检测到粒长QTL。另外, 在第10染色体定位到的粒宽QTL qkwid10 (143.4~143.7 Mb)和粒重QTLqkwei10 (147.0~147.3 Mb)同样在多环境下稳定表达。Peng等^[9]利用F_2:3群体在联合环境下将粒宽、粒重QTL定位在第10染色体物理位置143~147 Mb (bnlg1677~umc2172)之间。以上表明, 这两处区域遗传效应稳定, 是控制籽粒性状的重要位点, 可作为精细定位的重要候选区段。利用黄早四和Mo17构建的BC₃F₁群体对qklen1和qklw1遗传效应验证表明, 在BC₃F₁群体中qklen1和qklw1同样具有显著的遗传效应。由于BC₃F₁群体在207~212 Mb区段内发生重组交换的个体较少, 且2个QTL紧密连锁, 所以还不能判断该区段具有“ 一因多效” 的遗传效应, 还是两个独立基因在起作用, 需进一步自交利用大规模分离群体精细定位来验证。
4 结论控制不同性状的QTL可定位在相同或邻近的染色体区域; 部分主效位点可在多个环境下被检测到。其中, 粒长主效位点(qklen1)和粒长/粒宽主效位点(qklw1)在3个环境下均能被检测到, 且定位于染色体的相邻位置。qklen1和qklw1在BC₃F₁群体中同样具有显著的遗传学效应。以上结果为玉米籽粒性状主效QTL的精细定位和图位克隆提供了重要的研究基础。
The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。


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

[1]	Ray D K, Mueller N D, West P C, Foley J A. Yield trends are insufficient to double global crop production by 2050. PLoS One, 2013, 8(6): e66428[本文引用:1]
[2]	Sundaresan V. Control of seed size in plants. Proc Natl Acad Sci USA, 2005, 12: 17887-17888[本文引用:1]
[3]	Yang X, Ma H, Zhang P, Yan J, Guo Y, Song T, Li J. Characterization of QTL for oil content in maize kernel. Theor Appl Genet, 2012, 125: 1169-1179[本文引用:1]
[4]	Ribaut J M, Jiang C, Gonzalez D, Edmeades G O, Hoisington D A. Identification of quantitative trait loci under drought conditions in tropical maize: 2. Yield components and marker assisted selection strategies. Theor Appl Genet, 1997, 94: 887-896[本文引用:1]
[5]	Wen Y X, Zhu J. Multivariable conditional analysis for complex trait and its components. Acta Genet Sin, 2005, 32: 289-296[本文引用:1]
[6]	Borras L, Otegui M E. Maize kernel weight response to post flowering source-sink ratio. Crop Sci, 2001, 41: 1816-1822[本文引用:1]
[7]	Doebley J F, Gaut B S, Smith B D. The molecular genetics of crop domestication. Cell, 2006, 127: 1309-1321[本文引用:1]
[8]	李永祥, 王阳, 石云素, 宋燕春, 王天宇, 黎裕. 玉米籽粒构型与产量性状的关系及QTL作图. 中国农业科学, 2009, 42: 408-418 Li Y X, Wang Y, Shi Y S, Song Y C, Wang T Y, Li Y. Correlation analysis and QTL mapping for traits of kernel structure and yield components in maize. Sci Agric Sin, 2009, 42: 408-418 (in Chinese with English abstract)[本文引用:1]
[9]	Peng B, Li Y X, Wang Y, Liu C, Liu Z Z, Tan W W, Zhang Y, Wang D, Shi Y S, Sun B C, Song Y C, Wang T Y, Li Y. QTL analysis for yield components and kernel related traits in maize across multi-environments. Theor Appl Genet, 2011, 122: 1305-1320[本文引用:2]
[10]	黎裕, 王天宇, 石云素, 宋燕春. 基因组学方法在玉米种质资源研究中的应用. 植物遗传资源学报, 2003, 4: 256-260 Li Y, Wang T Y, Shi Y S, Song Y C. Applications of genomics approaches in studies on maize germplasm. J Plant Genet Resour, 2003, 4: 256-260 (in Chinese with English abstract)[本文引用:1]
[11]	Liu Y, Wang L W, Sun C L, Zhang Z X, Zheng Y L, Qiu F Z. Genetic analysis and major QTL detection for maize kernel size and weight in multi-environments. Theor Appl Genet, 2014, 127: 1019-1037[本文引用:1]
[12]	Zhang Z H, Liu Z H, Hu Y M, Li W H, Fu Z Y, Ding D, Li H C, Qiao M M, Tang J H. QTL analysis of kernel-related traits in maize using an immortalized F2 population. PLoS One, 2014, 9(2): e89645[本文引用:1]
[13]	Nikolic A, Andelkovic V, Dodig D, Drinic M S, Kravic N, Micic I D. Identification of QTLs for drought tolerance in maize: II. Yield and yield components. Genetica, 2013, 45: 341-350[本文引用:1]
[14]	Colasuonno P, Gadaleta A, Giancaspro A, Nigro D, Giove S, Incerti O, Mangini G, Signorile A, Simeone R, Blanco A. Development of a high-density SNP-based linkage map and detection of yellow pigment content QTLs in durum wheat. Mol Breed, 2014, 34: 1563-1578[本文引用:1]
[15]	Guo T T, Yang N, Tong H, Pan Q C, Yang X H, Tang J H, Wang J K, Li J S, Yan J B. Genetic basis of grain yield heterosis in an “immortalized F2” maize population. Theor Appl Genet, 2014, 127: 2149-2158[本文引用:1]
[16]	Zou G H, Zhai G W, Feng Q, Yan S, Wang A H, Zhao Q, Shao J F, Zhang Z P, Zou J Q, Han B, Tao Y Z. Identification of QTLs for eight agronomically important traits using an ultra-high-density map based on SNPs generated from high-throughput sequencing in sorghum under contrasting photoperiods. J Exp Bot, 2012, 63: 5451-5462[本文引用:1]
[17]	Chen D H, Ronald P. A rapid DNA minipreparation method suitable for AFLP and other PCR applications. Plant Mol Biol Rep, 1999, 17: 53-57[本文引用:2]
[18]	Elshire R J, Glaubitz J C, Sun Q, Poland J A, Kawamoto K, Buckler E S, Mitchell S E. A robust, simple genotyping- by-sequencing (GBS) approach for high diversity species. PLoS One, 2011, 6(5): e19379[本文引用:2]
[19]	Li C H, Li Y X, Shi Y S, Song Y C, Zhang D F, Buckler E S, Zhang Z W, Wang T Y, Li Y. Genetic control of the leaf angle and leaf orientation value as revealed by ultra-high density maps in three connected maize populations. PLoS One, 2015, 10(3): e0121624[本文引用:1]
[20]	Wan X Y, Wan J M, Jiang L, Wang J K, Zhai H Q, Weng J F, Wang H L, Lei C L, Wang J L, Zhang X. QTL analysis for rice grain length and fine mapping of an identified QTL with stable and major effects. Theor Appl Genet, 2006, 112: 1258-1270[本文引用:1]
[21]	Wang J K, Wan X Y, Crossa J, Crouch J, Weng J, Zhai H Q, Wan J M. QTL mapping of grain length in rice (Oryza sativa L. ) using chromosome segment substitution lines. Genet Res, 2006, 88: 93-104[本文引用:1]
[22]	Li H H, Ye G Y, Wang J K. A modified algorithm for the improvement of composite interval mapping. Genetics, 2007, 175: 361-374[本文引用:1]
[23]	Gupta P K, Rustgi S, Kumar N. Genetic and molecular basis of grain size and its relevance to grain productivity in higher plants. Genome, 2006, 49: 565-571[本文引用:1]
[24]	Peleg Z, Fahima T, Krugman T, Abbo S, Yakir D, Korol A B, Saranga Y. Genomic dissection of drought resistance in durum wheat × wild emmer wheat recombinant inbreed line population. Plant Cell Environ, 2009, 32: 758-779[本文引用:1]
[25]	Tuberosa R, Salvi S, Sanguineti M C, Land i P, Maccaferri M, Conti S. Mapping QTL regulating morpho-physiological traits and yield: case studies, shortcomings and perspectives in drought- stressed maize. Ann Bot, 2002, 89: 941-963[本文引用:1]
[26]	Li C H, Li Y X, Sun B C, Peng B, Liu C, Liu Z Z, Yang Z Z, Li Q C, Tan W W, Zhang Y, Wang D, Shi Y S, Song Y C, Wang T Y, Li Y. Quantitative trait loci mapping for yield components and kernel-related traits in multiple connected RIL populations in maize. Euphytica, 2013, 193: 303-316[本文引用:5]
[27]	Austin D F, Lee M. Comparative mapping in F2: 3 and F6: 7 generations of quantitative trait loci for grain yield and yield components in maize. Theor Appl Genet, 1996, 92: 817-826[本文引用:1]
[28]	张向歌, 王彬, 袁亮, 张晓祥, 时夏, 赵晓锋, 汤继华. 基于单片段代换系玉米子粒性状的QTL定位. 玉米科学, 2013, 21(6): 35-40 Zhang X G, Wang B, Yuan L, Zhang X X, Shi X, Zhao X F, Tang J H. QTL mapping for kernel related traits basing on the single segment substitution lines in maize. J Maize Sci, 2013, 21(6): 35-40 (in Chinese with English abstract)[本文引用:1]
[29]	张伟强, 库丽霞, 张君, 韩赞平, 陈彦惠. 玉米出籽率、籽粒深度和百粒重的QTL分析. 作物学报, 2013, 39: 455-463 Zhang W Q, Ku L X, Zhang J, Han Z P, Chen Y H. QTL analysis of kernel ratio, kernel depth, and 100-kernel weight in maize (Zea mays L. ). Acta Agron Sin, 2013, 39: 455-463 (in Chinese with English abstract)[本文引用:1]