小麦品种(系)抗叶锈病基因Lr10 位点基因型的多样性
韩淑晓1, 刘全兰2, 董洁1, 陈建省1, 田纪春1,*
1山东农业大学作物生物学国家重点实验室 / 山东省作物生物学重点实验室, 山东泰安271018
2青岛科技大学化工学院生物工程与技术系, 山东青岛266042
* 通讯作者(Corresponding author): 田纪春, E-mail:jctian@sdau.edu.cn, Te1: 0538-8242040
摘要

小麦抗叶锈病基因Lr10 的表达需与其紧密连锁的RGA2 基因调控, 按照这2个基因的连锁关系,Lr10 基因位点分成H1和H2两个古单倍型。为揭示我国小麦品种中Lr10 的遗传多样性, 对189个来自12省的小麦育成品种和58个品系的该基因位点变异进行了鉴定分析。绝大多数供试品种(系)为H2古单倍型, 其频率在育成品种和品系中分别为95.2% (180/189)和96.6% (56/58)。这2种古单倍型可进一步分为9种单倍型亚型, 其中H1-2、H2-4、H2-5、H2-6和H2-7亚型为首次报道。育成品种包含所有单倍型亚型, 以H2-1频率最高(69.3%), H2-3和H2-6频率最低(0.5%); 而在选育品系中仅检测到5种单倍型亚型, 其中H2-1频率最高(27.6%), H1-2频率最低(3.5%)。除H2-1与H2-2型外, 其他亚型与育成地显著相关(P< 0.05)。本研究结果支持Lr10 基因位点的古单倍型和单倍型亚型在育种中保持相对稳定的推断, 同时由于在育种过程中持续重组和变异而产生了新的单倍型亚型, 而且重组可以发生在H2和H1两种古单倍型之间。在现代育成品种和品系中, H1古单倍型所占比例已降至5%以下, 应加以保护。

关键词: 小麦; Lr10; 古单倍型; 单倍型亚型; 遗传多样性
Genotypic Diversity atLr10 Locus for Leaf Rust Resistance in Various Hexaploid Wheat Varieties (Lines)
HAN Shu-Xiao1, LIU Quan-Lan2, DONG Jie1, CHEN Jian-Sheng1, TIAN Ji-Chun1,*
1State Key Laboratory of Crop Biology / Shandong Provincial Key Laboratory of Crop Biology, Shandong Agricultural University, Tai’an 271018, China
2 Department of Bioengineering and Biotechnology, College of Chemical Engineering, Qingdao University of Science and Technology, Qingdao 266042, China
Abstract

Genotypes H1 and H2 are two ancient haplotypes onLr10 locus, which is a resistance gene to wheat leaf rust. They are identified based on the presence of bothLr10 andRGA2 in full length (H1), and the absence ofLr10 and chromatin rearrangement ofRGA2 (H2) due to chromosome 1AS reorganization during species evolution. Both haplotypes contain several subhaplotypes. To understand the genetic diversity onLr10 locus, we tested the frequencies of H1 and H2 haplotypes and their subhaplotypes in 189 wheat cultivars and 58 breeding lines from various provinces of China. The H2 haplotype was dominant in the developed cultivars and breeding lines with the frequencies of 95.2% (180/189) and 96.6% (56/58), respectively. Five novel subhaplotypes were detected in the developed cultivars, namely H1-2, H2-4, H2-5, H2-6, and H2-7. In cultivars, subhaplotype H2-1 was the most frequent (69.3%), whereas subhaplotypes H2-3 and H2-6 had the lowest frequencies (0.5%). Among the five subhaplotypes detected in breeding lines, subhaplotype H2-1 had the highest frequencies (27.6%), whereas subhaplotypes H1-2 was the least frequent (3.5%). Interestingly, the subhaplotypes except for H2-1 and H2-2 were highly breeding location-dependent (P< 0.05). These results support the hypothesis that both haplotypes have been maintained through a balanced polymorphism mechanism during wheat domestication. Furthermore, the five novel subhaplotypes suggest that genetic diversities resulting from recombination and other types of chromatin reorganizations (e.g., origin of the H2 haplotype and polymorphism for theLr10 locus) have been going on since the formation of the hexaploid wheat. Significantly, identification of the H1-2 subhaplotype is a clear evidence of occurrence of recombination between the two ancient haplotypes, which disagrees with the previous notion. Importantly, because the frequencies of H1 is lower than 5% in developed cultivars and breeding lines, actions for protecting germplasm with the H1 haplotype are highly recommended.

