解松峰，吉万全，王长有，胡卫国，李俊，张耀元，师晓曦，张俊杰，张宏，陈春环

小麦穗部性状的主基因+多基因混合遗传模型分析

解松峰1,2，吉万全1，王长有1，胡卫国3，李俊4，张耀元1，师晓曦1，张俊杰1，张宏1，陈春环1

（1西北农林科技大学农学院/旱区作物逆境生物学国家重点实验室/国家小麦改良中心杨陵分中心，陕西杨凌 712100；2中国富硒产业研究院/ 农业部富硒产品开发与质量控制重点实验室/富硒食品开发国家地方联合工程实验室/安康市富硒产品研发中心，陕西安康 725000；3河南省农业科学院小麦研究中心，郑州 450002；4四川省农业科学院作物研究所，成都 610066）

【目的】穗部性状是小麦重要的产量性状，在小麦产量构成中占据重要地位和作用。开展小麦穗部性状遗传研究、分析其遗传机制，对制定高产育种策略、提高育种效率提供理论和实践指导。【方法】以主茎穗长、小穗数、穗粒数、小穗粒数为指标，采用数量性状的主基因+多基因混合遗传模型方法，对不同生态环境条件下来自母本品冬34与父本BARRAN及其衍生的F7:8、F8:9代重组自交系群体（RIL）进行穗部性状的遗传模型分析与遗传参数估计，以确定控制各性状的基因数目，估计遗传效应值及遗传率。【结果】穗长和小穗数的最佳遗传模型均是B-2-1（PG-AI），符合2对连锁主基因+加性-上位性多基因遗传模型。穗长的多基因遗传率是90.64%，小穗数的多基因遗传率是89.52%，穗长的环境变异平均值占表型变异的比例为9.39%，小穗数的环境变异平均值占表型变异的比例为10.50%；穗粒数的最佳遗传模型是G-1（MX3-AI-A），符合3对加性-上位性主基因加多基因+加性混合遗传模型，主基因遗传率是69.39%，多基因遗传率是29.94%，环境变异平均值占表型变异的比例为2.18%。控制穗粒数的第1对主基因的加性效应值和第3对主基因的加性效应值数值相等，同是4.56，具有正向效应。第2对主基因的加性效应值与加性效应和第1对主基因×第2对主基因×第3对主基因的加性效应值相同，均是-1.64，且为负向效应。加性和加性×加性上位性互作效应值与加性和第2对主基因加性×第3对主基因加性上位性互作效应值相等，均是-6.02。加性和第1对主基因加性×第3对主基因加性上位性互作效应值是0.18，多基因的加性效应值是0.15，表现为较低的正向遗传效应；小穗粒数的最佳遗传模型是H-1（4MG-AI），符合4个主基因+加性-上位性遗传模型，主基因遗传率是81.50%。第1至4对主基因加性效应值分别为0.22、0.18、-0.20和0.24，加性和第1对主基因×第1对主基因的加性上位性互作效应值是-0.170，加性和第1对主基因×第3对主基因的加性效应值是0.240，加性和第1对主基因×第4对主基因的加性效应值是-0.20，加性和第2对主基因×第3对主基因的加性效应值与加性和第2对主基因加性×第4对主基因加性上位性互作效应值绝对值相同，效应相反，前者值是0.03，后者值是-0.03。加性和第3对主基因×第4对主基因×的加性效应值是0.06。【结论】小麦穗部性状以多基因遗传效应为主，符合数量遗传特征，易受环境影响。小穗粒数存在着主基因遗传特性，主基因遗传力较高，受环境影响小，小穗粒数可作为有效改良穗部性状早期选择的直接指标，实现单株定向选择，提高育种效率。

小麦；穗部性状；主基因+多基因；遗传效应

0 引言

【研究意义】普通小麦（L.，AABBDD，2=6=42）是世界上三大主粮作物之一，在中国，小麦总产量位于世界前列，小麦的持续增产、稳产对中国乃至世界粮食安全有着重要的影响[1]。随着国内外学者对影响小麦产量因素地不断研究，穗部性状的改善越来越成为提高小麦产量的关键因素[2-3]。因此，研究小麦穗部性状的遗传规律，对提高小麦育种效率及高产稳产具有重要的指导意义。【前人研究进展】自20世纪20年代以来，数量性状在数量遗传学的建立之后得到了广泛研究，与小麦产量高低有关的大多数性状属于数量性状。因此，数量遗传学与植物育种相结合，应用数理统计方法和遗传模型分析群体数量性状遗传规律具有重要的意义，也为现代植物育种提供了科学的依据[4-5]。同时，植物数量性状的基因分为主基因和微效多基因，它们共同存在并影响植物的表型[6]。对于数量性状遗传系统，建立了“主基因多基因遗传模型分离与分析方法”[7-8]，用于鉴定主要基因、多基因效应、基因互作[9]、上位性以及基因与环境的相互作用[9-10]。近年来，植物数量性状的主基因加多基因混合遗传模型在大豆[9]、花生[11]、水稻[12]、油菜[13]、紫薇[14]、黄瓜[15]、甜瓜[16]、菊花[17]、小麦[18]、卷心菜[19]等植物育种中得到了广泛的应用，该遗传模型在小麦穗部性状研究中也得到应用，有的研究表明穗长受加性-显性-上位性多基因混合遗传模型控制[20-21]或主基因＋多基因遗传模型控制[22-24]，有研究认为小穗数符合主基因+多基因模型[20-23]，也有研究发现穗粒数由主基因多基因混合遗传模型控制[20-23]，研究者采用的材料不同，得到的遗传规律也不完全一样。已有的研究通常在一个环境下研究穗部性状遗传特点，然后进行遗传学统计分析。【本研究切入点】迄今为止，对于多个环境下小麦穗部性状的数量遗传分析的研究鲜见报道。西北农林科技大学农学院吉万全教授课题组前期以品冬34为母本和BARRAN为父本配置杂交组合，用8年时间构建了高代重组自交系。【拟解决的关键问题】本研究在2年4个环境下，F7:8和F8:9重组自交系（共476份衍生的家系），通过主基因+多基因混合遗传模型在小麦穗部性状的应用，确定小麦群体穗部性状的遗传模型和基因作用方式，对是否存在有效的加性基因进行分离分析，以期揭示穗部性状的主要遗传成分及其表现，为小麦遗传基础研究及优良品系的选择提供理论依据。

1 材料与方法

1.1 试验材料

以品冬34为母本和BARRAN为父本配置杂交组合，采用单籽粒传递方法，构建高代重组自交系。品冬34具有矮杆、千粒重高、籽粒大等特点，是由中国农业科学院作物品种资源研究所选育。BARRAN是引进国外小粒种质、粒小、千粒重较低。亲本间以穗部性状为代表的田间农艺性状表型差异较大。2007年得到F0代，2009—2016年用7年的时间，直到2016年得到F7:8重组自交系（共476份衍生的家系）。在3个地点（陕西杨凌、河南原阳和四川广汉）、4个环境（即2016SY（环境E1）、2017SY（环境E2）、2017HY（环境E3）、2017SG（环境E4））种植，设置2个年度重复，即2015—2016年和2016—2017年。环境E1、环境E2为陕西杨凌西北农林科技大学干旱半干旱地区节水农业实验农场（34°17′30″N，108°3′53″E）。环境E3为河南省农业科学院现代农业示范试验基地，地点在新乡原阳县（35°0′27″N，113°4′42″E），环境E4为四川省农业科学院在德阳广汉（31°0′29″N，104°23′43″E）的育种试验基地。田间常规管理，每个生长季节在10月份播种，并在6月初收获，每个环境采用随机区组设计，2次重复，行长2 m，行距24 cm，株距10 cm，每行播种20粒。

