谷子NCED基因家族鉴定及其干旱胁迫响应表达模式分析.docx

    谷子NCED基因家族鉴定及其干旱胁迫响应表达模式分析    程鸿燕 郭昱 马芳芳摘要：基于转录组数据分析并结合谷子基因组数据库鉴定出10个候选NCED基因，对它们的结构特点、理化性质、启动子元件功能等进行分析，并以抗旱（GG）和干旱敏感谷子品种（JF16）为材料，对PEG胁迫前后差异基因的表达特点进行分析结果表明，谷子中编码NCED家族基因启动子中含较多与抗旱胁迫相关的功能元件;经PEG胁迫处理后，Seita.2G035400在GG和JF16中表达量均上调，但上调幅度有所差异该研究结果进一步加强了对植物NCED基因家族的了解，也为后续进行谷子抗旱机制和抗旱分子育种提供了理论借鉴依据关键词：谷子;NCED家族基因;干旱胁迫;表达;生物信息学： S515.01 文献标志码： A ：1002-1302（2019）01-0040-04谷子[Setaria italica （L.） P. Beauv.]是禾本科狗尾草属的粮草兼用作物，最早被驯化种植于中国北方，已有8 700多年的栽培史[1]，其营养价值高，富含人体所需的类胡萝卜素、氨基酸、叶酸和硒等成分[2]，至今仍被誉为“五谷之首”谷子叶面积小、自花授粉、蒸腾效率和水分利用效率高，具有耐干旱、耐贫瘠特性[3]。

但是近年来，随着全球气候变暖，可利用水资源日趋匮乏，干旱成为制约作物正常生长发育的关键因素之一发掘谷子耐旱基因，了解其基因家族在干旱胁迫条件下的表达模式，对谷子抗旱新品种选育具有极为重要的指导作用作物在干旱、盐碱、低温、冻害等非生物胁迫条件下，可通过自身调节来适应或抵制这些胁迫[4-5]脱落酸（abscisic acid，ABA）作为一种植物激素，是一种极为重要的逆境信号，不仅在植物生长发育过程不同阶段起重要作用，而且在干旱、低温、衰老、盐渍等多种逆境响应中起着重要调控作用，可通过介导叶片气孔关闭过程来抵御干旱胁迫[6-7]9-顺式-环氧类胡萝卜素双加氧酶（9-cis-epoxycarotenoid dioxygenase，NCED）属于类胡萝卡素双加氧裂解酶 （carotenoid cleavage dioxygenases，CCDs）家族成员，被认为是高等植物ABA生物合成过程中最重要的限速酶，有研究表明其表达量与植物体内ABA含量呈正相关[8-10]编码NCED的基因最初是从玉米突变体viviparousl4（vpl4）中被克隆[11]，随后在其他植物如番茄[12]、豇豆[13]、拟南芥[8，14]、柱花草[15]等中也克隆得到了NCED基因。

近年来，我国在NCED响应小麦、拟南芥、花生、烟草、水稻等作物干旱胁迫的研究中发现，干旱胁迫可以诱导NCED基因表达和内源ABA的积累，且在一定范围内，该基因表达与ABA积累随脱水延长而明显增强[16-20]，这说明NCED基因在提高植物抗旱能力方面具有重要研究价值谷子是重要的抗逆模式作物，目前，尚无对谷子中NCED家族进行系统研究的相关报道，干旱胁迫下谷子NCED基因家族情况及其结构特点和表达特征并不明晰本研究以耐旱品种勾勾母鸡咀（GG）和干旱敏感品种晋汾16（JF16）为试验材料，在人工气候室模拟昼夜条件对其进行培养，于苗期进行20%聚乙二醇6000（PEG-6000）模拟干旱处理后取材，用于表达谱测序从中找到抗旱相关基因NCED后，利用生物信息学方法从谷子基因组数据库中鉴定同源NCED基因，并系统分析其基因家族成员的结构特点、理化特性、启动子调控元件、进化关系以及该基因家族成员在干旱胁迫下的表达特点，旨在为进一步研究谷子抗旱分子育种提供有用信息1 材料与方法1.1 试验材料本试验所用植物材料为耐旱性品种勾勾母鸡咀（GG）和干旱敏感性品种晋汾16（JF16），均由山西农业大学农业生物工程研究所提供。

