寇程, 高欣, 李立群, 李扬, 王中华, 李学军^*
西北农林科技大学农学院 / 旱区作物逆境生物学国家重点实验室, 陕西杨凌 712100
^* 通讯作者(Corresponding author): 李学军, E-mail:xuejun@nwsuaf.edu.cn 第一作者联系方式: E-mail:koucheng7417@163.com
收稿日期:2015-04-20 接受日期:2015-07-20网络出版日期:2015-08-05基金:本研究由国家转基因生物新品种培育科技重大专项(2013ZX08002-003), 陕西省科技统筹创新工程计划项目(2014KTZB02-01-01)和陕西省重点科技创新团队计划项目(2014KCT-25)资助

摘要位于6A染色体的 TaGW2是控制小麦籽粒大小的关键基因, 已发现其第8外显子有一个T碱基插入等位变异, 其启动子区存在Hap-6A-A及Hap-6A-G等位变异。利用高分辨率熔解曲线分析技术(high resolution melting curve analysis, HRM)和 Hap-6A-P1/P2分子标记检测了316份小麦品种(系)的 TaGW2-6A基因在上述2个位点的等位变异, 分析了其不同等位变异与粒长、粒宽和千粒重的相关性, 并以大面积推广的大粒品种周麦22为例, 解析 TaGW2-6A基因优异等位变异在系谱选育中的遗传传递。共检测到61份T碱基插入等位变异(命名为977T基因型)和255份无T碱基插入的等位变异(977-基因型)材料; 在977T基因型中, Hap-6A-A (TA)和Hap-6A-G (TG)单倍型材料分别为29份和32份, 在977-基因型中, Hap-6A-A (-A)和Hap-6A-G (-G)单倍型材料分别为160份和95份。关联分析表明, 977T基因型与977-基因型的粒长( P<0.05)、粒宽( P<0.001)和千粒重( P<0.001)均有显著差异, Hap-6A-A单倍型与Hap-6A-G单倍型的粒长( P<0.05)、粒宽( P<0.05)和千粒重( P<0.001)也有显著差异。 TaGW2-6A基因编码区和启动子区等位变异之间存在相互作用, 共同调控小麦籽粒的大小, 其中TA单倍型比TG、-A、-G单倍型更能增加小麦的粒宽和粒重, 是优异的等位变异组合。周麦22为TA单倍型, 系谱分析表明, 该等位变异并非来源于亲本周8425B, 而是来源于亲本辉县红, 且TA单倍型能够稳定遗传, 但是在常规育种选择过程中可能会丢失。本研究筛选出的 TaGW2-6A优异等位变异TA单倍型材料及高通量分子检测方法为分子标记辅助育种提供材料和方法依据。

关键词:小麦籽粒性状; TaGW2基因; 等位变异; 高分辨率熔解曲线(HRM); 分子标记辅助选择
Composition and Selection of TaGW2-6A Alleles for Wheat Kernel Weight
KOU Cheng, GAO Xin, LI Li-Qun, LI Yang, WANG Zhong-Hua, LI Xue-Jun^*
College of Agronomy, Northwest A&F University / State Key Laboratory of Crop Stress Biology in Arid Areas, Yangling 712100, China

Abstract TaGW2 on chromosome 6A is a key gene governing kernel size of wheat ( Triticum aestivumL.). There are a single T-base insertion in the eighth exon of TaGW2 (977 bp) and two SNPs (Hap-6A-A and Hap-6A-G) in the promoter region. In this study, 316 wheat varieties (lines) were detected by high resolution melting curve (HRM) analysis and Hap-6A-P1/P2 molecular marker for TaGW2-6A allelic variations and their correlations with kernel length, kernel width, and thousand-kernel weight were analyzed. Furthermore, the TaGW2-6Aallelic variations were traced in the pedigree of the famous large-kernel variety Zhoumai 22. In the 977 bp position, 61 and 255 lines were detected with and without the T-base insertion, which were designated 977T and 977-genotypes, respectively. In the 977T genotypes, 29 lines were Hap-6A-A (TA) haplotype and 32 lines were Hap-6A-G (TG) haplotype. In the 977-genotypes, 160 lines were Hap-6A-A (-A) haplotype and 95 lines were Hap-6A-G (-G) haplotype. Significant difference was found in kernel length ( P < 0.05), kernel width ( P < 0.001) and thousand-kernel weight ( P < 0.001) between 977T and 977-genotypes. Similarly, significant difference was also found in kernel length ( P < 0.05), kernel width ( P < 0.05) and thousand-kernel weight ( P < 0.001) between Hap-6A-A and Hap-6A-G haplotypes. The allelic variation in TaGW2-6A encoding region and the promoter region jointly contributed to kernel size, and the TA haplotype was superior to -A, TG, and -G haplotypes in increasing kernel width and weight. According to pedigree analysis, the Zhoumai 22 inherited the TA haplotype from the parent Huixianhong, not from the popular parent Zhou 8425B. This haplotype is inheritable stably but tends to be lost in the process of wheat breeding. The results of this study provide not only a high-throughput molecular technique to detect wheat TA haplotype but also breeding materials in marker-assisted selection of wheat.

