,*南京农业大学大豆研究所 / 农业部大豆生物学与遗传育种重点实验室 / 国家大豆改良中心 / 作物遗传与种质创新国家重点实验室, 江苏南京 210095Cloning and analysis of candidate gene resistant to SC3 in soybean
XIANG Wen-Yang, YANG Yong-Qing, REN Qiu-Yan, JIN Tong-Tong, WANG Li-Qun, WANG Da-Gang, ZHI Hai-Jian,*Institute of Soybean, Nanjing Agricultural University / Key Laboratory of Soybean Biology and Genetic Breeding, Ministry of Agriculture / National Center for Soybean Improvement / State Key Laboratory of Crop Genetics & Germplasm Innovation, Nanjing 210095, Jiangsu, China通讯作者:
收稿日期:2019-04-2接受日期:2019-06-22网络出版日期:2019-07-15
| 基金资助: |
Received:2019-04-2Accepted:2019-06-22Online:2019-07-15
| Fund supported: |
作者简介 About authors
向文扬,E-mail:a15295572307@163.com。
摘要
关键词：
Abstract
Keywords：
PDF (3930KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
向文扬, 杨永庆, 任秋燕, 晋彤彤, 王丽群, 王大刚, 智海剑. 大豆抗SC3候选基因的克隆及分析[J]. 作物学报, 2019, 45(12): 1822-1831. doi:10.3724/SP.J.1006.2019.94054
XIANG Wen-Yang, YANG Yong-Qing, REN Qiu-Yan, JIN Tong-Tong, WANG Li-Qun, WANG Da-Gang, ZHI Hai-Jian.
大豆花叶病毒病是由大豆花叶病毒(SMV)引发的一种大豆病害, 在我国主要大豆产区都有分布, 给我国大豆的产量和品质带来严重影响[1]。目前大豆花叶病毒的防治主要从3个方面入手, 第一是控制SMV初侵染源, 带毒大豆种子是SMV初侵染来源[2], 优选大豆种子, 特别是不用前一年发病植株上收获的种子可以一定程度上防治SMV的危害; 第二是控制SMV的田间传播, 蚜虫是SMV田间传播的主要媒介[1,3], 及时发现并杀灭蚜虫能减少SMV在田间的传播; 第三是选育优良的大豆抗SMV品种, 这是最经济有效的防治方法。
国内外已经有许多关于抗SMV基因研究的报道。美国****从不同的抗SMV种质中分别发现了4个抗病位点Rsv1、Rsv3、Rsv4和Rsv5 [4,5,6,7], 并将其定位在D1b (第2染色体)、F (第13染色体)和B2 (第14染色体) 3个连锁群上, 并在随后的研究中发现这些抗病位点的多个等位位点; 我国针对不同的SMV株系从抗性种质中定位出了多个抗病基因, 分别命名为R3、RSC3Q、RSC4~RSC8、RSC10~RSC12、RSC14Q、RSC15、RSC17~RSC18、RSC20抗病位点等, 这些基因分布在D1b、F、C2 (第6染色体)和B2四个连锁群上。从物理和遗传图谱位置上来看, 美国和我国定位出来的抗病基因可能存在重叠或者处于相近区域。
目前关于R基因(resistance gene)的研究在国内外有很多的报道, 在不同的物种中已经有几十个抗性基因被成功克隆[8,9,10,11]。这些R基因介导对不同病原的抗性, 尽管其来源不同, 但是绝大多数的R基因具有高度保守的蛋白结构域, 其中典型的保守结构域包括核苷酸结合位点(nucleotide binding site, NBS)、亮氨酸拉链结构(leucine zipper, LZ)、丝氨酸/苏氨酸蛋白激酶结构域(serine-threonine kinase, STK)、toll蛋白/白细胞介素-1受体类似结构域(toll- interleukin-1 receptor, TIR)、亮氨酸富集重复结构(leucine-rich repeat, LRR)、跨膜结构域(transmembrane domain, TM)等, 其中TIR-NBS-LRR类结构基因是目前鉴定的R基因中最多的一类。
大豆抗SMV基因的分离和鉴定已经取得了一定的研究进展。Davis[12]在大豆品种L29中发现了抗SMV候选基因G533, 并通过烟草瞬时表达系统在烟草中共表达该基因与大豆花叶病毒蛋白CI, 结果表明该基因与SMV蛋白CI相互作用导致烟草植株出现了坏死症状, 说明G533可能参与对SMV的抗性。Tran等[13]在不含Rsv3的大豆品种瞬时表达Glyma.14G204700使其获得了对SMV株系G5H的系统抗性; 同时使用病毒诱导的基因沉默技术(virus-induced gene silencing, VIGS)沉默含有Rsv3大豆品种中的Glyma.14G204700使其丧失了对G5H的抗性。G533和Glyma.14G204700实际上为同一个基因, 该结果初步确认Glyma.14G204700基因与大豆对SMV的抗性有关。
大部分的植物基因含有内含子, 这些基因转录形成的前体mRNA需要经过剪接后才能形成成熟的mRNA, 前体mRNA的剪接方式可以分为组成性剪接(constitutive splicing)和选择性剪接(alternative splicing)[14,15], 其中选择性剪接是高等生物普遍存在的一种基因转录后调控机制, 能够增加基因产物的多样性, 提高植物对环境变化的适应能力, 已有研究表明, 可变剪接在植物抗病反应、抗胁迫逆境反应及生长发育过程中发挥着十分重要的作用[16,17,18]。按照剪接方式的不同, 可变剪接主要分为内含子保留(intron retention)、外显子丢失(exon skip)、可变供体(alternative donor site)、可变受体(alternative donor site)及外显子互斥(mutually exclusive exons)等[19], 内含子保留是最为常见的可变剪接形式[20]。研究植物的可变剪接可以提高对植物基因功能多样性的理解, 目前关于大豆抗SMV基因定位的研究已有很多, 但鉴定分离成功的抗病基因还很少, 候选基因可变剪接体的分析可以为大豆抗SMV基因的分离鉴定提供全新的视角。
杨永庆等[21]在PI96983品种中将大豆抗SC3基因Rscpm定位在了大豆第13染色体为345 kb的物理区间内, 本文在该研究的基础上对这个区间内的可能的抗病基因进行筛选, 搜索到了2个高度同源(99%同源性)的具有TIR-NBS-LRR典型抗病结构域的基因GmR47和GmR51, 克隆这2个基因并对其进行生物信息学分析, 同时检测了这2个基因及其剪接体在病毒诱导条件下的表达, 初步探讨了GmR47和GmR51的结构及其可能的功能, 为后续基因功能的研究奠定了基础。
1 材料与方法
1.1 材料
大豆品种南农1138-2、齐黄1号、Williams 82及PI96983均由国家大豆改良中心保存提供, SMV株系SC3由本课题组保存提供。其中齐黄1号和PI96983为抗SMV株系SC3的品种, Wiliams 82和南农1138-2是 SMV株系SC3的感病品种。1.2 材料种植及处理
本实验所使用的大豆材料均种植于温度为25℃、湿度为65℃、光照10 h/黑夜8 h的恒温培养室中; 在大豆第1对真叶展开时利用摩擦接种法[22]接种, 实验组处理为接种SMV株系SC3, 对照组为接种磷酸缓冲液(0.01 mol L-1磷酸钠, 0.01 mol L-1磷酸钾或0.05 mol L-1柠檬酸钠, pH 7)。1.3 植物总RNA的提取及第1链cDNA的反转录
提取植物总RNA的严格按照RNAprep pure Plant Kit (购于TIANGEN)说明书。利用琼脂糖凝胶电泳法对提取的RNA进行质量检测, 选取质量合格的RNA进行第1链cDNA的反转录, 按照PrimeScript 1st Strand cDNA Synthesis Kit (购于TAKARA)说明书进行反转录。1.4 基因克隆
由于GmR47及GmR51具有高度同源性, 故采用特异性较高的巢式PCR (Nested-PCR)进行扩增。以NCBI上公布的Williams 82基因组Wm82.a2.v1数据库为参考, 利用Primer premier 5软件设计引物(表1)。首先使用引物对GmR47-1和GmR51-1进行第1步PCR, 反应体系为50 μL, 包含酶25 μL 2×PrimerStar Master Mix (Takara)、上下游引物(浓度为10 μmol L-1)各1 μL、模板cDNA 1 μL、灭菌ddH2O 22 μL。PCR程序为95℃ 3 min; 98℃ 10 s, 58℃ 15 s, 72℃ 45 s, 30个循环; 72℃ 5 min。然后使用引物对GmR47-2和GmR51-2进行第2步PCR, 反应体系为50 μL, 包含酶25 μL 2×PrimerStar Master Mix (Takara)、上下游引物(浓度为10 μmol L-1)各1 μL、模板为第1步PCR产物稀释100倍后取1 μL、灭菌ddH2O 22 μL。反应程序同上。利用琼脂糖凝胶电泳法检测第2步PCR反应产物, 回收条带大小正确的产物。并将回收产物使用T4连接酶连接至T-easy载体pMD20-T, 将连接产物转化大肠杆菌感受态DH5α中过夜培养, 挑取单克隆, 将检测为阳性的单克隆送至通用生物公司测序。Table 1
表1
表1基因克隆引物及荧光定量引物
Table 1
| 引物名称 Primer name | 上游引物 Forward primer (5'-3') | 下游引物 Reverse primer (5'-3') |
|---|---|---|
| GmR47-1 | TTTCAGAATCTATTGATATTAGGG | TCTGTGAAACTATGCCTTGC |
| GmR51-1 | GGATTGGAGTGGGAGGA | GGGTTGAAGCAACGAAA |
| GmR47-2 | TGGGATTGGAGTGGGAGGAGT | CCATCAAAGTTGAATCCCAATTCAC |
| GmR51-2 | GGGAGGAGTTGAATTAGCG | CCAGTCTCCACCACTTTGTT |
| GmR47/51-qPCR | GTAGGGGTCTTTGTTGATGTTGATT | TCTGAGAAACTCCTATGTCCGCTA |
| Tubulin-qPCR | AGTATGAGGACGAGGAGGACGAT | TACGCATCACATAGCATAAGTAAGACAC |
新窗口打开|下载CSV
1.5 序列生物信息学分析
比对碱基及氨基酸序列使用BioXM软件; 通过NCBI BLASTP (1.6 荧光定量
按照荧光定量试剂盒ChamQ Universal SYBR qPCR Master Mix (购于VEZYME)说明书进行荧光定量反应, 内参基因选择Tubulin (登录号为AY907703), 每个样品和内参均设置3个平行实验, qPCR条件为95℃ 1 min; 98℃ 10 s, 60℃ 15 s, 72℃ 20 s, 40个循环。GmR47和GmR51相似度高达99%, 无法利用1对引物将其区分, 故将这2个基因的表达量统一讨论, GmR47、GmR51及Tubulin荧光定量引物见表1。分析荧光定量结果采用2-ΔΔCt法[23], 以接种磷酸缓冲液的对照组候选基因的表达量为基准品, 计算接种SC3实验组候选基因的相对表达量。1.7 可变剪接体的检测及诱导表达分析
GmR47和GmR51均存在3个含子, 分别在这3个内含子处设计了2对引物, 1对是跨内含子, 设计在2个相邻外显子上的常规引物, 另外1对是引物在内含子内的特异性引物, 内含子检测引物见表2。以PI96983、南农1138-2和齐黄1号3个大豆品种的cDNA为模板进行PCR, 检测可变剪接是否在这2个基因中存在, 此处使用的cDNA模板经过去基因组DNA处理, 排除cDNA模板中的基因组DNA的污染。按照HiScript III RT SuperMix for qPCR (+gDNA wiper)荧光反转录试剂盒(购于VEZYME)说明书进行去基因组DNA污染的反转录。可变剪接体的诱导表达分析以每个内含子处的总剪接体数为单位进行, 每个剪接体处的荧光定量引物见表2。Table 2
表2
表2可变剪接检测引物及荧光定量引物
Table 2
| 引物名称 Primer name | 引物 Primer (5'-3') |
|---|---|
| Intron1-Uni | F: TTCGGAGGCGATTGAACAT |
| R: TGCTCAATTTGATCAATTAATTAAACTCA | |
| Intron2-Uni | F: AAATTGTTCGACAGGAATCACC |
| R: TGGCAGAAACTCGACGTAATG | |
| Intron3-Uni | F: CTCTTGAATCCACCAAAGCAGG |
| R: AGTTCATATGATTAGAAGAACATGTAGAG | |
| Intron1-Spe | F: ATCTACTAGTCATATGGATTTGCTCAATTTGATCAATTAATTAAACTCA |
| R: CGGTACCCGGGGATCCGATTTTCGGAGGCGATTGAACAT | |
| Intron2-Spe | F: CGGTACCCGGGGATCCGATTAAATTGTTCGACAGGAATCACC |
| R: ATCTACTAGTCATATGGATTTGGCAGAAACTCGACGTAATG | |
| Intron3-Spe | F: CGGTACCCGGGGATCCGATTCTCTTGAATCCACCAAAGCAGG |
| R: ATCTACTAGTCATATGGATTAGTTCATATGATTAGAAGAACATGTAGAG | |
| Intron1-qPCR | F: GAAATGGAGGAACGCACTGC |
| R: AAGATCGACAAGAAGATATTACCTGAA | |
| Intron2-qPCR | F: GGCAAACGCAGTAGGTTATGG |
| R: ATACCATCAATAGAATGGCAGAAACT | |
| Intron3-qPCR | F: ATGTAGGAGAATGCTACTGAAACAGG |
| R: AGAATATCTCTTATTTCCAGCCTCAT |
新窗口打开|下载CSV
1.8 亚细胞定位
使用在引物5'端含有GateWay元件的引物扩增GmR47和GmR51完整的开放阅读框, 其中不包含终止子。利用GateWay系统构建载体, 经过BP反应及LR反应, 将目的基因构入pGWB5载体获得目的基因与GFP融合的载体pGWB5-GmR47-GFP和pGWB5-GmR51-GFP。使用农杆菌注射的方法侵染烟草[24], 侵染烟草48 h后, 利用Leica TCS SP2激光共聚焦荧光显微镜观察注射烟草细胞中GFP荧光信号的分布。2 结果与分析
2.1 候选基因的克隆
通过已公布的Williams 82基因组数据, 分别获得GmR47和GmR51的参考序列, 比对参考序列可知2个基因具有高度同源性。克隆大豆品种PI96983中的抗病候选基因GmR47和GmR51, 设计2对特异性引物分别对这2个基因进行巢式PCR扩增, 以大豆品种PI96983的cDNA为模板扩增获得大小为3500 bp左右的目标片段。测序结果显示, 大豆品种PI96983中GmR47和GmR51基因均含有一个大小为3288 bp的开放阅读框(ORF), 均编码1095个氨基酸残基。GmR47编码合成的蛋白质分子量为125.34 kD, 等电点为6.69; GmR51编码合成的蛋白质分子量为125.4 kD, 等电点为6.64。2.2 候选基因生物信息学分析
GmR47和GmR51编码的蛋白产物含有相同的蛋白结构域, 从N端到C端分别为TIR、NBS及LRR结构域, 是典型的TIR-NBS-LRR家族基因(图1), 预测的蛋白功能为烟草花叶病毒(TMV)抗性N蛋白。GmR47、GmR51在抗病品种PI96983及感病品种南农1138-2之间的氨基酸序列差异(图2)显示, GmR47在抗感品种之间存在8个SNP的差异, 导致5个氨基酸突变, 突变位点在TIR、NBS和LRR结构域中均有分布; GmR51在抗感品种之间存在8个SNP差异, 导致3个氨基酸突变, 突变位点分布在NBS和LRR结构域中。GmR47和GmR51在不同物种间的同源比对结果(图3)显示, GmR47和GmR51基因在野生大豆中的亲缘性最近, 其次是苜蓿。图1
新窗口打开|下载原图ZIP|生成PPT图1GmR47和GmR51编码蛋白质的保守结构域预测
A: GmR47编码蛋白质的保守结构域预测。B: GmR51编码蛋白质的保守结构域预测。
Fig. 1Conserved domain prediction of protein encoded by GmR47 and GmR51
A: conservative domain prediction of GmR47-encoded proteins. B: conservative domain prediction of GmR51-encoded proteins.
图2
新窗口打开|下载原图ZIP|生成PPT图2GmR47和GmR51编码的核苷酸序列在抗病品种(PI96983)及感病品种(南农1138-2)之间的比对结果
A: GmR47编码氨基酸序列的比对结果。B: GmR51编码氨基酸序列比对的结果。
Fig. 2Alignment result of the Amino acid sequence encoded by GmR47 and GmR51 in resistant (PI96983) and susceptible (Nannong 1138-2) cultivar
A: alignment result of amino acid sequence encoded by GmR47; B: alignment result of amino acid sequence encoded by GmR51.
图3
新窗口打开|下载原图ZIP|生成PPT图3GmR47和GmR51基因在不同物种间的同源进化树
Fig. 3Phylogenetic trees of GmR47 and GmR51 in different species
2.3 GmR47和GmR51基因编码蛋白的二级结构和三级结构预测
图4-A表明, 在GmR47二级结构中α螺旋结构占45.84%, 延伸链占15.07%, β转角结构占5.57%, 不规则卷曲占33.52%; GmR51编码的蛋白产物二级结构(图4-B)中α螺旋结构占45.02%, 延伸链占16.07%, β转角结构占5.57%, 不规则卷曲占33.33%。抗感病品种中GmR47和GmR51编码蛋白产物预测三级结构见图5, 显示GmR47和GmR51基因编码蛋白产物三级结构在抗感品种中存在差异。图4
新窗口打开|下载原图ZIP|生成PPT图4GmR47 (A)和GmR51 (B)基因编码的氨基酸序列二级结构预测结果
Fig. 4Prediction of secondary structure of amino acid sequence encoded by GmR47 (A) and GmR51 (B)
图5
新窗口打开|下载原图ZIP|生成PPT图5抗感品种中GmR47和GmR51基因编码合成的蛋白质三维结构预测结果
A: PI96983 GmR47; B: PI96983 GmR51; C: 南农1138-2 GmR47; D: 南农1138-2 GmR51。
Fig. 5Predicted three-dimensional structure of protein encoded by gene GmR47 and GmR 51
A: PI96983 GmR47; B: PI96983 GmR51; C: Nannong 1138-2 GmR47; D: Nannong 1138-2 GmR51.
2.4 候选基因在抗感品种之间的表达量差异
利用SC3接种抗感品种诱导候选基因的表达, 检测抗感品种之间表达量的差异。实验组采用SMV株系SC3接种, 对照组用磷酸缓冲液处理。分别在处理后的0、1、2、4、8、12、24、48、72 h以及7 d采集对照组和实验组的处理叶片。荧光定量的结果显示(图6), 候选基因在抗病品种中的表达量变化趋势是先升后降再升, 分别在2 h、7 d达到峰值, 最高值约为未接种时表达量的2倍; 在24 h达到最低值, 此时表达量与未接种时持平。候选基因在感病品种表达量的变化不大, 在12 h时取得最低值, 约为未接种时表达量的0.5倍, 在24 h达到峰值, 约为未接种时表达量的1.6倍; 由此看出, 感病品种中候选基因的表达量相对平稳, 而在抗病品种中表达量则有较大的起伏变化, 在抗病品种中的表达量要高于感病品种中。说明候选基因GmR47和GmR51能够积极的响应SMV的诱导提高表达量。图6
新窗口打开|下载原图ZIP|生成PPT图6PI96983和南农1138-2接种SMV后GmR47及GmR51基因的总表达量随时间的变化
Fig. 6Total expression of GmR47 and GmR51 genes at different time in PI96983 and Nannong 1138-2 after SMV inoculation
2.5 可变剪接体的分析
GmR47和GmR51基因均存在3个内含子, 为检测其内含子处的可变剪接现象, 在每个内含子处设计了1对常规引物和1对特异性引物, 以齐黄1号、PI96983、南农1138-2和Williams 82的去基因组污染的cDNA为模板进行PCR扩增。扩增的结果显示(图7), 在这3个内含子处, 跨内含子设计在相邻外显子上的常规引物均扩增出2条或2条以上的条带, 且设计在内含子上的特异性引物同样扩增出了条带。说明GmR47和GmR51基因在这3个内含子处均存在可变剪接的现象, 并且可变剪接现象在4个大豆品种中均存在。将常规引物PCR扩增产物回收后测序显示(图8), 3个内含子处均存在内含子保留可变剪接体, 由于GmR47和GmR51基因的高度同源性, 未能检测出每个基因的不同可变剪接体。图7
新窗口打开|下载原图ZIP|生成PPT图7可变剪接体的凝胶电泳检测结果
上排条带为常规引物扩增的结果, 下排条带为内含子特异性引物的扩增结果, 相邻的4个条带为4个品种间的扩增结果。A: 齐黄1号; B: PI96983; C: 南农1138-2; D: Williams 82; M: marker。
Fig. 7Gel electrophoresis detection results of alternative splicing
The upper bands are the result of amplification of theconventional primer, the lower bands are the amplification result of the intron-specific primer, adjacent four bands are amplification from four varieties. A: Qihuang 1; B: PI96983; C: Nannong 1138-2; D: Williams 82; M: marker.
图8
新窗口打开|下载原图ZIP|生成PPT图8可变剪接测序结果
A: GmR47基因可变剪接体示意图; B: IN1序列; C: IN2序列; D: IN3序列。
Fig. 8Alternative splicing sequencing results
A: schematic diagram of alternative splicing of GmR47 gene; B: IN1 sequence; C: IN2 sequence; D: IN3 sequence.
由于GmR47和GmR51基因具有高度同源性, 且基因组较大, 无法单独对某个剪接体进行分析, 为研究不同剪接体对SMV的响应情况, 本文以每个内含子处存在的总剪接体为单位进行表达量分析, 即在每个内含子内设计特异性引物进行荧光定量, 将3个内含子处存在的总剪接体分别命名为IN1、IN2和IN3。可变剪接诱导表达所采用的实验材料为PI9683和南农1138-2, 实验组处理为接种SC3病毒, 对照组为接种磷酸缓冲液, 抗病品种及感病品种均设置实验组和对照组。在处理后的0、1、2、4、8、12、24、48、72 h以及7 d内采集处理叶片, 并进行荧光定量分析。荧光定量结果(图9)显示, IN1的表达量(图9-A)在抗病品种中的趋势是先升再降再升, 在感病品种中的趋势是先降再升再降, 其中2、4、8、72 h在抗感品种之间存在显著差异; IN2的表达量(图9-B)在抗感病品种中的变化趋势是先升再降, 抗病品种中候选基因在24 h表达量达到峰值, 在感病品种中在12 h表达量达到峰值; IN3的表达量(图9-C) 在抗病品种中的变化趋势是先升后维持较高表达水平, 在感病品种中表达量变化不大。以上结果表明, IN1、IN2和IN3的的表达量都能响应病毒的诱导, 并且抗病品种中的表达量要高于感病品种。不同内含子间可变剪接体的表达量存在差异, IN1表达量随病毒侵染时间的变化较为明显, 而IN2和IN3的表达量则相对平稳。总的来说, GmR47和GmR51在3个内含子处均含有可变剪接现象, 而且能够响应病毒的诱导表达, 说明这些可变剪接体可能参与大豆对SMV的抗性。
图9
新窗口打开|下载原图ZIP|生成PPT图9PI96983和南农1138-2接种SMV后GmR47及GmR51基因在不同内含子处的剪切体的总表达量随时间的变化
A: 内含子1处的剪接体的总表达量; B: 内含子2处的剪接体的总表达量; C: 内含子3处的剪接体的总表达量。
Fig. 9Total expression of GmR47 and GmR51 genes at different time in PI96983 and Nannong 1138-2 after SMV inoculation
A: total expression of alternative splicing in Intron 1; B: total expression of alternative splicing in Intron 2; C: total expression of alternative splicing in Intron 3.
2.6 候选基因的亚细胞定位
利用Cell-PLoc (图10
新窗口打开|下载原图ZIP|生成PPT图10候选基因与GFP融合蛋白在烟草细胞中的亚细胞定位结果
Fig. 10Subcellular localization results of candidate genes and GFP fusion proteins in tobacco cells
3 讨论
本文发现候选基因能响应SMV的侵染而增加表达量, 说明候选基因能够参与对SMV的抗性反应, 但其表达量与CP含量的变化趋势并不一致, 推测其原因可能是植物抗病过程由多个基因共同作用, 单个基因的高表达并不一定能使抗性增强而直接降低病毒量。美国以不同抗性的种质定位出Rsv1、Rsv3、Rsv4、Rsv5等抗病基因, 目前已经从PI96983中初步鉴定了具有NBS-LRR结构的抗SMV基因3gG2[25], 该基因有可能为Rsv1, 从遗传距离和物理距离上来看, Rsv1和我国定位的Rsc3距离很近, 可能互为等位基因或处于相近区域; 韩国****鉴定L29中的Glyma.14G204700为抗病基因, 证明其有可能为Rsv3[8], 该基因与我国大白麻含有的Rsc4基因位点相同, 初步说明L29与大白麻可能含有相同抗病基因; 在我国, 目前确定为抗SMV基因的候选基因还较少, 黄赛花等[26]利用VIGS沉默齐黄1号中Rsc3Q候选基因GmZ15使抗SMV植株的抗病性发生改变, 初步证明GmZ15为抗病基因参与齐黄1号对SMV的抗性。本文根据生物信息学和表达量分析的结果, 初步推测GmR47和GmR51基因为抗SMV基因, 至于候选基因是否为SMV抗性基因, 还需要VIGS、转基因以及CRIPSR-CAS9基因靶向编辑等实验后续证明。
可变剪接可以使一个基因生成多种转录产物, 极大地提高基因产物的多样性, 是一种重要的转录后调控机制, 在植物的抗病、抗逆、生长发育中发挥重要作用。目前关于大豆抗病基因的研究主要是以基因转录生成的前体RNA组成型切割体为基础, 而较少基于可变剪接体, 但在其他作物上已有关于可变剪接体参与植物抗性反应的报道, 如水稻抗稻瘟病基因RAG5转录后存在2种可变剪接体, 其中RAG5-A介导对稻瘟病的抗性, 而另一种剪接形式RAG5-B则没有表现出对稻瘟病的抗性[27], 结果说明基因转录产生的不同剪接体可能具有截然不同的功能。本文在大豆抗病候选基因GmR47及GmR51中检测到可变剪接现象的存在, 发现这些可变剪接体能够响应SMV的侵染增加表达量, 并且在抗病品种和感病品种之间的表达量存在较为明显的差异, 说明这2个基因的可变剪接体可能参与抗病品种对SMV的抗性。GmR47和GmR51基因剪接体的发现, 为其基因功能的鉴定提供了新的视角。但由于GmR47和GmR51基因具有较高同源性, 且基因各个内含子处均存在可变剪接, 候选基因每个可变剪接体的存在形式还需要进一步的实验证明。
4 结论
从大豆品种PI96983中扩增了2个TIR-NBS- LRR结构基因, 分别命名为GmR47和GmR51。它们高度同源, 均含有一个3288 bp的开放阅读框, 编码1095个氨基酸。GmR47在抗病品种(PI96983)和感病品种(南农1138-2)中存在8个SNP差异, 最终导致5个氨基酸发生突变; GmR51在抗病品种和感病品种中也存在8个SNP差异, 最终导致3个氨基酸突变。GmR47和GmR51亲缘最近的物种是野生大豆。GmR47和GmR51基因编码的蛋白产物主要分布在细胞膜和核膜上。GmR47和GmR51基因在抗感病品种中均存在可变剪接现象, 2个基因及其可变剪接体均能够响应病毒的诱导提高表达量, 且在抗病品种中的表达量高于感病品种, 其中可变剪接体IN1和IN2的表达量变化较为明显, IN3的表达量则相对平稳。初步推测候选基因GmR47和GmR51为抗SMV基因, 其可变剪接体可能参与对SMV的抗性。参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]
DOI:10.1007/BF00222952URL [本文引用: 1]
[本文引用: 1]
DOI:10.1007/s001220100567URL [本文引用: 1]
In soybean, genes controlling resistance to numerous diseases have been shown to cluster to regions on several chromosomes. One such vital chromosomal region is on the soybean molecular linkage group (MLG) F flanked by the RFLP markers K644 and B212. Here, genes controlling resistance to bacterial blight, Phytophthora root rot, and several viral diseases, as well as QTLs conditioning resistance to corn earworm, root knot nematode, and white mold have been mapped. We have previously identified two classes (b and j) of disease resistance-related nucleotide binding site (NBS) sequences that localize to this cluster of genes. Using both cDNA and genomic analyses, we have studied one multi-gene family of sequences representing the previously reported class j NBS of soybean. This class of NBS resembles the RPS2-like NBS sequences. RPS2 and similar resistance genes are referred to as non-TIR because they do not encode motifs homologous to the Toll-Interleukin-1 region (TIR). By designing PCR primers that specifically target these non-TIR-NBS encoding sequences, we have amplified at least six class j sequence members from soybean. In addition, we have conducted genomic and cDNA library screenings to identify additional class j members. In all, we have characterized 12 class j NBS sequence members. These members have been mapped within a 2-cM region of the soybean F linkage group. We have also identified homoeologous chromosomal regions on linkage groups A2 and E that contain class j NBS sequences. A BLAST search of the GenBank database has shown that non-TIR NBS sequences are present across the legume family. We have compared these non-TIR sequences from other legumes with the soybean clones to assess the level of diversity within this class of disease resistance-related sequences.
DOI:10.1007/s11032-016-0595-3URL [本文引用: 1]
DOI:10.1126/science.8091210URL [本文引用: 2]
DOI:10.1126/science.7638602URL [本文引用: 1]
DOI:10.1105/tpc.9.6.879URL [本文引用: 1]
DOI:10.1073/pnas.95.4.1663URL [本文引用: 1]
[本文引用: 1]
DOI:10.1016/j.virol.2017.10.014URL [本文引用: 1]
URL [本文引用: 1]
选择性剪接是细胞生物体产生复杂蛋白的一个重要机制。迄今为止，对植物选择性剪接基因的研究远不如动物。尽管如此，仍在选择性剪接的调控机制、剪接特异性和影响因素及剪接功能等方面取得了一定的进展。本文就这方面的研究进展作一简要综述。
URL [本文引用: 1]
选择性剪接是细胞生物体产生复杂蛋白的一个重要机制。迄今为止，对植物选择性剪接基因的研究远不如动物。尽管如此，仍在选择性剪接的调控机制、剪接特异性和影响因素及剪接功能等方面取得了一定的进展。本文就这方面的研究进展作一简要综述。
DOI:10.1016/0092-8674(78)90078-8URL [本文引用: 1]
[本文引用: 1]
[本文引用: 1]
DOI:10.1007/s12374-009-9018-xURL [本文引用: 1]
The FCA gene, which is a component of the autonomous pathway that regulates flowering time, is an important example of how alternative processing can control plant development. We have previously characterized the FCA homolog, OsFCA, from a japonica-type rice cultivar and demonstrated that the polyadenylation site within intron 3, which can generate non-functional FCA-β, was conserved in rice. In this study, we detected five alternatively processed variants of OsFCA pre-mRNA, four of which were equivalents of FCA-α, -β, -γ, and -δ, in japonica-type Korean rice cultivars. The fifth transcript, referred to as OsFCA-ɛ, was similar to OsFCA-γ, except a part of the OsFCA intron 16 was retained. Unlike the FCA-γ protein, the OsFCA-γ protein contains a glycine-rich region at its N-terminus. We detected the OsFCA transcripts missing the region encoding the glycine-rich domain in the indica-type rice, but not in the japonica-type rice. We also found that the OsFCA-δ and OsFCA-ɛ transcripts were expressed in almost all of the different tissue types examined. Taken together, these results indicate that the alternative processing of the OsFCA transcript is more complex than its Arabidopsis counterpart.
DOI:10.1038/ng0102-13URL [本文引用: 1]
DOI:10.1104/pp.121.2.411URL [本文引用: 1]
DOI:10.1007/s00122-013-2092-yURL [本文引用: 1]
Soybean mosaic virus (SMV) is one of the most broadly distributed soybean (Glycine max (L.) Merr.) diseases and causes severe yield loss and seed quality deficiency. Multiple studies have proved that a single dominant gene can confer resistance to several SMV strains. Plant introduction (PI) 96983 has been reported to contain SMV resistance genes (e.g., Rsv1 and Rsc14) on chromosome 13. The objective of this study was to delineate the genetics of resistance to SMV in PI 96983 and determine whether one gene can control resistance to more than one Chinese SMV strain. In this study, PI 96983 was identified as resistant and Nannong 1138-2 was identified as susceptible to four SMV strains SC3, SC6, SC7, and SC17. Genetic maps based on 783 F2 individuals from the cross of PI 96983?×?Nannong 1138-2 showed that the gene(s) conferring resistance to SC3, SC6, and SC17 were between SSR markers BARCSOYSSR_13_1114 and BARCSOYSSR_13_1136, whereas SC7 was between markers BARCSOYSSR_13_1140 and BARCSOYSSR_13_1185. The physical map based on 58 recombinant lines confirmed these results. The resistance gene for SC7 was positioned between BARCSOYSSR_13_1140 and BARCSOYSSR_13_1155, while the resistance gene(s) for SC3, SC6, and SC17 were between BARCSOYSSR_13_1128 and BARCSOYSSR_13_1136. We concluded that, there were two dominant resistance genes flanking Rsv1 or one of them at the reported genomic location of Rsv1. One of them (designated as “Rsc-pm”) conditions resistance for SC3, SC6, and SC17 and another (designated as “Rsc-ps”) confers resistance for SC7. The two tightly linked genes identified in this study would be helpful to cloning of resistance genes and breeding of multiple resistances soybean cultivars to SMV through marker-assisted selection (MAS).
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
[本文引用: 1]
DOI:10.1534/genetics.166.1.493URL [本文引用: 1]
[本文引用: 1]
[本文引用: 1]
DOI:10.1105/tpc.112.107201URL [本文引用: 1]
Resistance (R) proteins recognize pathogen avirulence (Avr) proteins by direct or indirect binding and are multidomain proteins generally carrying a nucleotide binding (NB) and a leucine-rich repeat (LRR) domain. Two NB-LRR protein-coding genes from rice (Oryza sativa), RGA4 and RGA5, were found to be required for the recognition of the Magnaporthe oryzae effector AVR1-CO39. RGA4 and RGA5 also mediate recognition of the unrelated M. oryzae effector AVR-Pia, indicating that the corresponding R proteins possess dual recognition specificity. For RGA5, two alternative transcripts, RGA5-A and RGA5-B, were identified. Genetic analysis showed that only RGA5-A confers resistance, while RGA5-B is inactive. Yeast two-hybrid, coimmunoprecipitation, and fluorescence resonance energy transfer-fluorescence lifetime imaging experiments revealed direct binding of AVR-Pia and AVR1-CO39 to RGA5-A, providing evidence for the recognition of multiple Avr proteins by direct binding to a single R protein. Direct binding seems to be required for resistance as an inactive AVR-Pia allele did not bind RGA5-A. A small Avr interaction domain with homology to the Avr recognition domain in the rice R protein Pik-1 was identified in the C terminus of RGA5-A. This reveals a mode of Avr protein recognition through direct binding to a novel, non-LRR interaction domain.
