以陆地棉
Using upland cotton
病毒介导的基因沉默(virus-induced gene silencing, VIGS)是植物体内普遍存在的一种遗传免疫机制, 属于转录后基因沉默(post-transcriptional gene silencing, PTGS)[ 1]。当携带含有内源基因片段的病毒载体侵染寄主植物后, 能够激活植物自身的免疫系统; 在识别并降解病毒RNA的同时也会产生含内源目的基因的microRNA, 这些microRNA与靶基因的mRNA结合, 之后被Dicer酶降解, 从而使目的基因表达水平下降或功能丧失[ 2]。其中由烟草脆裂病毒(tobacco rattle virus, TRV)诱导的基因沉默(TRV- VIGS), 以其沉默效率高、持续时间长, 寄主植物病毒症状轻, 不会掩盖沉默表型, 且在各种组织均可产生基因沉默等优点, 成为目前应用最为广泛的一类基因沉默体系。TRV-VIGS已在茄属的番茄、烟草、辣椒[ 3, 4, 5]、拟南芥[ 6]、麻风树[ 7]、矮牵牛[ 8]、罂粟[ 9, 10, 11]等植物上成功应用。
植物 CLA 1基因( cloroplastos alterados 1 gene)编码1-deoxyxylulose 5-phosphate synthase, 参与叶绿体发育过程, 进化中高度保守, CLA1基因突变体 cla1-1有明显的白化表型[ 12], 是易于识别的标记性状。棉花是世界性重要的经济作物。近年来, 已有TRV-VIGS体系构建与应用的相关报道。以雷蒙德氏棉 CLA1基因序列设计目标沉默片段, Gao等[ 13]建立了基于农杆菌介导的棉花TRV-VIGS体系, 并在不同遗传背景的6个棉花品种材料中成功应用, 为棉花重要基因功能的鉴定奠定了基础。进一步通过TRV病毒抑制棉花抗黄萎病相关基因 GhNDR1和 GhMKK2表达, 使棉花受体更易感病, 快速鉴定了基因功能[ 14]。Qu等[ 15]发现TRV-VIGS不仅可用于沉默营养器官中基因, 还可用于沉默纤维等组织和生殖器官如花蕾中的基因。然而, TRV-VIGS能否广泛应用于棉花各品种, 即该体系是否受寄主基因型的影响还不清楚, 未见相关研究报道。近年来, 棉花基因组研究取得快速进展, 大量的棉花表达序列信息已在公共数据库释放(http://www.ncbi.nlm.nih.gov/dbEST/), 与四倍体棉种D基因组供体亲缘关系最近的雷蒙德氏棉, 已完成全基因组测序及序列分析[ 15, 16], 大量棉花重要基因的功能需进一步鉴定。建立基于目标性状的不同棉花受体基因瞬时干扰表达体系, 将加速棉花功能基因组学研究进程, 快速获得目标性状的候选功能基因。
本研究以拟南芥 AtCLA1基因为探针, 搜索雷蒙德氏棉( G. raimondiiUlbrich)全基因组序列[ 16], 获得 CLA1基因序列, 进一步设计引物, 通过同源克隆法从陆地棉遗传标准系TM-1中克隆出 CLA 1基因, 命名为 GhCLA1。构建了VIGS-GhCLA1载体, 以海7124为受体材料, 建立了棉花的TRV-VIGS体系。选择不同来源及不同生态区的34份棉花材料, 进行沉默体系重复试验。通过TRV-VIGS技术抑制受黄萎病诱导表达的 GhMAPKKK基因的表达, 发现海7124对黄萎病的感病性显著增强。与通过棉花遗传转化来研究基因功能相比, TRV-VIGS具有无基因型选择性、快速(当代可观察表型)和高通量等特点, 将在棉花基因功能研究中扮演重要角色。
34份棉花品种(系), 包括33份陆地棉材料和1份海岛棉抗黄萎病品种海7124 (表1)。33份陆地棉材料包括: 黄河流域棉区11个, 长江流域棉区8个, 西北内陆棉区9个, 北部特早熟棉区2个; 感黄萎病品系军棉1号, 抗黄萎病品系60182, 以及陆地棉遗传标准系TM-1。TRV-VIGS体系的构建以及 GhMAPKKK的功能研究均以海7124为研究材料。同时诱导34份材料 GhCLA1基因沉默, 观察叶片的白化情况, 以研究TRV病毒在不同遗传背景材料中的适用性。
pTRV1和pTRV2 VIGS载体由清华大学刘玉乐教授惠赠。农杆菌GV3101由本实验室保存。MES (Sigma原装)及卡那霉素、利福平、乙酰丁香酮等主要国产生化试剂, 均购自南京奥斌生物科技有限公司。
参考Velásquez等的制备方法[ 17]。以pTRV1和含有目的基因片段的pTRV2分别转化农杆菌GV3101。挑取新鲜培养的含pTRV1和目的基因片段的pTRV2的农杆菌GV3101单菌落, 分别接种到3 mL LB液体培养基(Kan, 50 μg mL-1; Rif, 50 μg mL-1)中, 28℃ 180转 min-1培养24 h; 以1%含量接入50 mL LB液体培养基 (Kan, 50 μg mL-1; Rif, 50 μg mL-1)中, 28℃ 180转 min-1培养12 h, 至菌液OD600为1.0左右; 1700× g离心5 min收集菌体细胞, 以适当体积的重悬液(10 mmol L-1 MgCl2, 10 mmol L-1 MES以及200 μmol L-1乙酰丁香酮)重悬至终浓度为0.8 (OD600); 将重悬液于室温下静置3 h以上, 含有pTRV1和目的基因片段的pTRV2的重悬液按体积比1∶1混匀, 用于注射棉花子叶。
以人工气候室培养到出苗后8 d, 子叶完全展平的棉花幼苗为试验接菌材料。先用注射器针头轻轻刺破子叶背面造成微伤口, 用去针头的注射器从伤口处注入1∶1混匀的含有pTRV1和含目的基因片段的pTRV2重悬液。于昼/夜温度23/22℃, 光/暗周期16 h/8 h的人工气候室内培养侵染后的棉花植株。在相同的环境下培养未经侵染处理和注射空载体(1∶1混和pTRV1和pTRV2)的对照植株, 2~5周观察不同处理棉花的表型, 并检测目的基因的表达情况。对每种材料处理10个单株。
2.1.1 棉花 CLA1基因的克隆及pTRV2-GhCLA1载体构建 以NCBI中释放的拟南芥 CLA1基因序列(GenBank登录号为AT4g15560)为探针, 搜索雷蒙德氏棉( G. raimondiiUlbrich)全基因组序列[ 15], 获得其序列信息(Gorai.011G184200.2), 进一步设计引物, 从陆地棉遗传标准系TM-1中克隆 GhCLA1基因 (GenBank登录号为KJ123647)。其ORF全长2163 bp, 编码720个氨基酸。序列相似性分析表明, 该基因与拟南芥 CLA1基因的氨基酸相似性为83.0%。而来自雷蒙德氏棉和陆地棉TM-1中的 CLA1同源基因其核苷酸和氨基酸序列相似性均达99%以上, 非常保守。进一步设计PCR引物, 并在引物的5'端分别加入限制性内切酶 EcoR I和 Kpn I的酶切位点和保护碱基(表2), 扩增位于棉花 GhCLA1基因ORF中部994~ 1414 bp基因片段, 目标扩增片段长度为421 bp。经测序验证后, 将该基因片段插入沉默载体pTRV2, 获得重组载体pTRV2-GhCLA1。
2.1.2 棉花TRV-VIGS体系的建立 以未注射农杆菌及注射含pTRV1/pTRV2空载体的农杆菌混合液的植株为对照, 用含有pTRV1/pTRV2- GhCLA1的农杆菌混合液侵染海岛棉海7124出苗8 d的子叶, 两周后观察真叶表型。结果表明注射pTRV1/pTRV2空载体的幼苗叶片未呈现白化, 表型与未经侵染的对照无异。而注射pTRV1/pTRV2-GhCLA1植株, 其真叶表现为叶脉及边缘部白化, 逐渐扩散到整个叶片; 从第2片真叶起, 白化现象明显, 新生叶几乎完全白化(图1)。
为检测TRV在棉花植株中的扩散情况, 提取注射pTRV1和pTRV2空载体3周后的棉株根、茎、叶组织RNA, 利用TRV-RNA2特异引物将病毒RNA2反转为cDNA, 并设计特异性引物扩增病毒RNA2序列片段, 以检测病毒在棉花不同组织器官的扩散。RT-PCR分析表明子叶接种病毒后, 病毒RNA在棉花根、茎、叶中均有积累(图2), TRV病毒侵染棉花后可扩散到棉花多种组织器官中, 为利用该体系高通量验证根、茎、叶等组织器官优势表达目标基因功能提供了依据。
对pTRV1/pTRV2-GhCLA1处理3周后, 白化表型明显的植株和仅空载体侵染处理的对照植株, 提取叶片总RNA, 通过半定量RT-PCR的方法检测 GhCLA1基因的表达水平。结果表明, 与对照相比, 侵染pTRV1/pTRV2-GhCLA1的植株 GhCLA1基因几乎不表达, 沉默效果明显(图3)。
2.1.3 棉花TRV-VIGS体系利用的广谱性 选择不同来源, 不同生态区的34份棉花材料, 进行沉默体系重复。以注射pTRV1/pTRV2空载体的植株为对照, 用含有pTRV1/pTRV2-GhCLA1的农杆菌混合液注射34份供试棉种子叶完全展平的幼苗, 诱导棉花 CLA1基因沉默。侵染3周后统一调查第2片真叶白化表型(图4)。发现在供试的不同棉种材料中, 其 CLA1抑制程度有差异, 如新陆早1号、朝阳棉1号、TM-1以及海7124等整个叶片白化, 中棉所12、晋棉25、豫棉20及华抗棉1号等叶脉明显白化, 叶肉轻度失绿。说明TRV病毒可应用于所有供试材料的 CLA1基因沉默, 但不同的棉花品种对TRV病毒的敏感性存在差异。利用TRV-VIGS技术诱导目的基因沉默受体遗传背景影响较小, 可广泛用于来源不同棉花资源材料中目标基因的功能鉴定。
在MAPK信号通路中包含至少3个顺序的激活蛋白激酶, 依次为MAPKKKs、MAPKKs及MAPKs[ 21]。MAPKKKs参与了病原物侵染植物后引起的过敏性反应[ 22]。 GhMAPKKK(GenBank登录号为FJ966894)是棉花MAPKKKs家族中的一个成员, 本研究表明该基因受黄萎病菌诱导表达。以生长到二叶一心期的海7124幼苗为受体, 接种黄萎病菌株Vd8, 分别检测接种0、24、48、96和144 h后根组织中 GhMAPKKK基因的表达水平。RT-PCR分析表明, 随着接种时间增加, GhMAPKKK基因表达明显增强, 接种后96 h表达量达到高峰, 接种后144 h表达量开始降低(图5)。说明 GhMAPKKK参与棉花抗黄萎病信号传导途径。
利用TRV-VIGS体系进一步验证了 GhMAPKKK在黄萎病抗性中的功能。选择 GhMAPKKK基因靠近3'区1519~2077 bp碱基位置, 基因专化性强, 扩增出556 bp的目标片段, 克隆到pTRV2载体中获得重组质粒pTRV2-GhMAPKKK。用含有pTRV1/pTRV2- GhMAPKKK以及pTRV1/pTRV2-GhCLA1和pTRV1/ pTRV2空载体的农杆菌重悬液, 分别侵染海岛棉中抗黄萎病栽培品种海7124出苗8 d、子叶展平的幼苗。3周后, 注射TRV2- GhCLA1的植株出现明显白化。提取pTRV2-GhMAPKKK处理以及注射空载体的对照植株叶片总RNA, 通过半定量RT-PCR的方法检测表明, 与对照植株相比, 侵染pTRV2- GhMAPKKK的植株 GhMAPKKK基因的表达水平显著降低, 沉默效果明显(图6)。
对pTRV2-GhMAPKKK处理以及注射空载体的对照海7124植株, 接种大丽轮枝菌菌株Vd8。一个月后观察黄萎病的发病情况。通过3次生物学重复分析, 注射空载体的对照海7124植株, 其平均病情指数为21.9, 而 GhMAPKKK基因沉默后的海7124植株感病性增强, 其平均病情指数达41.25。表明 GhMAPKKK基因参与了黄萎病诱导的过敏反应(图7)。
在病毒诱导的基因沉默中, 最早使用的沉默载体是以烟草花叶病毒(tobacco mosaic virus, TMV)和马铃薯病毒(potato virus X, PVX)为基础的, PVX比较稳定但是宿主范围小, TMV宿主范围广但不稳定。TMV和PVX系统的缺陷还在于二者均能表现出病毒本身产生的症状, 干扰那些能产生轻微表型的基因, 使结果不易区别[ 23]。同时这2种病毒在宿主植物的生长点和分生组织中不易扩散, 影响目标基因沉默的整体效果。而以烟草脆裂病毒(tobacco rattle virus, TRV)为基础的病毒载体具有明显的优势, TRV具有更广泛的宿主范围, 可以在包括植株生长点和分生组织的所有组织器官内扩散, 病毒侵染植株后的症状与其他病毒相比更加轻微。本研究以完全展平子叶的棉花为受体, 成功建立棉花病毒诱导基因沉默体系, 该体系具有稳定性和广谱性, 为快速高通量研究棉花基因功能奠定了基础。
构建pTRV-VIGS载体时需要考虑以下3个方面: (1)插入载体的目的基因片段的大小。要求适中, 以300~500 bp为宜。TRV病毒载体插入基因序列长度上限为1.5 kb, 若超过1.5 kb病毒不易在植物体内扩散或者丢失目的片段, 影响沉默效果[ 24]。(2)选择合适的目的基因区域插入载体。如果是单拷贝基因, 其ORF中的任何区域都可以用来作为插入片段, 若是多拷贝或家族基因, 则要根据研究目的选择基因区域。如果专化沉默某一家族中的一个基因, 应选择非保守区域或基因非翻译区(UTR)来区别高度保守的基因, 如果两基因的同源性达到85%以上就可能被同时沉默。如果同时沉默家族基因中功能相似的基因, 克服功能冗余, 则可选择高度保守区[ 25, 26]。(3)环境因素影响VIGS效果。温度是对沉默效果影响最大的因素之一, 适宜的温度有利于提高基因沉默的效率, TRV介导基因沉默的最佳温度为22℃左右[ 27]。本研究中, 克隆并选择棉花标记基因 GhCLA1插入载体的片段大小为421 bp, 目标验证基因 GhMAPKKK插入载体的片段大小为556 bp。选择播种后两周左右, 子叶完全展开而真叶尚未伸出时是接种病毒的最佳时期。同时农杆菌侵染棉花后, 选择23℃左右的生长温度, 沉默效果好。
陆地棉是异源四倍体, 基因组预测分析显示, 大约含7万个基因[ 28], 而2012年释放的四倍体棉种亲缘关系最近的D基因组供体种―雷蒙德氏棉, 其基因组中约含3.7万基因[ 16], 需要建立高通量、准确、易操作体系快速揭示基因功能。基于下胚轴的棉花组织培养转基因再生体系周期长, 受基因型影响大。建立的VIGS技术具有简单、快速、高通量、不受遗传背景影响等优点, 利用TRV-VIGS体系高通量研究棉花基因功能, 特别是快速查明表型明显, 作用大的重要质量性状基因, 与棉花生长发育、抗病抗逆等密切相关基因, 以及关键代谢通路基因, 将为目标性状基因的快速发掘与育种应用提供重要基因资源与依据。Gao等[ 14]通过TRV病毒抑制棉花抗黄萎病相关基因 GhNDR1和 GhMKK2表达, 阐明了其在棉花抗病研究中利用潜力。Gao等[ 29]通过棉花黄萎病菌侵染后蛋白组学差异表达分析, 结合差异表达基因的VIGS功能验证, 揭示棉酚、油菜素类酯、茉莉酸等对棉花抗黄萎病性具有重要作用。本研究利用建立的棉花TRV-VIGS体系, 证实沉默 GhMAPKKK基因后, 抗黄萎病的海岛棉海7124植株变得更易感病, 说明 GhMAPKKK基因参与了黄萎病诱导的过敏反应, 在棉花黄萎病抗病机制中发挥作用。
利用烟草脆裂病毒载体建立了基于病毒诱导的棉花基因沉默体系。利用该体系, 病毒可高效扩散到受体的根、茎、叶等器官, 受棉花不同遗传背景影响小。应用该体系, 成功抑制了受黄萎病菌显著上调的棉花 GhMAPKKK基因表达。与对照株相比, GhMAPKKK抑制表达的棉花植株接种黄萎病菌后, 更易感病。建立的棉花TRV-VIGS体系具有广谱性、灵敏性和高通量等特点, 将加速棉花重要基因鉴定和功能分析研究进程。
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|