1.2 性状测定

穗部性状调查和测定参考李立会的方法[25]，每个品系调查10株，如穗长（spike length，SL）：主穗基部小穗节至穗顶端（不含芒）的长度；主穗小穗数（spikelets of main ear，SME）：主穗小穗数总数包括不育小穗数，另计每穗有效（结实）小穗数；主穗粒数（grains of main ear，GME）：全穗粒数或单穗实际粒数或单株每穗平均结实粒数；小穗最多粒数（grains of small ear，GSE）：穗中部最多粒数的小穗。

1.3 数据分析

利用SPSS Statistics 17.0软件（SPSS，Chicago，IL，USA）进行表型数据、偏度及峰度系数估计值基本统计。RIL群体中各基因位点为纯合，其遗传方差等于加性方差，广义遗传力等于狭义遗传力，公式为：B2=a/(a+e)（a为加性方差，e为环境方差）。根据Choo等[26]提出的基因间互作方式分析原理，由各性状偏度系数（1）、峰度系数（2）的正负号及其相对大小可以估计基因的上位性互作及其互作方式。利用R语言包绘制小麦重组自交系穗部性状的频率分布直方图；利用章元明教授团队最新开发的R软件包SEA[27-28]进行主基因+多基因混合遗传分析。其遗传模型原理是根据盖钧镒等[6]和ZHANG等[29]提出的植物数量性状“主基因+多基因混合遗传模型中P1、P2和RIL群体的联合世代分析方法”，获得7类38个遗传模型，进一步计算极大似然函数值（maximum likelihood method，MLV）并转换为AIC（Akaike’s information criterion，AIC）值；选出最小或接近最小AIC值的备选遗传模型[30]，进行样本分布与模型理论分布的适合性检验，包括均匀性检验（12、22、32）、Smirnov检验（w2）和Kolmogorov检验（D）。根据适合性检验的结果及AIC值最小原则选出最优遗传模型[31]。通过SEA软件包中最小二乘法原理计算出最优遗传模型的遗传效应等一阶遗传参数[32]、遗传方差和遗传率等二阶遗传参数[6,29]，其中2p、2mg和2pg分别是群体表型方差、主基因和多基因遗传方差；2mg和2pg分别是主基因和多基因遗传率。

2 结果

2.1 穗部性状的表型特征及次数分布

对群体单株产量相关的穗部性状进行了测定，并

分析了该群体的主茎穗长、小穗数、小穗粒数、穗粒数在4个环境中的表现。由表1可知，亲本各性状在4个环境下平均极差变化为：穗长介于5.00—10.67，小穗数介于14.00—21.00，穗粒数介于33.75—50.00，小穗粒数介于1.95—2.79。群体各性状在4个环境下极差变化为9.25%—61.01%，平均极差分别为：穗长介于4.16—17.32，小穗数介于11.55—36.96，穗粒数介于19.15—80.78，小穗粒数介于1.39—3.45。同时，各性状平均变异系数从高到低依次为：小穗数（0.13%）、单穗粒数（0.19%）、穗长（0.23%）、小穗粒数（0.32%）。

表1 重组自交系及其亲本穗部性状最佳线性无偏预测描述性分析

SL：穗长；SME：小穗数；GME：单穗粒数；GSE：小穗粒数。*：差异达到显著水平（<0.05）；**：差异达到极显著水平（<0.01）。下同

SL: Spike length; SME: Spikelets of main ear; GME: Grains of main ear; GSE: Grains of small ear. *: difference at the 0.05 probability level; **: significant difference at the 0.01 level. The same as below

由此可以看出，穗部表型变异范围较大，有超亲分离现象，有一定程度的分散，适合做遗传分析。此外，各性状的偏度系数介于-0.92—5.72，峰度系数介于-0.99—17.04。性状的偏度及峰度介于0—1，表现近似正态分布的有：E1和E2环境下的小穗数，E2和E3环境下的单穗粒数；各性状的遗传力介于0.34—0.97，4个环境下平均遗传力由高到低依次是，穗长的遗传力最高，为0.92，其次是小穗粒数的遗传力是0.88，穗粒数的遗传力是0.76，小穗数的遗传力是0.70。从表1可以看出，E1、E3和E4环境下，穗长的偏度和峰度系数无显著差异，表明基因间无互作。E2环境下的穗长，E1、E2、E3和E4环境下的小穗数，E1、E2和E3环境下的单穗粒数，E1、E2和E3环境下的小穗粒数，这些性状的偏度和峰度系数都显著大于0，表明基因间存在互作作用。

图1显示了4种环境，分别是环境E1：2016SY、环境E2：2017SY、环境E3：2017HY和环境E4：2017SG下的穗部性状的的频率（柱），理论（实线）和组分（虚线）产量相关性状。具有单峰分布特点的有E2环境下穗长，E2、E3和E4环境下的小穗数，E3环境下的穗粒数，E1、E2、E3和E4环境下的小穗粒数，表明这些性状可能为多基因控制。E1、E2和E3环境下的穗长，E1环境下的小穗数，E1、E2和E4环境下的穗粒数，这些性状分布呈现连续变异，具有连续双峰或多峰分布特点，符合主基因加多基因的特征，表明可能存在主基因。

2.2 穗部性状遗传模型的适合性检测

根据遗传模型选取的AIC值最小准则，选取AIC值最小的一组遗传模型作为备选模型（表2），穗部性状备选模型的极大似然值和AIC值见表2。在4个环境下分别比较，B-2-1（PG-AI）、F-1（3MG-AI）、H-1（4MG-AI）和H-1（4MG-AI）的AIC值最低，分别为2016.217、2164.521、1804.726和1960.836，B-2-1（PG-AI）、F-1（3MG-AI）和H-1（4MG-AI）为穗长遗传的备选模型；B-2-1（PG-AI）和F-2（3MG-A）的AIC值较低，分别为2184.5、1855.7、2256.652和2451.369。B-2-1（PG-AI）和F-2（3MG-A）为小穗数遗传的备选模型；B-1-9（2MG-IE）、E-2-5（MX2-ER-A）、E-2-7（MX2-CE-A）和G-1（MX3-AI-A）的AIC值较低，分别为3068.409、3879.773、3886.283和3562.874，B-1-9（2MG-IE）、E-2-5（MX2-ER-A）、E-2-7（MX2-CE-A）和G-1（MX3-AI-A）为穗粒数的备选模型是。同理，F-2（3MG-A）和H-1（4MG-AI）为小穗粒数的备选模型。

