甘蔗应答黑穗病菌侵染的转录组与蛋白组研究及抗性相关基因挖掘

【摘要】甘蔗(Saccharumofficinarum)是最重要糖料作物,蔗糖占我国食糖总产92%。病害是影响甘蔗产量和质量的重要因素,甘蔗黑穗病(Sporisoriumscitamineum)是世界各植蔗区最主要的病害之一,也是我国甘蔗栽培上最重要的病害,种植抗病品种是控制该病最经济、最有效的措施。前人从相关细胞学、形态学、生理生化及抗病性遗传方面,探讨了甘蔗对黑穗病的抗性机制。目前,有关甘蔗响应黑穗病菌侵染的分子机制研究,基本还停留在以抗病甘蔗基因型接种病原后的材料进行个别几个差异表达基因的克隆与定量表达研究,这些研究初步反映了甘蔗抗黑穗病菌的分子机制受到多基因网络系统的调控,分子水平的机理与机制研究有待继续深入。利用高通量的技术手段,从转录组学和蛋白组学水平上,研究甘蔗对黑穗病菌胁迫的应答反应,有助于从整体上揭示甘蔗响应黑穗病菌侵染的代谢通路和分子调控网络,也有助于筛选与鉴定起关键作用的基因。本研究建立了一种基于实时荧光定量PCR技术检测甘蔗黑穗病菌的方法,在黑穗病症状未出现前快速、灵敏地检测病原菌,为甘蔗黑穗病菌的快速检测和早期预防提供依据。此外,以甘蔗黑穗病抗性基因型崖城05-179和感病主栽品种”ROC”22为植物材料,通过病原菌胁迫,建立了甘蔗—黑穗病菌之间的不亲和与亲和互作系统,利用高通量RNA-Seq测序技术与iTRAQ技术,通过多点时序比较,分析甘蔗响应黑穗病病原菌胁迫过程中,同一甘蔗基因型不同接种时间点和不同甘蔗基因型不同接种时间点的mRNA及蛋白表达水平的变化,构建甘蔗应答黑穗病菌侵染的表达谱,对于从分子水平上揭示甘蔗抗黑穗病性机制具有重要的理论指导和实践意义。实验还筛选鉴定了3类抗性相关基因,包括p-1,3-葡聚糖酶基因(β-1,3-glucanase)、几丁质酶基因(chitinase)和过氧化氢酶基因(catalase),旨在为甘蔗抗病分子育种提供优异基因资源。初步研究成果如下：1.甘蔗黑穗病是一种真菌性病害,在蔗区广泛流行,该病菌的早期诊断和准确检测将有助于该病害的治理。针对黑穗病菌bE基因(bEastmating-type,GenBank登录号为U61290.1)设计特异性的实时荧光定量PCR检测引物(bEQ-F/bEQ-R)与TaqMan探针,其检出限为10agbE质粒DNA和0.8ng甘蔗基因组DNA,比传统PCR(检出限为10fgbE质粒DNA和100ng甘蔗基因组DNA)更灵敏。实验通过检测人工接种甘蔗组培苗(FN40)进一步验证TaqMan技术的可靠性,结果显示,该技术在病原菌感染初期(12h)能够检测到黑穗病菌,适用于甘蔗无病种茎和幼苗的诊断。此外,通过不同甘蔗抗/感基因型(崖城05-179和”ROC”22)受黑穗病菌侵染后的定量检测实验,结果表明,该TaqMan技术具有进一步开发成为评价甘蔗基因型对黑穗病抗性技术的潜力。总的来说,所建立的技术体系可作为检测和诊断甘蔗黑穗病原菌的工具,实现灵敏、准确、快速和定量检测的目的。2.以甘蔗抗黑穗病无性系崖城05-179和感黑穗病品种”ROC”22蔗芽为材料,以清水处理24h时间点为对照,对接种甘蔗黑穗病菌24h、48h和120h时间点的蔗芽进行IlluminaRNA-Seq(测序)。借助维恩图(Vennydiagram),通过多点时序比较,分析甘蔗响应黑穗病病原菌侵染过程中,同一基因型不同时间点和不同基因型同一时间点的差异基因表达变化,构建了甘蔗响应黑穗病病原菌侵染过程的转录组表达谱。主要结果如下：(1)8个甘蔗样品转录组测序共获得36.68Gb数据,181,603,016对reads.分样品denovo组装后共获得148,605条Unigene,其中长度在500bp以上的有46,525条,长度在1.0kb以上的有20,798条；合并测序数据denovo组装后共获得99,824条Unigene,其中长度在500bp以上的有47,345条,长度在1.0kb以上的有22,091条。进行基于合并组装Unigene库的基因结构分析,包括ORF预测、SSR(simplesequencerepeat)分析以及样品间SNP(SingleNucleotidePolymorphism)分析,获得SSR标记5,095个。将Unigene与NT、NR,SwissProt、TrEMBL、COG、GO和KEGG数据库进行比对,共获得65,852条Unigene的生物信息学注释结果。另有33,972条Unigene(34.03%)未被注释到,推测该部分序列可能是新转录本,需进一步验证。本实验构建的转录组数据库可为甘蔗其他相关研究的基因注释提供参考。(2)进行各样品中基因表达量和差异表达基因的分析,以及对差异表达基因进行模式聚类、功能注释和富集性分析,并通过Q-PCR进行验证。结果显示：①黑穗病菌侵染后,两个甘蔗基因型被诱导的大多数差异表达基因和代谢通路都是常见的,几乎涵盖整个生命活动的各个方面。比对上的Unigene涉及细胞组件(cellularcomponent)、分子功能(molecularfunction)和生物过程(biologicalprocess)等功能。在此分析基础上,获得了与抗性相关的功能注释信息,如信号转导机制(Signaltransductionmechanisms)、能量产生与转导(Energyproductionandconversion)、离子转运机制(Inorganiciontransportandmetabolism)和防御机制(Defensemechanisms)等。Pathway聚集分析显示,差异基因参与了包括植物激素信号转导途径(Planthormonesignaltransductio)、类黄酮生物合成途径(Flavonoidbiosynthesis)、植物与病原菌互作途径(Plant-pathogeninteraction)和细胞壁防御途径(Cellwallfortificationpathway)等抗病相关代谢途径。可见,黑穗病菌激发了甘蔗多种抗病途径、差异表达基因涉及防御反应及信号转导等,提示甘蔗对黑穗病病原菌侵染响应的分子机制受到多基因网络系统的调控。②随着接种黑穗病菌时间的延长,激活的甘蔗差异表达基因数目增多,且上调表达的差异基因明显多于下调表达的差异基因。黑穗病菌在崖城05-179和”ROC”22中诱导的差异基因在很大程度上是相同,但总的来说,崖城05-179中分析到的显著性抗病相关代谢通路明显多于”ROC”22,且抗性相关基因的转录表达时间(24h～48h)较”ROC”22的(48h-120h)早,此外,抗病品种中差异转录本的上／下调表达情况相对较丰富,提示它们在非亲和互作中表现抗性特异性以及较早的时序性。以上应答黑穗病菌的甘蔗转录本潜在功能的分析,有助于我们更深入研究这些候选基因在甘蔗响应黑穗病病原菌侵染中的作用,同时也为开展甘蔗抗病功能基因组学研究奠定基础。3.在转录组学研究的基础上,以抗病和感病基因型接种黑穗病菌48h时间点的蔗芽为材料,利用iTRAQ技术,研究其蛋白组学的变化。结果表明：鉴定到的蛋白数量为4,251个,其中崖城05-179筛选到的差异表达蛋白有273个,“OC”22有341个。差异蛋白显著富集于苯丙氨酸代谢、苯丙烷类生物合成、次级代谢产物生物合成、苯并恶嗪生物合成、脂肪酸代谢和不饱和脂肪酸生物合成代谢通路中等,这些代谢通路直接或间接参与甘蔗对黑穗病菌的抗病反应。崖城05-179的蛋白组与转录组的关联性为0.1502,而”ROC”22的关联性为0.2466。甘蔗抗病和感病基因型中筛选出的与转录组关联的差异蛋白数分别为27个和10个,与差异基因表达趋势相同的差异蛋白的数量分别只有18个和10个,但均为关联上调表达,且其中近半数以上的关联差异蛋白与植物抗逆性密切相关,如过氧化物酶、病程相关蛋白PR-1、p-1,3-葡聚糖酶、热休克蛋白和植物凝集素等。研究结果为后续确定深入研究的代谢通路以及关键蛋白的鉴定提供依据。4.甘蔗响应黑穗病菌胁迫后,在转录组测序结果中,发现病程相关蛋白(pathogenesisrelatedproteins)的转录本差异表达,涉及10个葡聚糖酶基因(Glucanase)、26个几丁质酶基因(Chitinase)和1个过氧化氢酶基因(Catalase)。蛋白组与转录组关联分析结果也揭示了病程相关蛋白PR-1和p-1,3-葡聚糖酶等被甘蔗黑穗病菌诱导上调表达。植物病程相关蛋白是在某种病理或者与病理相关的环境条件下,植物体内诱导产生的的特异蛋白质,这些蛋白质在植物抗病反应过程中起重要作用,和植物系统获得性抗性(systemicacquiredresistance,SAR)紧密相关。本研究根据该3类典型PR蛋白的差异转录本信息,借助RT-PCR技术从甘蔗中分离出2个β-1,3-葡聚糖酶基因ScGluA1和ScGluD1、1个几丁质酶基因ScChi和1过氧化氢酶基因ScCATl。通过Q-PCR技术,分析差异表达基因在黑穗病菌侵染甘蔗过程中的时空表达特性,并应用基因工程手段,进行融合蛋白亚细胞定位及目的基因遗传转化研究,为后续甘蔗抗黑穗病遗传改良提供优异基因资源。(1)p-1,3-葡聚糖酶可以水解病原真菌细胞壁的主成分p-1,3-葡聚糖,已被证实在植物—病原菌互作中被诱导表达,且参与植物防御反应。本研究中,接种黑穗病菌后,甘蔗β-1,3-葡聚糖酶活性在抗病品种中较感病品种增加速度快且维持时间长,表明β-1,3-葡聚糖酶可能与甘蔗抗黑穗病性之间呈正相关。此外,实验从甘蔗接种黑穗病菌后的材料中,克隆到两个β-1,3-葡聚糖酶基因ScGluA1(GenBank登录号为KC848050)和ScGluD1(GenBank登录号为KC848051)。ScGluAl和ScGluD1推导的编码蛋白均存在糖苷水解酶第17家族的保守结构域,但两者之间的氨基酸序列一致性只有29.25%。系统发育树分析显示,ScGluA1和ScGluD1分别隶属于β-1,3-葡聚糖酶亚家族A和亚家族D。采用根癌农杆菌(Agrobacteriumtumefaciens)介导法,将携带目标基因和报告基因(绿色荧光蛋白,GFP)的重组菌株侵染洋葱表皮细胞,结果显示ScGluA1和ScGluD1定位于细胞外或质膜上。所构建的含pET32a-ScGluA1或pET32a-ScGluD1的Rosetta(DE3)菌株原核表达细胞,在添加NaCl、CdCl2、PEG、CuCl2或ZnS04的培养基中,表现出不同的耐受性。Q-PCR分析结果表明,ScGluA1在应答黑穗病菌中上调表达,ScGluD1基因小幅下调表达,暗示ScGluA1可能涉及甘蔗对黑穗病菌的防御反应,而ScGluD1可能被抑制。受其他环境因子(如水杨酸、茉莉酸甲酯、脱落酸、NaCl、CdCl2和自然干旱胁迫)处理后,ScGluA1和ScGluD1的基因表达模式与生物胁迫类似。在烟草中过表达ScGluA1,对烟草茄病镰刀菌蓝色变种(Fusariumsolanivar.coeruleum)和烟草灰葡萄孢菌(Botrytiscinerea)有抑制效果,且转ScGluA1烟草T0代植株叶片β-1,3-葡聚糖酶粗酶液,对烟草茄病镰刀菌蓝色变种的菌丝生长有一定的抑制作用。总的来说,β-1,3-葡聚糖酶可能作为甘蔗对黑穗病菌防御机制中的一个组分。ScGluA1的积极应答或许对甘蔗抗病性和抗逆反应有益,而ScGluD1的消极应答,显示其在甘蔗响应生物和非生物胁迫中,可能发挥不同的作用。(2)许多植物与病原菌互作会诱导产生几丁质酶(EC3.2.2.14),报道显示其与植物防御病原菌相关。本研究发现,几丁质酶活性与甘蔗黑穗病抗性之间呈正相关。实验克隆了一个ClassⅢ几丁质酶基因ScChi(GenBank登录号为KF664180),编码一个31.37kDa多肽。亚细胞定位显示,ScChi定位于细胞核、细胞质和质膜上。Q-PCR结果表明,ScChi在叶和蔗皮组织中高度表达。受黑穗病菌侵染后,ScChi转录水平在抗病品种中较感病品种高,且维持时间较长。ScChi同时受SA、H2O2、MeJA、ABA、NaCl、CuCl2、PEG和低温(4℃)胁迫的诱导表达。ScChi瞬时表达后的目标基因的转录水平及6个本氏烟(Nicotianabenthamiana)免疫相关标记基因均上调,而与对照相比,过表达pCAMBIA1301-ScChi的本氏烟叶片经组织化学分析,结果显示出更深的DAB染色效果及更高的导电率,表明存在高水平的H2O2积累。此外,瞬时表达ScChi有助于提高烟草对青枯菌(Pseudomonassolanacearum)、茄病镰刀菌蓝色变种和灰葡萄孢菌的防御作用,同时,转ScChi烟草T0代植株叶片的几丁质酶粗酶液,可以抑制茄病镰刀菌蓝色变种的菌丝生长。这些结果暗示,ScChi的表达与植物免疫之间关系密切。总的来说,ScChi能积极响应生物和非生物胁迫,表明该基因在甘蔗中是一类与逆境胁迫相关的抗病基因。(3)过氧化氢酶是铁卟啉酶,它作为一种有效的活性氧(ReactiveOxygenSpecies,ROS)清除剂,可以保护植物避免氧化损伤。过氧化氢酶的酶活性在甘蔗抗病品种(崖城05-179)应答黑穗病菌中总是高于感病品种(柳城03-182),表明过氧化氢酶活性可能与甘蔗黑穗病抗性呈正相关。为了从分子水平上了解过氧化氢酶的功能,实验从受黑穗病菌侵染的甘蔗材料中,分离到一个ScCAT1(GenBank登录号为KF664183)的cDNA序列。该基因预测编码492个氨基酸残基,其推导的蛋白序列与其他植物的过氧化氢酶同源性高。ScCAT1的重组大肠杆菌Rosetta细胞,在CuCl、CdCl2和NaCl的胁迫下长势更快,表明其耐受能力更强。Q-PCR分析结果显示,ScCAT1在蔗芽中的表达量相对较高,而在蔗皮和蔗髓的表达量中等。不同外源胁迫包括甘蔗黑穗病菌侵染、植物激素(SA、MeJA和ABA)处理、氧化胁迫(H2O2)、重金属(CuCl2)和高渗透(PEG和NaCl)胁迫,均诱导ScCAT1基因上调表达。亚细胞定位结果显示,ScCAT1定位在质膜和细胞质。此外,DAB、台盼蓝染色和电导率测定等组织化学分析表明,ScCAT1基因可能涉及甘蔗的免疫应答。总的来说,ScCAT1正响应生物和非生物胁迫的研究结果表明,ScCAT1参与保护甘蔗应答与氧化反应相关的环境刺激。
【作者】苏亚春；
【导师】许莉萍；
【作者基本信息】福建农林大学，农业生物技术，2014，博士
【关键词】甘蔗；甘蔗黑穗病菌；转录组；蛋白组；β=1,3-葡聚糖酶；几丁质酶；过氧化氢酶；

