删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

玉米穗行数主效位点 qKRN5.04精细定位与遗传效应解析

本站小编 Free考研考试/2021-12-26

白娜, 李永祥*, 焦付超, 陈林, 李春辉, 张登峰, 宋燕春, 王天宇, 黎裕, 石云素*
中国农业科学院作物科学研究所, 北京 100081
* 通讯作者(Corresponding authors): 李永祥, E-mail:liyongxiang@caas.cn; 石云素, E-mail:shiyunsu@caas.cn 第一作者联系方式: E-mail:baina0919@sina.cn
收稿日期:2016-06-20 基金:本研究由国家公益性行业(农业)科研专项, 作物种质资源保护与利用专项(201303007, 2015NWB030-04), 国家科技支撑计划项目 (2013BAD01B02-3), 国家自然科学基金项目(91335206)和中国农业科学院科技创新工程项目资助

摘要穗行数是影响玉米产量的重要因素之一, 其遗传机制解析和关键基因精细定位对开展分子育种具有重要的意义。本研究以穗行数仅有4行的“四路糯系选”和多穗行自交系“农531” (18~22行)为亲本, 构建了高代回交群体和次级定位群体(四路糯选系为供体亲本, 农531为轮回亲本)。通过对不同类型试验群体的多环境表型鉴定和基因型鉴定, 利用完备区间作图法(ICIM)进行穗行数主效QTL定位分析, 将穗行数主效位点 qKRN5.04定位到第5染色体136.3~140.0 Mb的区间之内; 遗传效应分析发现, 该位点在不同环境条件下最大可解释的表型变异为21.76%, 效应值为0.80~1.76行。通过次级分离群体重组事件分析可将其进一步定位到~300 kb区间内。本研究结果不仅为分子标记辅助选择提供了实用的InDel标记, 而且为玉米穗行数主效位点 qKRN5.04的图位克隆和候选基因挖掘奠定了重要的基础。

关键词:玉米; 穗行数; QTL; 精细定位
Fine Mapping and Genetic Effect Analysis of qKRN5.04, a Major QTL Associated with Kernel Row Number in Maize
BAI Na, LI Yong-Xiang*, JIAO Fu-Chao, CHEN Lin, LI Chun-Hui, ZHANG Deng-Feng, SONG Yan-Chun, WANG Tian-Yu, LI Yu, SHI Yun-Su*
Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing 100081, China
Fund:This study was supported by the Special Fund for Agro-scientific Research in the Public Interest, the Special Fund for Protection and Utilization of Crop Germplasm Resources (201303007, 2015NWB030-04), the National Key Technology Support Program of China (2013BAD01B02-3), the National Natural Science Foundation of China (91335206), and the Agricultural Science and Technology Innovation Program of Chinese Academy of Agricultural Sciences
AbstractKernel row number (KRN) is one of the important factors of maize yield. The genetic basis dissection and fine mapping of crucial gene would be greatly beneficial to maize molecular breeding. In present study, series of advanced backcross population and secondary segregation population were developed from a backcross of the derived line of Silunuo (four kernel rows, as the donor parent) and Nong 531 (18-22 kernel rows, as the recurrent parent). The QTL mapping for KRN was conducted by the inclusive complete interval mapping (ICIM) method in multiple environments. And a major KRN QTL, qKRN5.04, was mapped to the interval of 136.3-140.0 Mb on chromosome 5, with the largest phenotypic variation of 21.76% and the effect value of 0.80-1.76 row. Furthermore, according to recombinant analysis of secondary population, qKRN5.04 was fine mapped to the region of about 300 kb, which provided both practical InDel markers for marker-assisted selection and sufficient supports for the map-based cloning and candidate gene mining of the target locus.

Keyword:Maize; Kernel row number (KRN); Quantitative trait locus (QTL); Fine mapping
Show Figures
Show Figures