对备选模型进行一组适合性检验见表3，包括均匀性检验、Smirnov检验和Kolmogorov检验的5个统计量12、22、32、W2和D，选择具有最小AIC值和统计显著性水平个数最少的模型作为最优模型（表3）。由表3可知，穗长和小穗数的最佳遗传模型均B-2-1（PG-AI）模型，即2对连锁主基因+加性-上位多基因遗传模型。2对主基因连锁，具有明显的上位性作用，表示2对增效主基因在一起时，穗长有较大幅度的提高。表1中穗长最大的株系长于最长亲本5 cm，推测是2对增效主基因的聚合体。依次类推，G-1（MX3-AI-A）是穗粒数的最佳遗传模型，即3对加性-上位性主基因+加性多基因混合遗传模型。小穗粒数的最佳遗传模型是H-1（4MG-AI）模型，即4对主基因+加性上位性遗传模型。

2.3 穗部性状最适遗传模型的遗传参数估计

利用R软件包SEA对穗部性状的最适遗传模型估计其一阶参数和二阶参数（表4）。由表4中的一阶、二阶遗传参数可知：穗长的多基因遗传率是90.64%，小穗数的多基因遗传率是89.52%，穗长的环境变异平均值占表型变异的比例为9.39%，小穗数的环境变异平均值占表型变异的比例为10.50%；同样，穗粒数的主基因遗传率是69.39%，多基因遗传率是29.94%，环境变异平均值占表型变异的比例为2.18%。控制穗粒数的第1对主基因的加性效应值（da）和第3对主基因的加性效应值（dc）数值相等，同是4.56，具有正向效应。第2对主基因的加性效应值（db）与加性效应和第1对主基因×第2对主基因×第3对主基因的加性效应值（iabc）相同，均是-1.64，为负向效应。加性和加性×加性上位性互作效应值（iab）与加性和第2对主基因加性×第3对主基因加性上位性互作效应值（ibc）相等，均是-6.02，呈现较大的负向效应。加性和第1对主基因加性×第3对主基因加性上位性互作效应值（iac）值是0.18，多基因的加性效应值（d）是0.15，表现为较低的正向遗传效应；小穗粒数的主基因遗传率是81.50%，环境变异平均值占表型变异的比例为100.00%。第1至4对主基因加性效应值（da、db、dc和dd）分别为0.22、0.18、-0.20和0.24，加性和第1对主基因×第1对主基因的加性上位性互作效应值（iab）是-0.17，加性和第1对主基因×第3对主基因的加性效应值（iac）是0.24，加性和第1对主基因×第4对主基因的加性效应值（iad）是-0.20，加性和第2对主基因×第3对主基因×的加性效应值（ibc）与加性和第2对主基因加性×第4对主基因加性上位性互作效应值（ibd）绝对值相同，效应相反，前者值是0.03，后者值是-0.03。加性和第3对主基因×第4对主基因×的加性效应值（icd）是0.06。穗部性状的遗传参数见表4。

SL：穗长；SME：小穗数；GME：单穗粒数；GSE：小穗粒数；Frequency：频数；Density：密度

SL: Spike length; SME: Spikelets of main ear; GME: Grains of main ear; GSE: Grains of small ear

图1 RIL群体产量相关性状的次数（柱形）、理论（实线）与成分（虚线）分布

Fig. 1 Frequent (column), mixed (solid line, theoretical) and component (dotted line) distributions for wheat panicle-related traits in RILs

表2 小麦重组自交系品冬34×Barran F8、F9群体穗部性状最佳遗传模型分离分析的极大似然函数MLV值和Akaike信息准则AIC值

表续表2 Continued table 2

a：A模型表示1对主基因，无多基因遗传模型；B模型表示2对主基因，无多基因遗传模型；D模型表示1对主基因+多基因混合遗传模型；E模型表示2对主基因+多基因混合遗传模型；F模型表示3对主基因，无多基因遗传模型；G模型表示3对主基因+多基因混合遗传模型；H模型表示4对主基因+多基因遗传模型。b：MG：主基因模型；PG：多基因遗传模型；MX：主基因+多基因混合模型；A：加性效应；D：显性效应；I：互作；N：负向；E：相等；AI：加性上位性效应；EA：等加性；ED：显性上位；ER：隐性上位；AE：累加作用；CE：互补作用；DE：重叠作用；IE：抑制作用；CEA：全等加性；PEA：部分等加性；EEA：2个主基因等加性；EEEA：3个主基因等加性。例如： E-1模型MX2-ADI-AD，表示2对加性-显性-上位性主基因+加性-显性多基因混合遗传模型。MX3-CEA-A则表示3对等加性主基因+加性多基因混合遗传模型；粗体表示备选模型的AIC值（具有最低的2个AIC值）。下同

a: A: one major gene without polygene; B: two major genes without polygene; D: one major gene plus polygene mixed model; E: two major genes plus polygene mixed model; F: three major genes without polygene; G: three major genes plus polygene mixed model; H: four major genes plus polygene mixed model.b: MG: major gene model; PG: polygene model; MX: mixed major gene and polygene model; A: additive effect; D: dominance effect; I: interaction; N: negative; E: equal; AI: additive+epistasis effect; EA: equal additive effect; ED: epistasis dominance; ER: epistasis recessively; AE: accumulative effect; CE: complementary effect; DE: duplicate effect; IE: inhibition effect; CEA: congruent equal additive; PEA: partial equal additive; EEA: 2 major genes with equal additive effect; EEEA: 3 major genes with equal additive effect. Model E-1=MX2-ADI-AD, means mixed model with two major genes of additive-dominance-epistasis effects plus additive-dominance polygene. MX3-CEA-A: 3 major-genes with congruent equal additive effects plus polygenes mixed model. Bond latter are AIC for the selected optimal models. The same as below

表3 重组自交系群体（品冬34×Barran）穗部相关性状的最佳遗传模型适合性检验

12、22、32为均匀性检验统计量；W2为Smirnov检验统计量；D为Kolmogorov检验统计量；12、22、32检测统计量中括号内为相应的概率P；P1、P2：亲本；RILs：重组自交系

12,22,32are the statistic of Uniformity test, the numbers in brackets are the distribution values in theory;2is the statistic of Smirnov test;Dis the statistic of Kolmogorov test1; P1, P2: Parents; RILs: Recombinant inbred lines. SL: Spike length; SME: Spikelets of main ear; GME: Grains of main ear; GSE: Grains of small ear

表4 品冬34×Barran群体产量相关性状的部分模型的遗传参数

m：群体均方；d：主基因效应值；da：第一对主基因的加性效应；db：第二对主基因的加性效应值；dc：第三对主基因的加性效应值；dd：第四对主基因的加性效应值；iab(i*)：第1对主基因的加性×第2对主基因的加性上位性互作效应；iac：第1对主基因的加性×第3对主基因的加性效应；iad：第1对主基因的加性×第4对主基因的加性效；ibc：第2对主基因的加性×第3对主基因的加性效应；ibd：第2对主基因的加性×第4对主基因的加性效应；icd：第3对主基因的加性×第4对主基因的加性效应；iabc：3对主基因加性效应的互作值；[d]：多基因加性效应；2p：表型方差（群体方差）；2pg：多基因遗传方差；2mg：主基因遗传方差；2e：环境方差（误差方差）；2mg(%)：主基因遗传率；2pg(%)：多基因遗传率；—表示空缺

