朱玉磊, 王升星, 赵良侠, 张德新, 胡建帮, 曹雪连, 杨亚杰, 常成, 马传喜, 张海萍^*
安徽农业大学农学院 / 农业部黄淮南部小麦生物学与遗传育种重点实验室, 安徽合肥230036
^*通讯作者(Corresponding author): 张海萍, E-mail:zhhp20@163.com 收稿日期:2014-07-06 基金:本研究由国家现代农业产业技术体系建设专项, 安徽省小麦产业技术体系建设专项和教育部高等学校博士学科点专项(20113418120004)资助;

摘要利用分布于小麦全基因组的181对分子标记, 分析264份自然群体的基因型, 采用TASSLE软件的GLM和MLM模型检测与整穗发芽抗性紧密关联的标记位点, 发掘相关位点内的优异等位变异。在2012年和2013年室内整穗发芽率、2013年田间自然降雨整穗发芽率3个环境中, 共关联到20个显著位点(P<0.05), 分布于小麦染色体1AS、2DS、3AS、3BL、4AL、5AS、5BL、6BS、6DS、7AL和7BL上。分别位于2DS和7BL上的分子标记gwm102和barc340同时在3个环境下关联到, 属于稳定的抗性位点; 另有6个标记位点同时在2个环境下关联到; 其余12个标记位点仅在1个环境下关联到。位于7BL上的barc340标记位点为一新报道位点。从重复关联的8个标记位点内共检测出10种优异等位变异。barc28-229bp和barc28-217bp对提高整穗发芽抗性效应最显著, 主要分布在地方品种中(如遂宁坨坨麦等), 而gwm102-142bp和barc186-199bp效应虽然相对较小, 但多分布在推广品种中(如扬麦158等), 有利于穗发芽抗性分子育种的直接应用。

关键词:小麦; 穗发芽; 休眠; 分子标记; 关联分析
Exploring Molecular Markers of Preharvest Sprouting Resistance Gene Using Wheat Intact Spikes by Association Analysis
ZHU Yu-Lei, WANG Sheng-Xing, ZHAO Liang-Xia, ZHANG De-Xin, HU Jian-Bang, CAO Xue-Lian, YANG Ya-Jie, CHANG Cheng, MA Chuan-Xi, ZHANG Hai-Ping^*
Key Laboratory of Wheat Biology and Genetic Improvement in South Yellow & Huai River Valley, Ministry of Agriculture / College of Agronomy, Anhui Agricultural University, Hefei 230036, China
Fund:
AbstractTo improve pre-harvest sprouting (PHS) resistance in wheat breeding, it is important to explore marker loci and elite alleles associated with PHS resistance using intact spikes. In this study, a total of 181 markers were used to 264 genotype materials. General and mixed linear models (GLM and MLM) were used to analyze PHS phenotypic data in three environments (2012-in house, 2013-in house and 2013-in field). The results showed that twenty markers were identified by association analysis, and located on chromosomes 1AS, 2DS, 3AS, 3BL, 4AL, 5AS, 5BL, 6BS, 6DS, 7AL, and 7BL. The markersgwm102 on 2DS andbarc340 on 7BL were detected stably in three environments, among whichbarc340 was likely to be novel and needs to be further studied through biparental linkage mapping analysis. Six markers were detected in two environments, and the other loci linked with 12 markers were detected only in one environment. A total of ten elite alleles were further explored among the eight loci with repeated associations. The allelesbarc28-229bp andbarc28-217bpfor high PHS resistance were all distributed in local cultivars (e.g. Suiningtuotuomai). However, the allelesgwm102-142bp andbarc186-199bp with intermediate PHS resistance were mainly detected in released cultivars (e.g. Yangmai 158), which could be beneficial to wheat molecular breeding.

Keyword:Triticum aestivumL.; Preharvest sprouting; Dormancy; Molecular marker; Association analysis
Show Figures
Show Figures