Keyword: Wheat; Lr10; Ancient haplotype; Subhaplotype; Genetic diversity
引言

小麦是全球主要粮食作物之一, 种植面积仅次于玉米。叶锈病是小麦的主要病害之一, 发病严重时可减产50%[1]。小麦叶锈病的病原菌为 Puccinia triticina , 具有小种专化性, 对不同小麦品种具有不同的致病性[2]; 但病原小种易发生变异, 致使抗病品种丧失抗性, 特别是大面积推广品种都高抗当时流行的生理小种, 加速了病原菌的定向选择和新生理小种的形成。因此, 研究抗叶锈基因在小麦品种(系)中的多样性和分布规律, 对合理利用抗叶锈病基因资源, 延缓抗病品种的抗性丧失具有重要意义。

抗叶锈基因 Lr10 位于小麦1A染色体短臂上[3,4],与另一抗病基因同源序列 RGA2 紧密连锁, 在该位点存在H1和H2两种古单倍型[3]。H1古单倍型含有完整的 Lr10 基因及与其紧密连锁的 RGA2 基因, H2型仅含有 RGA2 的部分序列, 缺失 RGA2 的5′端部分序列及全部 Lr10 序列[3,4,5,6]。群体遗传学研究明, 这两种古单倍型原存在于一粒小麦( Triticum urartu )和野生二粒小麦( T. dicoccoides )中[3,4,5,6]。在小麦进化和选育过程中, Lr10 基因在普通小麦中被保留下来, 并对叶锈菌菌株( AvrLr10 ) 89-201 CBTB(TX)现抗性[5]。Stein等[5]检测了 Lr10 古单倍型在56份欧洲小麦品系中的分布规律, 发现H1古单倍型的比例为14.29% (8/56), H2古单倍型的比例为85.71% (48/56), 说明H2古单倍型是欧洲普通小麦的主要类型。Isidore等[6]设计了4对H2单倍型的特异引物(A、B、C和B_3kb)用于研究不同倍性小麦中H2单倍型的保守性, 共检测出3种H2亚型, 在欧洲普通小麦中仅检测出2种H2亚型。尽管前人对 Lr10 基因位点的两种古单倍型及其亚型在小麦中的分布规律已有较多研究, 但 Lr10 基因位点在中国小麦品种中的分布规律未见报道。

小麦育种过程中多选用少数综合性状好、配合力好的品种为亲本, 并以高产、优质、抗病虫和抗旱等为主要育种目标对后代进行更强有力的选择。这种强目的性的人工选择必将强化基因的瓶颈效应, 导致育成品种遗传基础变窄, 遗传多样性降低。与古老的地方品种和野生品种相比, 现代育成品种遗传多样性明显降低。为评价 Lr10 基因位点在中国小麦品种中的遗传多样性及可能经历的遗传瓶颈效应, 本研究选用247份中国小麦品种(系), 检测 Lr10 基因位点的古单倍型及其等位基因的遗传变异, 旨在阐明 Lr10 基因位点的H1和H2古单倍型在小麦驯化过程中的遗传瓶颈效应, Lr10 基因位点在中国小麦品种(系)中的等位变异及其遗传多样性, 及 Lr10 基因位点的多样性机制。

1 材料与方法
1.1 植物材料

189个小麦育成品种和58个高代品系(见), 均由国家小麦改良中心泰安分中心提供。189个育成品种包括山东103份、河北31份、河南22份、陕西10份、江苏7份、山西5份、北京4份、安徽2份、宁夏2份, 以及西藏、新疆和云南各1份; 58个高代品系均由国家小麦改良中心泰安分中心近年新育成。2011—2012年, 将所有品种(系)种植于山东农业大学试验站, 每个材料种植3行, 行长3.00 m, 行宽0.26 m, 四周各有1 m保护行。生育期间肥水按高产田管理, 未进行药剂防治病害, 小麦亦未发生倒伏。