【参考文献】
[1]王立慧.基于PFC的桩承式路堤土拱效应研究[D].西安建筑科技大学,岩土工程,2013,硕士.
[2]胡歆婷.商业银行理财产品监管的价值选择和制度构建[D].上海大学,法学理论,2014,硕士.
[3]张思思.京津冀地区交易效率与经济增长研究[D].兰州商学院,区域经济学,2013,硕士.
[4]李诚诚.印度的经济改革：动因、成果与前景[D].上海国际问题研究院,国际政治,2014,硕士.
[5]董春霞.凝血及纤溶指标的变化与冠心病的关系[D].山西医科大学,血液病学,2003,硕士.
[6]黄夏.论我国反垄断损害赔偿制度的构建[D].山东大学,经济法学,2013,硕士.
[7]张伯明.现代能量控制中心概念的扩展与前景展望[J].电力系统自动化,2003,15:1-6.
[8]孙阔.阿司匹林合成过程的拉曼光谱监控[D].广西科技大学,生物化工,2013,硕士.
[9]李方.明清时期湖南佛教造像的区域性比较研究[D].湖南工业大学,美术学,2012,硕士.
[10]张漫莉.分泌表达纤维素酶EGA的重组枯草芽孢杆菌构建及其高密度发酵[D].浙江理工大学,生物化学与分子生物学,2013,硕士.
[11]安静.城市化对区域景观格局和景观生态质量影响研究[D].南京师范大学,环境科学,2012,硕士.
[12]张玲琳.北京市交通运输能源消耗指数分解研究及与经济增长关系动态分析[D].北京交通大学,2014.
[13]李锦文.基于网络模式的复杂三维集装箱问题求解方法研究[D].兰州理工大学,计算机应用技术,2012,硕士.
[14]宋阳.新型嘧啶（1，3-二氧六环）并二氢嘧啶酮衍生物的合成研究[D].西北大学,化学生物学,2014,硕士.
[15]李利克.拒不支付劳动报酬罪案研究[D].黑龙江大学,法律,2013,硕士.
[16]李海桥.高速冲击问题的光滑粒子流体动力学的数值计算[D].中北大学,动力机械及工程,2014,硕士.
[17]张微.基于相场模型的软性磨粒流两相流场特性分析方法[D].浙江工业大学,机械电子工程,2012,硕士.
[18]赵春晖.光合细菌和重组大肠杆菌生物合成5-氨基乙酰丙酸的研究[D].浙江大学,生物化学工程,2004,硕士.
[19]王鑫.《朝花夕拾》汉英文本中语法衔接手段的对比研究[D].西南交通大学,外国语言学及应用语言学,2012,硕士.
[20]刘少成.叙大铁路聂家坡隧道爆破地震效应对地表建筑物的影响研究[D].西南交通大学,岩土工程,2013,硕士.
[21]胡霸桥.基于ZigBee技术的农业大棚环境温度监测系统设计[D].哈尔滨理工大学,通信与信息系统,2012,硕士.
[22]马文景.论“隐含作者”对凯特·肖邦短篇小说主题多样性的作用[D].鲁东大学,英语语言文学,2014,硕士.
[23]阳双梅.企业商业模式价值网络系统的复杂性研究[D].华侨大学,工商管理,2013,硕士.
[24]陈海焱,陈金富,段献忠.含风电场电力系统经济调度的模糊建模及优化算法[J].电力系统自动化,2006,02:22-26.
[25]马玉宝.基于数据仓库的知识发现模式的研究[D].合肥工业大学,计算机软件与理论,2003,硕士.
[26]薛丽萍.海南省碳基金的构建研究[D].海南大学,金融学,2014,硕士.
[27]蔺小军,杨艳,吴广,高源,陈悦,刘敏,刘明星.面向叶片型面的五轴联动柔性数控砂带抛光技术研究[J].航空学报.
[28]胡倩.城市化进程中武汉市花山中学发展战略研究[D].华中师范大学,教育经济与管理,2014,硕士.
[29]张春生.基于校园网的教务管理信息系统的设计与实现[D].大连理工大学,计算机技术（专业学位）,2013,硕士.
[30]赵珑.城市生活垃圾管理信息系统的开发[D].昆明理工大学,2004.
[31]顾莉洁.数字证书管理系统DCMS的研究与设计[D].苏州大学,计算机应用技术,2003,硕士.
[32]EvadneCowell.信息通信技术政策在中国和牙买加小学中实施情况的对比分析[D].东北师范大学,教育技术学,2012,硕士.
[33]常丽.情感营销和视觉表现的本土化研究[D].浙江理工大学,设计艺术学,2014,硕士.
[34]姜广宗.颈椎融合与非融合手术方式的三维有限元研究[D].山东大学,外科学,2013,硕士.
[35]张宁.西安都市农业发展战略研究[D].西北农林科技大学,农村与区域发展,2012,硕士.
[36]申幸娇.三个封闭群实验兔遗传多样性研究[D].中国食品药品检定研究院,免疫学,2013,硕士.
[37]方柯青.公司治理结构对无形资产影响的实证研究[D].浙江工商大学,会计学,2013,硕士.
[38]陈娜.海南荣丰花卉苗木高新技术产业园可行性研究报告[D].中南林业科技大学,农村与区域发展,2013,硕士.
[39]王萌萌.复方醇提物对预防小鼠急性酒精中毒的药效学研究[D].湖北中医药大学,药理学,2013,硕士.
[40]徐洋洋.一维磁铅石型铁氧体的制备、形成历程及性能研究[D].沈阳理工大学,应用化学,2013,硕士.
[41]费嘉祺.基于B/S模式的客户关系管理系统的设计与实现[D].西安电子科技大学,软件工程,2012,硕士.
[42]崔羽.中医药辨证治疗功能性腹胀的探索性研究[D].北京中医药大学,中西医结合临床,2013,硕士.
[43]刘韬,王宾.银行业数据仓库构建策略[J].华南金融电脑,2004,02:94-96.
[44]陈志华.高超声速飞行器建模与轨迹跟踪控制方法[D].哈尔滨工业大学,控制科学与工程,2014,硕士.
[45]牧骞.新型干法选煤系统的研究与设计[D].陕西科技大学,机械电子工程,2013,硕士.
[46]胡玥.资源型城市可持续发展的财政支持研究[D].安徽大学,财政学,2014,硕士.
[47]徐向南.YAP蛋白在肝细胞癌表达的临床意义及其与肝星状细胞的关系[D].复旦大学,临床医学,2013,博士.
[48]田田.基于实践共同体下英语新教师成长路径的案例研究[D].宁波大学,课程与教学论,2013,硕士.
[49]李晓伟.留地安置：现行法律框架下的困境与解脱[D].西南政法大学,民商法学,2012,硕士.
[50]曹超锋.陕西邮政商旅服务项目开发研究[D].西安科技大学,工商管理,2013,硕士.