利用普通六倍体小麦和西藏半野生小麦杂交衍生的重组自交系定位小麦芒长QTL
宫希1, 蒋云峰2, 徐彬杰2, 乔媛媛2, 华诗雨1, 吴旺1, 马建2, 周小鸿2, 祁鹏飞2,*, 兰秀锦2
1四川农业大学农学院, 四川成都 611130
2四川农业大学小麦研究所, 四川成都 611130
*通讯作者(Corresponding author): 祁鹏飞, E-mail: pengfeiqi@hotmail.com, Tel: 028-82650337

第一作者联系方式: E-mail: 360035030@qq.com

摘要

芒长是普通小麦的重要农艺性状, 受多个基因控制。本研究利用长芒的普通小麦郑麦9023与无芒的西藏半野生小麦Q1028构建一个重组自交系群体(186个株系); 采用SSR和DArT分子标记, 构建覆盖小麦全基因组的遗传图谱(2597 cM)。基于重组自交系群体两年芒长表型数据, 采用ICIM作图法对小麦芒长性状进行QTL定位分析。共检测到2个与芒长相关的QTL, 即 Qwa.sau-4AS Qwa.sau-5AL。它们分别位于4AS和5AL染色体上, 可分别解释7.4%和27.3%的表型变异。这2个QTL效应可能分别来源于钩芒基因 Hd与抑芒基因 B1。利用连锁标记进行基因型分析, 表明 Qwa.sau-4AS Qwa.sau-5AL对芒长的抑制效果具有累加效应, 且 Qwa.sau-5AL效应强于 Qwa.sau-4AS。本研究将为精细定位及克隆这2个QTL奠定基础。

关键词: 普通小麦; 芒长; QTL; SSR标记; DArT标记
Mapping QTLs for Awn Length in Recombinant Inbred Line Population Derived from the Cross between Common Wheat and Tibetan Semi-wild Wheat
GONG Xi1, JIANG Yun-Feng2, XU Bin-Jie2, QIAO Yuan-Yuan2, HUA Shi-Yu1, WU Wang1, MA Jian2, ZHOU Xiao-Hong2, QI Peng-Fei2,*, LAN Xiu-Jin2
1College of Agronomy, Sichuan Agricultural University, Chengdu 611130, China
2Triticeae Research Institute, Sichuan Agricultural University, Chengdu, 611130, China
Abstract

Awn length is an important agronomic trait in hexaploid wheat and controlled by multiple genes. A recombinant inbred line (RIL) population containing 186 lines were developed by crossing Tibetan semi-wild wheat accession Q1028 (female) and common wheat variety Zhengmai 9023, and a genome-wide genetic map (2597 cM) was constructed using SSR and DArT markers. QTLs controlling awn length were identified with two-year phenotypic data and the genetic map by using the Inclusive Composite Interval Mapping (ICIM) method. Two QTLs associated with awn length, designated Qwa.sau-4AS and Qwa.sau-5AL, were detected on chromosomes 4A and 5A, which explained 7.4% and 27.3% of phenotypic variations, respectively. According to genetic locations, effects of the two QTLs might be from Hdand B1 genes that could reduce awn length, respectively. Genetic analysis showed that Qwa.sau-4ASand Qwa.sau-5ALhad cumulative effects, with stronger inhibiting effect in Qwa.sau-5AL than in Qwa.sau-4AS. These results are valuable for fine mapping and cloning target genes controlling awn length in the future.

Keyword: Common wheat; Awn length; QTL; SSR marker; DArT marker

小麦有芒品种通常比无芒品种的产量高10%以上[1, 2], 因此芒也成为高产育种的考察性状之一。芒作为小麦穗器官的组成部分, 是植物长期进化、适应环境的结果[3], 是小麦光合作用和蒸腾作用的重要器官。芒的绿色细胞富含叶绿体与线粒体, 含量较颖片和外稃更高, 其形态和结构与叶片相似, 具有光合功能[4], 在小麦灌浆后期, 其光合活性高于旗叶[5]。叶片干枯后, 几天之内水分仍可运输到芒, 继续进行光合作用, 光合产物就近转运, 增加籽粒干物质积累[6, 7]