玉米产量是复杂的数量性状, 受微效多基因控制, 其主要构成要素包括单位面积穗数、穗粒数和百粒重, 而穗粒数由穗行数和行粒数组成。玉米穗行数具有较高的遗传力, 且与产量密切相关, 是重要的产量相关性状之一[1]。因此, 深入研究玉米穗行数关键调控基因及遗传效应, 对促进玉米高产育种理论与应用的发展具有重要的意义[2]
近年来, 在玉米穗行数遗传基础研究方面, 借助QTL作图, 已开展了大量工作。已报道的研究结果表明, 玉米穗行数遗传区域主要分布在bin2.03-2.04、bin4.05、bin4.08-4.09、bin5.04、bin7.02-7.03、bin8.05、bin9.03、bin10.03-10.04等染色体区段[3, 4, 5]。基因挖掘方面, 利用人工创制的突变体, 确定了一些重要的玉米穗行数调控基因。例如, 当barren inflorescence2 (bif2)、barren stalk1 (ba1)、sparse inflorescence1 (spi1)、Barren inflorescence1 (Bif1)、barren stalk1 (ba1)、Suppressor of sessile spikelets1 (Sos1)等基因突变后, 会引起小穗分枝减少, 从而导致穗行数减少[6, 7, 8]。而另一些基因的突变则会引起穗行数的增加, 这些基因包括FASCIATEDEAR2(FEA2)、td1、ramosa1 (ra1)、ramosa2 (ra2)、branched silkless1 (Bd1)、rgo1Zfl1Zfl2[9, 10, 11, 12, 13, 14]。但人工突变体与自然群体中存在的变异类型差别较大, 难以直接应用于育种实践。在此条件下, 充分挖掘自然群体中的既有变异类型, 基于图位克隆方法挖掘穗行数相关基因, 对解析玉米穗行数遗传调控机制和开展分子标记辅助选择具有重要参考价值。
基于分子标记辅助选择, 利用高代回交群体构建导入系进行QTL分析, 已经成为研究数量性状遗传机制的有效方法。导入系群体遗传背景基本一致, 极大程度地降低了遗传背景对表型的干扰, 从而有效提高目标表型性状鉴定的准确性; 同时, 由于导入系含有的外源片段数较少, 任何表型变化都可认为由导入片段而引起[15], 因此, 可实现主效位点遗传效应的准确估计。胡立宗等[16]基于玉米BC2F2群体, 采用复合区间作图法定位到了3个穗行数QTL; 李卫华等[17]利用单片段代换系群体结合显著性检验的方法检测到3个穗行数QTL; 齐欢欢等[18]利用导入系群体对玉米产量及产量相关性状定位, 两个环境下一致性地检测到2个穗行数QTL; Li等[19]利用130个染色体片段代换系(CSSLs)共检测到11个穗行数QTL; Tian等[20]利用BC5F2:3把一个穗行数的主效QTL qKRN10定位到标记ssr1430和umc1077之间。截至目前, 虽然利用自然变异的分离群体定位到一些穗行数主效QTL, 但真正实现克隆的穗行数主效位点仍然很少, 仅Liu等[21]利用近等基因系, 将一个穗行数主效位点KRN4定位到了1.2 kb的区间内, 并成功验证水稻的同源基因Unbranched3 (UB3)是KRN4的关键调控基因。
本研究基础材料为四路糯选系和农531。四路糯是起源于云南勐海、西双版纳等地的一个原始地方品种, 具有多分蘖、多果穗、少穗行、半有稃等一系列玉米原始性状特征; 农531是一个多穗行数玉米自交系, 最高可达22行。以上2份材料的穗行数性状差异极其显著。在本实验室前期的研究工作中, 利用穗行数只有4行的四路糯选系和多穗行亲本农531为亲本, 构建了F2:3家系分离群体, 借助多环境表型数据和基因型数据, 在基因组范围内共定位到5个穗行数主效QTL。其中, 位于第5染色体上的主效位点qKRN5.04单环境下最大可以解释13.55%的表型变异, 联合环境分析可以解释9.83%的表型变异[22], 且该位点位于bin5.04穗行数QTL富集区段之内。因此, 本研究以该位点为关键目标区段, 借助分子标记辅助选择, 通过构建一系列高代回交群体和大规模次级分离群体, 进行穗行数主效位点qKRN5.04的精细定位和遗传效应分析, 以期为该位点的图位克隆和候选基因挖掘提供参考和支持。
1 材料与方法1.1 试验材料以四路糯选系为供体亲本, 农531为轮回亲本, 从BC1代开始, 借助分子标记辅助选择, 筛选仅包含qKRN5.04导入的单株(去掉含有其他穗行数主效位点的基因型), 并通过连续的回交和自交, 获得了主效位点定位的相关试验群体(图1)。
图1
Fig. 1
Figure OptionViewDownloadNew Window
图1 qKRN5.04主效位点定位试验材料Fig. 1 Plant materials used for the mapping of qKRN5.04