m: The mean value of P1generation; d: main gene effect value; da: Addictive effect of the first pair major gene; db: Addictive effect of the second pair major gene; dc: Addictive effect of the third pair major gene; dd: Addictive effect of the fourth pair major gene; iab(i*): Addictive effect plus addictive effect of the 1stpair major genes×the 2ndpair major gene; iac: Additive effect plus additive effect of the 1stpair major gene × the third pair major gene; iad: Additive effect plus additive effect of the 1stpair major gene × the fourth pair major gene; ibc: Additive effect plus additive effect of the 2ndpair major gene × the third pair major gene; ibd: Additive effect plus additive effect of the 2ndpair major gene × the fourth pair major gene: icd: Additive effect plus additive effect of the third pair major gene × the fourth pair major gene; iabc: Additive effect plus additive effect of the 1stpair major gene × the 2ndpair major gene × the third pair major gene; [d]: the additive effects of polygenes; σ2e: environmental variance;2p: phenotypic variance;2pg: polygene variance;2mg: major gene variance;2e: environmental variance;2mg(%): heritability of major gene;2pg(%): heritability of polygene-Var;“—”in the cells mean the value is absent

3 讨论

3.1 穗部性状的遗传模型

近年来，数量性状主基因+多基因混合遗传模型在不同群体、不同类型植物中应用广泛[15,17,33-41]，以穗部为代表的农艺性状研究开始受到大家重视[41,20-24,42-52]，多数性状表现连续性较好，作为数量性状进行遗传分析，但随着遗传学的发展，研究发现部分农艺性状也呈现不连续的变异，如穗粒数中的多粒特性表现为主基因控制的遗传[53,44-46]。

本研究中在环境E1、E2、E3和E4环境下，穗长的最佳遗传模型分别是H-1（4MG-AI）、B-2-1（PG-AI）、H-1（4MG-AI）、F-1（3MG-AI），在环境E1和E3环境下均有相同的遗传模型。其中E1和E3环境下的遗传模型H-1（4MG-AI）表示4对主基因-加性上位性遗传模型，这与程洁等[46]的穗长符合2对主基因遗传模型结构相似，与毕晓静等[21]的多基因遗传模型不同，与李法计等[22]和魏艳丽等[24]的主基因+多基因遗传模型不同。E4环境下的遗传模型F-1（3MG-AI）3对主基因（等比例），加性-上位性主基因遗传模型。以各环境下的最佳模型作为候选模型进行适合性检验发现，E2环境下的遗传模型也就是穗长的最佳遗传模型：B-2-1（PG-AI）即2对连锁主基因+加性-上位多基因遗传模型，主要受多基因模型控制的遗传，这与前人的研究观点相似[20-21,54-56]，但也有通过模型研究表明穗长主基因控制存在[20-24]，这可能因为研究选育群体类型、环境条件、多环境胁迫条件不一样有关。主茎小穗数在4个环境下，前3个环境E1、E2和E3的最佳遗传模型相同，都是模型B-2-1（PG-AI），表示2对连锁主基因（等比例），加性-上位性多基因遗传模型。环境E4下最佳遗传模型是F-2（3MG-A），表示3对主基因（等比例），加性（无显性）主基因遗传模型。经过适合性检验表明小穗数的最佳遗传模型是B-2-1（PG-AI），这与程洁的主基因遗传模型结构相似，与魏艳丽等[24]、李法计等[22]的主基因+多基因遗传模型不同。与房敬业等[57]研究认为多小穗性主要受2对主效基因控制结果差异较大，与卢翔等[55]的２对主基因＋多基因模型以及杜希朋的微效多基因模型也不相同。本研究中穗长和小穗数在不同环境下的最佳遗传模型相同，这或许和基因的连锁性以及一因多效等有关，揭示出相关性状的遗传具有相似或相同性。表明在相近或相同的染色体区域或QTL位点含有控制有关表型的基因[58]。穗粒数在E1、E2、E3和E4每个环境下最佳的遗传模型依次是G-1（MX3-AI-A）、E-2-7（MX2-CE-A）、B-1-9（2MG-IE）和E-2-5（MX2-ER-A）。其中，E1环境下的G-1（MX3-AI-A）模型含义是3对加性-上位性主基因+加性多基因混合遗传模型，E2环境下的E-2-7（MX2-CE-A）模型含义是2对互补作用主基因（2对连锁主基因，等比例）加性多基因混合遗传模型，E3环境下的B-1-9（2MG-IE）模型的含义是2对抑制作用主基因（抑制B等位基因对基因座Aa的影响）模型，E4环境下的E-2-5（MX2-ER-A）模型的含义是2对隐性上位主基因（等比例）加性多基因混合遗传模型。穗粒数在4个环境下的最佳遗传模型是G-1（MX3-AI-A）即3对加性-上位性主基因加多基因+加性多基因混合遗传模型，本研究下的最佳多基因混合遗传模型结果与李法计等[22]、卢翔等[55]、杜希朋等[56]相似，与程洁等[46]和王健胜等[53]等的主基因遗传模型不同；小穗粒数在E1、E2、E3和E4每个环境下最佳的遗传模型依次是F-2（3MG-A）、H-1（4MG-AI）、H-1（4MG-AI）和H-1（4MG-AI），E2、E3和E4环境下模型均相同。E1环境下的模型是F-2（3MG-A），表示3对主基因(等比例)+加性（无显性）主基因遗传模型，E2、E3和E4环境下的模型均是H-1（4MG-AI）表示4对主基因+加性上位性遗传模型。小穗粒数在4个环境下的最佳遗传模型是H-1（4MG-AI）即4对主基因+加性-上位性遗传模型，这与程洁等[46]的主基因研究结果相似。

3.2 遗传模式对小麦穗部育种的指导意义