芒在逆境条件下有助于提高产量[8, 9]。小麦的芒有较多气孔和较高的蒸腾作用, 在干旱条件下可加快水分向麦穗的运输, 使其吸收更多养分, 保证籽粒发育。此外, 小麦芒在种子的传播, 抗御虫害以及防止脱粒等方面也具有重要的作用。白天空气干燥时芒保持弯曲, 而夜晚空气湿度增大引起芒的纤维排列变化, 使芒变直, 每日循环的空气湿度成为芒推动种子进入地面的能量来源, 证明芒在小麦繁殖进化的过程中起了重要的作用[10]。芒表面有坚硬的硅质毛, 不利于害虫的飞落, 降低害虫在麦穗上产卵的机会, 有利于抗御鸟雀和虫害[11]。芒的弹性又可防止穗间相互撞击, 避免籽粒脱落。

小麦芒长由多个基因控制, 其中效应较强且报道最多的是3个基因位点B1B2HdB1B2是芒抑制基因, 其中B1抑制作用最强, 可导致穗基部和中部芒极短或无芒[9, 12]B1位于5AL染色体末端, 与标记Xgwm291的遗传距离为4.99 cM [9]B2会造成整个穗部的短芒, 它被定位在中国春6BL-5和6BL-6缺失片段末端之间, 与分子标记Xwmc539Xgpw5130Xwmc748临近[13]; Hd基因控制钩芒的产生[14], 也会对芒长产生抑制作用, 被定位于4AS染色体Xfba78标记附近[15]

本研究构建了一个重组自交系群体, 利用两种分子标记构建了遗传图谱, 结合2年芒长表型数据, 定位了2个QTL, 为进一步精细定位芒长相关基因奠定了基础。

1 材料与方法
1.1 作图群体及田间试验

郑麦9023 (Triticum aestivum, AABBDD)是一个有芒普通小麦品种, 西藏半野生小麦Q1028 (Triticum aestivum, AABBDD)是一个无芒六倍体材料。它们的芒长差异显著(图1)。以郑麦9023为母本与Q1028杂交, F1代起连续自交, 通过单粒传法获得重组自交系(RIL)。

图1 郑麦9023(有芒)和Q1028(无芒)的穗芒形态Fig. 1 Awn feature of Zhengmai 9023 (with awn) and Q1028 (without awn)

将186个RIL株系及其亲本分别于2011— 2012 (F11)和2012— 2013 (F12)小麦生长季, 在四川农业大学小麦研究所温江试验田进行表型观测。试验田肥力中等、地力均匀。采用单粒播种, 每系种1行, 行长2 m, 行距0.3 m, 株距0.1 m。10月12日至15日播种, 分蘖前均苗, 使苗与苗之间实际间距大于0.08 m, 抽穗期单穗套袋自交。其他田间管理与大田生产一致。

1.2 芒长测定

小麦成熟后, 每株随机取5穗, 用游标卡尺测量侧芒的长度, 取平均值。

1.3 分子标记分析

取亲本及RIL系(F12)的幼苗叶片, 利用CTAB法提取小麦基因组DNA。根据Rö der等[16]和Somers等[17]发表的小麦遗传图谱, 选取1089对SSR引物用于亲本间全基因组多态性筛选, 其中107对引物在双亲间具有多样性, 用于检测186个RIL系。引物序列信息来源于http://wheat.pw.usda.gov/GG2/index. shtml, 引物由上海英俊生物有限公司合成。PCR扩增程序为94℃变性5 min; 94℃变性1 min; 60℃退火(根据引物不同而调整) 1 min; 72℃延伸1.5 min共40个循环; 72℃延伸7 min后降温到12℃。扩增产物经8%聚丙烯酰胺凝胶电泳检测, 银染显带。DArT芯片由Diversity Arrays Technology Pty. Ltd. (http://www. diversityarrays.com/)开发, 分析检测方法参照Wenzl等[18]的报道。

1.4 连锁图谱的构建和QTL定位

按Jiang等[19]描述的方法构建遗传连锁图, 采用QTL IciMapping 3.2软件对所调查的重组自交系以及芒长数据进行QTL分析[20]。基于逐步回归法(inclusive composite interval mapping, ICIM)作图[21], 采用1.0 cM 步移速度进行QTL分析, LOD阈值设为2.5。