1.2 Lr10 基因位点相关片段的PCR扩增

采用Doyle和Doyle[7]的方法从小麦幼苗提取的基因组DNA。采用引物Control[6]和rga1Pro[8], 通过PCR扩增检测H1和H2单倍型, 其中引物Control (1)用于扩增H1和H2两种单倍型的共同序列(1), 引物rga1Pro (1)特异性扩增 Lr10 基因的5′端非编码区及部分5′端编码区, 只在H1型个体中有扩增产物(1)。引物ThLr10_T/P、ThLr10_E/H和ThLr10_G/J[9]用于检测H1的亚型; 而引物A、B、C和B_3k[6]用于检测H2的亚型。所有引物均由生工生物工程(上海)股份有限公司合成。

在T-gradient Thermal Cycler (Bio-metra)上, 按已报道的方法[6,8,9]进行PCR扩增。扩增产物经琼脂糖电泳后, 1%~2%溴化乙锭染色, 用Bioshine GelX 1650凝胶成像系统(上海欧翔科学仪器有限公司)观察、照相、读带。

1.3 数据分析

Lr10 位点基因型按照PCR扩增片段的有或无命名。例如, 若两对引物Control和rga1Pro都有PCR产物, 则其古单倍型的类型为H1[8]; 若仅引物Control有PCR产物, 引物rga1Pro无扩增产物, 则其古单倍型的类型为H2[8]。H1和H2单倍型亚型的扩增带型如2所示。采用SAS 9.0软件[10]统计各古单倍型及其亚型的发生频率以及基因型。选用样品数较多的山东、河北、河南和陕西材料进行 Lr10 位点基因型与品种育成地的相关性分析。

1 引物扩增片段在 Lr10 位点的位置序列结构引自Isidore 等[6]2和5, 有修改。Fig . 1 Positions of amplified fragments used different primers at locus Lr10The sequence structures were modified from the original Figs. 2 and 5 in the article by Isidore et al.[6]

1 用于PCR扩增的引物序列 Tab le 1 Primers used in this study
2 按PCR产物命名的基因型 Table 2 Genotypes classified by different amplicons produced by different primers
2 结果与分析
2.1 中国小麦品种(系)中 Lr10 基因型的多样性

小麦品种或品系中均检测到H1和H2两种古单倍型(2), 但两种古单倍型的发生频率有显著差异(3)。在189份小麦品种中, 有180份材料的古单倍

型为H2, 占供试品种的95.2%, 仅有9份材料含有H1单倍型, 占供试品种的4.8%; 在58份高代品系中, H2和H1单倍型的频率分别为96.6% (56/58)和3.5% (2/58)。

2 Lr10 基因位点的引物在小麦品种(系)中的扩增谱M: 标准分子量DNA (Trans2K plus); 1: 濮兴02-16; 2: 红地963; 3: 石09-4276; 4: 高-W16; 5: 航麦0903; 6: LS6045; 7: 荷麦0302; 8: 西农9814; 9: 明麦1号; 10: 泰农984; 11: 22037; 12: 89100S。Fig. 2 Amplification profiles on locus Lr10 in wheat varieties (lines) and breeding linesM: tz DNA marker (Trans2K plus); 1: Puxing 02-16; 2: Hongdi 963;3: Shi 09-4276;4: Gao-W16;5: Hangmai 0903; 6: LS6045;7: Hemai 0302;8: Xinong 9814;9: Mingmai 1; 10: Tainong 984; 11: 22037; 12: 89100S.

3 小麦品种(系)中 Lr10 位点各单倍型及其亚型发生频率 Table 3 Frequencies of genotypes at locus Lr10 from common wheat varieties (lines) and breeding lines
2.2 Lr10 基因型遗传多样性与品种育成地的相关性