本研究对小麦穗部性状的育种实践具有重要指导意义。一方面，由于分析群体样本足够大，从而提高了分析结果的准确性。另一方面，从遗传率的大小来看，各性状的遗传率介于29.94—90.64，这与广义遗传率大小范围一致，同为穗长的遗传力最高。穗长的环境变异平均值占表型变异的比例为9.39%，小穗数的环境变异平均值占表型变异的比例为10.50%，穗粒数环境变异平均值占表型变异的比例为2.18%。说明穗部各性状的遗传贡献率较大，主要受到遗传因素影响，同样也受环境效应的影响。因此，小麦育种实践中要充分考虑主基因和微效多基因的影响，及时开展早期有效的单株选择，重点关注高代单株性状的后期稳定表现，排除环境影响，并通过适当的田间管理进行调节。穗部性状的穗长、主茎小穗数、穗粒数和小穗粒数等小麦产量形成的重要因素，也是提高小麦产量的重要途径[59]。通常，数量性状表现出个体间的表型差异，并且具有复杂的遗传基础，这涉及到整个基因组中的许多基因，对整体表型的贡献是可变的[60]，环境被定义为所有影响表型的非遗传变量的综合影响，增加了数量性状的复杂性[61]，主基因+多基因混合遗传模型方法与QTL模型检测结果连续性好[62]。只是前者应用更为方便、快捷、实用，并且只需要表型数据，研究结果还需要从分子机理上进一步通过QTL验证分析，最终为小麦分子标记辅助选择和图位克隆提高理论基础。为了满足不断增长的人口和不断增加的粮食消费的需求，需要不断增加小麦产量。高产小麦及其衍生品种在小麦新品种选育中发挥了重要作用。课题组自主构建的大重组自交系遗传分析群体虽然花费了大量时间，但其高代群体可以多年多点重复使用，其亲本农艺性状差别明显，含高产、抗病、优质等基因，遗传分析群体中表型优异的高代单株可以作为优良品系培育新品种使用。下一步将穗部性状结合高通量测序进行遗传定位分析，从基因序列等分子水平上进一步挖掘候选基因，对基因位点Meta及关联分析等，阐述其遗传调控网络，为小麦高产稳产的分子标记辅助育种提供理论依据。

4 结论

明确了穗长和小穗数的最佳遗传模型，发现了2对主基因连锁，具有明显的上位性作用。穗粒数由3对加性上位性主基因+加性多基因混合遗传模型。小穗粒数性状通过4对主基因+加性上位性遗传模型来控制它们的遗传。穗长和小穗数在不同环境下的最佳遗传模型相同，或许与基因的连锁性以及一因多效等有关。表明在相近或相同的染色体区域或QTL位点可能含有控制有关表型的基因。小麦穗部性状以多基因遗传效应为主，易受环境影响，表现为正向遗传，符合数量遗传特征。小穗粒数存在着主基因遗传特性，受正值的影响，主基因遗传力较高，受环境影响小，小穗粒数可作为有效改良穗部性状早期选择的直接指标，实现单株定向选择，提高育种效率。

[1] Hawkesford M J, Araus J L, Park R, Calderini D, Miralles D, Shen T, Zhang J, Parry M A J. Prospects of doubling global wheat yields., 2013, 2(1): 34-48.

[2] Ku M S, Agarie S. Nomura M, Fukayama H, Tsuchida H, Ono K, Hirose S, Toki S, Miyal M, Matsuoka M. High-level expression of maize phosphoenolpyruvate carboxylase in transgenic rice plants., 1999, 17(1): 76-80.

[3] Ficher R A, Ress D, Sarre K D, Lu Z M, Condon A G, Larque S A. Wheat yield associated with high stomatal conductance and photosynthetic rete and cooler canopies., 1998, 278(6): 1467-1475.

[4] KUMAR S. Quantitative genetics, molecular markers, and plant improvement., 2014,4(10): 502-511.

[5] Beche E, Benin G, Silva C L, Munaro L B, Marchese J A. Genetic gain in yield and changes associated with physiological traits in Brazilian wheat during the 20th century., 2014, 61: 49-59.

[6] 盖钧镒, 章元明, 王建康. 植物数量性状遗传体系. 北京: 科学出版社, 2003: 96-102.

Gai J Y, Zhang Y M, Wang J K. Genetic System of Quantitative Traits in Plants. Beijing:, 2003: 96-102. (in Chinese)

[7] 盖钧镒，章元明，王建康. QTL混合遗传模型扩展至2对主基因+多基因时的多世代联合分析. 作物学报, 2000, 26(4): 385-391.

Gai J Y, Zhang Y M, Wang J K. A joint analysis of multiple generations for QTL models extended to mixed two major genes plus polygene., 2000,26(4):385-391. (in Chinese)

[8] Wang J K, Podlich D W, Cooper M, DeLacy I H. Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits., 2001, 103: 804-816.

[9] Wang J K, Gai J Y. Mixed inheritance model for resistance to agromyzid beanfly (Zehntner) in soybean., 2001, 122(1): 9-18.

[10] Gai J Y. Segregation analysis on genetic system of quantitative traits in plants., 2006, 1(1): 85-92.

[11] 黄冰艳, 张新友, 苗利娟, 刘华, 秦利, 徐静, 张忠信, 汤丰收, 董文召, 韩锁义, 刘志勇. 花生油酸和亚油酸含量的遗传模式分析. 中国农业科学, 2012, 45(4): 617-624.

Huang B Y, Zhang X Y, Miao L J, Liu H, Qin L, Xu J, Zhang Z X, Tang F S, Dong W Z, Han S Y, Liu Z Y. Inheritance analysis of oleic acid and linoleic acid content ofL..,2012, 45(4): 617-624. (in Chinese)

[12] Li H, Ribaut J M. Inclusive composite interval mapping(ICIM) for digenic epistasis of quantitative traits in biparental populations., 2008, 116: 243-260.

[13] 汪文祥, 胡琼, 梅德圣, 李云昌, 周日金, 王会成, 洪涛, 付丽, 刘佳. 甘蓝型油菜分枝角度主基因+多基因混合遗传模型及遗传效应. 作物学报, 2016, 42(8): 1103-1111.

Wang W X, Hu Q, Mei D S, Li Y C, Zhou R J, Wang H C, Hong T, Fu L, Liu J. Genetic effects of branch angle using mixture model of major gene plus polygene inL.., 2016, 42(8): 1103-1111. (in Chinese)

[14] Ye Y J, Wu J Y, Feng L, Ju Y Q, Cai M, Cheng T R, Pan H T, Zhang Q X. Heritability and gene effects for plant architecture traits of crape myrtle using major gene plus polygene inheritance analysis., 2017, 225: 335-342.

[15] 曹齐卫,张允楠, 王永强, 杨桂兰, 孙小镭, 李利斌. 黄瓜节间长的主基因+多基因混合遗传模型分析. 农业生物技术学报, 2018, 26(2): 205-212.

Cao Q W, Zhang Y N, Wang Y Q, Yang G L, Sun X L, Li L B. Genetic analysis of internode length using mixed major- gene plus polygene inheritance model in., 2018, 26(2): 205-212. (in Chinese)

[16] Qi Z Y, Li J X, Raza M A, Zou X X, Cao L W, Rao L L, Chen L P. Inheritance of fruit cracking resistance of melon (L.) fitting E-0 genetic model using major gene plus polygene inheritance analysis., 2015, 189: 168-174.

[17] 赵倩茹, 钟兴华, 张飞, 房伟民, 陈发棣, 滕年军. 切花小菊绿心性状杂种优势与混合遗传分析. 中国农业科学, 2018, 51(5): 964-976.

Zhao Q R, Zhong X H, Zhang F, Fang W M, Chen F D, Teng N J. Heterosis and mixed genetic analysis of green-center trait of spray cut chrysanthemum., 2018, 51(5): 964-976. (in Chinese)

[18] Khan M I, Khattak G S S, Khan A J, Khan A J, Subhan F, Mohammad T, Ali A. Genetic control of flag leaf area in wheat () crosses., 2012,7(27): 3978-3990.