2 结果与分析
2.1 亲本与重组自交系的芒长

郑麦9023侧芒平均为6.2 cm, 顶芒为4.0 cm; Q1028无芒。RIL群体不同株系芒长差异巨大。两年的表型数据均呈现明显的双峰分布, 短芒株系的芒长为0.4~1.2 cm; 长芒株系的芒长为2.8~5.2 cm; 极端长芒株系较少, 这表明本群体中芒长性状应该由多基因控制并存在一个效应更强的位点(图2)。

图2 重组自交系群体两年芒长直方图Fig. 2 Frequency distribution of awn length in the RIL population for two years

2.2 芒长的QTL分析

共检测到2个芒长QTL, 分别位于4AS和5AL染色体, 两年平均可分别解释表型变异的7.4%和27.3% (表1)。位于4AS的QTL被命名为Qwa.sau-4AS。两年的表型数据均将其定位于Xgpw4448附近(图3-a), 相距1.3 cM。位于5AL染色体的QTL被命名为Qwa.sau-5AL, 其贡献率高于Qwa.sau-4AS (表1), 两年的数据均将其定位于wPt-1038附近, 相距1.1 cM (图3-b)。

图3 4A染色体(a)和5A染色体(b)遗传图谱及芒长性状分析的LOD曲线图
连锁图谱右边是标记名称, 左边是标记位置(cM)。红色和绿色曲线分别为2012和2013年数据; 虚线指示LOD阈值(LOD = 2.5)。
Fig. 3 Genetic maps of chromosomes 4A (a) and 5A (b) and the LOD plots of awn length
In the linkage map, marker names and positions (cM) are on the right and left, respectively. LOD scores in 2012 and 2013 are shown with the red and green curves, respectivley. The vertical dashed line is the threshold of significance (LOD = 2.5).

表1 郑麦9023 × Q1028群体的芒长QTL Table 1 QTLs for awn length identified in the population of ZM9023 × Q1028
2.3 不同基因型芒长效应

利用与2个QTL连锁的分子标记, 对186个RIL株系进行基因型筛选。去除83个在QTL区间发生重组的株系, 共得到26个同时含有2个QTL的株系, 其平均侧芒长为1.21 cm; 得到24个不含芒长QTL的株系, 其平均侧芒长为4.43 cm; 17个仅含Qwa.sau-5AL位点的株系, 其平均侧芒长为1.71 cm; 37个株系仅含Qwa.sau-4AS位点的株系, 其平均侧芒长为2.97 cm (图4)。

图4 RIL群体不同基因型的芒长箱型分布图Fig. 4 Box plot distribution of awn length in different genotypes of the RIL population
L1: Qwa.sau-5AL; L2: Qwa.sau-4AS ; L0: without QTL.

3 讨论

本研究基于郑麦9023和Q1028构建的RIL群体及其全基因组遗传图谱, 共检测到2个与芒长相关的QTL。其中Qwa.sau-4AS被定位于染色体4AS上的Xgpw4448附近, 与Sourdille等[15]的定位结果相比较, Xgpw4448Xfba78之间的遗传距离接近, Qwa.sau-4AS应与之前报道的Hd基因位于同一区域, 因此推测该QTL效应来源于钩芒基因HdQwa.sau-5ALwpt-1038连锁, 与杜斌等[9]的定位结果相比较, wpt-1038Xgwm291都位于5AL的末端区域, 与B1基因在同一区域, 并表现出同样强烈的芒抑制作用, 因此该QTL的效应应该来源于强抑芒基因B1

2个抑芒位点均来源于无芒亲本Q1028, 对芒长具有显著抑制, 并存在累加效应。通过比较不同基因型株系的芒长数据发现, 只携带Hd基因的株系平均芒长为2.97 cm, 而只携带B1基因的株系平均芒长为1.71 cm, 说明B1对芒长的抑制效果强于Hd基因。而当同一株系上同时存在这2个基因时, 平均芒长最小。由此可见, B1基因与Hd基因对芒长的抑制效果具有累加效应。这2个位点累计可解释33.9%~35.5%的表型变异, 推测存在未检测出的QTL。