在来自山东、河北、河南和陕西的育成品种 中, H2-1型( P =0.2701)和H2-2型( P =0.51715)与育成地没有相关性, H2-3 ( P =0.04701)、H2-4 ( P = 0.0464)、H2-5 ( P =0.0397)和H2-7型( P =0.0100)与育成地有显著相关。H2-1和H2-2型广泛分布于这4个省, H2-3和H2-4型仅分布在山东, H2-5型分布在山东、河南和陕西, H2-7型仅分布在河北(3); 而H2-6型在这4个省份中都没有分布, 仅在江苏的小麦品种中检测到(见)。总之, 除基因型H2-1和H2-2外, 其他基因型与其育成地存在显著相关性( P< 0.05)。

3 4省小麦品种中基因型的分布品种(系)序号同中。Series numbers of wheat varieties (lines) correspond to the numbers for variety names given in the Appendix.Fig. 3 Distribution of genotypes across common wheat varieties (lines) from four provinces

3 讨论
3.1 Lr10 基因位点的遗传瓶颈效应

长期的自然选择和人工驯化最终将导致小麦基因的遗传多样性降低[11]。尽管如此, 小麦基因组通过其可塑性可以弥补这种瓶颈效应。Dubcovsky和Dvorak[12]指出, 小麦可以通过其A、B和D基因组的基因补偿效应来弥补基因瓶颈效应造成的不良后果, 杂交育种可以使新育成的小麦获得优于亲本的基因型, 杂交后代通过基因突变、重组、缺失等方式也可以形成全新的基因型。Stein等[5]和Isidore等[6]的研究也明, 小麦的抗病基因同样经历了类似的遗传瓶颈效应, 但这些抗病基因经历瓶颈效应的程度和可能获得的新基因并没有得到深入研究。