1.2 田间试验2014年夏, 在北京昌平(BJ, 116.20° N, 40.22° E)和河北石家庄(SJZ, 114.26° N, 38.03° E)同时种植BC3F2和BC4F2群体, 进行单株分离群体的表型和基因型鉴定, 并全部自交以获得BC3F3和BC4F3群体的种子。当年冬天, 在海南三亚(HN, 109.30° N, 18.15° E), 根据北京和河北石家庄QTL定位结果, 一方面, 挑选主效区段纯合型的BC3F2和BC4F2自交后代的种子, 进行基于BC3F2:3和BC4F2:3家系的表型评价, 同时种植目标区段杂合的BC3F3和BC4F3分离群体, 进行表型和基因型鉴定, 并全部自交以获得BC3F4和BC4F4分离群体; 另一方面, 种植目标区段杂合的BC5F2群体, 进行单株表型和基因型鉴定, 并全部自交, 用于主效位点的进一步验证和精细定位群体的构建; 2015年, 在北京昌平种植BC3F4、BC4F4和BC5F3群体, 以获得更多的交换单株, 用于QTL效应检测和重组事件分析, 完成主效位点qKRN5.04精细定位。单行种植试验群体, 1次重复, 行长3.0 m, 行距0.6 m, 每行定植12株。各试验点田间管理与当地大田管理相同, 成熟后收获, 以玉米穗轴中部穗行数为准调查表型。
1.3 分子标记筛选及基因型鉴定基因型鉴定的分子标记, 一部分来自MaizeGDB (http://www.maizegdb.org/)公布的SSR标记, 该部分标记主要用于单片段主效QTL (qKRN5.04)导入片段的分子标记辅助选择。同时, 为了提高作图精度, 实现主效位点精细定位, 利用Illumina HiSeq 2000测序平台, 进行亲本四路糯和农531深度测序(20X), 利用Primer5.0设计开发针对qKRN5.04遗传区段的多态性分子标记, 以增加目标区段标记密度。
在幼苗长到6片叶取单株新鲜叶片, 采用CTAB法[23]提取植物基因组DNA。利用筛选到的多态性标记, 通过PCR扩增、8.0%非变性聚丙烯酰胺凝胶电泳、0.1% AgNO3染色、甲醛和NaOH显色, 鉴定所有群体基因型。
1.4 QTL定位及验证利用QTL IciMapping V3.3完备区间作图法(ICIM, inclusive composite interval mapping), 进行QTL检测, 共进行1000次的排列测验(permutation test)以确定QTL的LOD阈值, 同时计算QTL的遗传效应及其对表型的贡献率[24]。利用SAS9.2进行数据的方差分析(Proc ANOVA)、单标记分析(Proc GLM)及多重比较(Proc GLM/DUNCAN, α = 0.05)。利用t测验(Proc TTEST)检验不同交换类型与受体亲本间的表型差异。按Stuber等[25]标准判定基因作用方式, 当显性效应值与加性效应值比值的绝对值为0~0.2时认为是加性(additive, A), 0.2~0.8为部分显性(partial dominance, PD), 0.8~1.2为显性(dominance, D), 大于1.2为超显性(over dominance, OD)。利用一般线性模型(GLM)检测主效位点与环境间的互作效应, 模型为Y = G + E + G × E + error, 其中, Y代表家系穗行平均数, G代表家系基因型, E代表环境, error代表残差。

2 结果与分析2.1 不同高代回交群体及亲本的玉米穗行数多环境表型变异四路糯选系在不同环境下穗行数稳定为4行, 农531在14BJ、14SJZ、14HN、15BJ、15SJZ、15XX (河南新乡) 6个环境下穗行数均值分别为18.28、18.15、18.38、18.60、17.07和18.02, 变化范围均值在17~19行之间; 双亲之间存在极显著差异。不同环境下, 不同世代回交群体的穗行数均值为15.38~ 16.81行, 群体内部穗行数表型多样性高(表1)。
表1
Table 1
表1(Table 1)
表1 不同回交群体及其亲本的穗行数 Table 1 Kernel row number for parents and backcross populations with qKRN5.04
环境
Environment
亲本ParentqKRN5.04的回交群体 Backcross populations of qKRN5.04
四路糯
Silunuo
农531
Nong 531
世代
Generation
群体大小
Number
均值± 标准差
Mean ± SD
变化范围
Range
2014北京Beijing418.28± 1.88BC3F270516.41± 2.09* * 10-24
2014石家庄Shijiazhuang418.15± 1.21BC3F272116.20± 2.06* * 10-22
2014三亚Sanya418.38± 1.85BC3F311716.81± 2.13* * 12-24
2015北京Beijing418.60± 2.02BC3F429616.74± 2.12* * 10-24
2015石家庄Shijiazhuang417.07± 1.58BC3F432815.38± 1.97* * 8-24
2015新乡Xinxiang418.02± 1.87BC3F343316.44± 2.27* * 8-24
2015北京Beijing418.20± 1.02BC5F364016.41± 2.09* * 10-24
* * 代表P< 0.01水平下差异极显著。* * Significantly different at P< 0.01.

表1 不同回交群体及其亲本的穗行数 Table 1 Kernel row number for parents and backcross populations with qKRN5.04

2.2 多态性标记筛选依据穗行数主效位点qKRN5.04在染色体上的物理位置, 筛选多态性引物并开发标记。首先, 在MaizeGDB上查找SSR引物序列, 筛选多态性引物; 同时参考亲本四路糯和农531全基因组重测序结果, 设计InDel扩增引物, 共计230对, 其中筛选到24对引物扩增稳定、条带清晰, 用于高代回交群体的基因型鉴定(部分多态性引物见表2)。
表2
Table 2
表2(Table 2)
表2 目标区段内开发的多态性InDel标记 Table 2 Polymorphic InDel markers within the target QTL
标记
InDel
marker
物理位置
Physical position
(Mb, RefGen_v3)
退火温度
Annealing
temp. (° C)
正向引物
Forward primer (5° -3° )
反向引物
Reverse primer (5° -3° )
N5M35138.4157TCCAAATTCACGCGATCCCTGGTCTTTTCCCC
N5M2138.5655ACATTCTCTTTGGAAGCTATAGAAGATTCGTAATAGGGCATAAGACTAAC
N5M38139.7057ACGGCTAGCTCTCGTTTCCCCGCTCCTGCCTTCACTCCAA
N5M98139.7355TGTGGCAGTACTTGTCTTTTCACCATTTGGTTTATTTTCG
N5M136140.0255TCAGAGTCGCTCCTCATCCTCGCGTCCACGTCTTCCT
N5M5140.0253CCCAGAGGTTTGAGCTATCACTCAGATTTTGTTTTGCCTTTGAAAGACTCGAAGAC
N5M158140.4555ATTGTAGGTTCGATTGATGAGTTTCTTGCGAGTTTTGTTA
N5M91140.6255GGCTCTACCTCTTCGCACAACCCACACATAAGGCTCAAAA
N5M180140.8455GAGCACATGAAGTTGTAACATTAACTCATTGGAGGTATTT
N5M202141.3357TTTTAACGGAGGGCCTTTCTGGATTCCTATCTAGCCACTT
N5M60141.4755CATGGGTGGATATGGAGAATTGACTGGGGTGAAATGGAAT
N5M61141.8055GGTGTAGGTGGAGAGAACATGCAGCATGATCG

表2 目标区段内开发的多态性InDel标记 Table 2 Polymorphic InDel markers within the target QTL

2.3 玉米穗行数主效位点qKRN5.04的多群体验证与遗传效应分析 利用qKRN5.04区段内的最显著SSR分子标记(基于F2:3群体定位结果), 进行各回交群体的基因型鉴定, 然后结合表型数据, 利用单标记分析(single marker analysis, SMA)方法, 对BC3F2、BC4F2、BC3F3、BC4F3、BC5F2等分离群体, 进行不同基因型材料间穗行数的显著性测验(表3图2-A, B, C)。结果表明, 基因作用方式均呈现部分显性; 除2014年海南三亚的BC3F3群体外, 农531和四路糯基因纯合型的穗行数均达到显著性差异。为了减小偶然因素造成的影响, 还进行了基于家系表型和基因型的主效位点遗传效应分析。结果表明, 在2014年海南三亚, 农531纯合型BC3F2:3家系(N531/N531)和四路糯纯合家系(SLN/SLN), 穗行数均值分别为16.70、15.24, 遗传效应为1.46行, 且两种基因型间穗行数差异达到极显著水平(P = 5.40E-09) (图2-D)。
表3
Table 3
表3(Table 3)
表3 玉米穗行数主效位点qKRN5.04的多群体验证和遗传效应分析 Table 3 Effect validation and genetic analysis of multiple populations for major locusqKRN5.04 on kernel row number
环境
Environment
世代
Generation
基因型分析 Genotype analysis多重比较 Multiple comparison基因作用
方式
Gene action
农531类型
N531/N531
杂合类型
SLN/N531
四路糯类型
SLN/SLN
农531类型
N531/N531
杂合类型
SLN/N531
四路糯类型
SLN/SLN
2014北京
2014 Beijing
BC3F217.2415.6114.79abcPD
BC4F216.9116.1515.36abcPD
2014石家庄
2014 Shijiazhuang
BC3F216.7915.1114.53abcPD
BC4F216.5315.5814.70abcPD
2014三亚
2014 Sanya
BC3F316.5616.1415.78aaaPD
BC4F317.6616.3315.25aabbPD
BC5F217.2016.8715.90aabPD
PD refers to partial dominance effect; different letters in multiple comparison column indicate statistically significant differences among groups (Tukey’ s honestly significant difference (HSD) test, P < 0.05).
PD代表部分显性效应; 多重比较(P < 0.05)字母不同代表不同组之间存在显著性差异。

表3 玉米穗行数主效位点qKRN5.04的多群体验证和遗传效应分析 Table 3 Effect validation and genetic analysis of multiple populations for major locusqKRN5.04 on kernel row number

图2
Fig. 2
Figure OptionViewDownloadNew Window
图2 玉米穗行数主效位点qKRN5.04在不同环境条件下的遗传效应比较分析
A, B分别代表2014年北京和河北石家庄环境下的主效位点qKRN5.04在BC3F2和BC4F2群体内的遗传效应比较; C代表2014年北京和河北石家庄2个环境下(群体合并计算) qKRN5.04的遗传效应比较; D代表2014年海南三亚环境下不同基因型BC3F2:3家系间的穗行数差异比较, * * 在代表水平P< 0.01差异极显著; N531、SLN分别是农531和四路糯选系的缩写。Fig. 2 Genetic effect comparisons for kernel row number related major QTL (qKRN5.04) under different environments
A, B represent the genetic effects analysis ofqKRN5.04(BC3F2 and BC4F2 population) in 2014 Beijing and 2014 Shijiazhuang, respectively; C represents the genetic effects analysis ofqKRN5.04in 2014 Beijing and 2014 Shijiazhuang (total populations); D represents the KRN comparison of different genotypes of BC3F2:3 family in 2014 Hainan, * * represents significant difference at P< 0.01; N531 and SLN are the abbreviation of Nong 531 and Silunuo, respectively.


2.4 基于多个高代回交群体的玉米穗行数主效位点qKRN5.04定位分析 在2014北京, 用BC3F2和BC4F2群体可同时将qKRN5.04定位在136.3~138.5 Mb的染色体区段之内, LOD值分别达37.56和16.68, 可以分别解释21.76%和8.68%的表型变异(表4图3-A)。在2014石家庄, 利用同样的BC3F2和BC4F2定位群体, 可以在紧密相邻的区域范围之内, 将该QTL定位于138.5~140.0 Mb之内, 可以解释的表型变异分别为15.29%和6.44% (表4图3-B)。遗传效应分析表明, 减效等位变异同样来自供体亲本四路糯, 效应值为0.80~1.76行。2015年, 在北京还种植了BC3F4、BC4F4、BC5F3分离群体, 其中用BC3F4和BC5F3群体可同时将qKRN5.04定位至相同染色体的139.7~140.0 Mb区段之内。同时发现, 在不同环境下该主效位点的基因作用方式以加性和部分显性为主。总之, 本研究定位结果与前期基于F2:3分离群体结果高度一致, 表明该位点能够在多个高代回交群体稳定地表达, 为重要的玉米穗行数主效位点。
利用一般线性模型分析主效位点与环境的互作效应(71个BC3F2:3家系、30个BC3F3:4家系和56个BC4F3:4家系, 同时种植于北京、河北石家庄和河南新乡3个环境), 没有检测到显著的互作效应, 表明qKRN5.04能够在不同环境中稳定表达。
表4
Table 4
表4(Table 4)
表4 不同回交群体玉米穗行数主效位点qKRN5.04的多环境定位分析 Table 4 QTL mapping for qKRN5.04using multiple backcross population in different environments
环境
Environment
世代
Generation
群体大小
Number
of plant
标记区间
Marker interval
物理位置
Position (Mb)
LOD加性效应
Add
显性效应
Dom
基因作用方式
Gene action
贡献率
PVE (%)
2014北京
2014 Beijing
BC3F2705umc1060-N5M2136.3-138.537.56-1.23-0.41PD21.76
BC4F2858umc1060-N5M2136.3-138.516.68-0.800.04A8.68
2014石家庄
2014 Shijiazhuang
BC3F2721N5M2-N5M5138.5-140.027.85-1.120.15A15.29
BC4F21046N5M2-N5M5138.5-140.014.25-1.00-0.12A6.44
2014北京
2014 Beijing
BC3F4296N5M98-N5M136139.7-140.02.59-0.470.34PD4.14
BC4F4189N5M35-N5M38138.4-139.63.09-0.180.11PD14.52
BC5F3640N5M98-N5M136139.7-140.08.91-1.76-0.58PD18.73
A and PD represent additive effect and partial dominance effect, respectively; PVE: phenotypic variation explained.
A和PD分别代表加性效应和部分显性效应。

表4 不同回交群体玉米穗行数主效位点qKRN5.04的多环境定位分析 Table 4 QTL mapping for qKRN5.04using multiple backcross population in different environments

图3
Fig. 3
Figure OptionViewDownloadNew Window
图3 不同环境条件下BC3F2和BC4F2群体的玉米穗行数主效位点qKRN5.04定位
图A和B分别代表qKRN5.04在2014年北京和2014年石家庄的定位结果。Fig. 3 QTL mapping of qKRN5.04using BC3F2 and BC4F2 populations in different environments
A and B represent the QTL mapping for qKRN5.04under Beijing and Shijiazhuang in 2014, respectively.


2.5 玉米穗行数主效位点qKRN5.04的精细定位 在利用多个回交群体进行主效QTL验证定位的工作基础之上, 还进一步构建了包含398个单株的BC5F3群体, 结合标记加密, 可将该主效位点定位至3.8 Mb的染色体区段之内(umc1060-N5M5) (图4-C)。在此基础上, 通过目标区段内多态性标记的进一步加密, 进行了不同重组类型与轮回亲本间的表型差异分析, 以实现穗行数主效位点qKRN5.04的精细定位(图4)。在umc1060-N5M5区段内, 可进一步筛选出3个多态性标记, 可将BC5F3群体分为12种重组类型(图4-D)。以受体亲本农531为对照(P2), 进行了各重组类型与对照组的t检验分析, 结果表明, 标记区段N5M98-N5M136内含四路糯等位基因的重组类型(R1、R6-R8), 其穗行数显著高于受体亲本类型(P2), 标记N5M98-N5M136之间含农531等位基因的交换类型(R2-R5, R9-R12), 穗行数与受体亲本类型之间无显著差异(图4-D)。因此, 推测调控玉米穗行数的关键基因位点位于标记N5M98-N5M136之间约300 kb的区间之内。
图4
Fig. 4
Figure OptionViewDownloadNew Window
图4 穗行数主效位点qKRN5.04的精细定位
a表示主效位点初定位区段, b和c分别表示BC5F2和BC5F3群体的定位结果, d表示亲本和12种重组类型, 表格表示每种重组类型相应的穗行数和单株个数, 黑色和白色分别代表来自“ 四路糯” (P1)和来自“ 农531” (P2)的染色体片段, * * 表示在0.01水平显著, * 表示在0.05水平显著。Fig. 4 Fine mapping of major QTL qKRN5.04 for kernel row number
a shows bin5.04 of Chr.5; b shows qKRN5.04 in BC5F2, and c shows qKRN5.04 in BC5F3, d indicates nine recombinant types (R1-R12) and two parents, and the table shows their corresponding kernel row number, with numbers of individuals for each recombinant on the right. Black and white bars indicate chromosomal segments from SLN (P1) and N531 (P2), respectively. * * and * mean significant difference at the P < 0.01 and 0.05, respectively.


3 讨论本研究的目标区段qKRN5.04初定位于81~150 Mb约69 Mb的区间, 最大可以解释13.55%的表型变异[22]。该定位结果与之前很多研究得出的结果区段重叠。例如, Yan等[5]、谭巍巍等[26]利用F2:3家系在第5染色体上分别定位到穗行数QTL rn5(80.8~ 151.6 Mb)、QTL qKRE1-5-2 (80.4~142.3 Mb); 吕学高等[27]利用大小为400的F2群体, 也在bin5.04检测到1个穗行数主效QTL RN3, 位于标记bnlg1660- phil09188之间, 贡献率为15.1%。周强等[3]通过meta-QTL分析, 在第5染色体定位到4个穗行数一致性QTL区间, 其中包括本研究的qKRN5.04区段。以上证据表明, qKRN5.04所在区段对玉米穗行数调控具有重要作用, 可作为穗行数遗传基础研究和候选基因挖掘的重要热点区域。
基因型与环境的互作对数量性状遗传定位具有重要的影响[28], 甚至可导致定位结果的偏差。但整体上, 主效QTL位点更容易在不同环境中被稳定地检测到[29]。本研究中, 在2014年北京, 利用BC3F2和BC4F2群体, 可同时将qKRN5.04定位在136.3~ 138.5 Mb之间; 在石家庄, 利用相同群体, 可将该位点定位在138.5~140.0 Mb之间, 与北京试验点定位区间紧密相邻。另外, 无法检测到显著的主效位点与环境间的互作效应。进一步表明, qKRN5.04能够在不同环境中稳定表达, 为环境不敏感型玉米穗行数主效位点。
初定位的区间较大, 不利于候选基因的挖掘和图位克隆, 利用回交构建导入系等次级作图群体是数量性状主效位点挖掘基因的有效策略之一[30]。导入系的独特优势在于, 除特定的染色体片段来源于供体外, 其他遗传组成与受体完全相同, 消除了遗传背景的干扰, 固此可通过对导入系与表型性状的显著性分析把目标性状基因或QTL定位于较小的区段内, 实现精细定位[31]。Zhang等[32]利用导入系BC4F3将控制穗粒数的主效位点qKN缩小到480 kb的范围内。Teng等[33]用单片段导入系(SIL)将株高主效位点qPH3.1定位到12.6 kb区间内。Liu等[21]利用近等基因系将控制穗行数的主效位点KRN4缩小到1.2 kb的区间内。本研究在已有初定位基础上, 利用分子标记辅助选择构建了针对主效区段的多个高代回交群体和次级分离群体, 在准确验证原有定位结果的同时, 利用重组事件分析将玉米穗行数主效位点qKRN5.04定位到约300 kb的染色体区段之内, 为下一步的候选基因挖掘提供了重要的工作支持。
在本研究定位到的~300 kb区间之内, 以玉米B73基因组为参考, 发现该区段内共存在10个基因, 其中7个能够在玉米雌穗中得以表达, 但仅有2个基因有功能注释, 分别为GRMZM2G300841和GRMZM2G060253。序列分析发现, GRMZM2G 300841在拟南芥中的同源基因是At3g24315, 即编码磷脂酰肌醇转运蛋白(PITPs)的成员之一SEC20转运蛋白。研究发现AtSec20在拟南芥花器官中高强度地特异表达[34, 35]。GRMZM2G060253在拟南芥中存在同源基因At4g23800, 编码HMG转录因子。研究表明, HMG转录因子具有保守的HMG-box结构域, 该结构域与生物生殖发育关系密切, 且有证据表明该基因与玉米雌雄花序发育具有一定的关系[36]。另外一个基因GRMZM5G897364功能尚未验证, 但是表达谱分析发现, 该基因在发育的玉米雌穗中表达量最高[35]。总之, 依据目前定位结果尚不足以确定候选基因, 仍需依据交换事件分析进一步缩小目标区段, 在此基础上克隆基因, 并通过转基因及互补试验验证功能, 挖掘并深入分析玉米穗行数主效位点qKRN5.04的关键调控基因。
4 结论以“ 四路糯” 选系和“ 农531” 分别为供体亲本和轮回亲本, 针对玉米穗行数主效位点qKRN5.04, 构建多个高代回交群体, 将其定位到第5染色体的136.3~140.0 Mb区间之内。同时, 通过次级分离群体重组事件分析, 将其进一步缩小到约300 kb的染色体区间范围内, 从而为目标位点的图位克隆提供了重要的工作基础。
The authors have declared that no competing interests exist.

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


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

[1]Li J Z, Zhang Z W, Li Y L, Wang Q L, Zhou Y G. QTL consistency and meta-analysis for grain yield components in three generations in maize. Theor Appl Genet, 2011, 122: 771-782[本文引用:1]
[2]李成璞, 白苇, 翟立红, 陶勇生, 张祖新. 玉米穗行数QTL及其互作分析. 植物遗传资源学报, 2011, 12: 965-970
Li C P, Bai W, Zhai L H, Tao Y S, Zhang Z X. Identifying and interaction assay of QTL for row number per ear of maize. J Plant Genet Resour, 2011, 12: 965-970 (in Chinese with English abstract)[本文引用:1]
[3]周强, 王平喜, 程备久, 朱苏文, 谢传晓. 玉米穗行数性状QTL的元分析. 玉米科学, 2014, 22(2): 35-40
Zhou Q, Wang P X, Cheng B J, Zhu S W, Xie C X. Meta-analysis of QTL for ear row number in maize. J Maize Sci, 2014, 22(2): 35-40 (in Chinese with English abstract)[本文引用:2]
[4]谭巍巍, 李永祥, 王阳, 刘成, 刘志斋, 彭勃, 王迪, 张岩, 孙宝成, 石云素, 宋燕春, 杨德光, 王天宇, 黎裕. 在干旱和正常水分条件下玉米穗部性状QTL分析. 作物学报, 2011, 37: 235-248
Tan W W, Li Y X, Wang Y, Liu C, Liu Z Z, Peng B, Wang D, Zhang Y, Sun B C, Shi Y S, Song Y C, Yang D G, Wang T Y, Li Y. QTL mapping of ear traits of maize under different water regimes. Acta Agron Sin, 2011, 37: 235-248 (in Chinese with English abstract)[本文引用:1]
[5]Yan J B, Tang H, Huang Y Q, Zheng Y L, Li J S. Quantitative trait loci mapping and epistatic analysis for grain yield and yield components using molecular markers with an elite maize hybrid. Euphytica, 2006, 149: 121-131[本文引用:2]
[6]Barazesh S, McSteen P. Barren inflorescence1 functions in organogenesis during vegetative and inflorescence development in maize. Genetics, 2008, 179: 389-401[本文引用:1]
[7]Barazesh S, Nowbakht C, McSteen P. Sparse inflorescence1, Barren inflorescence1 and Barren stalk1 promote cell elongation in maize inflorescence development. Genetics, 2009, 182: 403-406[本文引用:1]
[8]Wu X T, Skirpan A, McSteen P. Suppressor of sessile spikelets1 functions in the ramosa pathway controlling meristem determinacy in maize. Plant Physiol, 2009, 149: 205-219[本文引用:1]
[9]Taguchi-Shiobara F, Yuan Z, Hake S, Jackson D. The fasciated ear2 gene encodes a leucine-rich repeat receptor-like protein that regulates shoot meristem proliferation in maize. Genes Dev, 2001, 15: 2755-2766[本文引用:1]
[10]Vollbrecht E, Springer P S, Goh L, Buckler E S, Martienssen R. Architecture of floral branch systems in maize and related grasses. Nature, 2005, 436: 1119-1126[本文引用:1]
[11]Clermont Y, Bortiri E. Ramosa2 encodes a lateral organ boundary domain protein that determines the fate of stem cells in branch meristems of maize. Plant Cell, 2006, 18: 574-585[本文引用:1]
[12]Chuck G, Muszynski M, Kellogg E, Hake S, Schmidt R J. The control of spikelet meristem identity by the branched silkless1 gene in maize. Science, 2002, 298: 1238-1241[本文引用:1]
[13]Kaplinsky N J, Freeling M. Combinatorial control of meristem identity in maize inflorescences. Development, 2003, 130: 1149-1158[本文引用:1]
[14]Bomblies K, Wang R L, Ambrose B A, Schmidt R J, Meeley R B, Doebley J. Duplicate FLORICAULA/LEAFY homologs Zfl1 and Zfl2 control inflorescence architecture and flower patterning in maize. Development, 2003, 130: 2385-2395[本文引用:1]
[15]周红菊, 穆俊祥, 赵胜杰, 余四斌. 水稻高世代回交导入系耐盐性的遗传研究. 分子植物育种, 2005, 3: 716-720
Zhou H J, Mu J X, Zhao S J, Yu S B. Genetic analyses of salt tolerance in an advanced backcross population of rice. Mol Plant Breed, 2005, 3: 716-720 (in Chinese with English abstract)[本文引用:1]
[16]胡利宗, 刘均革, 郭晋杰, 赵永锋, 祝丽英, 宋占权, 陈景堂. 基于玉米BC2F2群体的穗部性状QTL分析. 华北农学报, 2010, 25(4): 73-77
Hu L Z, Liu J G, Guo J J, Zhao Y F, Zhu L Y, Song Z Q, Chen J T. QTL analysis of ear traits based on BC2F2 population in maize (Zea may L. ). Acta Agr Boreali-Sin, 2010, 25(4): 73-77 (in Chinese with English abstract)[本文引用:1]
[17]李卫华, 王洪秋, 袁亮, 张向歌, 谢慧玲, 胡彦民, 汤继华. 利用单片段代换系群体定位玉米穗部性状的QTL. 河南农业大学学报, 2013, 47: 143-146
Li W H, Wang H Q, Yuan L, Zhang X G, Xie H L, Hu Y M, Tang J H. Mapping of the QTL for ear related traits using a series of single segment substitution lines in maize. J Henan Agric Univ, 2013, 47: 143-146 (in Chinese with English abstract)[本文引用:1]
[18]齐欢欢, 段利超, 胡伟, 黄娟, 冯阳, 黄亚群, 祝丽英, 张祖新, 岳兵. 利用导入系群体对玉米产量及产量相关性状进行定位分析. 玉米科学, 2013, 21(4): 24-27
Qi H H, Duan L C, Hu W, Huang J, Feng Y, Huang Y Q, Zhu L Y, Zhang Z X, Yue B. Identification of significant loci for yield and yield-related traits in maize with introgression lines. J Maize Sci, 2013, 21(4): 24-27 (in Chinese with English abstract)[本文引用:1]
[19]Li F, Jia H T, Liu L, Zhang C X, Liu Z J, Zhang Z X. Quantitative trait loci mapping for kernel row number using chromosome segment substitution lines in maize. Genet Mol Res, 2014, 13: 1707-1716[本文引用:1]
[20]Tian B H, Wang J H, Wang G Y. Confirmation of a major QTL on chromosome 10 for maize kernel row number in different environments. Plant Breed, 2014, 133: 184-188[本文引用:1]
[21]Liu L, Du Y F, Shen X M, Li M F, Sun W, Huang J, Liu Z L, Tao Y S Zheng Y L, Yan J B, Zhang Z X. KRN4 controls quantitative variation in maize kernel row number. PLoS Genet, 2015, 11(11): e1005670[本文引用:2]
[22]焦付超, 李永祥, 陈林, 刘志斋, 石云素, 宋燕春, 张登峰, 黎裕, 王天宇. 特异玉米种质四路糯的穗行数遗传解析. 中国农业科学, 2014, 47: 1256-1264
Jiao F C, Li Y X, Chen L, Liu Z Z, Shi Y S, Song Y C, Zhang D F, Li Y, Wang T Y. Genetic dissection for kernel row number in the specific maize germplasm four-rowed waxy corn. Sci Agric Sin, 2014, 47: 1256-1264 (in Chinese with English abstract)[本文引用:2]
[23]Chen D H, Ronald P C. A rapid DNA min preparation method suitable for AFLP and other PCR applications. Plant Mol Biol Rep, 1999, 17: 53-57[本文引用:1]
[24]Li H, Ye G J, Wang J K. A modified algorithm for the improvement of composite interval mapping. Genetics, 2007, 175: 361-374[本文引用:1]
[25]Stuber C W, Edwards M D. Wendel J F. Molecular marker facilitated investigations of quantitative trait loci in maize: II. Factors influencing yield and its component traits. Crop Sci, 1987, 27: 639-648[本文引用:1]
[26]谭巍巍, 王阳, 李永祥, 刘成, 刘志斋, 彭勃, 王迪, 张岩, 孙宝成, 石云素, 宋燕春, 杨德光, 王天宇, 黎裕. 不同环境下多个玉米穗部性状的QTL分析. 中国农业科学, 2011, 44: 233-244
Tan W W, Wang Y, Li Y X, Liu C, Liu Z Z, Peng B, Wang D, Zhang Y, Sun B C, Shi Y S, Song Y C, Yang D G, Wang T Y, Li Y. QTL analysis of ear traits in maize across multiple environments. Sci Agric Sin, 2011, 44: 233-244 (in Chinese with English abstract)[本文引用:1]
[27]吕学高, 蔡一林, 陈天青, 徐德林, 王伟林, 刘志斋, 王久光. 玉米穗部性状QTL定位. 西南大学学报(自然科学版). 2008, 30(2): 64-70
Lyu X G, Cai Y L, Chen T Q, Xu D L, Wang W L, Liu Z Z, Wang J G. QTL mapping for ear traits in maize (Zea mays L. ). J Southwest Univ (Nat Sci Edn), 2008, 30(2): 64-70 (in Chinese with English abstract)[本文引用:1]
[28]Zhuang J Y, Lin H X, Lu J, Qian H R, Hittalmani S, Huang N, Zheng K L. Analysis of QTL × environment interaction for yield components and plant height in rice. Theor Appl Genet, 1997, 95: 799-808[本文引用:1]
[29]Fulton T. RFLP mapping of the rice genome. In: Rice Genetics: II. Proceedings of the Second International Rice Genetics Symposium, 14-18 May 1990. 2015. pp 435-442[本文引用:1]
[30]毛传澡, 程式华. 水稻农艺性状QTL定位精确性及其影响因素的分析. 农业生物技术学报, 1999, 7: 386-394
Mao C Z, Cheng S H. Analysis of accuracy and influence factor in QTL mapping about agronomic traits in rice (Oryza sativa L. ). J Agric Biol, 1999, 7: 386-394 (in Chinese with English abstract)[本文引用:1]
[31]Howell P M, Lydiate D J, Marshall D F. Towards developing intervarietal substitution lines in Brassica napus using marker-assisted selection. Genome, 1996, 39: 348-358[本文引用:1]
[32]Zhang G D, Wang X P, Wang B, Tian Y C, Li M, Nie Y X, Peng Q C, Wang Z L. Fine mapping a major QTL for kernel number per row under different phosphorus regimes in maize (Zea mays L. ). Theor Appl Genet, 2013, 126: 1545-1553[本文引用:1]
[33]Teng F, Zhai L H, Liu R X, Bai W, Wang L Q, Huo D G, Tao Y S, Zheng Y L, Zhang Z X. ZmGA3ox2, a cand idate gene for a major QTL, qPH3. 1, for plant height in maize. Plant J, 2013, 73: 405-416[本文引用:1]
[34]Heyndrickx K S, Klaas V. Systematic identification of functional plant modules through the integration of complementary data sources. Plant Physiol, 2012, 159: 884-901[本文引用:1]
[35]Bolduc N, Yilmaz A, Mejia-Guerra M K, Morohashi K, O’Connor D, Grotewold E, Hake S. Unraveling the KNOTTED1 regulatory network in maize meristems. Genes Dev, 2012, 26: 1685-1690[本文引用:2]
[36]牛林海, 杨国栋. 玉米HMG全长cDNA的克隆及其分析. 山东农业大学学报(自然科学版), 2002, 33: 239-241
Niu L H, Yang G D. Cloning and expression analysis of a full-length cDNA clone in maize. J Shand ong Agric Univ (Nat Sci), 2002, 33: 239-241 (in Chinese with English abstract)[本文引用:1]
相关话题/基因 环境 遗传 北京 鉴定