[19] Cao X W, Cui H M, Li J, Xiong A S, Hou X L, Li Y. Heritability and gene effects for tiller number and leaf number in non-heading Chinese cabbage using joint segregation analysis., 2016, 203: 199-206.

[20] 李树华, 张文杰, 白海波, 吕学莲, 董建力, 惠建, 魏亦勤, 康学兵. 春小麦穗部性状的主基因+多基因遗传分析. 中国农学通报, 2017, 33(6): 20-26.

Li S H, Zhang W J, Bai H B, Lü X L, Dong J L, Hui J, Wei Y Q, Kang X B. Genetic analysis of major gene plus polygene of spike traits of spring wheat.. 2017, 33(6): 20-26. (in Chinese)

[21] 毕晓静, 史秀秀, 马守才, 韩芳, 亓佳佳, 李清峰, 王志军, 张改生, 牛娜. 小麦农艺性状的主基因+多基因遗传分析. 麦类作物学报, 2013, 33(4): 630-634.

Bi X J, Shi X X, Ma S C, Han F, Qi J J, Li Q F, Wang Z J, Zhang G S, Niu N. Genetic analysis of agronomic traits related to yield based on majou gene plus polygene model in wheat., 2013, 33(4): 630-634. (in Chinese)

[22] 李法计, 常鑫, 王宇娟, 宋全昊, 田芳慧, 孙道杰. 小麦重组自交系群体9个重要农艺性状的遗传分析. 麦类作物学报, 2013, 33(1): 23-28.

Li F J, Chang X, Wang Y J, Song Q H, Tian F H, Sun D J. Genetics analysis of nine important agronomic traits in wheat population of recombinant inbred lines., 2013, 33(1): 23-28. (in Chinese)

[23] 朱欣果, 万洪深, 李俊, 郑建敏, 唐宗祥, 杨武云. 人工合成小麦育种优势的主基因+多基因混合遗传分析. 南京农业大学学报, 2018, 41(4): 625-632.

Zhu X G, Wan H S, Li J, Zheng J M, Tang Z X, Yang W Y. Mixed major-genes plus polygenes inheritance analysis for breeding superiority in synthetic hexaploid wheat., 2018, 41(4): 625-632. (in Chinese)

[24] 魏艳丽, 王彬龙, 李瑞国, 蒋会利, 张安静. 大穗小麦穗部性状的遗传分析. 麦类作物报, 2015, 35(10): 1366-1371.

Wei Y L, Wang B L, Li R G, Jiang H L, Zhang A J. Genetic analysis on spike characteristics of wheat variety with large spike., 2015, 35(10): 1366-1371. (in Chinese)

[25] 李立会, 李秀全. 小麦种质资源描述规程和数据标准. 北京: 中国农业出版社, 2006.

Li L H, Li X Q.. Beijing: China agriculture press, 2006. (in Chinese)

[26] Choo T M, Reinbergs E. Estimation of the number of genes in doubled haploid populations of barley()., 1982, 24(3): 337-341.

[27] 章元明，盖钧镒，王永军. 利用P1、P2和DH或RIL群体联合分离分析的拓展.遗传，2001，23(5): 467-470.

Zhang Y M, Gai J Y, Wang Y J. An expansion of joint segregation analysis of quantitative trait for using P1, P2 and DH or RIL populations., 2001, 23(5): 467-470. (in Chinese)

[28] 曹锡文, 刘兵, 章元明. 植物数量性状分离分析windows软件包SEA的研制. 南京农业大学学报, 2013, 36: 1-6.

Cao X W, Liu B, Zhang Y M.SEA: a software package of segregation analysis of quantitative traits in plants., 2013, 36: 1-6. (in Chinese)

[29] Zhang Y M, Gai J Y, Yang Y H. The EIM algorithm in the joint segregation analysis of quantitative traits., 2003, 81: 157-163.

[30] Akaike H. On entropy maximization principle.. The Netherland:1977: 27-41.

[31] Cai C C, Tu J X, Fu T D. The genetic basis of flowering time and photoperiod sensitivity in rapeseed (L.)., 2008, 44: 326-333.

[32] Gamble E E. Gene effects in corn (L.): I. Separation and relative importance of gene effects for yield., 1962, 42: 339-348.

[33] 李英双, 胡丹, 聂蛟, 黄科慧, 张玉珂, 张园莉, 佘恒志, 方小梅, 阮仁武, 易泽林. 甜荞株高和茎粗的遗传分析. 作物学报, 2018, 44(8): 1185-1195.

Li Y S, Hu D, Nie J, Huang K H, Zhang Y K, Zhang Y L, She H Z, Fang X M, Ruan R W, Yi Z L. Genetic analysis of plant height and stem diameter in common buckwheat., 2018, 44(8): 1185-1195. (in Chinese)

[34] 郝贤伟, 徐秀红, 许家来, 崔胜利, 王传义, 张兴伟, 任夏, 朱佩, 张忠锋. 烤烟耐烤性的遗传效应. 中国农业科学, 2012, 45(23): 4939-4946.

Hao X W, Xu X H, Xu J L, Cui S L, Wang C Y, Zhang X W, Ren X, Zhu P, Zhang Z F. Genetic effects of holding curing potential in flue-cured tobacco., 2012, 45(23): 4939-4946. (in Chinese)

[35] 张保雷, 张卫东, 高庆荣, 王茂婷, 李楠楠, 张艳玉, 王慧娜, 高建华, 赵兰飞，茹振刚. 温光敏雄性不育小麦BNS育性的遗传效应分析. 中国农业科学, 2013, 46(8): 1533-1542.

Zhang B L, Zhang W D, Gao Q R, Wang M T, Li N N,Zhang Y Y, Wang H N, Gao J H, Zhao L F, Ru Z G. Genetic analysis on male sterility of thermo-photo-sensitive male sterile line BNS in wheat., 2013, 46(8): 1533-1542. (in Chinese)

[36] 王金社, 李海旺, 赵团结, 盖钧镒. 重组自交家系群体4对主基因加多基因混合遗传模型分离分析方法的建立. 作物学报, 2010, 36(2): 191-201.

Wang J S, Li H W, Zhao T J, Gai J Y. Establishment of segregation analysis of mixed inheritance model with four major genes plus polygenes in recombinant inbred lines population., 2010, 36(2): 191-201. (in Chinese)

[37] 张晓芬, 陈晓慧, 陈斌, 韩华丽, 耿三省. 农业生物技术学报, 2013, 21(4): 407-412.

Zhang X F, Chen X H, Chen B, Han H L, Geng S S. Genetic analysis of trichome density on the main stem and leaves in a recombinant inbred lines population derived from wild pepper (L.)., 2013, 21(4): 407-412. (in Chinese)

[38] 吴浪, 刘婧仪, 梁燕. 番茄绿果与红果颜色性状遗传的研究. 园艺学报, 2016, 43(4): 674-682.

Wu L, Liu J Y, Liang Y. Inheritance on fruit color character between green and red of Tomato., 2016, 43(4): 674-682. (in Chinese)

