冬小麦低温下基因表达谱分析及海藻糖基因家族研究

【摘要】温度是限制作物地理分布和生长发育的主要因素之一,低温会给作物的产量造成巨大损失。小麦是世界上主要的粮食作物,但也是受低温影响严重的作物,低温极大的限制了小麦在高纬度地区的分布。实践证明,选育抗寒小麦品种不仅可以降低低温造成的减产,而且是发展高寒地区小麦种植的有效手段。因此,研究小麦的抗寒机制以及对抗寒基因转录的整体水平进行探索,对于全面了解小麦的抗寒分子基础以及培育抗寒新品种具有重要意义。本研究以强抗寒品种东农冬麦1号和弱抗寒品种济麦22为试验材料,分别在4℃6h、4℃30d和-12℃6h三个低温处理下进行分蘖节取材,采用Solexa高通量测序技术进行转录组分析,从中筛选出大量的差异表达基因及重要的代谢途径。为了进一步验证表达谱中两个冬小麦品种的差异代谢途径抗寒性的影响,在室外自然降温过程中对两个品种分蘖节的生理指标及膜质脂肪酸含量进行了比较分析,对外源施加海藻糖对冬小麦抗寒性的影响进行了研究,并在-20℃处理条件下对两个品种的海藻糖-6-磷酸合成酶（TPS）家族基因进行了荧光定量PCR分析。本研究为小麦抗寒遗传改良提供了基因资源,也为进一步了解冬小麦越冬分子机制奠定基础。主要研究结果如下：1.利用Solexa/Illumina测序平台,共获得224886contigs,平均长度为272nt,N50长为378。经过聚类,contigs被进一步组装成unigenes,数量为80704,平均长度为686nt,N50长为1036nt。大约有64%的unigenes被比对到了Nr数据库中。而16551个unigenes没有与任何序列比对上。AlA、A2B、A3B和B3A的cleanreads数目达到2.5M时,B1A和B2达到2M时,检测到的基因数目几乎趋于饱和,基本不在上升,说明这些样本的测序量已经基本覆盖到了取样部位细胞中正在表达的全部基因。2.为了验证数字表达谱测序结果中基因表达量的准确性,运用RT-PCR对随机选取的差异表达基因进行了表达量的分析。8个随机差异unigenes的24个数据点中有19个表达趋势及方向与表达谱的结果一致,一致度达到79%。由此证明Solexa测序的结果达到了进一步试验分析的要求。3.通过对三个低温处理下两个冬小麦品种差异表达基因的分析,发现符合FDR<0.001且|log2Ratio|≥2的差异表达基因有6471个。三个低温处理共有的差异表达基因439个,且在东农冬麦1始终上调表达且被注释功能的基因共有27个,其中被注释为碳水化合物代谢过程的unigene33313B2和NBS-LRR抗病蛋白功能的unigene39910B2以及细胞质囊功能的unigene49064B2等三个基因在三个低温处理下的表达量在东农冬麦1号中均比济麦22上调了10倍以上,另外还有17个unigenes在驯化时期在东农冬麦1号中的表达量比济麦22上调了10倍以上,而进入冷冻期这些基因的表达量开始下降,但是仍高于济麦22。这些基因多被注释为结合功能、激酶活性、逆境响应、氨基酸蛋白激酶和抗病蛋白等功能4.经过GO分析,有22724个序列在GO数据库中具有注释,而这些序列又进一步被分配到生物过程、细胞组件和分子功能等子集中,其中大多数unigenes功能集中在生物过程的细胞过程和代谢过程,细胞组件的细胞、细胞部分和细胞器,分子功能的结合和催化活性等生物学过程。在此基础上分别对三个低温下的分子功能部分进行了分析,并且以correctedp-value<0.05为阈值,两个冬小麦品种在经历4℃驯化6h和4℃驯化30d两个处理时GO注释得到的基因功能分类明显少于-12℃处理6h,且在前两个处理（4℃）下结合活性较-12℃处理显著富集。而在植株经历-12℃处理后,抗氧化活性和离子结合等能力明显增强。并且抗氧化活性、氧化还原酶活性和糖基转移酶活性等三个与低温抗性密切相关的子集仅在-12℃处理下被注释到。5.18764个序列经过COG分类后被聚成25类,注释为一般功能的unigenes有7151个,在25类中比例最高达到38.11%,而细胞外结构和细胞核结构比例最低分别为0.04%和0.01%,unigenes数量分别为24和7个。其中被注释为与脯氨酸代谢有关的氨基糖转移和代谢的unigenes有2128个;与可溶性糖代谢有关的碳水化合物转移和代谢的unigenes有3408个；与细胞膜脂不饱和脂肪酸代谢有关的脂质转移和代谢和细胞膜/外膜的unigenes分别有1713和4207个。6.对三个低温处理组进行pathway分析后,分别鉴定出113、103和100条pathway,且分别注释到1519、995和732个差异基因。通过比较三个低温处理下基因显著富集pathway的结果,发现氨基糖和核苷酸糖代谢、淀粉和蔗糖代谢、氧化磷酸化、半乳糖代谢、丙氨酸,天冬氨酸和谷氨酸代谢及甘氨酸,丝氨酸和苏氨酸代谢等6条pathway在东农冬麦1号经历-12℃处理6h后显著富集,这些pathway与植物抗逆密切相关。7.对两个冬小麦品种表达谱中的差异基因进行功能分类,发现被注释为转录因子功能的unigenes共有51个,其中CBF的19个,WRKY有16个,ERF有6个,MYB有5个,NAC有3个还有2个Zinc-finger;被注释为抗病蛋白TBS-LRR的有24个unigenes,RGA的有9个,抗病蛋白有11个,病程相关蛋白有2个；被注释为逆境响应蛋白膜蛋白有8个unigenes,冰再结晶阻遏蛋白有9个,WCOR有6个,热激蛋白有2个,诱导响应蛋白有2个以及1个脱水蛋白；被注释为谷胱甘肽-抗坏血酸循环关键酶共有21个unigenes;被注释为糖代谢关键酶的糖基转移酶有49个,葡聚糖酶有9个,葡萄糖-6-磷酸脱氢酶有3个；注释为脂肪酸合成及降解的unigenes共有19个,其中脂肪氧合酶有10个,磷脂酶A有4个,FabF有2个,FabG有2个,FabI有1个。8.因为在前期的数字表达谱分析中发现两个冬小麦品种在低温下的很多代谢途径都存在差异,所以为了更进一步的洞察东农冬麦1号和济麦22两个抗寒性不同的品种的遗传差异。对两个品种在低温下的一些抗寒生理指标进行了测定,在室外自然降温的过程中可溶性糖和脯氨酸在东农冬麦1中含量明显高与济麦22,说明它们在东农冬麦1中起到了更强的渗透调剂作用,从而增加了分蘖节的保水能力,降低了质膜的过氧化程度进而降低丙二醛含量。9.对冬小麦低温下数字表达谱分析发现,某些脂肪酸代谢途径相关基因的表达量在两个品种中表现出明显差异。所以对低温下冬小麦的膜脂脂肪酸含量及组分配比进行了测定分析,东农冬麦1号的不饱和脂肪酸含量高于济麦22,这进一步说明了极端低温下东农冬麦1号的膜脂稳定性高于济麦22。10.通过对低温下冬小麦幼苗外源施加海藻糖,发现适当浓度的外源海藻糖可以在一定程度上提高冬小麦的抗寒性。鉴定了9个小麦海藻糖-6-磷酸合成酶（TPS）基因并对它们的结构、进化及功能进行了分析。9个TaTPS基因中均包含HAD-likedomain、UDP-forming、20Glycotransf20和TrehalosePpase4个保守结构域。11.经过比较发现小麦和水稻的TPS在进化树中关系较近,其中OsTPS2vs.TaTPS8,OsTPS4vs.TaTPS4、5,OsTPS5vs.TaTPS9在进化树中分别被分在同一个亚组之中,相对应的两个物种的TPS基因可能具有相同或相似的生物学功能。12.在-20℃冷冻条件下分析了两个冬小麦品种中9个TaTPS基因的表达模式,TaTPS1、3、7、8、9在东农冬麦1号的根、茎中各时间点均表现出上调表达,对东农冬麦1号的抗寒性起到重要的作用。TaTPS8在济麦22的根、叶中保持上调表达；TaTPS4和TaTPS7在济麦22的茎中各时间点也都表现出上调表达,但是上调并不显著。
【作者】谢冬微；
【导师】李卓夫；
【作者基本信息】东北农业大学，作物遗传育种，2014，博士
【关键词】冬小麦；低温；表达谱；海藻糖-6-磷酸合成酶；