小麦收获前遇连阴雨天气或处于湿热环境时间较长, 易发生穗发芽(pre-harvest sprouting, PHS)。穗发芽导致种子内部储藏物质消耗, 千粒重下降和减产, 同时营养和加工品质劣化^{[ 1]}。采用准确、可靠的评价方法筛选抗穗发芽资源, 从中挖掘抗性基因, 开发功能标记, 最终培育抗穗发芽品种是解决这一问题的根本途径。
小麦收获前穗发芽抗性受多种因素共同影响, 如种皮颜色、种子的休眠特性、小麦的穗部形态(有芒或无芒等)、颖壳内发芽抑制物质、籽粒吸水特性、α-淀粉酶活性、脱落酸(ABA)和赤霉素(GA)等因素, 以及外界环境因素如光照、温度和湿度等, 其中种子的休眠性被认为是穗发芽抗性的主要遗传因素^{[ 2, 3, 4, 5]}。目前, 对小麦穗发芽抗性的评价主要有籽粒发芽法、整穗发芽法和α-淀粉酶活性测定法等。籽粒发芽法仅能反映小麦种子本身的休眠状况, 却忽略了颖壳等外部抑制物质的作用。α-淀粉酶活性测定法操作相对较复杂, 不适于批量测定。相对于前两者, 整穗发芽法操作简单, 适合较多试验材料的批量测定, 且更能全面衡量小麦穗子的整体抗性水平, 主要有发芽纸保湿法、塑料袋保湿法、纱布保湿法、模拟降雨法等^{[ 6]}。
通过解析小麦穗发芽抗性分子机制, 发现小麦穗发芽抗性受多基因控制, 存在主效基因。利用DH或RIL等人工作图群体, 控制穗发芽抗性的主效基因被定位在2AL ( Xwmc1045和 Xgwm296), 3AS ( Xfbb370、 Xbarc310、 Xbcd907、 Xbarc57、 Xbarc321、 Xbarc12和 Xbarc310), 3AL ( Xgwm155和 Xwmc153), 4AL ( Xbarc170、 Xgwm397、 Xgwm269、 Xgwm937、 Xgwm894、 Xgwm637、 Xzxq118和 Xhbe03), 5D ( Xcfd40和 Xbarc1097), 6BS ( Xwmc104)及7DL ( XMST101)上^{[ 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21]}。而基于自然群体的关联分析是定位和挖掘数量性状位点的另一条有效途径。利用这一方法, Arif等^{[ 22]}检测出34个与小麦整穗发芽性状相关的位点, 分别位于1A、2A、2D、4A、5B、5D、6A、6B、7A和7D上, 其中2D、5B和7A是携带穗发芽抗性位点的密集区域; Kulwal等^{[ 23]}除发现位于1BS、2BS、2BL、2DS、4AL、6DL和7DS上已报道的QTL外, 还鉴定出位于7BS的新位点。
本试验采用发芽纸保湿法对264份小麦材料进行整穗发芽抗性测定, 采用6%变性聚丙烯酰胺凝胶电泳对其进行基因型分析, 利用TASSLE软件的一般线性模型(GLM)和混合线性模型(MLM)对基因型和表型值进行关联分析, 检测与整穗发芽抗性紧密关联的分子标记位点, 并进一步挖掘优异等位变异及其相应载体材料, 为小麦抗穗发芽分子育种提供参考。
1 材料与方法1.1 试验材料本课题组前期收集并保存410份育种亲本材料, 利用均匀分布于小麦21条染色体上的80对多态性丰富的SSR引物对其进行遗传多样性分析, 从中筛选出亲缘关系相对较远, 且穗发芽抗性不同的264份材料作为本试验材料, 包括推广品种204份、国外引进品种4份、地方品种18份、中间品系和高代材料38份。2011年和2012年, 所有材料均秋播于安徽农业大学合肥大杨店试验站, 每份材料种2行, 行长2 m, 行宽20 cm, 每行40粒, 常规田间管理。
1.2 穗发芽抗性评价采用整穗发芽率(germination percentage, GP)来评价穗发芽抗性水平, 分为抗穗发芽、中间型和感穗发芽3个抗性等级, 其整穗发芽率分别为0~30%、31%~60%和61%~100%^{[ 6]}。
2012年和2013年分别进行室内鉴定。于小麦蜡熟后期(茎秆、叶片和麦穗全部变黄)收获每份材料10个主茎穗, 5穗为一个重复, 室内风干1 d, 立即存放于-20℃冰箱以维持种子的休眠性, 待收获全部材料后一并进行整穗发芽试验。先把整穗在灭菌水中充分浸泡10 min, 用发芽纸卷裹, 保鲜袋保湿, 置人工气候培养箱中(温度20℃, 光周期14 h昼/10 h夜, 相对湿度80%)培养7 d。于第7天取出并除去保鲜袋, 但保留发芽纸, 于电热恒温干燥箱(150℃)快速烘干, 阻止其继续发芽; 再于阳光下暴晒至籽粒易剥取; 记录每个穗子穗发芽率后取两重复平均值, 以胚部出现中破裂迹象为鉴定标准^{[ 1]}。
2013年小麦收获期遇到连阴雨天气, 于田间进行了重复鉴定。每份材料留10个生理成熟一致主茎穗直接在田间自然发芽, 1周后收回室内, 置于电热恒温干燥箱中(150℃)快速烘干, 手工脱粒, 统计每份材料发芽种子数和总粒数, 计算平均整穗发芽率。
采用SPSS19.0软件分别对3个环境下小麦整穗发芽率进行统计描述和方差分析。
1.3 分子标记筛选及其扩增条件选取15份亲缘关系较远的材料进行SSR引物筛选, 筛选标准为带型清晰、多态性丰富。引物信息来自Graingenes (http://wheat.pw.usda.gov/GG2/ index.shtm/), 由生工生物工程(上海)有限公司合成引物。最终筛选出180对SSR引物, 另外补充小麦穗发芽基因功能标记 Vp1B3^{[ 6]}, 总共181个分子标记用于关联分析。
采用SDS饱和酚法提取小麦基因组DNA。在NanoVue Plus微量分光光度计上检测DNA 浓度, 于-20℃冰箱保存。PCR扩增体系为10 μL, 包括10× buffer (含有2.0 mmol L^-1 Mg²⁺) 1.0 μL, 2.5 mmol L^-1 dNTPs 0.8 μL, 5 U μL^-1 Taq DNA polymerase 0.11 μL, 10 μmol L^-1引物各0.4 μL, 2.0 μL模板DNA (50~ 60 ng μL^-1), ddH₂O 5.29 μL。反应程序为94℃预变性5 min; 35个循环(95℃变性30 s, 48~62℃退火30 s, 72℃延伸30 s); 72℃延伸10 min; 4℃保存。用6%变性聚丙烯酰胺凝胶电泳检测扩增产物。利用软件Gelpro32 (http://www.mediacy.com/)分析带型。
1.4 关联分析参考Mackay等^{[ 24]}的方法, 从180对SSR引物中选取随机分布于小麦全基因组的105对独立SSR引物对群体遗传结构进行评估。采用STRUCTURE软件^{[ 25]}的混合模型来确定群体结构, K值取1~11, 每个 K值单独运行6次。当 K=4时曲线出现最大拐点, 因此将264份材料分为4个类群。
利用TASSLE软件^{[ 26]}的一般线性模型(general linear model, GLM)和混合线性模型(mixed linear model, MLM), 将群体结构分析所得的 Q值和亲缘关系 K作为协变量, 进行标记与性状的关联分析, 选取两种模型共同关联到的显著( P<0.05)标记位点, 并计算这些标记位点对表型变异的解释率。
参考文自翔等^{[ 27]}描述的方法, 对整穗发芽抗性显著关联的标记位点进一步挖掘其优异等位变异, 明确各等位变异的表型效应及其载体材料。
a_i =
式中, a_i代表第 i个等位变异的表型效应值, x_ij为携带第 i个等位变异在第 j份材料中表型值, n_i为具有第 i等位变异的材料数, N_k为携带无效等位变异(某一位点中, 若其中一个等位变异被视为“有效”, 则其他等位变异相应地为“无效”)的第 k份料的表型值, n_k为具有无效等位变异的材料数。若 a_i > 0, 则认为该等位变异为增效等位变异, 反之为减效等位变异。

2 结果与分析2.1 不同环境下整穗发芽抗性表型分析264份材料在3个环境中均表现出丰富的表型变异, 整穗发芽率变幅分别为0~92%、0~100%和0~96%, 其中2012年室内和2013年田间自然环境下的变异系数均在0.80以上。2013年室内总体平均值最高, 为69%, 2012年室内和2013年田间自然环境下总体平均值较低, 分别为25%和16% (表1)。
表1
Table 1
表1(Table 1)

表1 264份材料在3个环境下的整穗发芽率变异 Table 1 Variation of germination percentage of intact spike in 264 wheat genotypes detected in three environments 环境 Environment 最小值 Minimum (%) 最大值 Maximum (%) 平均值±标准差 Mean ± SD (%) 变异系数 CV 2012室内 2012-in house 0 92 25±22 0.85 2013室内 2013-in house 0 100 69±29 0.42 2013田间 2013-in field 0 96 16±13 0.81

表1 264份材料在3个环境下的整穗发芽率变异 Table 1 Variation of germination percentage of intact spike in 264 wheat genotypes detected in three environments

不同穗发芽抗性等级材料的分布在不同环境中表现不一, 2012室内和2013田间环境下分布趋势相似, 均有60%以上的材料分布在0~30%抗性等级中, 分布到61%~100%等级中的材料不满10%。2013室内频率分布不同于前两者, 其分布在0~30%和31%~60%等级中的材料均不满20%, 超过60%的材料分布在61%~100%等级中(图1)。究其原因可能是2013年小麦种子成熟期间的温度以及收获前降雨的影响, 造成整体供试材料穗发芽抗性水平下降。
方差分析(表2)显示, 3种环境下测得供试材料间整穗发芽率差异极显著( P=0.0001), 说明小麦穗发芽抗性既受基因型的控制, 也易受环境影响。供试材料的整穗发芽率在3种环境间的相关性均达极大, 为0.469; 2013室内和2013田间相关系数为0.334。综合分析上述结果说明, 整穗发芽抗性易受环境影响, 但在不同年份间, 以及同一年份内不同环境下, 不同供试材料整穗发芽率表现出相对稳定的趋势。
图1
Fig. 1

	Figure OptionViewDownloadNew Window
	图1 3个环境下不同级别穗发芽抗性材料的分布Resistant、Intermediate和Susceptible表示抗穗发芽、中间型和感穗发芽, 发芽率分别为0~30%、31%~60%和61%~100%。Fig. 1 Frequencies of materials with different resistances to preharvest spouting evaluated in three environmentsThe germination percentage of resistant, intermediate and susceptible cultivars separately was 0-30%, 31%-60%, and 61%-100%.

表2
Table 2
表2(Table 2)

表2 供试264份材料整穗发芽抗性的方差分析 Table 2 ANOVA of sprouting resistance in 264 genotypes used in this study 变异来源Source of variance 自由度 df 平方和 SS 均方 MS F值 F-value P值 P-value 基因型Genotype 263 23.24 0.09 2.41 0.0001 环境Environment 2 42.61 21.31 580.91 0.0001 误差Error 526 19.29 0.04 总变异Total 791 85.14

表2 供试264份材料整穗发芽抗性的方差分析 Table 2 ANOVA of sprouting resistance in 264 genotypes used in this study

显著水平( P<0.01)。其中, 2012室内和 2013室内相关系数最小, 为0.286; 2012室内与2013田间相关系数最
2.2 整穗发芽抗性基因定位与分子标记发掘两种模型共同检测到20个标记位点与整穗发芽抗性显著相关( P<0.05)。表3为MLM模型的检测结果。其中与2012室内关联的分子标记位点有10个, 分别位于2DS ( gwm102和 barc1146)、3AS ( barc321和 gwm369)、4AL ( wmc617)、5AS ( barc186)、5BL ( barc156)、6BS ( barc14)和7BL ( barc1181和 barc340), 各标记位点对表型贡献为2.2%~10.2%, 其中6BS上的 barc14对表型贡献最大(10.2%)。与2013室内关联的分子标记位点有11个, 分别位于染色体1AS ( barc28)、2DS ( gwm102)、3BL ( barc344和 Vp1B3)、4AL ( gwm610和 barc78)、5BL ( barc109和 barc156)、6DS ( barc173)、7AL ( barc281)和7BL ( barc340)上, 各标记位点对表型贡献为3.1%~8.0%, 其中4AL上的 barc78对表型贡献最大(8.0%)。与2013田间关联的分子标记位点有9个, 分别位于染色体1AS ( barc28)、2DS ( gwm102和 barc1146)、3AS ( barc321和 barc57)、5AS ( barc186)、5DL ( barc347)和7BL ( barc1181和 barc340)上, 各标记位点对表型贡献为3.5%~15.0%, 其中 barc321和 barc57对表型贡献相对较大(10.6%和15.0%), 均位于3AS染色体区域, 且距离较近。
表3
Table 3
表3(Table 3)

表3 混合线性模型检测到的与穗发芽抗性相关的标记位点( P<0.05)及其对表型贡献 Table 3 Marker loci associated with preharvest sprouting resistance ( P<0.05) detected by mixed linear model and their phenotypic contributions 标记 Marker 染色体 Chromosome 2012室内 2012-in house 2013室内 2013-in house 2013田间 2013-in field P R2(%) P R2(%) P R2(%) barc28 1AS — — 0.0005 7.4 0.0050 5.3 gwm102 2DS 0.0244 3.2 0.0094 4.1 0.0154 3.5 barc1146 2DS 0.0129 4.4 — — 0.0134 4.3 barc321 3AS 0.0048 5.7 — — 0.0000 10.6 barc57 3AS — — — — 0.0015 15.0 gwm369 3AS 0.0006 5.1 — — — — barc344 3BL — — 0.0048 5.3 — — Vp1B3 3BL — — 0.0187 3.5 — — wmc617 4AL 0.0443 3.4 — — gwm610 4AL — — 0.0018 5.3 — — barc78 4AL — — 0.0302 8.0 — — barc186 5AS 0.0108 2.2 — — 0.0005 4.0 barc109 5BL — — 0.0184 5.1 — — barc156 5BL 0.0345 2.3 0.0127 3.1 — — barc347 5DL — — — — 0.0121 4.4 barc14 6BS 0.0206 10.2 — — — — barc173 6DS — — 0.0020 6.0 — — barc281 7AL — — 0.0052 5.2 — — barc1181 7BL 0.0268 3.2 — — 0.0054 4.2 barc340 7BL 0.0273 3.8 0.0043 5.3 0.0107 4.4

表3 混合线性模型检测到的与穗发芽抗性相关的标记位点( P<0.05)及其对表型贡献 Table 3 Marker loci associated with preharvest sprouting resistance ( P<0.05) detected by mixed linear model and their phenotypic contributions

3种环境下均检测到与整穗发芽抗性相关的分子标记位点有2个, 分别为 gwm102和 barc340; 2种环境下均检测到的分子标记位点有6个, 分别为 barc28、 barc1146、 barc321、 barc186、 barc156和 barc1181; 其余12个分子标记位点仅在一种环境下被检测到。说明小麦整穗发芽抗性虽然存在效应稳定的QTL, 但受基因型与环境互作的显著影响, 与方差分析结果吻合。
2.3 优异等位变异及载体材料的鉴定对8个重复检测到的标记作等位变异分析, 共鉴定出10个降低整穗发芽率的优异等位变异(表4), 其载体主要为6份材料, 即茶淀红麦、遂宁坨坨麦、扬麦158、梓潼女儿麦、川麦42和扬麦20。从各等位变异的分布频率来看, barc186-199bp在供试材料中所占比例最大(62.1%, 164/264), 其次是 barc1181- 256bp (38.6%, 102/264), barc28-229bp最小(0.7%, 2/264)。从各等位变异的表型效应来看, barc28- 229bp和 barc28-217bp在3种环境下的穗发芽程度显著低于其他类型, 其来源主要是地方品种, 如遂宁坨坨麦、涪陵须须麦和梓潼女儿麦等, 在小麦育种中难以直接利用; 而广泛分布在推广品种中的 barc186-199bp和 gwm102-142bp虽然效应值相对较小, 但易于直接利用, 应该引起重视。
不同载体材料携带的等位变异及其组合不同, 如遂宁坨坨麦, 其同时携带 barc28-217bp、 barc321- 182/172bp和 barc1181-256bp三种等位变异组合; 扬麦20同时携带 barc186-199bp和 barc340-207bp两种组合; 梓潼女儿麦同时携带 barc1146-197bp和 barc340-194bp两种组合(表4)。这可能对表型起到一种聚合的效应, 从而导致其整穗发芽抗性水平较高。
进一步分析上述10种优异等位变异在不同抗性等级材料中的分布, 发现携带优异等位变异的材料大都表现抗穗发芽特性, 说明它们对整穗发芽抗性有不同程度的贡献(表4)。如携带 barc28-229bp和 barc28-217bp的材料全部表现抗穗发芽; 携带 gwm102-142bp、 barc1146-197bp和 barc321-182/ 172bp的材料中表现抗穗发芽的比率分别为65.2% (45/69)、73.3% (11/15)和66.7% (16/24)。

3 讨论业已报道, 与小麦穗发芽相关的QTL涉及多条染色体和众多区段^{[ 12, 13, 15, 16, 19, 22, 23]}。本研究采用GLM和MLM模型对264份材料不同环境下的穗发芽表型性状及其基因型进行关联分析, 20个标记位点被共同鉴定出与小麦整穗发芽抗性关联, 分别位于1AS、2DS、3AS、3BL、4AL、5AS、5BL、5DL、6BS、6DS、7AL和7BL上, 这一结果支持上述结论。但本研究鉴定的部分标记位点与相同染色体上已报道的位点距离较远, 如 barc28(1AS)与Mohan等^{[ 5]}定位的QTL(连锁标记为 Xbarc24和 Xbarc119), barc340(7BL)与Liu等^{[ 17]}鉴定的2个QTL(标记区间分别为 Xbarc255- Xgwm43和 Xwmc517- Xbarc278)都较远。本研究是否发现了新位点, 有待进一步验证。今后将着手开发多种标记(SNP、EST-SSR、EST、SCAR和CAPS等)用于加密标记图谱, 同时构建适当的作图群体, 以准确定位这些关联位点。
另有一些已报道的位点在本研究中被重复鉴定, 但其表型的贡献率在不同试验中不尽相同。如3AS上的 barc321和 barc57, 在本研究中其对表型解释效应为5.7%~15.0%, 而在张海萍等^{[ 14]}的试验中为25.6%~48.3%。本研究检测到的 wmc617 (4AL)、 gwm610 (4AL)和 barc347(5DL)与Chen等^{[ 16]}、Torada等^{[ 21]}和Munkvold等^{[ 28]}定位的QTL具有相同标记, 但表型效应较小。说明小麦抗穗发芽性状受群体、遗传背景、环境条件等影响很大, 此外, 基因型与环境互作也很明显, 本研究多数位点(12个)只能在单一环境下被检测到, 方差分析也证实这一点。
本研究检测到若干整穗发芽抗性关联位点在不同环境中较为稳定, 如 gwm102和 barc340在3种环境下均显著关联目标性状, 且与之较近的标记 barc1146和 barc1181也能在2种环境下被检测到, 与Kulwal等^{[ 8, 23]}和Liu等^{[ 13]}的研究结果相似。说明生产上有可能通过分子标记辅助选择技术改良小麦穗发芽抗性。在8个重复检测位点上, barc28-229bp、 barc28-217bp、 gwm102-142bp和 barc340-207bp具有显著提高整穗发芽抗性的作用; 前二者可能是提高整穗发芽抗性的主要优异等位变异, 但其主要由地方品种携带, 无法直接用于小麦育种, 但应重视利用其创制中间材料。另外, 推广品种中存在一些优异等位变异, 便于直接利用, 如本研究发现扬麦158、扬麦20、川麦42均携带控制低穗发芽率(0~30%)等位变异。我们认为, 育种中可以尝试利用分布在育成品种中、效应适中且稳定的等位变异或变异组合(如 gwm102-142bp、 barc156-204bp、 barc340-207bp和 barc340-194bp)来提高穗发芽抗性水平; 在多种环境下表现稳定的QTL, 虽然其效应可能不是最大, 但比在特定环境下表达的QTL更有育种价值^{[ 29]}。
表4
Table 4
表4(Table 4)

表4 与整穗发芽抗性显著关联的分子标记位点优异等位变异的表型效应和相应载体材料 Table 4 Phenotypic effects and carrier materials of elite alleles at loci strongly associated with PHS resistance 等位变异 Allele 表型效应 ai 材料数 No. of genotypes 比例 Ratio (%) 抗穗发芽 Resistant 中间型 Intermediate 感穗发芽 Susceptible 典型载体 Typical carrier E1 E2 E3 barc28-229bp -0.25 -0.42 -0.13 2 0.7 2 0 0 茶淀红麦 Chadianhongmai barc28-217bp -0.15 -0.49 -0.14 4 1.5 4 0 0 遂宁坨坨麦 Suiningtuotuomai gwm102-142bp -0.10 -0.13 -0.08 69 26.1 45 11 13 扬麦158 Yangmai 158 barc1146-197bp -0.12 -0.08 -0.11 15 5.7 11 1 3 梓潼女儿麦 Zitongnü’ermai barc321-182/172bp* -0.10 -0.04 -0.07 24 9.0 16 3 5 遂宁坨坨麦 Suiningtuotuomai barc156-204bp -0.12 -0.19 -0.06 38 14.4 26 6 6 川麦42 Chuanmai 42 barc186-199bp -0.07 -0.01 -0.08 164 62.1 94 30 40 扬麦20 Yangmai 20 barc340-207bp -0.06 -0.04 -0.05 41 15.5 24 8 9 扬麦20 Yangmai 20 barc340-194bp -0.05 -0.06 -0.03 73 27.6 43 15 15 梓潼女儿麦 Zitongnü’ermai barc1181-256bp -0.03 -0.03 -0.05 102 38.6 58 18 26 遂宁坨坨麦 Suiningtuotuomai E1: 2012-in house; E2: 2013-in house; E3: 2013-field; * indicates two alleles on one locus in a genotype. E1: 2012室内; E2: 2013室内; E3: 2013田间; * 表示同一引物在同一材料中有两种等位变异。

表4 与整穗发芽抗性显著关联的分子标记位点优异等位变异的表型效应和相应载体材料 Table 4 Phenotypic effects and carrier materials of elite alleles at loci strongly associated with PHS resistance

4 结论鉴定出与小麦整穗发芽抗性显著关联的分子标记, 证实2DS、3AS和7BL染色体区域含有稳定的抗穗发芽位点, 1AS和7BL等染色体上可能存在新位点。小麦整穗发芽抗性受基因型和环境共同影响。
The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.

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

[1]	Groos C, Gay G, Perretant R M, Bernard G M, Charmet D G. Study of the relationship between pre-harvest sprouting and grain color by quantitative trait loci analysis in a white red grain bread-wheat cross. Theor Appl Genet, 2002, 104: 39-47[本文引用:2][JCR: 3.658]
[2]	Tan M K, Sharp P J, Lu M Q, Howes N. Genetics of grain dormancy in a white wheat. Aust J Agric Res, 2006, 57: 1157-1165[本文引用:1][JCR: 1.328]
[3]	Mares D J, Mrva K. Mapping quantitative trait loci associated with variation in grain dormancy in Australian wheat. Aust J Agric Res, 2001, 52: 1257-1265[本文引用:1][JCR: 1.328]
[4]	Kottearachchi N S, Uchino N, Kato K, Miura H. Increased grain dormancy in white-grained wheat by introgression of preharvest sprouting tolerance QTLs. Euphytica, 2006, 152: 421-428[本文引用:1][JCR: 1.643]
[5]	Gatford K T, Eastwood R F, Halloran G M. Germination inhibitors in bracts surrounding the grain of Triticum tauschii. Funct Plant Biol, 2002, 29: 881-890[本文引用:2][JCR: 2.471]
[6]	Yang Y, Zhao X L, Xia L Q, Chen X M, Xia X C, Yu Z, He Z H, Roder M. Development and validation of a Viviparous-1 STS marker for preharvest sprouting tolerance in Chinese wheats. Theor Appl Genet, 2007, 115: 971-980[本文引用:3][JCR: 3.658]
[7]	Roy J K, Prasad M, Varshney R K, Balyan H S, Blake T K, Dhaliwal H S, Singh H, Edwards K J, Gupta P K. Identification of a microsatellite on chromosomes 6B and a STS on 7D of bread wheat showing an association with preharvest sprouting tolerance. Theor Appl Genet, 1999, 99: 336-340[本文引用:1][JCR: 3.658]
[8]	Kulwal P L, Kumar N, Gaur A, Khurana P, Khurana J P, Tyagi A K, Balyan H S, Gupta P K. Mapping of a major QTL for preharvest sprouting tolerance on chromosome 3A in bread wheat. Theor Appl Genet, 2005, 111: 1052-1059[本文引用:2][JCR: 3.658]
[9]	Mohan A, Kulwal P, Singh R, Kumar V, Mir R R, Kumar J, Prasad M, Balyan H S, Gupta P K. Genome-wide QTL analysis for pre-harvest sprouting tolerance in bread wheat. Euphytica, 2009, 168: 319-329[本文引用:1][JCR: 1.643]
[10]	Osa M, Kato K, Mori M, Shindo C, Torada A, Miura H. Mapping QTLs for seed dormancy and the Vp1 homologue on chromosome 3A in wheat. Theor Appl Genet, 2003, 106: 1491-1496[本文引用:1][JCR: 3.658]
[11]	Mori M, Uchino N, Chono M, Kato K, Miura H. Mapping QTLs for grain dormancy on wheat chromosome 3A and group 4 chromosomes, and their combined effect. Theor Appl Genet, 2005, 110: 1315-1323[本文引用:1][JCR: 3.658]
[12]	Liu S B, Cai S B, Robert G, Chen C X, Bai G H. Quantitative trait loci for resistance to pre-harvest sprouting in US hard white winter wheat Rio Blanco. Theor Appl Genet, 2008, 117: 691-699[本文引用:2][JCR: 3.658]
[13]	Liu S B, Bai G H. Dissection and fine mapping of a major QTL for preharvest sprouting resistance in white wheat Rio Blanco. Theor Appl Genet, 2010, 121: 1395-1404[本文引用:3][JCR: 3.658]
[14]	张海萍, 冯继明, 常成, 马传喜, 张秀英, 闫长生, 游光霞, 肖世和. 中国小麦地方品种籽粒强休眠特性的主效基因鉴定. 农业生物技术学报, 2011, 19: 270-277 Zhang H P, Feng J M, Chang C, Ma C X, Zhang X Y, Yan C S, You G X, Xiao S H. Investigation of main loci contributing to strong seed dormancy of Chinese wheat land race. J Agric Biotechnol, 2011, 19: 270-277 (in Chinese with English abstract)[本文引用:2]
[15]	朱占玲, 田宾, 刘宾, 谢全刚, 田纪春. 小麦整穗发芽的QTL定位分析. 山东农业科学, 2010, (6): 19-23 Zhu Z L, Tian B, Liu B, Xie Q G, Tian J C. Quantitative trait loci analysis for pre-harvest sprouting using intact spikes in wheat (Triticum aestivum L. ). Shand ong Agric Sci, 2010, (6): 19-23 (in Chinese with English abstract)[本文引用:2]
[16]	Chen C X, Cai S B, Bai G H. A major QTL controlling seed dormancy and pre-harvest sprouting resistance on chromosome 4A in a Chinese wheat land race. Mol Breed, 2008, 21: 351-358[本文引用:3][JCR: 3.251]
[17]	Liu S B, Bai G H, Cai S B, Chen C X. Dissection of genetic components of preharvest sprouting resistance in white wheat. Mol Breed, 2011, 27: 511-523[本文引用:2][JCR: 3.251]
[18]	Mares D J, Mrva K, Cheong J, Williams K, Watson B, Storlie E, Sutherland M, Zou Y. A QTL located on chromosome 4A associated with dormancy in white- and red grained wheats of diverse origin. Theor Appl Genet, 2005, 111: 1357-1364[本文引用:1][JCR: 3.658]
[19]	Ogbonnaya F C, Imtiaz M, Ye G, Hearnden P R, Hernand ez E, Eastwood R F, Van Ginkel M, Shorter S C, Winchester J M. Genetic and QTL analysis of seed dormancy and preharvest sprouting resistance in the wheat germplasm CN10955. Theor Appl Genet, 2008, 116: 891-902[本文引用:2][JCR: 3.658]
[20]	Zhang X Q, Li C D, Tay A, Lance R, Mares D, Cheong J, Cakir M, Ma J H, Appels R. A new PCR-based marker on chromosome 4AL for resistance to pre-harvest sprouting in wheat (Triticum aestivum L. ). Mol Breed, 2008, 22: 227-236[本文引用:1][JCR: 3.251]
[21]	Torada A, Ikeguchi S, Koike M. Mapping and validation of PCR-based markers associated with a major QTL for seed dormancy in wheat. Euphytica, 2005, 143: 251-255[本文引用:2][JCR: 1.643]
[22]	Rehman Arif M A, Neumann K, Nagel M, Kobiljski B, Lohwasser U, Börner A. An association mapping analysis of dormancy and pre-harvest sprouting in wheat. Euphytica, 2012, 188: 409-417[本文引用:2][JCR: 1.643]
[23]	Kulwal P, Ishikawa G, Benscher D, Feng Z Y, Yu L X, Jadhav A, Mehetre S, Sorrells M E. Association mapping for pre-harvest sprouting resistance in white winter wheat. Theor Appl Genet, 2012, 125: 793-805[本文引用:3][JCR: 3.658]
[24]	Mackay I, Powell W. Methods for linkage disequilibrium mapping in crops. Trends Plant Sci, 2007, 12: 57-63[本文引用:1][JCR: 11.808]
[25]	Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation on study. Mol Ecol, 2005, 14: 2611-2620[本文引用:1][JCR: 6.275]
[26]	Bradbury P J, Zhang Z W, Kroon D E, Casstevens T M, Ramdoss Y, Buckler E S. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics, 2007, 23: 2633-2635[本文引用:1][JCR: 5.323]
[27]	文自翔, 赵团结, 郑永战, 刘顺湖, 王春娥, 王芳, 盖钧镒. 中国栽培和野生大豆农艺及品质性状与SSR标记的关联分析: II. 优异等位变异的发掘. 作物学报, 2008, 34: 1339-1349 Wen Z X, Zhao T J, Zheng Y Z, Liu S H, Wang C E, Wang F, Gai J Y. Association analysis of a agronomic and quality traits with SSR markers in Glycine max and Glycine soja in China: II. Exploration of elite alleles. Acta Agron Sin, 2008, 34: 1339-1449 (in Chinese with English abstract)[本文引用:1][CJCR: 1.667]
[28]	Munkvold J D, Tanaka J, Benscher D, Sorrells M E. Mapping quantitative trait loci for preharvest sprouting resistance in white wheat. Theor Appl Genet, 2009, 119: 1223-1235[本文引用:1][JCR: 3.658]
[29]	Fulton T M, Beck-Bunn T, Emmatty D, Eshed Y, Lopez J, Petiard V, Uhlig J, Zamir D, Tanksley S D. QTL analysis of an advanced backcross of Lycopersicon peruvianum to the cultivated tomato and comparisons with QTLs found in other wild species. Theor Appl Genet, 1997, 95: 881-894[本文引用:1][JCR: 3.658]