删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

烟草表达抗病基因同源物(RGA)的鉴定及RGA-SSR标记的开发

本站小编 Free考研考试/2021-12-26

袁清华, 谢锐鸿, 张振臣, 马柱文, 李集勤, 李淑玲, 陈俊标*
广东省农业科学院作物研究所, 广东广州 510640
* 通讯作者(Corresponding author): 陈俊标, E-mail:jbchen2007@126.com, Tel: 86-20-87596595 收稿日期:2013-07-25 基金:本研究由广东省烟草专卖局科技计划项目(201106), 广东省科技计划项目(2010B020302004)和广东省烟草专卖局(公司)科技项目(200905)资助。

摘要烟草是研究植物与病原菌互作的理想材料。鉴定烟草抗病基因及其同源物对揭示抗病机制具有重要意义。近年来公共数据库不断增长的EST序列为烟草表达RGA的鉴定提供丰富的数据。本研究通过拼接GenBank收录的412 325条烟草EST序列, 获得149 606条Uni-EST序列。随后利用已克隆的112个植物R基因蛋白序列对其扫描, 检测出1113个NtRGA, 其中有273、546、53、102和30个分别包含NBS-LRR、LRR-PK、LRR、PK和Mlo结构域, 另有109个未检测到结构域。通过序列比对将1071个NtRGA定位于N. benthamiana基因组712个位点上。经搜索, 从72个NtRGA中检测出78个SSR, 根据其侧翼序列设计64对引物。54对成功从烟草基因组DNA中扩增出清晰条带, 9对在24个普通烟草品种间检测出多态性, 检出等位基因数2~4个, 平均2.56个; 41对在6个烟草种间检测出多态性, 检出等位基因数2~4个, 平均2.61个。

关键词:烟草; 表达序列标签; 抗病基因同源物; SSR
Identification of Expressed Resistance Gene Analogues (RGAs) and Development of RGA-SSR Markers inNicotiana
YUAN Qing-Hua, XIE Rui-Hong, ZHANG Zhen-Chen, MA Zhu-Wen, LI Ji-Qin, LI Shu-Ling, CHEN Jun-Biao*
Crops Research Institute, Guangdong Academy of Agricultural Sciences, Guangzhou 510640, China

AbstractTobacco is an ideal experimental system for studing on plant-pathogen interaction. Identification of tobaccoR gene and resistance gene analogs is propitious to elucidating the underlying resistant mechanisms. In recent years, the growing public tobacco EST data provide rich source for identifying expressed RGA. In this study, 149 606 Uni-EST were assembled from 412 325 ESTs of tobacco in GenBank. By scanning the Uni-EST with 112 plant R gene protein sequences 1113 NtRGAs were identified. These expressed RGAs comprised 273, 546, 53, 102, and 30 of NBS-LRR, LRR-PK, LRR, PK, and Mlo domains encodingR genes, respectively. No domain was detected in the rest of 109 RGAs. By aligning sequence 1079 NtRGAs were allocated on 712 loci inN. benthamianagenome. A total of 78 simple sequence repeats (SSRs) were identified from 72 NtRGAs. Sixty-four primer pairs were designed base on the flanking sequence of SSR. Among them, 54 primer pairs were amplified with clear bands from tobacco genomic DNA. Nine primer pairs were detected to have polymorphism among 24 varieties ofNicotiana tabacum with two to four alleles (on average 2.56 alleles). Forty-one primer pairs were detected to have polymorphism among six species inNicotiana with two to four alleles (on average 2.61 alleles).

Keyword:Nicotiana; Expressed sequence tags; Resistance gene analogs; SSR
Show Figures
Show Figures