[39] 彭辉, 陈发棣, 房伟民, 蒋甲福, 陈素梅, 管志勇, 廖园. 切花小菊分枝性状杂种优势表现与遗传分析, 园艺学报, 2013, 40(7): 1327-1336.

Peng H, Chen F D, Fang W M, Jiang J F, Chen S M, Guan Z Y, Liao Y. Heterosis and mixed genetic analysis of branch traits of cut chrysanthemum., 2013, 40(7): 1327-1336. (in Chinese)

[40] 江建华, 洪德林, 郭媛, 张启武. 粳稻穗角与谷粒性状的相关性及谷粒性状遗传分析. 植物学报, 2009, 44(2): 167-177.

Jiang J H, Hong D L, Guo Y, Zhang Q W. Correlation between panicle angle and grain traits, and genetic analysis of grain traits in Japonica rice ()., 2009, 44(2): 167-177. (in Chinese)

[41] 江建华, 张启武, 洪德林. 粳稻穗部性状遗传分析. 植物学报, 2010, 45(2): 182-188.

Jiang J H, Zhang Q W, Hong D L. Genetic analysis of panicle traits inssp. japonica., 2010, 45(2): 182-188. (in Chinese)

[42] 匡勇, 罗丽华, 周倩倩, 何云礼, 范西林, 肖颖慧. 水稻籼粳交重组自交系群体穗部性状的相关和遗传分析. 华北农学报, 2011, 26(3): 72-78.

Kuang Y, Luo L H, Zhou Q Q, He Y L, Fan X L, Xiao Y H. Genetic and correlation analysis of the panicle traits of Recombinant Inbred Lines derived from an Indica/Japonica rice cross., 2011, 26(3): 72-78. (in Chinese)

[43] 郑建敏, 蒲宗君, 李式昭, 李俊杨. 人工合成小麦CI-LD抗穗发芽遗传特性分析. 麦类作物学报, 2015, 35(4): 464-470.

Zheng J M, Pu Z J, Li S Z, Li J y. Genetic analysis of pre-harvest sprouting resistance in synthetic wheat CI-LD., 2015, 35(4): 464-470. (in Chinese)

[44] 杨兴圣, 梁子英, 李华, 沈玮囡, 李美霞, 奚亚军, 王竹林, 刘曙东. 普通小麦籽粒性状的主基因＋多基因遗传模型分析. 麦类作物学报, 2013, 33(6): 1119-1127.

Yang X S, Liang Z Y, Li H, Shen W N, Li M X, Xi Y J, Wang Z L, Liu S D. Analysis on genetic model of grain characteristics in common wheat by mixed inheritance model of major genes plus polygenes., 2013, 33(6): 1119-1127. (in Chinese)

[45] 张安静, 张俊祖, 刘凤琴, 罗洪溪, 王彬龙, 赵会利.应用极大似然法分析小麦穗长的遗传. 中国农学通报, 2006, 9: 182-185.

Zhang A J, Zhang J Z, liu F Q, Luo H X, Wang B L, Zhao H L. The genetics law on the main stem ear length of long spike wheat lines using maximum likelihood methods., 2006, 22(9): 182-185. (in Chinese)

[46] 程洁, 周荣全, 吴玉川, 宋新颖, 林琪, 穆平. 不同水分条件下小麦穗部性状的遗传分析. 华北农学报, 2015, 30(增刊): 146-151.

Cheng j, Zhou R Q, Wu Y C, Song X Y, Lin Q, Mu P. Genetic analysis of spike traits in wheat cultivated in contrasted water conditions in wheat., 2015, 30(suppl.):146-151. (in Chinese)

[47] Cui F, Li J, Ding A. Conditional QTL mapping for plant height with respect to the length of the spike and internode in two mapping populations of wheat., 2011, 122(8): 1517-1536.

[48] Cui F, Ding A M, Li J. QTL detection of seven spike related traits and their genetic correlations in wheat using two related RIL populations,, 2012, 186(1): 177-192.

[49] 任哓波, 兰秀锦, 汪加丽. 人工合成小麦穗部特异性状的遗传分析. 四川农业大学学报, 2006, 24(4): 375-380.

Ren X B, Lan X J, Wang J L. Genetic analysis on special character in spike of synthetic hexaploid wheat., 2006, 24(4): 375-380. (in Chinese)

[50] 王新宁, 杜旭烨, 李斌, 王振林, 贺明荣, 李安飞, 贾继增, 孔令让. 察雅折达29×偃展 1号重组自交系群体主要农艺性状遗传分析. 山东农业科学, 2010, 7: 17-19.

Wang X N, Du X Y, Li B, Wang Z L, He M R, Li A F, Jia J Z, Kong L R. Genetic analysis of main agronomic traits in recombinant inbred lines of chayazheda 29×Yanzhan., 2010, 7: 17-19. (in Chinese)

[51] 林志强, 郑燕, 蔡英杰, 黄姗, 李志勇, 沈伟伟, 郑秀娟, 梁康迳. 水稻长穗大粒 RIL 群体产量、穗部和谷粒性状的遗传分析. 福建农林大学学报(自然科学版), 2011, 40(5): 449-454.

Lin Z Q, Zheng Y, Cai Y J, Huang S, Li Z R, Shen W W, Zheng X J, Liang K J. Genetic analysis on yield，panicle and grain traits in rice RIL population of long panicle and big grain., 2011, 40(5): 449-454. (in Chinese)

[52] 向道权, 黄烈健, 曹永国, 戴景瑞. 玉米产量性状主基因-多基因遗传效应的初步研究. 华北农学报, 2001, 16(3): 1-5.

Xiang D Q, Huang L J, Cao Y G, Dai J R. A preliminary study on genetic effect of maize yield component traits based on major gene and polygene mixed inheritance., 2001, 16(3): 1-5. (in Chinese)

[53] 王健胜, 王辉, 刘伟华, 武军, 李立会. 小麦-冰草多粒新种质及其多粒性遗传分析. 中国农业科学, 2009, 42(6): 1889-1895.

Wang J S, Wang H, Liu W H, Wu J, Li L H. The large kernel number in the novel wheat-agropyron germplasm 3228 and its inheritance analysis., 2009, 42(6): 1889-1895. (in Chinese)

[54] 闫林. 大穗小麦西农9814主要性状遗传分析及性状改良研究[D]. 杨凌: 西北农林科技大学, 2009.

Yan L. Genetic analysis of main traits and research on traits improvement for big ears wheat xinong 9814[D]. Yangling: Northwest A &F University, 2009. (in Chinese)

[55] 卢翔, 张锦鹏, 王化俊, 杨欣明, 李秀全, 李立会．小麦-冰草衍生后代3558-2穗部相关性状的遗传分析和QTL定位. 植物遗传资源学报, 2011, 12(1): 86-91．

Lu X, Zhang J P, Wang H J, Yang X M, Li X Q, Li L H. Genetic analysis and QTL mapping of wheat spike traits in a derivative line 3558-2 from wheat agropyron cristatum Offspring., 2011, 12(1): 86-91. (in Chinese)