将这2个品种种植在混合均匀的基质（营养土 ∶ 蛭石=3 ∶ 1）上，然后将其放于人工气候室，模拟外界昼夜条件进行14 h光培养（28 ℃）和10 h暗培养（23 ℃），出苗培养21 d后，将2个品种谷子幼苗用一定量的20% PEG-6000模拟干旱胁迫处理0.5 h，而对照组幼苗用同等量的蒸馏水处理相同时间之后在对照组和处理组中的每个品种中分别随机取出3株幼苗，进行RNA提取和表达谱测序1.2 试验方法根据表达谱测序结果查询到NCED基因名，并通过查阅相关文献和代谢途径，明确其注释功能;利用Phytozome v 12.1（https：//phytozome.jgi.doe.gov/pz/portal.html）基因组数据库进行搜索，查找NCED基因家族信息利用GSDS 2.0（http：// 500 bp序列）相关元件并进行统计;采用MAGA 7.0软件对谷子NCED蛋白及拟南芥、水稻、小麦和玉米等近源物种中NCED蛋白构建进化树，进行亲缘关系比对2 结果与分析2.1 谷子NCED基因的基本信息及结构分析在Phytozome v12.1的谷子数据库中找到10个NCED家族基因，并對其基本信息（染色体位置、CDS序列和氨基酸序列）进行统计与分析，结果发现，谷子NCED家族基因分布在不同的染色体上，其中1号、4号和9号染色体上各包含1个成员，2号和3号染色体上各包含2个成员，5号染色体上的NCED家族基因最多，为3个（表1）。

谷子NCED家族基因CDS序列的平均长度为1 773 bp，不同成员的CDS序列长度之间也存在差异，其中Seita. 9G156500的CDS序列最长，为1 965 bp;Seita. 3G408300的CDS序列最短，为1 638 bp（表1）为进一步分析谷子NCED家族基因的结构特征，使用GSDS 2.0软件对NCED家族成员的基因序列和CDS序列进行结构分析，结果显示，NCED家族不同成员之间的结构存在较大差异，其中Seita. 2G035400和Seita. 3G391000具有相同的结构，都没有内含子和上下游，只有编码区序列;其余8个基因都含有内含子和编码区序列，但数目差异较大，并且Seita. 1G288400和Seita. 4G062800与其余6个基因结构差异更大，未发现上下游结构（图1）2.2 谷子NCED理化性质根据豫谷1号NCED的cDNA序列预测到的氨基酸序列可知，谷子NCED家族基因氨基酸序列的平均长度为590 aa，而不同成员的氨基酸序列长度之间存在差异，其中Seita. 9G156500的氨基酸序列最长，为654 aa;Seita. 3G408300的氨基酸序列最短，为545 aa（表1）。