Keyword:Grain traits of wheat; TaGW2; Allelic variation; High resolution melting curve (HRM); Marker-assisted selection
Show Figures
Show Figures










小麦是全球主要的粮食作物, 占世界总耕地面积的四分之一^[1]。小麦产量由单位面积穗数、穗粒数和千粒重三因素构成, 其中千粒重在产量构成因素中受遗传特性的影响最大, 广义遗传力高达59%~80%^[2]。研究表明, 千粒重和穗粒数的增加对陕西关中地区小麦品种产量的提高有较大贡献, 其中千粒重的作用更为突出^[3]。胡延吉等^[4]提出, 随着品种产量水平的提高, 单株穗数或单位面积穗数一般没有相应增加或略有减少, 穗粒重成为进一步提高品种产量水平的主导因素。Wu等^[5]研究表明, 小麦粒重的显著增加对1945年到2010年间中国多个小麦主产区产量的持续增长有较大贡献。可见, 在品种选育过程中注重对千粒重的选择是提高小麦产量的有效方法。
GW2基因是对作物粒重有较大贡献的一个功能基因。在水稻中, OsGW2基因编码一个新型E3泛素连接酶, 通过26S蛋白酶体系统负调控种子细胞的分裂^[6]。在拟南芥中, 转AtGW2基因拟南芥籽粒千粒重明显高于野生型^[7]。在玉米中, 连锁不平衡、关联分析及基因表达分析均表明ZmGW2基因控制玉米粒重^[8]。在小麦中, Hong等^[9]通过构建TaGW2RNAi表达载体并转化到小麦品种石4185中, 发现转化后石4185籽粒的粒宽、粒厚和粒重都显著增加, 证明TaGW2基因负调控小麦籽粒大小; Yang等^[10]通过同源克隆得到小麦的粒重基因TaGW2, 发现兰考大粒中TaGW2-6A基因第8外显子977 bp处有一“ T” 碱基的插入, 关联分析表明该等位变异对小麦粒重、粒宽有显著影响; Su等^[11]克隆了TaGW2-6A基因的启动子区, 序列比对分析发现存在Hap-6A-A (-593A和-739G)和Hap-6A-G (-593G和-739A)两种等位变异型, 相关性分析表明Hap-6A-A是小麦大粒的优异等位变异。
综上, TaGW2-6A基因的编码区“ T” 插入等位变异与启动子区Hap-6A-A和Hap-6A-G等位变异均对小麦种子大小有重要影响。目前, 关于该基因编码区和启动子区等位变异的相互作用及两个区域优异等位变异基因型系谱分析的相关研究还未见报道。
本研究以316份小麦品种(系)构成的自然群体为材料, 利用HRM分析技术和Hap-6A-P1/P2分子标记检测TaGW2-6A基因编码区和启动子区的等位变异, 研究4种等位变异与小麦粒长、粒宽和粒重的关系, 同时以大面积推广的大粒品种周麦22为例分析了TaGW2-6A基因优异等位变异在品种系谱选育中的遗传传递。旨在为粒重基因TaGW2-6A优异等位变异的分子标记辅助选择育种提供材料和方法依据。
1 材料与方法1.1 供试材料及其性状测定由108份我国小麦微核心种质、174份黄淮麦区育成的小麦品种(系)及34个具有代表性的国外品种, 共计316个小麦品种(系)构成的自然群体为研究材料。在2012— 2014连续2年在西北农林科技大学北校区小麦育种试验田种植, 田间顺序排列, 人工播种; 每个品种种植1行, 行长2 m, 行距0.23 m, 田间管理同大田。
在小麦收获晾干后, 从每份材料随机取30粒完整饱满的籽粒, 按其长度方向紧密排成一行测定其长度, 再按籽粒宽度方向放置测定其宽度, 然后换算成单个籽粒的粒长和粒宽(mm), 每份材料重复测量3次, 取平均值。随机取200粒种子称重, 每份材料2次重复, 换算为千粒重(g)。
用SPSS 19.0和Microsoft Excel 2007计算和统计分析调查数据, 用LSR (least significant ranges)方法比较基因型间的差异显著性。
1.2 TaGW2-6A基因等位变异的分子检测每个品种取0.2~0.3 g新鲜叶片, 液氮预冷后粉碎放入1.5 mL离心管, 按照Gawel等^[12]的方法提取小麦基因组DNA。
1.2.1 TaGW2-6A基因编码区“ T” 碱基插入等位变异
根据兰考大粒TaGW2-6A基因T插入位点附近序列设计引物HRM-P (F: 5′ -ATGGGTGCTGCGGA AAGT-3′ ; R: 5′ -CGCTCCAGCTATCTGGTGAA-3′ ), 扩增片段长度为38 bp。引物由宝生物工程(大连)有限公司合成, 纯化级别为PAGE。用Light Cycler 480 (Roche)进行PCR扩增及HRM分析, 反应体系为20 μ L, 含80 ng模板DNA, 上下游引物各0.4 µ mol L^-1, 10 μ L Light Cycler HRM Master Mix (Roche), 1.5 mol L^-1 Mg²⁺ (Roche); 反应程序, 首先进行普通PCR (95℃ 5 min; 95℃ 30 s, 58℃ 30 s, 72℃ 2 s, 40个循环; 72℃ 10 min), 再对其产物进行HRM分析, 即以95℃ 30 s使荧光染料结合到碱基上, 再降温至59℃ 30 s使解开的DNA单链重新结合成双链结构, 随后逐渐升温到95℃, 每升高1℃收集25次荧光信号, 最后迅速降到37℃, 保持10 s, 完成整个HRM分析反应。
1.2.2 HRM-P引物扩增产物序列分析 中国春和兰考大粒的HRM反应产物经3%琼脂糖凝胶电泳分离, 将目标片段挖胶回收, 用pUCm-T载体克隆, 由生工生物工程(上海)股份有限公司测序。每个品种至少测定10个克隆。采用DNAMAN软件分析引物的扩增产物序列等位变异, 并在NCBI数据库(http://www.ncbi.nlm.nih.gov/)上进行比对。
1.2.3 TaGW2-6A基因启动子区等位变异 利用Su等^[12]设计的Hap-6A-P1/P2分子标记检测316份供试材料中TaGW2-6A基因启动子区Hap-6A-A和Hap-6A-G等位变异。
1.3 周麦22的系谱分析以同时含有两种优异等位变异(编码区为T插入型; 启动子区为Hap-6A-A单倍型)的大粒品种周麦22为例, 通过HRM分析技术及Hap-6A-P1/P2分子标记检测周麦22系谱中各亲本TaGW2-6A基因编码区及启动子区等位变异的组成, 进而分析TaGW2- 6A基因T碱基插入等位变异以及启动子区Hap-6A- A和Hap-6A-G等位变异的遗传传递规律。周麦22于2007年通过国家审定, 截至2014年夏收, 在黄淮麦区累计推广面积460多万公顷, 是大粒高产品种的代表。

2 结果与分析2.1 供试材料TaGW2-6A基因等位变异的组成分析 HRM分析结果显示HRM-P引物能够明显的区分出兰考大粒与中国春高分辨率熔解曲线的差异。目的片段扩增产物测序显示, 中国春扩增产物只有一种序列, 而兰考大粒扩增产物存在两种序列, 一种序列与中国春的序列完全相同, 推测其为6B、6D基因组TaGW2扩增的产物; 另一种序列比中国春多一个T碱基插入(图1), 将其在NCBI数据库进行比对后发现与TaGW2-6A-LK目的序列完全相同, 由此判断HRM分析技术显示出的熔解曲线的差异是由TaGW2-6A基因T碱基插入引起的。
图1
Fig. 1

	Figure OptionViewDownloadNew Window
	图1 HRM-P引物在中国春和兰考大粒中扩增产物的序列比对Fig. 1 Sequence alignment of PCR products in Chinese Spring and Lankaodali amplified with HRM-P primer

316份材料中有61份材料的熔解曲线呈现与兰考大粒的熔解曲线相近, 255份材料的熔解曲线呈现与中国春的熔解曲线相近(图2), 由此得出供试材料中T碱基插入材料为61份(977T基因型), 无T碱基插入的材料为255份(977-基因型)。
图2
Fig. 2

	Figure OptionViewDownloadNew Window
	图2 977T和977-基因型的高分辨率熔解曲线分析 红色、蓝色、绿色、黄色曲线分别表示兰考大粒、中国春、977T基因型和977-基因型熔解曲线。Fig. 2 High-resolution melting curve analysis of 977T and 977-genotypes The melting curves of Lankaodali, Chinese Spring, 977T genotype, and 977-genotype are in red, blue, green, and yellow, respectively.

用Hap-6A-P1/P2引物进行PCR扩增, 产物在1%琼脂糖凝胶上电泳检测, 区分启动子区Hap-6A-A (167 bp)和Hap-6A-G (218 bp)两种等位变异材料(图3), 结果Hap-6A-A等位变异材料189份, Hap-6A-G等位变异材料127份。
图3
Fig. 3

Figure OptionViewDownloadNew Window

图3 TaGW2-6A基因启动子区的多态性分析 Hap-6A-A和Hap-6A-G等位变异的扩增片段长度分别为167 bp和218 bp。M: DL2000 marker; 1: 周麦22; 2: 西农2208; 3: 陕229; 4: 蚂蚱麦; 5: 中麦349; 6: 小偃22; 7: 藁优9618。Fig. 3 Polymorphism analysis on promoter region of TaGW2-6A The sizes of amplified fragments were 167 bp for genotype Hap-6A-A and 218 bp for genotype Hap-6A-G. M: DL2000 marker; 1: Zhoumai 22; 2: Xinong 2208; 3: Shaan 229; 4: Mazhamai; 5: Zhongmai 349; 6: Xiaoyan 22; 7: Gaoyou 9618.

根据HRM和Hap-6A-P1/P2分子标记分析结果, 将316份材料分成TA、TG、-A和-G四种单倍型, 分别为29、32、160和95份。表1列出了主要品种的相关信息。
表1
Table 1
表1(Table 1)

表1 部分重要品种名称及基因型表型数据 Table 1 Part of important varieties and their genotype-phenotype data 品种 Variety 单倍型 Haplotype 粒长 KL (mm) 粒宽 KW (mm) 粒重 TKW (g) 来源 Origin 育成年份 Year released 兰考大粒 Lankaodali TA 8.05 3.62 51.25 中国河南 Henan, China — 西农6028 Xinong 6028 TA 6.50 3.46 52.00 中国陕西 Shaanxi, China 1948 豫麦49 Yumai 49 TA 6.74 3.52 48.94 中国河南 Henan, China 1998 水原86 Suweon 86 TA 6.76 3.58 48.07 韩国 Public of Korea — 辉县红 Huixianhong TA 6.53 3.63 42.15 中国河南 Henan, China — 陕213 Shaan 213 TA 7.30 3.22 40.98 中国陕西 Shaanxi, China 1986 加麦5号Jiamai 5 TA 6.41 3.54 48.33 加拿大 Canada — 周麦20 Zhoumai 20 TA 6.61 3.63 56.29 中国河南 Henan, China 2003 绵7198-8 Mian 7198-8 TA 6.64 3.50 49.66 中国四川 Sichuan, China — 西农1376 Xinong 1376 TA 6.84 3.44 40.61 中国陕西 Shaanxi, China 1988 西农2208 Xinong 2208 TA 7.57 3.16 40.87 中国陕西 Shaanxi, China 1997 周8425B Zhou 8425B TA 7.45 3.76 52.68 中国河南 Henan, China 1988 周麦22 Zhoumai 22 TA 6.54 3.46 51.15 中国河南 Henan, China 2007 蚂蚱麦 Mazhamai TA 6.51 3.72 36.83 中国陕西 Shaanxi, China 1930 豫麦2号 Yumai 2 TA 6.12 3.23 37.31 中国河南 Henan, China 1977 晋麦2418 Jinmai 2418 TA 6.69 3.07 43.87 中国福建 Fujian, China — 百农3217 Bainong 3217 TA 6.70 3.02 35.16 中国河南 Henan, China 1981 冀麦2号 Jimai 2 TA 7.00 3.28 46.73 中国河北 Hebei, China — 陕农7859 Shaannong 7859 TA 7.57 3.31 47.58 中国陕西 Shaanxi, China 1984 康定小麦 Kangding Xiaomai TA 6.21 3.18 38.45 中国四川 Sichuan, China — 白蚂蚱 Baimazha TA 5.47 2.93 31.95 中国宁夏 Ningxia, China — 青春28 Qingchun 28 TA 6.98 3.46 49.29 中国青海 Qinghai, China 1970 毕麦26 Bimai 26 TA 7.13 3.27 40.62 中国贵州 Guizhou, China 1979 去麦34 Qumai 34 TA 6.59 3.13 44.70 中国云南 Yunnan, China — 同家坝小麦 Tongjiaba Xiaomai TA 6.21 2.99 31.25 中国四川 Sichuan, China — 红花麦 Honghuamai TA 6.31 3.23 34.98 中国四川 Sichuan, China — 成都光头 Chengduguangtou TA 6.34 2.98 37.30 中国四川 Sichuan, China 1952 酱麦 Jiangmai TA 6.50 3.26 36.81 中国贵州 Guizhou, China — 汉中白 Hanzhongbai TA 6.21 2.93 32.30 中国陕西 Shaanxi, China — 陕253 Shaan 253 TG 7.43 3.70 42.95 中国陕西 Shaanxi, China 1994 西农2611 Xinong 2611 TG 7.21 3.35 43.76 中国陕西 Shaanxi, China 1991 武农886 Wunong 886 TG 6.38 3.33 51.75 中国陕西 Shaanxi, China 2009 邯6172 Han 6172 TG 6.69 3.25 47.77 中国河北 Hebei, China 1995 农林10号 Norin 10 TG 7.47 3.73 49.54 日本 Japan — 临旱957 Linhan 957 TG 7.37 3.60 51.65 中国山西 Shanxi, China — 藁优9618 Gaoyou 9618 TG 6.74 3.10 36.27 中国河北 Hebei, China 2005 扬麦5号 Yangmai 5 TG 6.94 3.00 36.69 中国江苏 Jiangsu, China 1991 西农9814 Xinong 9814 TG 6.65 3.41 52.70 中国陕西 Shaanxi, China — 胜利麦 Triumph TG 6.21 2.77 37.83 美国 United States — 中农28 Villa Glori TG 6.64 2.97 33.10 意大利 Italy — 高原506 Gaoyuan 506 TG 6.93 3.34 46.05 中国青海 Qinghai, China 1971 宁春4号 Ningchun 4 -A 7.12 3.46 40.79 中国宁夏 Ningxia, China 1981 中麦895 Zhongmai 895 -A 7.17 3.11 47.87 中国河北 Hebei, China 2012 天麦989 Tianmai 989 -A 6.92 3.31 44.23 中国陕西 Shaanxi, China 2012 周麦12 Zhoumai 12 -A 5.99 2.99 34.50 中国河南 Henan, China 1999 周麦13 Zhoumai 13 -A 6.39 3.26 47.93 中国河南 Henan, China — 周麦9号 Zhoumai 9 -A 5.53 2.90 38.93 中国河南 Henan, China 1980 温394 Wen 394 -A 6.96 3.15 32.40 中国河南 Henan, China — 碧蚂6号 Bima 6 -A 5.24 3.14 36.32 中国陕西 Shaanxi, China 1948 矮丰3号 Aifeng 3 -A 6.45 3.03 33.70 中国陕西 Shaanxi, China 1970 咸农39 Xiannong 39 -A 7.01 3.09 37.44 中国陕西 Shaanxi, China 1964 小偃6号 Xiaoyan 6 -A 6.63 3.38 43.09 中国陕西 Shaanxi, China 1981 小偃22 Xiaoyan 22 -A 6.97 3.28 38.86 中国陕西 Shaanxi, China 2003 郑麦9023 Zhengmai 9023 -A 7.12 3.15 46.36 中国河南 Henan, China 2009 中麦349 Zhongmai 349 -A 6.62 2.79 42.85 中国河北 Hebei, China 2009 新麦9号 Xinmai 9 -A 6.49 3.26 42.99 中国河南 Henan, China 1997 存麦1号 Cunmai 1 -A 6.76 3.38 44.25 中国河南 Henan, China 2011 浚麦55 Xunmai 55 -A 6.27 3.06 46.70 中国河南 Henan, China — 淮麦20 Huaimai 20 -A 7.27 3.46 44.07 中国江苏 Jiangsu, China 2003 丰产3号 Fengchan 3 -G 6.34 3.13 36.86 中国陕西 Shaanxi, China 1964 鲁麦1号 Lumai 1 -G 6.19 3.21 38.13 中国山东 Shandong, China 1975 Data are the means over two years (2012-2013 and 2013-2014 wheat seasons). “ — ” indicates released year unknown. 数据为2012-2013和2013-2014两年度的平均值。— 表示育成时间未知。

表1 部分重要品种名称及基因型表型数据 Table 1 Part of important varieties and their genotype-phenotype data

2.2 TaGW2-6A基因等位变异与籽粒表型相关性分析供试材料粒重表型数据与TaGW2-6A等位变异关联分析结果表明, 977T基因型与977-基因型材料的粒长(P< 0.05)、粒宽(P< 0.001)和千粒重(P< 0.001)差异显著, 而Hap-6A-A与Hap-6A-G基因型之间粒长(P< 0.05)、粒宽(P< 0.05)和千粒重(P< 0.001)也有显著差异。TA单倍型与TG单倍型材料的粒长、粒宽和千粒重差异不显著, 其中TA单倍型与-A单倍型(P< 0.05)和-G单倍型(P< 0.001)存在显著差异; TG单倍型与-A单倍型的3个籽粒性状无显著差异, 但二者都与-G单倍型(P< 0.05)呈显著差异(表2)。可见, TA单倍型比TG、-A、-G单倍型更能增加小麦的粒宽和粒重, 是优异的等位变异组合。
2.3 周麦22的系谱分析根据分子标记检测结果, 周麦22为TA单倍型, 而其系谱中相关亲本在TaGW2-6A基因编码区和启动子区的等位变异呈多样性(图4)。周麦22是以周麦12/豫麦49为母本、周麦13为父本杂交选育而成的, 其中周麦12和周麦13是-A单倍型材料; 周麦13又是由周麦骨干亲本周8425B和周麦9号杂交选育而成, 其中周8425B是TA单倍型; 豫麦49遗传了周麦22的1/4血统, 而豫麦49是TA单倍型。可见, 周麦22的TA单倍型并非来源于周麦13的亲本周8425B, 而是来源于豫麦49, 并通过系谱向上追溯到豫麦2号。豫麦2号的一个亲本65(14) 3是矮丰3号(-A单倍型)的姊妹系, 另一个亲本是辉县红(TA单倍型), 因此推定, 豫麦2号的TA单倍型来源于小麦地方品种辉县红。可见, TaGW2-6A基因的TA单倍型能够在品种选育过程中稳定遗传, 但如果选择不当也可能丢失, 如蚂蚱麦的TA单倍型没有传递到碧蚂6号, 周8425B的TA单倍型没有传递到周麦13 (图4)。
表2
Table 2
表2(Table 2)

表2 TaGW2-6A基因等位变异的表型关联分析 Table 2 Correlation analysis of TaGW2-6A polymorphism and phenotype 基因型/单倍型 Genotype/haplotype 粒长 Kernel length (mm) 粒宽 Kernel width (mm) 千粒重 1000-kernel weight (g) 2012-2013 977T 6.67± 0.07 aA 3.30± 0.03 aA 42.40± 0.82 aA 977- 6.41± 0.03 bA 3.13± 0.02 bB 38.07± 0.40 bB Hap-6A-A 6.52± 0.04 aA 3.19± 0.02 aA 40.46± 0.47 aA Hap-6A-G 6.39± 0.05 bA 3.11± 0.02 bA 36.60± 0.57 bB TA 6.73± 0.10 aA 3.36± 0.05 aA 43.43± 1.19 aA -A 6.48± 0.04 bA 3.17± 0.02 bB 39.96± 0.51 bB TG 6.62± 0.10 abA 3.25± 0.05 abAB 41.46± 1.13 abAB -G 6.31± 0.06 cC 3.07± 0.03 cB 34.96± 0.66 cC 2013-2014 977T 6.61± 0.07 aA 3.26± 0.03 aA 41.66± 0.80 aA 977- 6.35± 0.03 bA 3.10± 0.02 bB 37.72± 0.39 bB Hap-6A-A 6.45± 0.04 aA 3.16± 0.02 aA 39.95± 0.45 aA Hap-6A-G 6.32± 0.05 bA 3.09± 0.02 bA 36.29± 0.55 bB TA 6.68± 0.10 aA 3.29± 0.05 aA 42.64± 1.16 aA -A 6.41± 0.04 bA 3.14± 0.02 bA 39.46± 0.49 bA TG 6.54± 0.10 abA 3.23± 0.05 abA 40.77± 1.10 abA -G 6.24± 0.06 cB 3.04± 0.03 cB 34.78± 0.64 cB In each growing year, values followed by different capitals and lowercases are significantly different between TaGW2-6A genotypes or among TaGW2-6A haplotypes at P < 0.001 and P < 0.05, respectively. 同一年度中, 不同大写和小写字母分别表示TaGW2-6A位点2种基因型之间或启动子区4种单倍型之间在0.001和0.05概率水平下差异显著。

表2 TaGW2-6A基因等位变异的表型关联分析 Table 2 Correlation analysis of TaGW2-6A polymorphism and phenotype

图4
Fig. 4

	Figure OptionViewDownloadNew Window
	图4 周麦22系谱 ●: TA单倍型; ○: -A单倍型; △ : -G单倍型; ?: 没有进行分子检测。Fig. 4 Pedigree of Zhoumai 22 ●: TA genotype; ○: -A genotype; △ : -G genotype; ?: not detected.

3 讨论在TaGW2-6A基因T插入等位变异的品种筛选过程中, Yang等^[10]利用AS-PCR技术来筛选T基因型材料, Liu等^[13]针对特异等位位点设计了多对AS-PCR引物, 发现该类引物在扩增过程中无法达到稳定检测特异等位位点的目的。我们在利用AS-PCR技术检测的过程也发现其重复性差的问题。目前为止, 小麦粒重性状相关的分子标记被开发的较多, 但能够直接用于辅助育种的很少^[14]。高分辨率熔解曲线分析(high resolution melting curve analysis, HRM)是Gundry等^[15]于2003年提出并最初应用于SNP检测分析的一项新技术, 具有灵敏度高、特异性强、速度快、通量高、操作简便等优点, 近些年已经用于标记辅助选择、突变检测及种质鉴定等研究^[16]。例如, 利用该技术, Ganopoulos等^[17]分析菜豆SSR标记检测结果, 实现了菜豆品种的快速准确鉴定; Jeong等^[18]利用COSII标记和Waxy基因的SNP标记数据, 对辣椒种质资源进行了鉴定; Han等^[19]利用四倍体苜蓿的多个SNP位点数据, 对其亲本及后代群体进行SNP分型; Hofinger等^[20]将HRM技术用于大麦eIF4E等位变异筛选。本研究将HRM技术引入小麦TaGW2-6A基因等位变异分析, 我们设计的HRM-P引物组合能有效区分两种基因型, 准确率达到100%且重复性良好, 为TaGW2-6A基因编码区优异等位变异的分子标记辅助选择提供了高通量的分子检测方法。
Yang等^[10]以兰考大粒与中国春为亲本构建F_2:3群体, 通过TaGW2-6A基因型与表型的关联分析, 证实该基因位点T插入突变与籽粒的表型相关, 能够显著增加粒宽和粒重。Su等^[11]利用自然群体进行关联分析, 发现Hap-6A-A相对于Hap-6A-G是优异等位变异类型。本研究进一步表明, TaGW2-6A基因TA单倍型的粒宽和粒重大于TG、-A和-G单倍型, TG和-A单倍型的粒宽和粒重显著大于-G单倍型, 而TG单倍型与-A单倍型对粒宽、粒重的作用不显著, 说明TaGW2-6A基因TA单倍型是优异的等位变异组合。然而, 在977T基因型中也出现如藁优9618、中农28、定西24等粒重较小的材料, 其TaGW2-6A基因启动子区为Hap-6A-G等位变异; 同时, 在977-基因型中也存在周麦18、中麦895等籽粒较大的材料, 其TaGW2-6A基因启动子区为Hap-6A-A等位变异。根据本研究结果推测, TaGW2-6A基因编码区和启动子区等位变异之间存在相互作用, 共同调控TaGW2-6A基因表达, 从而控制小麦籽粒的大小和粒重。
本研究利用HRM分析技术以及Hap-6A-P1/P2分子标记, 在蚂蚱麦、西农6028、辉县红、周8425B、豫麦49、西农1376、豫麦2号、周麦22等黄淮麦区的骨干亲本中均检测到TA单倍型, 说明本研究采用的技术策略可用于小麦TaGW2-6A优异等位变异组合的检测和筛选。另一方面, 利用该方法我们进行了周麦22 TA单倍型的系谱追踪, 发现TA单倍型是从辉县红、豫麦2号、豫麦49传递到周麦22中, 而周麦22的其他亲本, 如蚂蚱麦、周8425B, 其TA单倍型没能传递给后代, 说明育种过程中如果选择不当, 将会导致优异基因丢失。本研究采用的分子标记检测技术和方法, 可有效监测小麦遗传改良中TaGW2-6A基因优异等位变异及其单倍型的传递过程, 提高传统育种的目标性和育种速度, 为加速培育突破性新品种提供技术支撑。
4 结论利用HRM分析技术和Hap-6A-P1/P2分子标记对316份小麦品种(系)进行了检测, 发现约一半品种(160份, 50.6%)为-A单倍型, 其次为-G单倍型(95份, 30.1%), TA和TG单倍型品种较少, 分别为29份(9.2%)和32份(10.1%), 尤其是优异等位变异组合TA单倍型的频率最低。TaGW2-6A基因编码区和启动子区等位变异之间存在相互作用, 共同调控TaGW2-6A基因从而控制小麦籽粒的大小。
The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.


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

[1]	Fox J L. Whatever happened to GM wheat?Nat Biotechnol, 2009, 27: 974-976[本文引用:1]
[2]	庄巧生. 中国小麦品种改良及系谱分析. 北京: 中国农业出版社, 2003. pp 497-518 Zhuang Q S. Analysis of Chinese Wheat Impronement and Pedigree. Beijing: China Agriculture Press, 2003. pp 497-518(in Chinese)[本文引用:1]
[3]	许为刚, 胡琳, 吴兆苏, 盖钧镒. 关中地区小麦品种产量与产量结构遗传改良的研究. 作物学报, 2000, 26: 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. Acta Agron Sin, 2000, 26: 352-358 (in Chinese with English abstract)[本文引用:1]
[4]	胡延吉, 赵檀方. 小麦高产育种中粒重作用的研究. 作物学报, 1995, 21: 671-678 Hu Y J, Zhao T F. Studies on the effect of grain weight in breeding of high-yield wheat. Acta Agron Sin, 1995, 21: 671-678 (in Chinese with English abstract)[本文引用:1]
[5]	Wu W, Li C J, Ma B L, Shah F, Liu Y, Liao Y C. Genetic progress in wheat yield and associated traits in China since 1945 and future prospects. Euphytica, 2013, 196: 155-168[本文引用:1]
[6]	Song X J, Huang W, Shi M, Zhu M Z, Lin H X. A QTL for rice grain width and weight encodes a previously unknown RING-type E3 ubiquitin ligase. Nat Genet, 2007, 39: 623-630[本文引用:1]
[7]	姜颖, 曹言勇, 路运才, 唐保军, 王利锋, 李会勇. 拟南芥中RING型E3泛素连接酶基因AtGW2的克隆和功能分析. 植物遗传资源学报, 2011, 12: 448-454 Jiang Y, Cao Y Y, Lu Y C, Tang B J, Wang L F, Li H Y. Cloning and functional analysis of Arabidopsis thaliana AtGW2, a RING-type E3 ubiquitin ligase protein. J Plant Genet Resours, 2011, 12: 448-454 (in Chinese with English abstract)[本文引用:1]
[8]	Li Q, Li L, Yang X H, Warburton M L, Bai G H, Dai J R, Li J S, Yan J B. Relationship, evolutionary fate and function of two maize co-orthologs of rice GW2 associated with kernel size and weight. BMC Plant Biol, 2010, 10: 1-15[本文引用:1]
[9]	Hong Y T, Chen L F, Du L P, Su Z Q, Wang J F, Ye X G, Li Q, Zhang Z Y. Transcript suppression of TaGW2 increased grain width and weight in bread wheat. Funct Integr Genomics, 2014, 14: 341-349[本文引用:1]
[10]	Yang Z B, Bai Z Y, Li X L, Wang P, Wu Q X, Yang L, Li L Q, Li X J. SNP identification and allelic-specific PCR markers development for TaGW2, a gene linked to wheat kernel weight. Theor Appl Genet, 2012, 125: 1057-1068[本文引用:3]
[11]	Su Z Q, Hao C Y, Wang L F, Dong Y C, Zhang X Y. Identification and development of a functional marker of TaGW2 associated with grain weight in bread wheat (Triticum aestivum L. ). Theor Appl Genet, 2011, 122: 211-223[本文引用:2]
[12]	Gawel N J, Jarret R L. A modified CTAB DNA extraction procedure for Musa and Ipomoea. Plant Mol Biol Rep, 1991, 9: 262-266[本文引用:2]
[13]	Liu J, Huang S M, Sun M Y, Liu S Y, Liu Y M, Wang W X, Zhang X R, Wang H Z, Hua W. An improved allele-specific PCR primer design method for SNP marker analysis and its application. Plant Methods, 2012, 8: 1-9[本文引用:1]
[14]	Xu Y B, Crouch J H. Marker-assisted selection in plant breeding: from publications to practice. Crop Sci, 2008, 48: 391-407[本文引用:1]
[15]	Gundry C N, Vand ersteen J G, Reed G H, Pryor R J, Chen J, Wittwer C T. Amplicon melting analysis with labeled primers: a closed-tube method for differentiating homozygotes and heterozygotes. Mol Diagnostics and Genet, 2003, 49: 396-406[本文引用:1]
[16]	甘四明, 李发根, 翁启杰, 李梅, 周长品, 于晓丽. 高分辨率熔解及其在植物DNA变异检测中的应用. 基因组学与应用生物学(网络版), 2010, 29: 5. DOI: DOI: 105376/gab. cn. 201029. 0005. http://zwxb.chinacrops.org/article/2015/%20//editors.%20wanfangdata.%20com.%20cn/openaccessperiodical/Details/6ef1e733-8363-4a26-b115-5d1482d20353 Gan S M, Li F G, Wong Q J, Li M, Zhou C P, Yu X L. High-resolution melting and its applications in plant dna variation detection. Genomics Appl Biol (Online), 2010, 29: 5. DOI: DOI:10.5376/gab.cn.2010.29.0005. http://editors.wanfangdata.com.cn/openaccessperiodical/Details/6ef1e733-8363-4a26-b115-5d1482d20353(in Chinese with English abstract)[本文引用:1]
[17]	Ganopoulos I, Bosmali I, Madesis P, Tsaftaris A. Microsatellite genotyping with HRM (high resolution melting) analysis for identification of the PGI common bean variety Plake Megalosperma Prespon. Eur Food Res Technol, 2012, 234: 501-508[本文引用:1]
[18]	Jeong H J, Jo Y D, Park S W, Kang B C. Identification of Capsicum species using SNP markers based on high resolution melting analysis. Genome, 2010, 53: 1029-1040[本文引用:1]
[19]	Han Y H, Khu D M, Monteros M J. High-resolution melting analysis for SNP genotyping and mapping in tetraploid alfalfa (Medicago sativa L. ). Mol Breed, 2011, 29: 489-501[本文引用:1]
[20]	Hofinger B J, Jing H C, Hammond-Kosack K E, Kanyuka K. High-resolution melting analysis of cDNA-derived PCR amplicons for rapid and cost-effective identification of novel alleles in barley. Theor Appl Genet, 2009, 119: 851-865[本文引用:1]