[56] 杜希朋, 闫媛媛, 刘伟华, 高爱农, 张锦鹏, 李秀全, 杨欣明, 车永和, 郭小敏.蚂蚱麦×碧玉麦杂交F2代部分农艺性状的遗传分析. 麦类作物学报, 2011, 31(4): 624-629.

Du X P, Yan Y Y, Liu w H, Gao A N, Zhang J P, Li X Q, Yang X M, Che Y H, Guo X M. Genetic analysis on several important agronomic traits in F2generation of mazhamai ×quality., 2011, 31(4): 624-629. (in Chinese)

[57] 房敬业, 孙东发. 多小穗小麦51885的多小穗性遗传及与剑叶关系的初步研究. 华中农业大学学报, 2004, 34: 61-67.

Fang J Y, Sun D F. Inheritance of supernumerary spikelets and its relations with flag leaf characters of supernumerary spikelets line 51885 in wheat., 2004, 34: 61-67. (in Chinese)

[58] 范平, 詹克慧, 孙建英, 王淑凤, 赵国山. 小麦主要性状的遗传模型分析. 河南农业大学学报, 1999, 33(3): 231-234.

Fan P, Zhan k h, Sun J Y, Wang S F, Zhao G S. Analysis on the genetic models for main characters in wheat., 1999, 33(3): 231-234. (in Chinese)

[59] 许为刚, 胡琳, 吴兆苏, 盖钧镒. 关中小麦品种产量与产量结构遗产改良的研究. 作物学报, 2000, 26(3): 352-358.

Xu W G, Hu L, Wu Z S, Gai J Y. Studies on genetic improvement of yield and yield components of wheat cultivars in Mid-Shaanxi area., 2000, 26(3): 352-358. (in Chinese)

[60] Holland J B. Genetic architecture of complex traits in plants., 2007, 10: 156-161.

[61] Xu H, Zhu J. Statistical approaches in QTL mapping and molecular breeding for complex traits., 2012, 57: 2637-2644.

[62] 王春娥, 盖钧镒, 傅三雄, 喻德跃, 陈受宜. 大豆豆腐和豆乳得率的遗传分析与 QTL 定位. 中国农业科学, 2008, 41: 1274-1282.

Wang C E, Gai J Y, Fu S X, Yu D Y, Chen S Y. Inheritance and QTL mapping of tofu and soymilk output in soybean., 2008, 41: 1274-1282. (in Chinese)

Genetic analysis of panicle related traits in wheat with major gene plus polygenes mixed model

XIE SongFeng1,2, Ji WanQuan1, WANG ChangYou1, Hu WeiGuo3, Li Jun4, ZHANG YaoYuan1, SHI XiaoXi1, ZHANG JunJie1, ZHANG Hong1, Chen ChunHuan1

(1College of Agronomy, Northwest A&F University/State Key Laboratory of Crop Stress Biology in Arid Areas/Yangling Sub-centre, National Wheat Improvement Centre, Yangling 712100, Shaanxi;2Key Laboratory of Se-enriched Food Development, Ankang R&D Center for Se-enriched Prducts, Ankang 725000, Shaanxi;3Wheat Research Centre, Henan Academy of Agricultural Sciences, Zhengzhou 450002;4Crop Research Institute, Sichuan Academy of Agricultural Sciences, Chengdu 610066)

【Objective】Panicle traits are important yield traits of wheat, occupying an important position and role in wheat yield composition. Carrying out genetic research on wheat panicle traits and analyzing its genetic mechanism provide theoretical and practical guidance for formulating high-yield breeding strategies and improving breeding efficiency. 【Method】Based on the length of the main stem, the number of spikelets, the number of grains per spike, and the number of spikelets, the main gene + polygene mixed genetic model of quantitative traits was used to obtain the parental product 34 and the male parent under different ecological conditions. BARRAN and its derived F7:8, F8:9generation recombinant inbred line population (RIL) were used for genetic model analysis and genetic parameter estimation of panicle traits to determine the number of genes controlling various traits, and to estimate genetic effect values and heritability.【Result】The best genetic model for panicle length and spikelet number were B-2-1 (PG-AI), which was consistent with two pairs of linked major genes + additive-epistasis polygene genetic model. The polygenic heritability of spike length was 90.64%, the polygenic heritability of spikelet number was 89.52%, the average of environmental variation of spike length accounted for 9.39% in phenotypic variation, and the average of environmental variation of spikelet number accounted for 10.50% in phenotypic variation; Major gene heritability was 69.39%, Polygenes heritability rate was 29.94%, and the average environmental variation accounted for 2.18% in phenotypic variation. Additive effect value of the first pair of main genes controlling the number of spikes and the additive effect value of the third pair of major genes are equal, and the same was 4.56, which has a positive effect. The additive effect value of the second pair of major genes was the same as the additive effect of the first pair of major genes × the second pair of major genes × the third pair of major genes, both of which were -1.44, and are negative effects. The additive and additive × additive epistasis interaction values were equal to the additive and the second pair of major gene additions × the third pair of major gene additive epistatic interactions, both of which were -6.02. Additive and the first pair of major gene additive × the third pair of main gene additive epistatic interaction effect value is 0.18, the multi-gene additive effect value is 0.15, showing a lower positive genetic effect; H-1(4MG-AI) was best-fitting genetic model for the spikelet number traits, which showed that their inheritance was controlled by incorporating four major genes additive-epistasis genetic model. The heritability of the main gene was 81.50%. The additive effect values of the main genes in the first to fourth pairs were 0.22, 0.18, -0.20, and 0.24, respectively, the additive and epistatic interactions of the first pair of major genes × the first pair of major genes were -0.170, the additive effect value of the additive and the first pair of major genes × the third pair of major genes was 0.240. the additive effect value of the additive and the first pair of major genes × the fourth pair of major genes was -0.200, additive and the second pair of major genes × the third pair of major genes × additive effect value and additive and the second pair of major gene additive × fourth pair of major gene additive epistatic interaction value absolute value, the effect in contrast, the former value was 0.030, and the latter value was -0.030. The additive effect value of the additive and the third pair of major genes × the fourth pair of major genes was 0.060.【Conclusion】The panicle traits of wheat are mainly polygenic genetic effects, which are in line with quantitative genetic characteristics and are susceptible to environmental influences. The number of spikelet grains has the genetic characteristics of the main gene. The main gene has high heritability and is affected by the environment. The number of spikelets can be used as a direct indicator to effectively improve the early selection of panicle traits, achieving single plant directional selection and improving breeding efficiency.

wheat; panicle traits; major gene + polygene; genetic effect

2019-06-19；

接受日期：2019-08-06

国家重点研发计划（2016YFD0102004）、农业部作物基因资源与种质创制陕西科学观测实验站项目；农业部富硒产品开发与质量控制重点实验室（试运行)、陕西省农业科技创新转化项目（NYKJ-2015-037)、陕西省创新能力支撑计划（2018TD-021，2018PT-31）、富硒食品开发国家地方联合工程实验室（陕西）项目

解松峰，E-mail：xiesongfengboheng@163.com。通信作者吉万全，E-mail：jiwanquan2003@126.com

（责任编辑 李莉）