【参考文献】
[1]陈作才.我国公职人员财产申报制度研究[D].吉林大学,法律,2014,硕士.
[2]马威.小城镇污水处理适用技术工艺的分析与评价[D].西安建筑科技大学,市政工程,2013,硕士.
[3]檀建群.二氧化碳加氢合成低碳烯烃催化剂研究[D].长春工业大学,工业催化,2013,硕士.
[4]孔淑兰,张焕水,张召生,张承慧.无线通讯网络功率和流量的预联合控制(英文)[J].自动化学报,2007,07:761-764.
[5]张颖.论破产公司董事的义务[D].中国政法大学,民商法学,2014,硕士.
[6]刘珊.初期到中期北原白秋诗风的转变[D].河北大学,日语语言文学,2014,硕士.
[7]张宽.蒙特卡罗粒子输运模拟中反应堆物理关键参数计算方法研究与实现[D].中国科学技术大学,核能科学与工程,2014,硕士.
[8]谢继勇.基于神经网络的港口清淤工程预警系统[D].天津大学,2004.
[9]张霄力,赵军.任意切换下不确定线性切换系统的鲁棒镇定[J].自动化学报,2002,05:859-861.
[10]汪海峰.水热法制备多铁性材料BiFeO_3及性能研究[D].东北大学,材料学,2010,硕士.
[11]蒋玥妍.服务贸易对中国经济增长影响的实证研究[D].华东师范大学,国际经济与贸易,2013,硕士.
[12]彭洁.《查泰莱夫人的情人》中的凝视[D].湖南科技大学,外国语言文学,2013,硕士.
[13]曹江卫.矿用自卸车车厢的力学特征及优化设计研究[D].辽宁工程技术大学,车辆工程,2012,硕士.
[14]王月.北京城郊土壤镉对植物、微生物的毒害及其影响因素研究[D].首都师范大学,环境科学,2014,硕士.
[15]张敏.印度城市化问题探析[D].河北师范大学,世界史,2012,硕士.
[16]洪国城.城中村私有住房管理机制研究[D].华侨大学,建筑学（专业学位）,2013,硕士.
[17]本报记者安仲文.紧盯稀有资源基金潜伏“类稀土”概念[N].上海证券报,2011-03-21S03.
[18]吴煌鑫.广州市北部山区扶贫开发工作调查研究[D].电子科技大学,公共管理（专业学位）,2012,硕士.
[19]潘燕.基于GPS/GPSOne混合定位嵌人式车载导航系统的设计与实现[D].电子科技大学,软件工程（专业学位）,2012,硕士.
[20]蒋云龙.前置反硝化生物滤池处理生活污水的脱氮试验研究[D].哈尔滨工业大学,土木工程,2014,硕士.
[21]胡宝珍.唐代诗与画之相关性研究[D].西北大学,中国古代文学,2014,硕士.
[22]张干.三价铬体系超声—脉冲电沉积Cr/SiC复合镀层的研究[D].湖南工业大学,材料加工工程,2012,硕士.
[23]金兰.媒体栏目收视率影响因素研究[D].浙江工业大学,2013.
[24]李文升.分布式资源与电网相互作用的机理及其协同调度技术的研究[D].山东大学,电气工程（专业学位）,2012,硕士.
[25]杨意.日本学生汉语声调习得难度等级考察[D].南京大学,汉语国际教育,2013,硕士.
[26]高燕.纳米结构Co/ZnO颗粒膜的结构、磁性和磁电阻效应[D].山西师范大学,化学,2014,硕士.
[27]张德林.二连盆地阿尔善浅气藏的地震特征[J].石油地球物理勘探,1999,01:107-118+122.
[28]巫洪章.能源材料Al/Ni复合物快速放热和g-C_3N_4光解水性能的理论模拟研究[D].上海大学,2013.
[29]顾珏弘.透过法国媒体分析法语中的英语借词现象[D].上海外国语大学,法语语言文学,2012,硕士.
[30]邬霞.成骨细胞玻璃化冷冻保护剂的研究[D].大连理工大学,化学工程,2004,硕士.
[31]张涛.大鼠Ⅲ度房室传导阻滞模型中慢病毒介导的RNA干扰抑制心室肌KCNJ2基因表达对心室率影响的研究[D].山东大学,外科学,2014,博士.
[32]钱荣钧.文北三维偏移速度的选取和偏移效果分析[J].石油地球物理勘探,1986,05:494-502+539.
[33]何海霞.论公司成立前交易的责任承担[D].中国社会科学院研究生院,商法（专业学位）,2014,硕士.
[34]周洁琼,王立德,王涛,申萍.基于虚拟链路交换式以太网的列车通信网络可靠性分析[J].铁道学报,2014,02:55-61.
[35]史雪飞,胡纪五.LONWORKS总线式有毒气体智能仪表的研制[J].冶金自动化,2000,03:27-30.
[36]苏红键,赵坚.产业专业化、职能专业化与城市经济增长——基于中国地级单位面板数据的研究[J].中国工业经济,2011,04:25-34.
[37]黄娇.Cu-ligand基团修饰的多酸簇合物的合成与性质研究[D].河北师范大学,无机化学,2012,硕士.
[38]曹志超.清水河县城镇体系结构优化与可持续发展研究[D].内蒙古师范大学,人文地理学,2013,硕士.
[39]杨曦.认知语言学视域中英汉仿拟的运作机理及文化内涵探究[D].西安外国语大学,语言学及应用语言学,2013,硕士.
[40]何石轩.考虑回路控制的间歇过程优化与控制方法研究[D].东北大学,控制理论与控制工程,2010,硕士.
[41]孙瑞玲.基于石墨烯复合材料及石墨烯量子点的制备并用于一些酚类环境污染物的分析[D].南昌大学,2014.
[42]马世伟.激光跟踪与欺骗干扰演示系统研究[D].哈尔滨工业大学,微电子学与固体电子学,2014,硕士.
[43]明慧芳.视频信号的数字加密系统研究[D].武汉理工大学,信号与信息处理,2004,硕士.
[44]宋瑞荣.T型有机器件自旋注入性质的理论研究[D].山东师范大学,凝聚态物理,2013,硕士.
[45]李纪恩.上市公司吸收合并中股东权益保护问题之研究——从TCL吸收合并案谈起[D].对外经济贸易大学,国际法学,2004,硕士.
[46]隋思慧.子品牌与“小群分众广告”[D].吉林大学,传播学,2013,硕士.
[47]李永芳.基金会资产运营与资金增值研究[D].浙江大学,法律（专业学位）,2012,硕士.
[48]康华靖,陶月良,权伟,王伟,欧阳竹.植物光合CO_2响应模型对光下(暗)呼吸速率拟合的探讨[J].植物生态学报,2014,12:1356-1363.
[49]秦丽.装饰性语言在摄影艺术中的应用与研究[D].哈尔滨理工大学,设计艺术学,2013,硕士.
[50]彭博.基于对象特征的高分辨率遥感影像特定目标提取[D].辽宁工程技术大学,地图制图学与地理信息工程,2012,硕士.