摘要
褐飞虱是水稻生产危害中最严重的虫害之一,发掘水稻中抗褐飞虱的新基因并培育抗性品种是现今控制该虫害最有效的举措。广西农家种质——籼稻47-1连续3年的苗期、成株期对褐飞虱都表现一致高抗,遗传稳定性良好,具有明显的趋避性和抗生性,定位其抗性基因并探究其机理,对丰富抗褐飞虱基因的多样性和培育抗性水稻品种具有重大意义。本研究通过利用抗褐飞虱材料47-1与感褐飞虱材料9311杂交构建了包含91个单株的F2群体,采用虫体增重法进行褐飞虱抗性表型鉴定。结果显示F2群体抗性单株(虫体增重倍数≤1.1)与感性单株(虫体增重倍数>1.1)的数量比为61∶30,卡方检验(
虫害对水稻的生长发育和产量的影响极大,而褐飞虱是水稻生产危害最严重的虫害之一。褐飞虱是专一性寄主水稻的一种同翅目飞虱科(Homoptera: Delphacidae)昆
近年来,随着分子标记的发展和水稻高密度连锁图谱的构建,分别在水稻第2、3、4、6、10、11和12号染色体定位了42个水稻抗褐飞虱主效基因,并克隆了10多个水稻抗褐飞虱基
水稻褐飞虱抗性鉴定方法众多,通用的有苗期集团法、强制取食法、成株期抗虫法等。其中邓钊
本研究所用的抗褐飞虱水稻材料—籼稻品系47-1由广西农业科学院水稻所提供,是一份从广西省柳州市融水苗族自治县收集的农家品种,其在褐飞虱蔓延的田间依然生长健康、分蘖旺盛、结实正常,对褐飞虱表现出高抗特性。本研究于连续3年(2020、2021和2022年)的6月份在广西大学试验基地对材料进行了抗褐飞虱鉴定。在同一基地于2021年5月份以47-1为抗虫父本,感性材料9311为母本进行杂交获得F1,2021年10月份F1自花授粉构建包含91个单株的F2分离群体,其分离群体抗褐飞虱差异表现明显。2022年5月份将F2群体种植于广西大学试验基地。
参考Qiu
苗龄为12~14d的两个亲本植株移栽到同一个塑料桶(d=29 cm, h=20 cm),每桶4株苗(抗、感各2株,同一亲本对角线移栽),同时使用纱网罩将塑料桶罩住,防止褐飞虱入侵。植株长至30 d时,将100头褐飞虱投入每个纱网罩中,并在接虫后1 d、2 d、3 d、4 d、5 d、6 d、7 d、8 d统计植株上着虫数,试验设置5个重复。数据处理同上。
苗龄为12~14 d的两个亲本植株移栽到同一个塑料杯中(d=10 cm, h=16 cm),将小一号的塑料杯(d=8 cm, h=10 cm)底部打上圆孔,将水稻苗穿过圆孔罩住。植株长至30 d时,从圆孔处接入10头2~3龄的褐飞虱若虫,并用脱脂棉堵住圆孔防止褐飞虱飞出。连续8 d统计水稻苗上存活的褐飞虱数目。每个亲本设置8个重复。数据处理同上。
采集新鲜的水稻叶片200 mg磨成细粉,参照Murray
引物名称 Primer name | 正向引物(5'—3') Forward primers (5'—3') | 反向引物(5'—3') Reverse primers (5'—3') | 起始物理位置(bp) Initial location | 结束物理位置(bp) End location |
---|---|---|---|---|
04M0.453 | CATTTCTCCTCTGTCTTCCA | AATCTCCATTAAGCTCCTCC | 441997 | 442098 |
04M3.688 | CATTGGACAGTCACATGTTT | TGGTCCAATGACTAAACACA | 4643332 | 4643475 |
04M6.387 | GTAGAGGAGGGTGACGGC | CATCCTCCTCCCACCATC | 6156077 | 6156205 |
04M04514 | TGATTTAACCGTAGATTGGA | TCCAACATATTCCCAAAAAC | 8215220 | 8215085 |
04M10.677 | GTCATGGAGTAGGTCGAACA | TCATTGGGTAAGCAAACTCT | 10701910 | 10702060 |
Bph40-CDS | ATGCTTTGGGAGAGCTTTTATG | TCACCAGAGAAGAGAGCATCTG | 4490101 | 4484997 |
物理位置参考日本晴基因组
Physical location based on Nipponbare genome sequence
基于F2单株表型值构建的抗、感池间已获得有多态性的引物,在其位置附近筛选更多亲本间存在多态性的引物,使用获得的多态性引物对作图群体中91个F2单株的基因型进行检测。利用软件Join Map 3.0构建F2群体的局部遗传连锁图谱并获得分子标记间的遗传距离。基于抗虫鉴定获得的F2群体表型数据以及基因型数据,使用软件ICIM对F2群体中的抗性QTL位点进行检测。
温室采集新鲜亲本叶片,剪碎置于锡纸中,短期保存于液氮,长期保存-80 ℃冰箱。水稻RNA提取采用TRIzol法。TRIzol购买自全式金公司(TransZol Up货号:ET111)。提取的RNA置于冰上,用1.5%的琼脂糖凝胶对RNA进行质量鉴定。观测到明亮28 s带后使用超微量分光光度计进行RNA浓度及质量鉴定,RNA鉴定合格后使用诺唯赞试剂盒HiScript III RT SuperMix for qPCR(+gDNA wiper)(货号:R323)反转录为cDNA。
参考已知Bph40的CDS序列设计扩增引物(
通过连续3年对籼稻47-1进行褐飞虱抗性的鉴定,发现其遗传稳定性较好。3年苗期接虫结果表明籼稻47-1对广西南宁田间的褐飞虱自然群体具有高抗性,植株表现健壮并呈正常绿色,2020年抗性值为2.17~2.76,均值为2.49;2021年抗性值为2.37~3.17,均值为2.80;2022年抗性值为2.43~3.03,均值为2.80(
年份 Year | 时期 Period | 编号 Number | 平均值±标准差 Mean±SD | 极大值 Max. | 极小值 Min. | 标准误 SE |
---|---|---|---|---|---|---|
2020 | 苗期 | 47-1 |
2.49±0.3 | 2.76 | 2.17 | 0.17 |
9311 | 8.16±0.17 | 8.33 | 8.00 | 0.10 | ||
成株期 | 47-1 |
2.25±0.2 | 2.50 | 2.00 | 0.14 | |
9311 | 8.50±0.25 | 8.75 | 8.25 | 0.14 | ||
2021 | 苗期 | 47-1 |
2.80±0.4 | 3.17 | 2.37 | 0.23 |
9311 | 8.12±0.20 | 8.30 | 7.90 | 0.12 | ||
成株期 | 47-1 |
2.33±0.3 | 2.75 | 2.00 | 0.22 | |
9311 | 8.50±0.43 | 8.75 | 8.00 | 0.25 | ||
2022 | 苗期 | 47-1 |
2.80±0.3 | 3.03 | 2.43 | 0.19 |
9311 | 8.18±0.22 | 8.43 | 8.03 | 0.13 | ||
成株期 | 47-1 |
1.92±0.2 | 2.25 | 1.75 | 0.17 | |
9311 | 8.33±0.14 | 8.50 | 8.25 | 0.08 |
**表示差异极显著(P<0.01);下同
** indicates extremely significant difference (P<0.01);The same as below

图1 抗感材料苗期和成株期抗虫表型
Fig.1 BPH resistant phenotype of resistant and susceptible materials at seedling and adult stage
A:苗期接虫前;B:苗期接虫后; C: 47-1成株期接虫前; D: 47-1成株期接虫后; E: 9311成株期接虫前; F: 9311成株期接虫后
A: Before BPH infestation in seeding stage; B: After BPH infestation in seeding stage; C: Before BPH infestation in 47-1 adult stage; D: After BPH infestation in 47-1 adult stage; E: Before BPH infestation in 9311 adult stage; F: After BPH infestation in 9311 adult stage
为了进一步探讨抗性材料47-1对褐飞虱的趋避性和抗生性特征,本研究通过对褐飞虱虫体的宿主选择和存活试验来进行衡量。宿主选择试验结果表明,接虫后8 d内,抗性材料47-1植株上的褐飞虱数量均低于感性材料9311,在第3天后一直达到显著水平,并且相对着虫率呈下降趋势,感性材料9311反之;第8天时47-1植株上褐飞虱平均相对着虫率仅为16%,9311植株上为84%,表明抗性材料47-1对褐飞虱具有显著的趋避性(

图2 籼稻47-1对褐飞虱的生理抗性特征
Fig.2 Physiological resistance characterization of indica rice 47-1 resistance to BPH
A:宿主选择试验;B:虫体存活试验
A: BPH host choice test; B: Number of BPH survival
本研究采用褐飞虱取食48 h后的虫体增重倍数作为表型值,对亲本和91个F2单株进行褐飞虱增重试验,结果表明,亲本47-1和9311对褐飞虱群体分别表现出高抗和高感表型,平均虫体增重倍数为0.83和1.32。F2群体单株中抗虫性变异较广,虫体增重倍数最小值为0.67,最大值为1.75,平均值为1.07。褐飞虱虫体增重鉴定结果表明,定位群体中91个F2单株的抗感分离呈连续性分布(

图3 47-1/9311 F2群体的虫体增重倍数分布
Fig.3 Distribution of BPH weight increased multiple of 47-1/9311 F2 population
F2表型 F2 phenotype | F2个体数 Number of F2 individuals | F2家系虫体增重倍数 BPH weight increased multiple of F2 line |
---|---|---|
高抗 HR | 13 | ≤0.9 |
抗 R | 48 | 0.9~1.1 |
高感 HS | 30 | >1.1 |
HR: High resistance; R: resistance; HS: High susceptibility
使用亲本47-1和9311筛选水稻12条染色体上分布的1599对分子标记,结果两亲本间共获得236个具有多态性的标记。基于褐飞虱虫体增重法从F2群体中选取10个高抗单株和10个高感单株构建抗感DNA池,标记筛选发现在第4号染色体上发现3个多态性标记。在3个标记附近增加2个多态性标记,获取91个F2单株在这5个多态性标记上的基因型,构建局部染色体连锁图谱。结合F2群体的表型和基因型数据,利用软件ICIM分析抗性位点所在染色体区段的具体位置。结果表明,在第4号染色体的分子标记04M0.453与04M3.688间检测到1个明显的QTL位点,最大LOD值为5.9,在群体中对抗褐飞虱表型贡献率达到25.9%(

图4 构建的局部遗传连锁图及抗性QTL位置
Fig.4 Construction of local genetic linkage map and locationof resistance QTL
虚线LOD值=2.5。PEV表示基因位点解释了F2群体中抗褐飞虱表型变异的25.9%,红色三角形表示抗性位点
The dotted line indicates LOD score =2.5. PEV indicated that the gene locus explained 25.9% of the phenotypic variation of BPH resistance in F2 population. The red triangle represents the resistance gene locus
根据已知的日本晴和9311基因组的信息,其分子标记04M0.453与04M3.688间有1个已经克隆的抗褐飞虱基因Bph4

图5 Bph40(t)和Bph40氨基酸序列分析
Fig.5 Amino acid sequence analyse of Bph40(t) and Bph40
A:Bph40(t)和Bph40氨基酸序列比较;数字代表氨基酸位置;LRR为亮氨酸富集的重复序列缩写;P为脯氨酸,T为苏氨酸,S为丝氨酸,L为亮氨酸;B:基于Bph40(t)氨基酸序列及其同源序列构建的系统发育树;方框内为Bph40(t)及已克隆Bph40的氨基酸序列
A: Amino acid sequence comparison of Bph40(t) and Bph40; Numbers represent amino acid positions; LRR is the abbreviation of leucine rich repeat; P is proline; T is threonine; S is serine; L is leucine; B: Phylogenetic tree developed with amino acid sequence of Bph40(t) and its homologous sequences; Box means amino acid sequence of Bph40(t) and cloned Bph40
随着农药防治的大规模推广,褐飞虱“生物型”也出现更多的变异,通过大量人力物力培育的抗性品种逐渐失去抗性,例如携带Bph1的IR26、携带Bph2的IR36和IR4
目前水稻褐飞虱抗性鉴定常用的方法为苗期集团法,存在制备F3群体的时间长,耗费种子量大、需要褐飞虱虫量大等缺点。邓钊
Bph40与Bph30均定位在第4号染色体上,是同一类编码具有2个LRR结构域的抗性基因。目前水稻第4号染色体上定位克隆的抗褐飞虱基因最多。例如在第4号染色体上的2个区段内定位了7个具有褐飞虱抗性的主效基因:其中一个区段在分子标记C946和RM5953之间,物理距离约4.9 Mb,覆盖了Bph12(t)、Bph15、Bph17和Bph20(t) 4个基因;另一个区段在分子标记G271和RM273之间,物理距离约4.1 Mb,覆盖了Bph12、Bph18(t)和Bph
本研究通过褐飞虱虫体增重试验,利用籼稻47-1和9311构建的F2分离群体,采用褐飞虱取食48 h后的虫体增重倍数作为表型值发掘了1个抗褐飞虱的主效基因,通过氨基酸序列分析表明该基因是Bph40的1个新的等位基因。本研究通过褐飞虱虫体增重法鉴定水稻F2遗传群体的表型值,具有时间短、见效快、操作简单等特点,为快速发掘抗褐飞虱基因提供一定参考基础,同时定位的新主效基因Bph40(t),为多基因聚合育种选育新品种提供了抗源多样性,一定程度促进了防治褐飞虱危害工作的进步。
参考文献
梁云涛,王春连,刘丕庆,傅强,Kshirod K J,赵开军. 水稻抗褐飞虱基因研究和利用现状. 植物遗传资源学报, 2008,9(1): 119-124,129 [百度学术]
Liang Y T,Wang C L,Liu P Q,Fu Q,Kshirod K J,Zhao K J. Research and application status of rice genes resistance to brown planthopper. Journal of Plant Genetic Resources, 2008, 9(1): 119-124,129 [百度学术]
黄得润, 陈洁, 赖凤香, 刘光杰, 庄杰云. 东乡野生稻抗褐飞虱QTL分析. 作物学报, 2012, 38(2): 210-214 [百度学术]
Huang D R,Chen J,Lai F X,Liu G J,Zhuang J Y. QTL analysis of resistance to brown planthopper in Dongxiang wild rice. Acta Agronomica Sinica, 2012, 38(2): 210-214 [百度学术]
Backus E A,Serrano M S,Ranger C M. Mechanisms of hopperburn: An overview of insect taxonomy, behavior, and physiology. Annual Review of Entomology,2005, 50: 125-151 [百度学术]
Li Z H,Xue Y X,Zhou H L,Li Y,Usman B,Jiao X Z,Wang X Y,Liu F,Qin B X, Li R B, Qiu Y F. High resolution m apping and breeding application of a novel brown planthopper resistance gene derived from wild rice (Oryza rufipogon Griff). Rice, 2019, 12(1): 41 [百度学术]
Jairin J,Phengrat K,Teangdeerith S,Vanavichit A,Toojinda T. Mapping of a broad-spectrum brown planthopper resistance gene, Bph3, on rice chromosome 6. Molecular Breeding, 2007, 19(1): 35-44 [百度学术]
邓一文,刘裕强,王静,陈学伟,何祖华. 农作物抗病虫研究的战略思考. 中国科学:生命科学, 2021,51(10):1435-1446 [百度学术]
Deng Y W,Liu Y Q,Wang J,Chen X W,He Z H. Strategic thinking and research on crop disease and pest resistance in China. Scientia Sinica(Vitae), 2021,51(10):1435-1446 [百度学术]
邓钊,石少阶,王卉颖,上官欣欣,刘丙芳,荆胜利,杜波,陈荣智,祝莉莉,何光存. 籼稻资源WD15515中抗褐飞虱QTL的定位研究. 作物学报, 2016,42(3):353-360 [百度学术]
Deng Z,Shi S J,Wang H Y,Shang-Guan X X,Liu B F,Jing S L,Du B,Chen R Z,Zhu L L, He G C. Analysis of QTLs for brown planthopper in indica rice WD15515. Acta Agronomica Sinica, 2016,42(3):353-360 [百度学术]
Qiu Y F,Guo J P,Jing S L,Zhu L L,He G C. High-resolution mapping of the brown planthopper resistance gene BPH6 in rice and characterizing its resistance in the 9311 and Nipponbare near isogenic backgrounds. Theoretical and Applied Genetics, 2010, 121: 1601-1611 [百度学术]
Lin J B,Wang X Y,Li Y,Bi F G,Cheng L,Huang F K,Li R B,Qiu Y F. Fine mapping, candidate genes analysis, and characterization of a brown planthopper (Nilaparvata lugens Stål) resistance gene in the rice variety ARC5984. Euphytica, 2020, 216:13 [百度学术]
Murray M G, Thompson W F. Rapid isolation of high molecular-weight plant DNA. Nucleic Acids Research, 1980, 8 (19):4321-4325 [百度学术]
Shi S J,Wang H Y,Nie L Y,Tan D,Zhou C,Zhang Q,Li Y,Du B,Guo J P,Huang J,Wu D,Zheng X H,Guan W,Shan J S,Zhu L L,Chen R Z,Xue L J,Walling L L,He G C. Bph30 confers resistance to brown planthopper by fortifying sclerenchyma in rice leaf sheaths. Molecular Plant, 2021, 14: 1714-1732 [百度学术]
陶蓉,潘根,江玲,程遐年,刘裕强,万建民. 籼稻品种'BP360e'抗褐飞虱遗传分析及基因定位. 南京农业大学学报, 2019, 42(1): 14-20 [百度学术]
Tao R,Pan G,Jiang L,Cheng D N,Liu Y Q,Wan J M. Genetic analysis and gene mapping of brown planthopper resistance in Indica rice variety 'BP360e'. Journal of Nanjing Agricultural University, 2019, 42(1): 14-20 [百度学术]
Zhao Y,Huang J,Wang Z Z,Jing S L,Wang Y,Ouyang Y D,Cai B D,Xin X F,Liu X,Zhang C X,Pan Y F,Ma R,Li Q F,Jiang W H,Zeng Y,Shangguan X X,Wang H Y,Du B,Zhu L L,Xu X,Feng S Q,He S Y,Chen R Z,Zhang Q F and He G C. Allelic diversity in an NLR gene BPH9 enables rice to combat planthopper variation. Proceedings of the National Academy of Sciences of the United States of America, 2016, 113(45): 12850-12855 [百度学术]
Guo J P,Xu C X,Wu D,Zhao Y,Qiu Y F,Wang X X,Ouyang Y D,Cai B D,Liu X,Jing S L,Shangguan X X,Wan H Y,Ma Y H,Hu L,Wu Y,Shi S J,Wang W L,Zhu L L,Xu X,Chen R Z,Feng Y Q,Du B, He G C. Bph6 encodes an exocyst-localized protein and confers broad resistance to planthoppers in rice. Nature Genetics, 2018, 50(2): 297-306 [百度学术]
刘光杰,郑宜才,桂丽琴,沈君辉. 水稻品种抗稻飞虱鉴定方法的比较研究. 浙江农业学报, 1999(11): 306-310 [百度学术]
Liu G J,Zheng Y C,Gui L Q,Shen J H. A comparative study of the rice planthopper resistance identification method. Acta Agriculturae Zhejiangensis, 1999(11): 306-310 [百度学术]
Ramalingam J,Vera Cruz C M,Kukreja K,Chittoor J M,Wu J L,Lee S W,Baraoidan M,George M L, Cohen M B,Hulbert S H. Candidate defense genes from rice, barley, land maize and their association with qualitative and quantitative resistance in rice. Molecular Plant-Microbe Interactions, 2003, 16(1): 14-24 [百度学术]
Chen J W,Wang L,Pang X F,Pan Q H. Genetic analysis and fine mapping of a rice brown planthopper (Nilaparvata lugens Stål) resistance gene bph19(t). Molecular Genetics and Genomics , 2006, 275: 321-329 [百度学术]
Monna L,Lin H X,Kojima S,Sasaki T,Yano M. Genetic dissection of a genomic region for quantitative trait locus Hd3, into two loci, Hd3a and Hd3b, controlling heading date in rice. Theoretical and Applied Genetics, 2002, 104: 772-778 [百度学术]
Wang G L,Mackill D J,Bonman J M,McCouch S R,Champoux M C,Nelson R J. RFLP mapping of genes conferring complete and partial resistance to blast in a durably resistant rice cultivar. Genetics, 1994, 136: 1421-1434 [百度学术]