Lr10 基因是小麦中4个已测序抗叶锈病基因之一[4](http://www.ncbi.nlm.nih.gov/genbank/), 全长为3935 bp, 与其他抗叶锈病基因序列的相似度较低[4]。该基因编码919个氨基酸长度的CC-NBS-LRR类型蛋白质[9], 在蛋白CC结构域的N端经历多样性选择(diversifying selection)[13]。另一抗病基因同源序列 RGA2 Lr10 基因紧密连锁, 是 Lr10 基因现出抗性必需的调控基因[13]。陈万权和王剑雄[14]采用基因推导法在来自世界各地的76个小麦种质资源中鉴定出4个材料含 Lr10 基因; 丁艳红等[15]采用STS分子标记在28个小麦微核心种质中鉴定出2个材料含 Lr10 基因。本研究在189份小麦品种中仅鉴定出9个同时含有 Lr10 基因和 RGA2 基因的品种, 其中2个品种含有 Lr10 基因和缺失的 RGA2 基因(3); 在58份小麦品系中没有鉴定出同时含有完整 Lr10 基因和 RGA2 基因的材料, 仅有2个品系同时含有 Lr10 基因和缺失的 RGA2 基因(3)。上述研究结果明, Lr10 基因经历了强烈的瓶颈效应, 推测导致这种结果的原因有2个: 其一, 当前流行的叶锈病菌已经对 Lr10 基因产生抗性, 因而, 小麦抗病育种工作者已不选育含有此基因的材料[2]; 其二, Lr10 基因可能对某些稀有致病小种或未知种有抗性, 但对当前流行的叶锈病病原菌没有效果, 且该基因的达对产量或其他农艺性状有负面效应, 最终导致育种工作者不选育含有此基因的材料[4]。我们认为这两种解释都有其合理性, 同时也认为国内外学者采用研究方法的不同也可能导致上述解释的差异。

Isidore等[6]发现, 欧洲小麦中 Lr10 位点的H2-1型占优势。本研究检测结果明, 该亚型在中国小麦中也有最高的发生频率, 其次是是H2-2 (3), 同时也说明, H2-1基因型可能是世界性的优势基因型, 但导致H2-1基因型占优势的原因还需进一步探究。另外, Lr10 位点的瓶颈效应说明需要保护那些基因型发生频率低的材料, 如含H1-1基因型的材料, 这类材料含有较完整的 Lr10 基因和 RGA2 基因。

3.2 Lr10 基因位点的进化

致病菌的高度变异为植物抗病基因的进化提供了很强的自然选择压, 关于两者的相关性主要有2种假说, 即trench-warfare模型[16]和arm-race模型[17]。前者认为抗病基因进化形成的多样性依赖于致病菌的变异频率, 频率高时抗病基因的多样性高, 频率低时抗病基因的多样性低, 在此模型下, 抗病基因的多样性较高, 可以形成新基因和保留已有基因, 如拟南芥中的基因 Rpm1[16] RPP8[17] RPP13[17]; 后者是birth-death模型[18]的补充完善, 强调抗病基因的多样性有明显的方向性, 即选择更具抗性的基因, 在此模型下, 抗病基因的多样性很低, 有明显的方向性, 如番茄中对尖孢镰孢菌( Fusarium oxysporum )产生抗性的抗病基因[19] Lr10 基因的进化模式为trench-warfare模型, 在野生群体中具有较高的遗传多样性[4,20]。Sela等[20]推测重组使 Lr10 基因具有较高的多样性; Isidore等[6]推测逆转录子的插入、转座子的缺失以及 Lr10 基因与 RGA2 基因间元件的不等重组等事件形成H2古单倍型; Isidore等[6]还推测H2和H1古单倍型之间非常稳定, 不会发生重组。本研究的结果支持Sela等[20]的结论, 即H2和H1古单倍型确实存在于小麦品种中, 而且H2-1和H2-2基因型分布广泛并相对稳定, 同时又修正了Isidore等[6]关于重组不能发生在H2和H1古单倍型之间的推断。

本研究鉴定出了5个新的单倍型亚型(H1-2、H2-4、H2-5、H2-6和H2-7, 见), 证明了重组机制在小麦育种过程中持续存在, 杂交育种可以创造出新的 Lr10 基因型。本研究结果证明, 在小麦育种过程中确实有些材料丢失了 Lr10 基因、 RGA2 基因及两基因间的元件, 但小麦通过重组机制可以产生新基因型, 以获得基因的可塑性来创造新的抗病性状。

4 结论

247份中国小麦育成品种和高代品系中H2型为优势基因型, 其中H2-1型出现频率最高。在育成品种中发现5个新单倍型亚型, 分别是H1-2、H2-4、H2-5、H2-6和H2-7, 其频率分别为1.1%、1.1%、2.1%、0.5%和1.6%。在现代小麦品种中, H1古单倍型所占比例低于5%, 亟需进行遗传资源保护。

材料信息及其 Lr10 位点基因型 Appendix Information of materials and their genotypes at locus Lr10
() 材料信息及其 Lr10 位点基因型 Appendix Information of materials and their genotypes at locus Lr10
() 材料信息及其 Lr10 位点基因型 Appendix Information of materials and their genotypes at locus Lr10
The authors have declared that no competing interests exist.

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

参考文献
[1] Komer J A. Genetics of resistance to wheat leaf rust. Annu Rev Phytopathol, 1996, 34: 435-455 [本文引用:1] [JCR: 9.875]
[2] Hu Y-Y(胡亚亚), Zhang N(张娜), Li L-M(李林懋), Yang W-X(杨文香), Liu D-Q(刘大群). Analysis of wheat leaf rust resistance genes in 14 wheat cultivars or lines. Acta Agron Sin(作物学报), 2011, 37(12): 2158-2166 (in Chinese with English abstract) [本文引用:2] [CJCR: 1.8267]
[3] Feuillet C, Gabriele S, Keller B. Molecular cloning of a new receptor-like kinase gene encoded at the Lr10 disease resistance locus of wheat. Plant J, 1997, 11: 45-52 [本文引用:4] [JCR: 6.16]
[4] Sela H, SpiridonL N, Petrescu A J, Akerman M, Mand el-Gutfreund Y, Nevo E, Loutre C, Keller B, Schulman A H, Fahima T. Ancient diversity of splicing motifs and protein surfaces in the wild emmer wheat (Triticum dicoccoides) Lr10 coiled coil (CC) and leucine-rich repeat (LRR) domains. Mol Plant Pathol, 2012, 13: 276-287 [本文引用:7] [JCR: 3.899]
[5] Stein N, Feuillet C, Wicker T, Schlagenhauf E, Keller B. Subgenome chromosome walking in wheat: a 450-kb physical contig in Triticum monococcum L. spans the Lr10 resistance locus in hexaploids wheat (Triticum aestivum L. ). Proc Natl Acad Sci USA, 2000, 97: 13436-13441 [本文引用:5] [JCR: 9.737]
[6] Isidore E, Scherrer B, Chalhoub B, Feuillet C, Keller B. Ancient haplotypes resulting from extensive molecular rearrangements in the wheat A genome have been maintained in species of three different ploidy levels. Genome Res, 2005, 15: 526-536 [本文引用:11] [JCR: 13.608]
[7] Doyle J J, Doyle J L. Isolation of plant DNA from fresh tissue. Focus, 1990, 12: 13-15 [本文引用:1]
[8] Scherrer B, Keller B, Feuillet C. Two haplotypes of resistance gene analogs have been conserved during evolution at the leaf rust resistance locus Lr10 in wild and cultivated wheat. Funct Integrat Genomic, 2002, 2: 40-50 [本文引用:4] [JCR: 2.842]
[9] Feuillet C, Travella S, Stein N, Albar L, Nublat A, Keller B. Map-based isolation of the leaf rust disease resistance gene Lr10 from the hexaploid wheat (Triticum aestivum L. ) genome. Proc Natl Acad Sci USA, 2003, 100: 15253-15258 [本文引用:3] [JCR: 9.737]
[10] Cody R. Learning SAS by Example: A Programmer’s Guide, SAS Institute Inc. , Cary, NC, USA, 2007 [本文引用:1]
[11] Nalam V J, Vales M I, Watson C J W, Kianian S F, Riera-Lizarazu O. Map-based analysis of genes affecting the brittle rachis character in tetraploid wheat (Triticum turgidum L. ). Theor Appl Genet, 2006, 112: 373-381 [本文引用:1] [JCR: 3.297]
[12] Dubcovsky J, Dvorak J. Genome plasticity a key factor in the success of polyploid wheat under domestication. Science, 2007, 316: 1862-1866 [本文引用:1]
[13] Loutre C, Wicker T, Travella S, Galli P, Scofield S, Fahima T, Feuillt C, Keller B. Two different CC-NBS-LRR genes are required for Lr10-mediated leaf rust resistance in tetrsploid and hexaploid wheat. Plant J, 2009, 60: 1043-1054 [本文引用:2] [JCR: 6.16]
[14] Chen W-Q(陈万权), Wang J-X(王剑雄). Genes for leaf and stem rust resistance in 76 wheat genetic resources. Acta Agron Sin(作物学报), 1997, 23(6): 655-663 (in Chinese with English abstract) [本文引用:1] [CJCR: 1.8267]
[15] Stahl E A, Dwyer G, Mauricio R, Kreitman M, Bergelson J. Dynamics of disease resistance polymorphism at the Rpm1 locus of Arabidopsis. Nature, 1999, 400: 667-671 [本文引用:1] [JCR: 36.28]
[16] Holub E B. The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet, 2001, 2: 516-527 [本文引用:2] [JCR: 38.075]
[17] Michelmore R W, Meyers B C. Clusters of resistance genes in plants evolve by divergent selection and a birth and death process. Genome Res, 1998, 8: 1113-1130 [本文引用:3] [JCR: 13.608]
[18] Takken F, Rep M. The arms race between tomato and Fusarium oxysporum. Mol Plant Pathol, 2010, 11: 309-314 [本文引用:1] [JCR: 3.899]
[19] Sela H, Loutre C, Keller B, Schulman A, Nevo E, Korol A, Fahima T. Rapid linkage disequilibrium decay in the Lr10 gene in wild emmer wheat (Triticum dicoccoides) populations. Theor Appl Genet, 2011, 122: 175-187 [本文引用:1] [JCR: 3.297]