无芒小麦中国春作为一个普通小麦的模式植物, 对其芒长性状已经有较深入研究。普遍认为其无芒特性是由HdB1B2联合作用所致[15]。本试验所采用的无芒亲本为西藏半野生小麦, 许多系统分类及聚类研究已经证明其与中国春具有密切的亲缘关系[22, 23, 24]。未被检测的QTL效应可能源于6BL染色体的B2, 我们采用连锁图谱的6B染色体包含着与B2基因连锁的标记Xwmc182, 但该位点及其所在区域与芒长无显著连锁, 因此未检测到预期的QTL。其原因可能是所采用的2个亲本之间B2基因型相同。

本研究结果为HdB1基因后续精细定位与克隆提供了分子标记基础。我们已经利用与钩芒基因Hd连锁的分子标记, 从F8的自交系群体中筛选到杂合单株构建近等基因系群体, 正在进行精细定位相关工作。

4 结论

共定位2个与芒长性状相关的QTL, 分别位于4AS和5AL染色体, 共解释34.7%的表型变异。这2个QTL应是钩芒基因Hd和抑芒基因B1。2个基因对芒长的抑制作用具有累加效应, 且B1的抑制效果强于Hd。本研究为后续筛选近等基因系及基因精细定位奠定了基础。

The authors have declared that no competing interests exist.

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