蛋白质分子量在 60 585.91～70 862.56 ku之间，平均值为64 656.88 ku，其中Seita. 9G156500的蛋白质分子量最大（70 862.56 ku），Seita. 5G315700的蛋白质分子量最小（60 585.91 ku）;等电点介于5.81（Seita. 4G062800）～7.66（Seita. 5G315700）之间，均值为6.37，并且除Seita. 5G315700的蛋白质等电点呈弱碱性外，其余9个均呈酸性，可见它们发挥其最大催化活性的最适pH值有所差异（表1）2.3 谷子NCED基因家族系统进化树为进一步分析谷子NCED基因家族之间及其与其他植物NCED之间亲缘关系的远近，了解其进化过程，将所得到的谷子NCED基因编码的蛋白序列与其他物种（水稻、玉米、高粱、拟南芥、大麦、小麦和狗尾草）进行BLAST蛋白同源性检索，并下载其他物种的NCED蛋白序列，之后将其与谷子NCED基因家族成员进行比对，并构建进化树（图2）结果发现，谷子NCED基因家族中，Seita. 5G315800、Seita. 5G315600和Seita. 5G315700单独聚为一大类，暗示这3个基因虽然与谷子中一些基因亲缘关系较近，但在进化过程中它们所编码的蛋白序列存在较多氨基酸变异位点;而其余7个基因与其他7个物种的NCED聚为一大类，说明谷子的这7个基因亲缘关系较近，并且在这一大类中，谷子Seita. 2G035400与狗尾草中NCED基因 （Sevir. 2G040100）是相同的，它们与高粱的NCED亲缘关系最近，与水稻、小麦的亲缘关系次之;谷子Seita. 9G156500与大麦NCED的亲缘关系最近;谷子Seita. 3G391000与玉米NCED的亲缘关系最近;而Seita. 9G156500和Seita. 3G391000与拟南芥NCED的亲缘关系的远近各次于大麦和玉米。

这一结果与植物经典分类结果相一致2.4 谷子NCED家族基因启动子调控元件功能分析通过PlantCARE软件对谷子NCED家族基因启动子（基因上游1 500 bp）顺式作用元件进行功能预测，有利于进一步鉴定该家族基因在谷子抗旱机制中所起的作用本研究主要对激素响应元件和逆境胁迫响应元件进行了统计，结果（表2）表明，每个基因启动子中均含有厌氧诱导响应元件，有部分基因启动子中含有低温响应和热响应元件;除基因Seita. 5G315600外，其余9个基因启动子中均含有ABA响应元件;除基因Seita. 9G156500、Seita. 1G288400和Seita. 3G408300外，其余7个基因均含有MeJA响应元件，并且Seita. 5G315800中最多（为12个）;有5个基因启动子中含有SA响应元件，它们分别是Seita. 3G391000、Seita. 3G408300、Seita. 2G174500、Seita. 5G315600和Seita. 5G315700;除基因Seita. 2G035400、Seita. 9G156500和Seita. 2G174500外，其余7个基因启动子中均含有GA响应元件，并且Seita. 4G062800中最多（为3个）;有5个基因启动子均含有MYB响应元件，它们分别是Seita. 9G156500、Seita. 3G391000、Seita. 2G174500、Seita. 5G315800和Seita. 5G315600。

2.5 谷子NCED家族基因成员在干旱胁迫下的表达模式对谷子幼苗进行干旱胁迫处理后，根据表达谱测序结果，对勾勾母鸡咀和晋汾16中NCED基因表达量进行分析，这2个品种中均含有Seita.2G035400基因从图3中可以看出，GG和JF16经不同处理后其体内Seita.2G035400基因均有表达，并且经干旱胁迫处理后该基因的表达量与对照相比均升高，但该基因在干旱敏感性品种JF16中表达量增幅明显高于在耐旱品种GG中表达量增幅（图3）3 讨论与结论脱落酸ABA与植物生长发育调控、气孔开闭以及对逆境适应和抵御之间关系密切，是一种尤为重要的激素调节物质[21-22]因此，有必要分析ABA生物合成关键基因，明确其家族基因结构与功能，掌握其在ABA生物合成中所起作用，进而通过调控植物体内ABA含量来达到植物自身生理调控目的NCED基因就是参与ABA生物合成调控的一类关键基因[9]从谷子NCED家族基因启动子元件功能预测发现，有很多调控元件与植物干旱胁迫响应物质有关，如ABRE和motif IIb（CCGCCGCGCT）是ABA响应元件，而ABA可通过调节气孔开闭或者促进根系对水分吸收等途径来调控植物对干旱胁迫的适应性，Neill等研究认为，ABA是植物的。