烟草是我国乃至世界重要的经济作物。农业生产中, 烟草青枯病、花叶病毒病和黑胫病等严重影响烟叶产量和品质。据全国烟草病虫害预测预报及综合防治网统计, 2010年和2011年, 全国16个烤烟主产省烟草病虫害发生面积达80万公顷, 产量损失6万吨, 产值损失7亿元。因此, 防治烟草病虫害对烟叶生产具有重要作用, 而烟草抗病基因( R基因)和抗病基因同源物(RGA)的克隆则对揭示烟草抗病机制以及制定正确的病虫害防治措施具有重要的意义。
植物 R基因对病原菌无毒基因编码蛋白的识别起关键作用[ 1]。近年来, 通过图位克隆或转座子标签法克隆的植物 R基因超过100个[ 2, 3, 4, 5](http://prgdb. crg.eu/wiki/Species_with_ R-genes)。尽管植物 R基因克服的病原物种类不同, 其编码的氨基酸序列却存在类似结构域, 如核苷酸结合位点(nucleotide binding site, NBS)、富含亮氨酸重复序列(leucine-rich repeats, LRR)、丝氨酸/苏氨酸激酶(serine-threonine kinase, STK)、亮氨酸拉链结构(leucine zippers, LZ)、跨膜结构域(transmembrane domain, TM)、Toll白介素-1区域(toll-interleukin-1 region, TIR)等[ 6, 7, 8, 9]。这些保守结构域为 R基因和RGA的快速鉴定克隆提供捷径。
根据氨基酸序列保守结构域可将抗病基因分为6大类[ 2]。第1类含有NBS和LRR结构域(NBS- LRR)。基于N端是否缺失TIR, 此类 R基因又可分为2个亚类, 即TIR-NBS-LRR和non-TIR-BS-LRR, 例如抗烟草花叶病毒的N基因含有TIR-NBS-LRR结构域[ 7, 10], 而抗丁香假单胞杆菌( Pseudomonas syringae)的拟南芥 Rps2基因则含有CC-NBS-LRR结构域[ 11]。第2类含有LRR和PK结构域(LRR-PK), 例如拟南芥 Fls2和水稻 Xa21基因[ 12, 13]。第3类含有一个胞外LRR结构域, 例如拟南芥 RPP27基因[ 14]。第4类只含有PK结构域, 例如番茄 Pto和香瓜 At1基因[ 15, 16]。第5类包含具有不同抗性机制的其他 R基因, 如玉米 Hm1和大麦 Mlo基因[ 3, 17]
早期RGA的分离主要通过PCR扩增 R基因保守结构域, 此法已成功从拟南芥[ 18, 19]、大豆[ 20]、水稻[ 21]、玉米[ 22]、小麦[ 23]、烟草[ 24, 25]和其他植物[ 26, 27, 28]中分离出RGA。Bertioli等[ 29]根据NBS设计简并引物, 从花生栽培种和野生种中分离到78个RGA。Gao等[ 25]基于NBS和PK结构域从烟草品种 Nicotiana repanda分离出100个RGA。
与PCR扩增相比, 采用数据挖掘方法从基因组中鉴定RGA效率更高, 基于已公布的植物全基因组序列, 迄今已成功从拟南芥[ 10]、水稻[ 30]、蒺藜苜蓿[ 31]和百脉根[ 32]基因组中鉴定大量RGA。Meyers等[ 10]从拟南芥基因组中鉴定出149个NBS-LRR编码基因和58个其他相关基因。Ameline-Torregrosa等[ 31]从蒺藜苜蓿全基因组草图中鉴定出333个非冗余NBS- LRR, 并预测全基因组含有400~500个NBS-LRR基因。近期, Li等[ 32]从百脉根基因组中成功分离出158个NBS序列。
由于植物基因组中存在大量丧失生物学功能的假基因, 而表达基因占植物基因组比例不到1%, 根据基因组DNA序列开发的RGA多为不表达的假基因[ 32], 严重干扰 R基因的克隆。因此, 鉴定具有表达活性的RGA对 R基因的克隆显得尤为重要。近期, 通过数据挖掘已从植物EST序列中鉴定出系列RGAs。Liu等[ 33]利用454技术测序获得169万条菜豆EST, 随后从中鉴定出364个RGA。Liu等[ 34]利用54个植物R基因氨基酸序列扫描GenBank收录的花生EST序列, 成功鉴定出385个表达花生RGAs。
RGA鉴定后常被用于开发RGA分子标记。基于RGA开发的分子标记主要有RFLP[ 35]、STS[ 36]、SSCP[ 37]、CAPS[ 38]和SSR[ 34]标记等。Sanz等[ 35]根据燕麦RGA序列设计31个RFLP探针, 并成功将53个RGA-RFLP标记定位到燕麦六倍体遗传图谱上( Avena byzantina cv. Kanota × Avena sativa cv. Rollo)。近期, Liu等[ 34]从25个花生RGA中开发出28个SSR标记, 并将其中一个标记RGA121定位到连锁群AhIV上。由于SSR标记具有共显性、多态性高、操作简单且重复性好等优点, 与其他标记相比, 从RGA中开发SSR更有利于抗病基因的定位与克隆。
目前公共核酸数据库GenBank收录的烟草EST序列已达412 325条(截至2012年7月20日), 几乎涵盖不同时期, 不同组织表达的基因, 为烟草表达RGA鉴定提供丰富的资源。鉴于此, 本研究拟利用数据挖掘方法从烟草EST序列中鉴定表达RGA, 继而开发RGA-SSR标记。这将为烟草 R基因的定位与克隆奠定基础。
1 材料与方法1.1 植物材料普通烟草品种24个(表1), 其他烟草种6个, 种植于广东省农业科学院试验农场, 于生长期取叶片提取基因组DNA。
1.2 烟草EST序列拼接采用TIGR基因指数聚类工具(TIGR gene indices clustering tools, 简称TGICL,http://compbio.dfci.harvard.edu/tgi/software/), 对GenBank收录的EST序列进行拼接。当EST序列重叠超过50个核苷酸且相似度达到90%, 同时序列不匹配长度不超过20个核苷酸时被认为符合拼接要求。
表1
Table 1
表1(Table 1)
表1 供试烟草栽培品种 Table 1 Varieties and wild species of tobacco used in this study
编号 No.材料
Material
类型
Type
1ATNARELLO白肋烟 Burley
2KY26白肋烟 Burley
3大白筋599 Big white burley 599白肋烟 Burley
4白肋9号 White burley 9白肋烟 Burley
5阿波烟 Aboyan白肋烟 Burley
6Margland 609白肋烟 Burley
7广黄5号 Guanghuang 5烤烟 Flue-cured tobacco
8丰字1号 Fengzi 1烤烟 Flue-cured tobacco
9革新6号 Gexin 6烤烟 Flue-cured tobacco
10革杂基 Gezaji烤烟 Flue-cured tobacco
11红花大金元 Honghuadajinyuan烤烟 Flue-cured tobacco
12金星 Jinxing烤烟 Flue-cured tobacco
13白骨牛利 Baiguniuli晒烟 Sun-cured tobacco
14大秋根2 Daqiugen 2晒烟 Sun-cured tobacco
15大山沟 Dashangou晒烟 Sun-cured tobacco
16大同 Datong晒烟 Sun-cured tobacco
17大叶密合 Dayemihe晒烟 Sun-cured tobacco
18大种白毛 Dazhongbaimao晒烟 Sun-cured tobacco
192040雪茄烟 Cigar
2074-16雪茄烟 Cigar
21夏湾娜 Xiawanna雪茄烟 Cigar
22Black sea samsun香料烟 Oriental tobacco
23Cuban Havana雪茄烟 Cigar
24Dark Virginia烤烟 Flue-cured tobacco
25克利夫兰氏烟草
Nicotiana clevelandii A. Gray
多室烟草组
Nicotianasect. Polydicliae
26残波烟草
Nicotiana repandaWilld
残波烟草组
Nicotiana sect. Repandae
27迪勃纳氏烟草
Nicotiana debneyi Domin
香甜烟草组
Nicotiana sect . Suaveolentes
28黄花烟草 Nicotiana rustica L.黄花烟草组
Nicotiana sect . Rusticae
29绒毛状烟草
Nicotiana tomentosiformisGood sp.
绒毛烟草组
Nicotiana sect. Tomentosae
30林烟草 Nicotiana sylvestris Speg. & Comes林烟草组
Nicotianasect. Sylvestres

表1 供试烟草栽培品种 Table 1 Varieties and wild species of tobacco used in this study

1.3 数据挖掘烟草RGA利用112个已克隆的植物 R基因(表2)氨基酸序列扫描Uni-EST, 从中鉴定表达RGAs。序列比对采用tBlastn工具, 比对分值大于或等于100, 且E值小于或等于1E-10的Uni-EST被认为是候选烟草RGAs (NtRGAs)。
1.4 NtRGA在 N. benthamiana基因组中的定位 利用Blast工具将NtRGA通过序列比对定位在本塞姆氏烟草( Nicotiana benthamiana)全基因组草图上(http://solgenomics.net/organism/Nicotiana_benthamiana/ nome)。根据比对结果选择比对分值最高、E值最小, 且比对分值大于50, E值小于1E-10的基因组序列作为RGA定位区段。
1.5 烟草RGA-SSR标记的开发采用Perl脚本MIcroSAtellite (MISA, http://pgrc. ipk-gatersleben.de/misa/)从RGA序列搜索SSR位点。设置以下SSR搜索限制因子: 最短15 bp, 二核苷酸基序最少重复8次, 三核苷酸基序至少重复5次, 四五六核苷酸基序至少重复4次, 不包含单核苷酸类型。基于SSR侧翼序列, 采用Primer Premier 5软件设计引物。新引物随后用于检测24个烟草栽培种和6个野生种的多态性。

2 结果与分析2.1 烟草EST序列拼接截至2012年6月7日, GenBank收录的烟草属EST序列共412 325条, 其中 N. tabacum334 384条, N. benthamiana 56 102条, N. langsdorffii× N. sanderae 12 448条, N. sylvestris 8583条, N. attenuata 355条, 其他烟草种453条。通过登陆GenBank, 将上述EST序列以FASTA格式下载保存, 作为烟草RGAs开发的原始数据。由于原始EST数据存在大量冗余序列, 为提高EST序列质量, 获得比原始EST序列更长, 且来自同一位点的共有序列, 采用TGICL工具拼接GenBank收录的烟草EST序列。结果由412 325条EST序列拼接出149 606条Uni-EST序列, 其中contigs 45 137条, singletons 101 169条, 序列最长为2312 bp, 最短为431 bp, 平均长度为874 bp。
2.2 NtRGAs鉴定利用tBlastn软件将112个已知的 R基因氨基酸序列与烟草Uni-EST序列比对分析表明, 比对E值小于或等于1E-10时有3个 R基因(RPW8.1, RPW8.2, xa27)无法检测到与之匹配的烟草Uni-EST, 其余109个 R基因共检测到6963条Uni-EST。由于具有相同结构域的 R基因序列存在相似性, 序列比对时具有相同结构域的 R基因常与同一Uni-EST序列匹配, 导致Uni-EST序列被重复计数。去除重复计数后, 共有1113个Uni-EST与109个 R基因序列匹配(数据未发表)。其中273个含NBS-LRR结构域, 546个含LRR-PK结构域, 53个含胞外LRR结构域, 102个只含PK结构域, 30个含Mlo结构域, 其余109个
表2
Table 2
表2(Table 2)
表2 112个用于本研究的植物 R基因 Table 2 One hundred and twelve known Rgenes from plants used in this study
物种
Species
R基因
Rgene
蛋白质ID
Protein ID
结构域
Structure
Aegilops tauschiiCre1AAM94164NBS, LRR, others
Arabidopsis thalianaEFRNP_197548LRR, PK
Arabidopsis thalianaER-ErectaNP_180201LRR, PK
Arabidopsis thalianaFLS2NP_199445LRR, PK
Arabidopsis thalianaHRTAAF36987NBS, LRR, PK, others
Arabidopsis thalianaPEPR1NP_177451LRR, PK
Arabidopsis thalianaRAC1AAS01763TIR, NBS, LRR, others
Arabidopsis thalianaRCY1BAC67706NBS, LRR, others
Arabidopsis thalianaRFO1NP_178085PK
Arabidopsis thalianaRLM3NP_001031652NBS, TIR, others
Arabidopsis thalianaRPM1NP_187360NBS, LRR, others
Arabidopsis thalianaRPP1NP_190034NBS, LRR, others
Arabidopsis thalianaRPP13NP_190237NBS, LRR, others
Arabidopsis thalianaRPP27CAE51863LRR
Arabidopsis thalianaRPP4NP_193420TIR, NBS, LRR, others
Arabidopsis thalianaRPP5NP_193428TIR, NBS, LRR, others
Arabidopsis thalianaRPP8NP_199160NBS, LRR, others
Arabidopsis thalianaRps2NP_194339NBS, LRR, others
Arabidopsis thalianaRps4NP_199338TIR, NBS, LRR, others
Arabidopsis thalianaRPS5NP_17268NBS, LRR, others
Arabidopsis thalianaRPW8.1AAK09266RPW8, others
Arabidopsis thalianaRPW8.2AAK09267RPW8, others
Arabidopsis thalianaRRS1NP_001078715NBS, LRR, others
Arabidopsis thalianaRTM1NP_172067
Arabidopsis thalianaRTM2NP_568144
Arabidopsis thalianaSSI4AAN86124TIR, NBS, LRR, others
Beta vulgarisHs1AAB48305
Capsicum annuumBs3ABW82012NBS, LRR, others
Capsicum annuumBs3-EABW82011NBS
Capsicum chacoenseBs2AAF09256NBS, LRR, others
Cucumis meloAt1AAL47679PK
Cucumis meloAt2AAL62332PK
Cucumis meloFOM-2ABB91438NBS, LRR, others
Cucumis meloVATNBS, LRR, others
Glycine maxKR1AAL56987TIR, NBS, LRR, others
Glycine maxRps1-k-1AAX89382NBS, LRR, others
Glycine maxRps1-k-2AAX89383NBS, LRR, others
Helianthus annuusPl8AAT08955NBS, LRR, others
Hordeum vulgareMLA1ACZ65507NBS, LRR, others
Hordeum vulgareMLA10AAQ55541NBS, LRR, others
Hordeum vulgareMLA13AAO16014NBS, LRR, others
Hordeum vulgareMloCAB06083Mlo
Hordeum vulgareRPG1ABK51311PK
Hordeum vulgareMla12AAO43441NBS, LRR, others

表2 112个用于本研究的植物 R基因 Table 2 One hundred and twelve known Rgenes from plants used in this study

(表2)
Table 2 Continued
(表2)(Table 2 Continued)
(表2) 112个用于本研究的植物 R基因 Table 2 Continued One hundred and twelve known Rgenes from plants used in this study
物种
Species
R基因
R gene
蛋白质ID
Protein ID
结构域
Structure
Hordeum vulgareMla6CAC29242NBS, LRR, others
Hordeum vulgareRdg2aADK47521NBS, LRR, others
Lactuca sativaDm3AAD03156NBS, LRR, others
Linum usitatissimumL6AAA91022TIR, NBS, LRR, others
Linum usitatissimumMAAB47618TIR, NBS, LRR, others
Linum usitatissimumP2AAK28805TIR, NBS, LRR, others
Nicotiana benthamianaSerk3AADO86982LRR, PK
Nicotiana benthamianaSerk3BADO86983LRR, PK
Nicotiana glutinosaNAAA50763TIR, NBS, LRR, others
Nicotiana tabacumIVRCAA08776
Oryza sativaPIBBAA76282NBS, LRR, others
Oryza sativaPi-taAAO45178NBS, LRR, others
Oryza sativaXA1BAA25068NBS, LRR, others
Oryza sativaxa21BAE93934LRR, PK
Oryza sativaPi2ABC94598NBS, LRR, others
Oryza sativaPi36ABI64281NBS, LRR, others
Oryza sativaPi9ABC18336NBS, LRR, others
Oryza sativaPid2ACR15163PK
Oryza sativaXa13ABD78944
Oryza sativaXa27AAY54163
Oryza sativaXa5AAV53715
Oryza sativaPi5-1ACJ54697NBS, LRR, others
Oryza sativaPi5-2ACJ54698NBS, LRR, others
Oryza sativaPid3ACN79514NBS, LRR, others
Oryza sativaPikm1-TSBAG72135NBS, LRR, others
Oryza sativaPikm2-TSBAG72136NBS, LRR, others
Oryza sativaPikp-2ADV58351NBS, LRR, others
Oryza sativaPitBAH20862NBS, LRR, others
Oryza sativaPiz-tABC73398NBS, LRR, others
Oryza sativaXa26ABD84047LRR, PK
Phaseolus vulgarisPGIPCAA46016LRR, PK
Solanum acauleRx2CAB56299NBS, LRR, others
Solanum bulbocastanumRpi-blb1AAP86601NBS, LRR, others
Solanum bulbocastanumRpi-blb2AAZ95005NBS, LRR, others
Solanum demissumR1AAL39063NBS, LRR, others
Solanum habrochaitesCf-4CAA05270LRR, PK
Solanum habrochaitesCf4ACAA05265LRR, PK
Solanum lycopersicumI-2AAD27815NBS, LRR, others
Solanum lycopersicumAsc-1AAF67518
Solanum lycopersicumBs4AAR21295TIR, NBS, LRR, others
Solanum lycopersicumHeroCAD29729NBS, LRR, others
Solanum lycopersicumLeEIX1AAR28377LRR, PK
Solanum lycopersicumLeEIX2AAR28378LRR, PK
Solanum lycopersicumMi1.2AAC67238NBS, LRR, others

(表2) 112个用于本研究的植物 R基因 Table 2 Continued One hundred and twelve known Rgenes from plants used in this study

(表2)
Table 2 Continued
(表2)(Table 2 Continued)
(表2) 112个用于本研究的植物 R基因 Table 2 Continued One hundred and twelve known Rgenes from plants used in this study
物种
Species
R基因
R gene
蛋白质ID
Protein ID
结构域
Structure
Solanum lycopersicumSw-5AAG31013NBS, LRR, others
Solanum lycopersicumTm-2AAQ10735NBS, LRR, others
Solanum lycopersicumTm-2aAAQ10736NBS, LRR, others
Solanum lycopersicumVe1AAK58682LRR, PK
Solanum lycopersicumVe2AAK58011LRR, PK
Solanum lycopersicumCf-5AAC78591LRR, PK
Solanum pimpinellifoliumCf-2AAC15779LRR, PK
Solanum pimpinellifoliumPrfAAF7630NBS, LRR, others
Solanum pimpinellifoliumPtoAAC48914PK
Solanum pimpinellifoliumCf-9CAA05277LRR, PK
Solanum pimpinellifoliumCf9BCAA05273LRR, PK
Solanum tuberosumGpa2AAF04603NBS, LRR, others
Solanum tuberosumGro1.4AAP44390TIR, NBS, LRR, others
Solanum tuberosumR3aAAW48299NBS, LRR, others
Solanum tuberosumRxCAB50786NBS, LRR, others
Solanum tuberosumRY-1CAC82811TIR, NBS, LRR, others
Triticum aestivumLr1ABS29034NBS, LRR, others
Triticum aestivumLr10AAQ01784NBS, LRR, others
Triticum aestivumLr21ACO53397NBS, LRR, others
Triticum aestivumLr34ADK62371NBS, LRR, others
Triticum aestivumPm3AAQ96158NBS, LRR, others
Zea maysHm1NP_001105920
Zea maysHm2ABY68564
Zea maysRp1-DAAD47197NBS, LRR, others

(表2) 112个用于本研究的植物 R基因 Table 2 Continued One hundred and twelve known Rgenes from plants used in this study

未发现结构域。将与 R基因序列匹配的Uni-EST初步确定为烟草RGA, 并暂命名为NtRGA。
2.3 NtRGA在本塞姆氏烟草( Nicotiana benthamiana)基因组中的定位 将NtRGA定位到烟草基因组上对分离特定候选抗病基因/QTL具有重要意义。随着 N. benthamiana全基因组草图的公布, 通过序列比对可直接将NtRGA定位到 N. benthamiana基因组上。用Blastn工具在 N. benthamiana全基因组草图中搜索与NtRGAs序列匹配的位点表明, 当分值大于50且E值小于1E-10时, 1113个NtRGA有1071个在 N. benthamiana基因组中找到相似的序列, 其中965个NtRGAs (90.7%)同时与基因组中多个片段匹配, 平均每个NtRGA与9.67个基因组片段匹配, 最多的一个(CL7158Contig2)与529个基因组片段匹配。鉴于植物基因在长期的进化过程中普遍存在多拷贝现象, 本研究取分值最高且E值最低的基因组序列作为NtRGA的基因组位点。最终1071个NtRGA被定位到712个基因组位点上, 共有218个基因组位点存在一个位点与多个NtRGAs配对(数据未发表)。进一步分析发现, 与同一个基因组位点配对的NtRGAs存在两种类型: (1)序列高度同源, 与基因组同一区段配对(图1); (2)序列无同源性, 与基因组同一位点不同区段配对(图2)。由于真核生物的外显子存在选择性剪切, 同一段基因序列可能会转录成不同的mRNA。属于类型1的NtRGA可能来自同一基因的不同转录本, 同时也有可能是不同烟草种间的同源基因。而类型2的NtRGA则可能来自同一基因的不同部位。EST通常只是完整基因的一部分, 因此, 来自相同基因的EST序列可能由于缺乏重叠区而无法被拼接在一起。
观察NtRGA在 N. benthamiana基因组中的位置表明, NtRGA在基因组中分布不均匀, 共有17处基因组区段存在多个NtRGA串联。例如全长63 kb的基因组scaffold (序列ID: Niben.v0.3.Scf25265845)包含4个NtRGA。
图1
Fig. 1
Figure OptionViewDownloadNew Window
图1 CL1421Contig1、CL1421Contig2与 N. benthamiana基因组序列比对结果A: CL1421Contig1与基因组scaffold (序列ID: Niben.v0.3.Scf24878592)的1078~1448 bp和1547~2749 bp区段匹配;B: CL1421Contig2与基因组支架Scf24878592的1276~1449 bp和1549~2786 bp区段匹配。Fig. 1 Alignment of CL1421Contig1 and CL1421Contig2 to N. benthamianagenomeA: CL1421Contig1 matches with 1078-1448 bp and 1547-2749 bp in genomic scaffold (sequence ID: Niben.v0.3.Scf24878592);B: CL1421Contig2 matches with 1276-1449 bp and 1549-2786 bp in genomic scaffold (sequence ID: Niben.v0.3.Scf24878592).


2.4 SSR标记的开发利用Perl脚本MISA从1113条NtRGA序列搜索SSR, 共发现78个SSR位点, 分布于72条序列, 其中含有SSR一个以上的NtRGA序列6条, 平均每 939 348 bp含1个SSR。根据SSR侧翼序列, 利用Primer Premier 5软件设计SSR引物, 共设计引物64对, 8条NtRGA序列由于SSR侧翼序列太短或者结构复杂而未能成功设计引物。新引物随后在24个烟草栽培品种和6个野生种中扩增, 并通过聚丙烯酰胺凝胶电泳检测其多态性(表3)。结果发现, 供试引物有54对在栽培种中能扩增出清晰条带, 其余10对无扩增产物或产生类似RAPD的非特异性条带。在可扩增引物中, 46对扩增产物长度与预期相近, 7对扩增产物长度大于预期, 另有1对扩增片段比预期短。可扩增的54对引物中有9对在烟草栽培品种间检测出多态性, 占可扩增引物的16.7%。检测等位基因数2~4个, 总共23个, 平均2.56个。54对在烟草栽培种中可扩增SSR引物在6个烟草野生种中均成功扩增。与栽培种相反, SSR引物在野生种中检测出较高的多态性, 共有41对引物检测出多态性, 占可扩增引物的75.9%, 检测等位基因数为2~4个, 总共92个, 平均2.61个。图3为引物RGA-63在烟草中的扩增带型, RGA-63在24个普通烟草资源和6个其他烟草种中均检测出4个等位基因。
图2
Fig. 2
Figure OptionViewDownloadNew Window
图2 FS423566和CL17188Contig1与烟草基因组序列比对结果A: FS423566与基因组scaffold (序列ID: Niben.v0.3.Scf24993858)的4904~5308 bp区段匹配; B: CL17188Contig1与基因组scaffold (序列ID: Niben.v0.3.Scf24993858)的5494~6282 bp区段匹配。Fig. 2 Alignment of FS423566 and CL17188Contig1 to N. benthamianagenomeA: FS423566 matches with 4904-5308 bp in genomic scaffold (sequence ID: Niben.v0.3.Scf24993858); B: CL17188Contig1 matches with 5494-6282 bp in genomic scaffold (sequence ID: Niben.v0.3.Scf24993858).


3 讨论作为“基因对基因”抗病机制的组成部分, 植物 R基因在识别病原物特异无毒基因编码产物中起重要作用[ 39]。本研究中从GenBank收录的烟草EST数据中成功鉴定出1113个RGA, 并通过序列比对将其定位到 N. benthamiana基因组上, 说明利用EST数据能够高效鉴定RGA。
过去几年, 通过数据挖掘已成功从甘蔗、小麦、玉米等作物鉴定RGA。Dilbirligi等[ 40]采用4种不同搜索策略在小麦中查找RGA序列, 包括结构域搜索, 单个或多个基序搜索, 共有序列搜索以及单一全长序列搜索。结果表明单一全长序列搜索效果最好。当E值≤ E-10, 采用单一全长序列共检测到243个含NBS-LRRs和101个其他类型的RGA序列。Xiao等[ 41]利用改良AFLP、RACE和数据挖掘方法鉴定玉米RGA和类 R基因EST时, 发现数据挖掘方法最为有效。Rossi等[ 42]通过严格的比对条件(E值表3个主要 R基因簇, 即NBS-LRR, LRR-TM和PK。上述研究报告表明不同E值对RGA搜索结果影响很大。本研究在E值≤E-10时, 共鉴定出1113个RGA, 其数量是小麦的3倍[ 40]。两者搜索结果差异大的原因有二, 首先是用于挖掘NtRGA的烟草EST数量多于小麦, 本研究共使用412 325条烟草EST序列, 而Dilbirligi等[ 40]的研究仅使用78 221条小麦EST序列。其次是用于比对分析的 R基因数量不同, 本研究采用112个 R基因进行比对, 远超过用于小麦比对的 R基因数量(22个)。
通过EST挖掘RGA的一个优点是所检测到的RGA均为表达的基因, 相反, 通过基因组序列挖掘的RGA存在不表达的假基因。Li等[ 32]从百脉根基因组鉴定的RGA中发现有65个编码不完整的蛋白序列, 最终被确定为假基因。Meyers等[ 10]从拟南芥基因组中搜索含NBS-LRR结构域的RGA, 发现至少有12个NBS-LRR基因由于框架移动和无义突变而演变为假基因。
表3
Table. 3
Figure OptionViewDownloadNew Window
表3 烟草RGA-SSR的基本特征Table. 3 Characteristic of RGA-SSR in Nicotiana

表3
Table. 3 Continued
Figure OptionViewDownloadNew Window
表3 烟草RGA-SSR的基本特征Table. 3 Continued Characteristic of RGA-SSR in Nicotiana

图3
Fig. 3
Figure OptionViewDownloadNew Window
图3 引物RGA-63在烟草中的扩增带型1: ATNARELLO; 2: KY26; 3: 大白筋599; 4: 白肋9号; 5: 阿波烟; 6: Margland 609; 7: 广黄5号; 8: 丰字1号; 9: 革新6号; 10: 革杂基; 11: 红花大金元; 12: 金星; 13: 白骨牛利; 14: 大秋根2; 15: 大山沟; 16:大同; 17: 大叶密合; 18: 大种白毛; 19: 2040; 20:74-16; 21: 夏湾娜; 22: Black Sea Samsun; 23: Burley Hampton; 24: Dark Virginia; 25: N. clevelandii A. Gray; 26: N. repandaWilld; 27: N. debneyi Domin; 28: N. rustica L.; 29: N. tomentosiformis; 30: N. sylvestris Speg. & Comes。Fig. 3 Amplification results of RGA-63 in Nicotiana1: ATNARELLO; 2: KY26; 3: Big White Burley 599; 4: White Burley 9; 5: Aboyan; 6: Margland 609; 7: Guanghuang 5; 8: Fengzi 1; 9: Gexin 6;10: Gezaji; 11: Honghuadajinyuan; 12: Jinxing; 13: Baiguniuli; 14: Daqiugen 2; 15: Dashangou; 16: Datong; 17: Dayemihe; 18: Dazhongbaimao; 19: 2040; 20: 74-16; 21: Xiawangna; 22: Black Sea Samsun; 23: Burley Hampton; 24: Dark Virginia; 25: N. clevelandii A. Gray; 26: N. repandaWilld; 27: N. debneyi Domin; 28: N. rustica L.; 29: N. tomentosiformis; 30: N. sylvestris Speg. & Comes.

由于目前多数植物 R基因是通过图位克隆和转座子标签法克隆得到[ 43], 因此, 将RGA定位到基因组上显得十分重要。本文通过序列比对将其中的1071个NtRGA定位到 N. benthamiana基因组中712个位点上, 为后续烟草 R基因的克隆奠定了基础。本研究尚有42个NtRGA无法定位于 N. benthamiana基因组, 推测有以下原因: (1)目前公布的 N. benthamiana全基因组序列仅为草图, 尚有部分基因组区段未测序; (2)本研究鉴定的NtRGA仅部分来自 N. benthamiana, 来自其他烟草种的NtRGA由于种间基因组差异而无法定位。本研究通过分析NtRGA在基因组中的位置发现, RGA在基因组中的分布并不均匀, 部分RGA成簇存在。这与前人在其他植物中的研究结果相似[ 44, 45] R基因的串联方便 R基因的遗传变异, 有利于新 R基因的进化。Bertioli等[ 46]通过分析花生、百脉根和苜蓿发现反转座子与一些抗病基因簇关联。其他几个假设例如复制、基因转换和不等位交换已被用于解释 R基因成簇和进化[ 1]
如何开发利用RGA0序列是RGA鉴定后需要解决的问题。本文通过SSR搜索从RGA序列中开发出62个RGA-SSR标记, 为后续抗病基因的定位克隆提供了有利的帮助。近期的研究表明, SSR在普通烟草具有较高的多态性, Bindler等[ 47]利用SSR构建普通烟草遗传图谱, 5119对SSR引物中有2415对(47%)在亲本间检测出多态性。本研究仅有16.7%的引物在普通烟草中检测到多态性, 且检测的等位基因数仅有2~4个。本研究开发的SSR标记多态性低的原因可能与标记来自表达基因序列有关。由于表达基因维持其生物学功能, 在进化过程中承受的选择压力远高于全基因的平均水平, 导致其序列高度保守。
4 结论利用植物 R基因蛋白质序列扫描EST数据是鉴定烟草表达RGA的有效途径。烟草表达RGA在基因组中分布不均匀, 其结构域以LRR-PK和NBS- LRR为主。烟草RGA-SSR在种内缺乏多态性, 但在种间具有较高的多态性。
The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.

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

1Ellis J, Dodds P, Pryor T. The generation of plant disease resistance gene specificities. Trends Plant Sci, 2000, 5: 373-379[本文引用:2][JCR: 11.808]
2Sanseverino W, Roma G, Simone M D, Faino L, Melito S, Stupka E, Frusciante L, Ercolano M R. PRGdb: a bioinformatics platform for plant resistance gene analysis. Nucl Acids Res, 2009, 38: D815[本文引用:2]
3Johal G S, Briggs S P. Reductase activity encoded by the Hm1 disease resistance gene in maize. Science, 1992, 258: 985-987[本文引用:2]
4Whitham S, McCormick S, Baker B. The N gene of tobacco confers resistance to Tobacco mosaic virus in transgenic tomato. Proc Natl Acad Sci USA, 1996, 93: 8776-8781[本文引用:1][JCR: 9.737]
5Dixon M S, Jones D A, Keddie J S, Thomas C M, Harrison K, Jones J D G, Lane C. The tomato Cf-2 disease resistance locus comprises two functional genes encoding leucine-rich repeat protein. Cell, 1996, 84: 451-459[本文引用:1][JCR: 31.957]
6Bent A F. Plant disease resistance genes: function meets structure. Plant Cell, 1996, 8: 1757-1771[本文引用:1][JCR: 9.251]
7Meyers B C, Dickerman A W, Michelmore R W, Sivaramakrishnan S, Sobral B W, Young N D. Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide-binding superfamily. Plant J, 1999, 20: 317-332[本文引用:2][JCR: 6.582]
8Hulbert S H, Webb C A, Smith S M, Sun Q. Resistance gene complexes: evolution and utilization. Annu Rev Phytopathol, 2001, 39: 285-312[本文引用:1][JCR: 10.229]
9Dangl J L, Jones J D. Plant pathogens and integrated defence responses to infection. Nature, 2001, 411: 826-833[本文引用:1][JCR: 38.597]
10Meyers B C, Kozik A, Griego A, Kuang H, Michelmore R W. Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis. Plant Cell, 2003, 15: 809-834[本文引用:4][JCR: 9.251]
11Bent A F, Kunkel B N, Dahlbeck D, Brown K L, Schmidt R, Giraudat J, Leung J, Staskawicz B J. RPS2 of Arabidopsis thaliana: a leucine-rich repeat class of plant disease resistance genes. Science, 1994, 265: 1856-1860[本文引用:1]
12Dunning F M, Sun W, Jansen K L, Helft L, Bent A F. Identification and mutational analysis of Arabidopsis FLS2 leucine-rich repeat domain residues that contribute to flagellin perception. Plant Cell, 2007, 19: 3297-3313[本文引用:1][JCR: 9.251]
13Song W Y, Wang G L, Chen L L, Kim H S, Pi L Y, Holsten T, Gardner J, Wang B, Zhai W X, Zhu L H, Fauquet C, Ronald P. A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. Science, 1995, 270: 1804-1806[本文引用:1]
14Tor M, Brown D, Cooper A, Woods-Tor A, Sjoland er K, Jones J D, Holub E B. Arabidopsis downy mildew resistance gene RPP27 encodes a receptor-like protein similar to CLAVATA2 and tomato Cf-9. Plant Physiol, 2004, 135: 1100-1112[本文引用:1][JCR: 6.555]
15Martin G B, Brommonschenkel S H, Chunwongse J, Frary A, Ganal M W, Spivey R, Wu T, Earle E D, Tanksley S D. Map-based cloning of a protein kinase gene conferring disease resistance in tomato. Science, 1993, 262: 1432-1436[本文引用:1]
16Taler D, Galperin M, Benjamin I, CohenY, Kenigsbuch D. Plant eR genes that encode photorespiratory enzymes confer resistance against disease. Plant Cell, 2004, 16: 172-184[本文引用:1][JCR: 9.251]
17Buschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, van Daelen R, van der Lee T, DieRGLarde P, Groenendijk J, Topsch S, Vos P, Salamini F, Schulze-Lefert P. The barley Mlo gene: a novel control element of plant pathogen resistance. Cell, 1997, 88: 695-705[本文引用:1][JCR: 31.957]
18Botella M A, Coleman M J, Hughes D E, Nishimura M T, Jones J D G, Somerville S C. Map positions of 47 Arabidopsis sequences with sequence similarity to disease resistance genes. Plant J, 1997, 12: 1197-1211[本文引用:1][JCR: 6.582]
19Aarts M G M, Hekkert B L, Holub E B, Beynon J L, Stiekema W J, Pereira A. Identification of R-gene homologous DNA fragments genetically linked to disease resistance loci in Arabidopsis thaliana. Mol Plant Microbe Interact, 1998, 11: 251-258[本文引用:1][JCR: 4.307]
20Graham M A, Marek L F, Lohnes D, Cregan P, Shoemaker R C. Expression and genome organization of resistance gene analogs in soybean. Genome, 2000, 43: 86-93[本文引用:1][JCR: 1.668]
21Mago R, Nair S, Mohan M. Resistance gene analogues from rice: cloning, sequencing and mapping. Theor Appl Genet, 1999, 99: 50-57[本文引用:1][JCR: 3.658]
22Collins N C, Webb C A, Seah S, Ellis J G, Hulbert S H, Pryor A. The isolation and mapping of disease resistance gene analogs in maize. Mol Plant Microbe Interact, 1998, 11: 968-978[本文引用:1][JCR: 4.307]
23Seah S, Sivasithamparam K, Karakousis A, Lagudah E S. Cloning and characterization of a family of disease resistance gene analogs from wheat and barley. Theor Appl Genet, 1998, 97: 937-945[本文引用:1][JCR: 3.658]
24Leng X, Xiao B, Wang S, Gui Y, Wang Y, Lu X, Xie J, Li Y, Fan L. Identification of NBS-type resistance gene homologs in tobacco genome. Plant Mol Biol Rep, 2010, 28: 152-161[本文引用:1][JCR: 5.319]
25Gao Y L, Xu Z L, Jiao F C, Yu H Q, Xiao B G, Li Y P, Lu X P. Cloning, structural features, and expression analysis of resistance gene analogs in tobacco. Mol Biol Rep, 2010, 37: 345-354[本文引用:2][JCR: 2.506]
26Wan H J, Zhao Z G, Malik A, Qian C T, Chen J F. Identification and characterization of potential NBS-encoding resistance genes and induction kinetics of a putative cand idate gene associated with downy mildew resistance in Cucumis. BMC Plant Biol, 2010, 10: 186[本文引用:1][JCR: 4.354]
27Huettel B, Santra D, Muehlbauer J, Kahl G. Resistance gene analogues of chickpea (Cicer arietinum L. ): isolation, genetic mapping and association with a Fusarium resistance gene cluster. Theor Appl Genet, 2002, 105: 479-490[本文引用:1][JCR: 3.658]
28Nair R A, Thomas G. Isolation, characterization and expression studies of resistance gene cand idates (RGCs) from Zingiber spp. Theor Appl Genet, 2007, 116: 123-134[本文引用:1][JCR: 3.658]
29Bertioli D J, Leal-Bertioli S C, Lion M B, Santos V L, Pappas G, Cannon S B, Guimaraes P M. A large scale analysis of resistance gene homologues in Arachis. Mol Gen Genomics, 2003, 270: 35-45[本文引用:1]
30Monsi B, Wisser R J, Pennill L, Hulbert S H. Full-genome analysis of resistance gene homologues in rice. Theor Appl Genet, 2004, 109: 1434-1447[本文引用:1][JCR: 3.658]
31Ameline-Torregrosa C, Wang B B, O’Bleness M S, Deshpand e S, Zhu H Y, Roe B, Young N D, Cannon S B. Identification and characterization of nucleotide-binding site-leucine-rich repeat genes in the model plant Medicago truncatula. Plant Physiol, 2008, 146: 5-21[本文引用:2][JCR: 6.555]
32Li X Y, Cheng Y, Ma W, Zhao Y, Jiang H Y, Zhang M. Identification and characterization of NBS-encoding disease resistance genes in Lotus japonicus. Plant Syst Evol, 2010, 289: 101-110[本文引用:4][JCR: 1.312]
33Liu Z, Crampton M, Todd A, Kalavacharla V. Identification of expressed resistance gene-like sequences by data mining in 454-derived transcriptomic sequences of common bean (Phaseolus vulgaris L. ). BMC Plant Biol, 2012, 12: 42[本文引用:1][JCR: 4.354]
34Liu Z, Feng S, Pand ey M K, Chen X, Culbreath A K, Varshney R K, Guo B. Identification of expressed resistance gene analogs from peanut (Arachis hypogaea L. ) expressed sequence tags. J Integr Plant Biol, 2013, 55: 453-461[本文引用:3][JCR: 3.75][CJCR: 0.7555]
35Sanz M J, Loarce Y, Fominaya A, Vossen J H, Ferrer E. Identification of RFLP and NBS/PK profiling markers for disease resistance loci in genetic maps of oats. Theor Appl Genet, 2013, 126: 203-218[本文引用:2][JCR: 3.658]
36Loarce Y, Sanz M J, Irigoyen M L, Fominaya A, Ferrer E. Mapping of STS markers obtained from oat resistance gene analog sequences. Genome, 2009, 52: 608-619[本文引用:1][JCR: 1.668]
37Tantasawat P A, Poolsawat O, Prajongjai T, Chaowiset W, Tharapreuksapong A. Association of RGA-SSCP markers with resistance to downy mildew and anthracnose in grapevines. Genet Mol Res, 2012, 11: 1799-1809[本文引用:1][JCR: 0.994]
38Palomino C, Fernández-Romero M D, Rubio J, Torres A, Moreno M T, Millán T. Integration of new CAPS and dCAPS-RGA markers into a composite chickpea genetic map and their association with disease resistance. Theor Appl Genet, 2009, 118: 671-682[本文引用:1][JCR: 3.658]
39Flor H H. The complementary genic systems in flax and flax rust. Adv Genet, 1956, 8: 29-54[本文引用:1][JCR: 4.85]
40Dilbirligi M, Gill K S. Identification and analysis of expressed resistance gene sequences in wheat. Plant Mol Biol, 2003, 53: 771-787[本文引用:3][JCR: 3.518]
41Xiao W K, Xu M L, Zhao J R, Wang F G, Li J S, Dai J R. Genome-wide isolation of resistance gene analogs in maize (Zea mays L. ). Theor Appl Genet, 2006, 113: 63-72[本文引用:1][JCR: 3.658]
42Rossi M, Araujo P G, Paulet F, Garsmeur O, Dias V M, Chen H, Van Sluys M A, D’Hont A. Genomic distribution and characterization of EST-derived resistance gene analogs (RGAs) in sugarcane. Mol Genet Genomics, 2003, 269: 406-419[本文引用:1][JCR: 2.881]
43Xiao W K, Zhao J, Fan S C, Li L, Dai J R, Xu M L. Mapping of genome-wide resistance gene analogs (RGAs) in maize (Zea mays L. ). Theor Appl Genet, 2007, 115: 501-508[本文引用:1][JCR: 3.658]
44He L, Du C, Covaleda L, Xu Z, Robinson A F, Yu J Z, Kohel R J, Zhang H B. Cloning, characterization, and evolution of the NBS-LRR-encoding resistance gene analogue family in polyploid cotton (Gossypium hirsutum L. ). Mol Plant Microbe Interact, 2004, 17: 1234-1241[本文引用:1][JCR: 4.307]
45Peñuela S, Danesh D, Young N D. Targeted isolation, sequence analysis, and physical mapping of nonTIR NBS-LRR genes in soybean. Theor Appl Genet, 2002, 104: 261-272[本文引用:1][JCR: 3.658]
46Bertioli D J, Moretzsohn M C, Madsen L H, Sand al N, Leal-Bertioli S C, Guimaraes P M, Hougaard B K, Fredslund J, Schauser L, Nielsen A M, Sato S, Tabata S, Cannon S B, Stougaard J. An analyses of synteny of Arachis with Lotus and Medicago sheds new light on the structure, stability and evolution of legume genomes. BMC Genomics, 2009, 10: 45[本文引用:1][JCR: 4.397]
47Bindler G, Plieske J, Bakaher N, Gunduz I, Ivanov N, Van der Hoeven R, Ganal M, Donini P. A high density genetic map of tobacco (Nicotiana tabacum L. ) obtained from large scale microsatellite marker development. Theor Appl Genet, 2011, 23: 219-230[本文引用:1][JCR: 3.658]
相关话题/基因 序列 烟草 结构 鉴定