参考文献
[1] Biscoe P V, Littleton E J, Scott R K. Stomatal control of gas exchange in barley awns. Ann Appl Biol, 1973, 75: 285-297 [本文引用:1]
[2] Blum A . Photosynthesis and transpiration in leaves and ears of wheat and barley varieties. J Exp Bot, 1985, 36: 432-440 [本文引用:1]
[3] 黄瑾, 骆惠生, 张勃, 贾秋珍, 金明安, 曹世勤, 金社林. 普通小麦芒的遗传分析. 甘肃农业科技, 2011, (2): 11-12
Huang J, Luo H S, Zhang B, Jia Q Z, Jin M A, Cao S Q, Jin S L. Genetic analysis of mount of common wheat. Gansu Agric Sci Tech, 2011, (2): 11-12 (in Chinese with English abstract) [本文引用:1]
[4] 巴青松, 傅兆麟, 白凡杰. 小麦芒的研究. 淮北煤炭师范学院学报(自然科学版), 2010, 31(1): 29-33
Ba Q S, Fu Z L, Bai F J. The research of wheat awns. J Huaibei Coal Ind Teachers Coll (Nat Sci Edn), 2010, 31(1): 29-33 (in Chinese with English abstract) [本文引用:1]
[5] Maydup M L, Antonietta M, Graciano C, Guiamet J J, Tambussi E A. The contribution of the awns of bread wheat (Triticum aestivum L. ) to grain filling: responses to water deficit and the effects of awns on ear temperature and hydraulic conductance. Field Crops Res, 2014, 167: 102-111 [本文引用:1]
[6] 王忠, 顾蕴洁, 高煜珠. 麦芒的结构及其光合特性. 植物学报, 1993, 35: 921-928
Wang Z, Gu Y J, Gao Y Z. Structure and photosynthetic characteristics of awns of wheat and barley. Acta Bot Sin, 1993, 35: 921-928 (in Chinese with English abstract) [本文引用:1]
[7] Li X J, Wang H G, Li H B, Zhang L Y, Teng N J, Lin Q Q, Wang J, Kuang T Y, Li Z S, Li B, Zhang A M, Lin J X. Awns play a dominant role in carbohydrate production during the grain-filling stages in wheat (Triticum aestivum). Physiol Plant, 2006, 127: 701-709 [本文引用:1]
[8] Rebetzke G J, Bonnett D G, Reynolds M P. Awns reduce grain number to increase grain size and harvestable yield in irrigated and rainfed spring wheat. J Exp Bot, 2016, 67: 2573-2586 [本文引用:1]
[9] 杜斌, 崔法, 王洪刚, 李兴锋. 小麦芒长抑制基因B1近等基因系的鉴定及遗传分析. 分子植物育种, 2010, 8: 259-264
Du B, Cui F, Wang H G, Li X F. Characterization and genetic analysis of near-isogenic lines for awn-inhibitor gene B1 of common wheat. Mol Plant Breed, 2010, 8: 259-264 (in Chinese with English abstract) [本文引用:4]
[10] Elbaum R, Zaltzman L, Burgert I, Fratzl P. The Role of wheat awns in the seed dispersal unit. Science, 2007, 316: 884-886 [本文引用:1]
[11] 姚国新, 张强, 吴建涛, 胡广隆, 李自超. 利用近等基因系对水稻芒基因AWN3-1的遗传定位. 中国农业大学学报, 2010, 15(5): 1-5
Yao G X, Zhang Q, Wu J T, Hu G L, Li Z C. Mapping awn gene AWN3-1 with near-isogenic line of rice. J China Agric Univ, 2010, 15(5): 1-5 (in Chinese with English abstract) [本文引用:1]
[12] Kosuge K, Watanabe N, Kuboyama T, Melnik V M, Yanchenko V I, Rosova M A, Goncharov N P. Cytological and microsatellite mapping of mutant genes for spherical grain and compact spikes in durum wheat. Euphytica, 2008, 159: 289-296 [本文引用:1]
[13] Li H, Han Y C, Guo X X, Xue F, Wang C Y, Ji W Q. Genetic effect of locus B2 inhibiting awning in double-ditelosomic 6B of Triticum durum DR147. Genet Resour Crop Evol, 2015, 62: 407-418 [本文引用:1]
[14] 金善宝. 中国小麦学. 北京: 中国农业出版社, 1996. p 13
Jin S B. Wheat Science in China. Beijing: China Agriculture Press, 1996. p 13 (in Chinese) [本文引用:1]
[15] Sourdille P, Cadalen T, Gay G, Gill B, Bernard M. Molecular and physical mapping of genes affecting awning in wheat. Plant Breed, 2002, 121: 320-324 [本文引用:3]
[16] Röder M S, Korzun V, Wendehake K, Plaschke J, Tixier M H, Leroy P, Ganal M W. A microsatellite map of wheat. Genetics, 1998, 149: 2007-2023 [本文引用:1]
[17] Somers D J, Isaac P, Edwards K. A high-density microsatellite consensus map for bread wheat (Triticum aestivum L. ). Theor Appl Genet, 2004, 109: 1105-1114 [本文引用:]
[18] Wenzl P, Carling J, Kudrna D, Jaccoud D, Huttner E, Kleinhofs A, Kilian A. Diversity arrays technology (DArT) for whole-genome profiling of barley. Proc Natl Acad Sci USA, 2004, 101: 9915-9920 [本文引用:1]
[19] Jiang Y F, Lan X J, Luo W, Kong X C, Qi P F, Wang J R, Wei Y M, Jiang Q T, Liu Y X, Peng Y Y, Chen G Y, Dai S F, Zheng Y L. Genome-wide quantitative trait locus mapping identifies multiple major loci for brittle rachis and threshability in Tibetan semi-wild wheat (Triticum aestivum ssp. tibetanum Shao). PloS One, 2014, 9: e114066 [本文引用:1]
[20] Meng L, Li H, Zhang L, Wang J. QTL IciMapping: Integrated software for genetic linkage map construction and quantitative trait locus mapping in biparental populations. Crop J, 2015, 3: 269-283 [本文引用:1]
[21] Li H H, Ye G Y, Wang J K. A modified algorithm for the Improvement of Composite Interval Mapping. Genetics, 2007, 175: 361-374 [本文引用:1]
[22] Wang Y J, Wang C Y, Zhang H, Yue Z N, Liu X L, Ji W Q. Genetic analysis of wheat (Triticum aestivum L. ) and related species with SSR markers. Genet Resour Crop Evol, 2013, 60: 1105-1117 [本文引用:1]
[23] Wang H Y, Wang X E, Liu D J. Assessment of genetic diversity of Yunnan, Tibetan, and Xinjiang wheat using SSR markers. J Genet Genomics, 2007, 34: 623-633 [本文引用:1]
[24] Kosuge K, Watanabe N, Kuboyama T, Melnik V M, Yanchenko V I, Rosova M A, Goncharov N P. Comparison of the genetic diversity between Triticum aestivum ssp. tibetanum Shao and Tibetan wheat land races (Triticum aestivum L. ) by using intron-splice junction primers. Euphytica, 2008, 159: 289-296 [本文引用:1]