不同降温模式下高加索三叶草的转录组比较分析

为探讨高加索三叶草（Trifolium ambiguum Bieb.）对不同降温模式低温胁迫的响应机制,本试验以缓慢降温（Gradual cooling,GC）和骤然降温（Sudden cooling,SC）处理对其进行了胁迫处理。通过转录组测序,比较了2个低温处理组与常温对照组（CK）间的转录组差异,并筛选了抗寒相关的基因和转录因子。结果表明：与CK相比,GC组和SC组中分别有8 020个和6 289个差异表达基因（DEGs）;GC组的DEGs主要在光合作用、光合作用天线蛋白、淀粉和蔗糖代谢等通路显著富集,SC组的DEGs主要在淀粉和蔗糖代谢、倍半萜和三萜生物合成、类黄酮生物合成等通路显著富集。此外,2处理组共有的DEGs仅在淀粉和蔗糖代谢通路显著富集;一些涉及淀粉和蔗糖代谢与植物激素信号转导通路的基因在2种降温模式下均上调表达,可作为潜在的抗寒基因,如BAM3,BFRUCT1,SUS,GH3.1,SAPK2等。进一步分析发现响应低温胁迫的转录因子主要分布在AP2/ERF,ZFP,MYB等家族。     

沙鞭全长转录组测序及生物信息学分析

为了明晰沙鞭(Psammochloa villosa)的转录组特征，本研究利用PacBio Sequel测序平台首次对其进行全长转录组测序和数据分析，结果共获得323 309个clean reads,自我矫正产生环形一致性序列(CCS)673 540个，预测得到蛋白编码区(CDS)序列28 447个；MISA软件共搜索得到沙鞭93 563个简单重复序列(SSR),分布于56 824条unigene上；转录本基因功能注释，共有166 541条序列得到NR注释，结果显示沙鞭与二穗短柄草(Brachypodium distachyon)亲缘最近；KOG数据库比对将97 892个unigene分为25个功能类别，其中注释较多的功能为翻译后修饰、蛋白质周转和分子伴侣；GO数据库比对共注释到87 930个unigene,分为生物过程、细胞组成和分子功能3大类及62个亚类分支；KEGG结果表明碳水化合物代谢、信号转导、能量代谢通路中注释基因较多。本研究结果丰富了沙鞭的遗传信息，为今后沙鞭关键耐旱基因的挖掘提供了理论依据。     

菠萝蜜响应低温胁迫的转录组分析

通过研究抗寒性不同的菠萝蜜品系,筛选菠萝蜜抗寒关键基因,为菠萝蜜抗寒改良和耐寒良种选育提供理论依据。以自然低温胁迫后广西本土的菠萝蜜抗寒品系、泰国引进的低温敏感型菠萝蜜品系为研究材料,结合田间表型持续观测,并运用转录组测序技术从分子水平上对2个不同抗寒特性菠萝蜜进行生物信息学分析。两个品种共得到30585个差异表达基因,上调17083个,下调13502个。GO功能富集分析显示,差异表达基因显著富集在光合作用的光捕获、光反应、光合作用的调节,对光强度的反应、光系统、光系统I、光系统II,叶绿素、NADH、色素结合功能,以及叶绿素、类黄酮、类胡萝卜素、萜类等次生代谢物的代谢过程。KEGG途径富集分析结果显示,光合作用相关途径、淀粉与蔗糖代谢途径、植物激素信号转导途径及MAPK信号通路为菠萝蜜响应低温胁迫的重要代谢途径。差异表达基因注释到WRKY、ERF、bHLH、NAC、MYB等家族转录因子,与植物应答低温胁迫密切相关。光合作用相关途径、淀粉与蔗糖代谢途径、植物激素信号转导途径及MAPK信号通路是菠萝蜜应答低温胁迫的重要代谢途径,WRKY、ERF、bHLH、NAC、MYB等家族转录因子可能参与菠萝蜜应答低温胁迫。     

紫花苜蓿响应低温胁迫转录因子的鉴定及表达分析

研究证明,转录因子广泛参与植物的生长发育过程并响应多种生物/非生物胁迫信号途径。低温胁迫可导致紫花苜蓿(Medicao sativa)减产、越冬率降低以及生产年限缩短。利用新一代高通量转录组测序技术,本研究对4℃低温胁迫下紫花苜蓿中响应低温胁迫的转录因子基因进行了鉴定,并对其表达情况进行分析,结果表明,转录组测序共获得78 925条Unigene序列和3 448个差异表达基因。其中,从3 448个差异表达基因中共鉴定出43个转录因子家族,共251个基因被显著地诱导表达。不同转录因子家族基因受低温胁迫诱导表达不同。本研究有助于在整体水平加深了解紫花苜蓿转录因子表达特性,为进一步研究这些响应低温胁迫的转录因子的功能提供参考。     

木豆响应低温胁迫差异表达基因分析

为探讨木豆[Cajanus cajan (L.) Millsp.]响应低温胁迫的机制，本研究以1份耐低温(MD)和1份低温敏感(QZ)木豆为试验材料，经4℃低温胁迫处理后进行RNA-seq测序分析，结果表明：2份木豆RNA-seq获得了丰富的转录本信息，对转录因子和差异表达基因的基因功能进行注释。木豆响应响应低温胁迫的转录因子主要包括MYB，AP2-EREBP等；在两份材料中，GO分析发现DEGs注释在细胞组成中“膜”和生物学过程中“细胞过程”，KEGG分析均显著富集在“核糖体”通路中，而MD中还显著富集在植物激素信号转导、昼夜节律-植物等通路。通过两份材料对比，在MD材料中发现特有与耐低温相关的150个DEGs，主要为α-1，4-葡萄糖苷酶，参与碳水化合物、脂质代谢等途径。还筛选37个与耐寒相关重要差异表达基因和14条通路，并对10个DEGs进行qRT-PCR验证。本文从分子水平阐述了木豆低温胁迫下的响应机制，为未来木豆的抗寒育种奠定了理论基础。     

千穗谷TALE转录因子家族的生物信息学分析

TALE转录因子广泛存在于植物中，对植物的生长发育和对外界环境的响应发挥着重要作用。本研究在千穗谷(Amaranthus hypochondriacus)基因组中共鉴定出11个TALE基因。系统发育分析将千穗谷TALE家族蛋白分为BELL和KNOX两个亚家族，同一亚家族中的AhTALEs基因具有相似的基因结构，编码蛋白也具有类似的结构域。蛋白二级结构和三级结构分析都表明AhTALEs蛋白由螺旋-角-螺旋的结构构成。此外，对千穗谷与其他物种的TALE基因共线性关系、组织表达特性和AhTALEs基因启动子区的顺式作用元件分析发现，不同AhTALEs亚家族在各组织中的表达模式不同，并且发现启动子区包含大量非生物胁迫和激素响应元件，这为AhTALE家族的进化分析提供了依据。qRT-PCR分析发现AhTALEs基因在低温和干旱胁迫下的表达水平有不同程度的变化，表明AhTALEs基因能够响应低温和干旱胁迫应答。本研究为开展千穗谷中TALE转录因子在非生物胁迫中的生物学功能研究提供了重要参考。     

高加索三叶草响应不同降温模式低温胁迫的代谢组学分析

为探究高加索三叶草（Trifolium ambiguum Bieb.）对不同降温模式低温胁迫的代谢响应机制，本试验对其分别进行缓慢降温（Gradual cooling，GC）和骤然降温（Sudden cooling，SC）处理，利用广泛靶向代谢组学技术检测叶片代谢物，并比较2个处理组与常温对照组（CK）间的代谢物差异。结果表明：在高加索三叶草叶片中共检测出894种代谢物；与CK相比，在GC和SC处理组中分别筛选出124个和139个差异代谢物（Differential metabolites，DMs）；2处理组的DMs均主要富集在氨基酸代谢和次生代谢产物的生物合成；紫铆因含量仅在GC组中显著增加（P < 0.05），酪氨酸和新橙皮苷含量仅在SC组中显著增加（P < 0.05）；γ-氨基丁酸、亮氨酸、组氨酸、赖氨酸、柠檬酸、矢车菊素和木犀草素含量在2处理组均显著增加（P < 0.05）。这些代谢物含量的增加可能是高加索三叶草适应低温胁迫的主要机制。     

詹纳斯激酶/信号转导蛋白转录激活因子信号通路在动物骨骼肌发育中作用的研究进展

詹纳斯激酶/信号转导蛋白转录激活因子（JAK/STAT）信号通路级联反应已被确定为肌肉发育的关键因素之一,其激活后影响骨骼肌成肌细胞的增殖、分化等。这条通路主要由3个部分组成,即詹纳斯激酶（JAK）、信号转导蛋白转录激活因子（STAT）及詹纳斯激酶受体。文章对JAK/STAT信号转导通路的基本途径及特点,以及它在动物骨骼肌修复生长过程中作用的研究进展进行总结。     

白羊草NAC转录因子家族生物信息学分析

NAC转录因子家族是植物特有的一类重要转录因子,广泛参与植物生长、发育以及器官的建成、逆境胁迫应答等过程。为进一步研究白羊草（Bothriochloa ischaemum） NAC转录因子家族,本研究利用生物信息学方法对白羊草NAC基因所编码的蛋白理化性质、二级结构、保守基序、亚细胞定位、潜在磷酸化位点及系统进化关系进行了分析。结果显示,具有NAC结构域的21条白羊草氨基酸序列共分为10个亚家族,蛋白二级结构以无规则卷曲为主要组成部分,含有5个保守基序,大多数NAC蛋白定位于细胞核,均含有潜在的丝氨酸（Ser）、苏氨酸（Thr）、和酪氨酸（Tyr）磷酸化位点。本研究将为白羊草NAC基因家族的进一步功能分析奠定重要的研究基础。     

基于转录组数据库的高羊茅HD-ZipⅠ转录因子的鉴定及表达模式解析

高等植物特有的同源异型-亮氨酸拉链蛋白(HD-Zip)家族转录因子在调控植物发育与逆境响应中起着重要作用。以生产中广泛应用的冷季型草高羊茅为研究材料,利用其目前已有的转录组序列,在Bioedit软件构建本地数据库。以水稻14个HD-ZIP I蛋白序列作为搜索请求,搜索本地数据库获得高羊茅中对应序列。结合序列分析及PCR克隆技术得到高羊茅13个HD-Zip I类成员的全长序列。系统进化树及MEME结构域分析发现,高羊茅13个HD-Zip I成员与水稻HD-Zip I家族成员进化关系接近,蛋白结构域的数目及排布方式也基本一致。采用荧光定量PCR技术分析短期及长期聚乙二醇6000模拟的干旱处理下11个FaHD-Zip I基因在高羊茅根颈中的mRNA相对表达水平。结果表明,短期干旱胁迫下,除FaHOX4和FaHOX6外,其余基因相对表达水平均发生一定程度的上调或者下调,且FaHOX22的响应程度最大;长期干旱胁迫下,所有基因的表达水平均在7或者14 d上调,FaHOX22在14 d时的表达水平比对照上调76.3倍。为基因组序列缺乏的非模式植物快速有效地进行基于转录组数据库的基因家族克隆及其功能研究提供了一种很好的方法。     

油莎豆CeWRKY转录因子响应非生物胁迫的功能表征

油莎豆亦称油莎草,地上部分是优质牧草,地下块茎可积累大量油脂,抗逆性强、适应性广,是一种新型饲草和油料作物。WRKY是一类数量庞大的转录因子超家族,广泛参与调控植物生长发育、物质代谢和胁迫应答等生理过程。然而,有关于油莎豆WRKY家族成员及其功能知之甚少。本研究基于转录组分析,鉴定出67条油莎豆WRKY转录因子。系统发育树结果显示,油莎豆CeWRKY与拟南芥AtWRKY类似,分为I、Ⅱ、Ⅲ3个亚家族。CeWRKYs均具有典型的WRKY结构域和锌指结构域,其中7条蛋白序列的WRKY结构域存在变异,可能与WRKY成员功能歧化相关。油莎豆CeWRKY蛋白的序列长度、相对分子量、等电点等存在较大差异,均为亲水性蛋白,不含跨膜结构,亚细胞定位预测均位于细胞核。CeWRKY基因在块茎发育时期的表达模式具有时空特异性。不同CeWRKY基因参与不同逆境胁迫响应,其中CeWRKY8基因分别在干旱、盐、碱胁迫处理的油莎豆幼苗中表达谱变化显著,可能是参与油莎豆逆境胁迫响应调控的一个主效转录因子。本研究为深入解析油莎豆CeWRKY转录因子的多样性生物学功能以及介导逆境胁迫应答的分子机制奠定了基础,亦为优质高抗...     

紫花苜蓿bZIP转录因子MsbZIP1的克隆、表达特性和遗传转化分析

bZIP（Basic leucine zipper）是植物中一类非常重要的转录因子,参与植物从生长发育到抗性调控的多个生物学过程。为了了解紫花苜蓿（Medicago sativa）中bZIP转录因子的特性,解析紫花苜蓿MsbZIP1基因的生物学功能,从而阐述MsbZIP1基因响应紫花苜蓿抗逆调控机制,本研究利用RACE技术从紫花苜蓿中获得1个MsbZIP1基因的全长cDNA,该序列全长1 176 bp,编码361个氨基酸,预测分子量为42.3 kD,等电点为6.5。分析发现,该蛋白含有bZIP家族典型的BRLZ碱性结构域和亮氨酸拉链,属于bZIP家族蛋白。进化树分析表明,该蛋白属于bZIP转录因子C亚族,与拟南芥（Arabidopsis thaliana）的AtbZIP63具有很高的同源性,推测可能具有与该类蛋白相似的功能。qRT-PCR分析表明,MsbZIP1基因对干旱、高盐、高温、低温,以及脱落酸（Abscisic acid,ABA）和生长素（Auxin,IAA）处理都有不同程度的响应,推测该基因可能参与调控紫花苜蓿多种非生物胁迫。本试验通过构建表达载体PCAMBIA3301-MsbZIP1,以农杆菌介导的花序浸染法转化拟南芥,经后代筛选、扩繁和分子检测,得到7株超表达的转基因拟南芥。本研究第一次分离了紫花苜蓿C亚族bZIP转录因子,初步确定紫花苜蓿MsbZIP1基因响应多种逆境胁迫的反应,并获得了阳性转基因材料,并为进一步探索该类转录因子在紫花苜蓿抗逆性调控中的作用奠定了基础。     

羊草DREB转录因子的系统发育和功能研究

羊草是兼具经济价值和生态价值的重要牧草,具有耐寒、耐旱和耐盐碱的特点。在植物应对非生物胁迫的过程中,DREB转录因子起到关键作用。然而,目前在GenBank的核酸数据库和EST数据库中羊草仅有2条DREB序列,其中只有1个DREB基因的功能得到实验验证。本研究从羊草转录组测序数据中1次挖掘了26条DREB EST。用羊草EST和拟南芥基因组中DREB基因的蛋白序列构建了系统发育树, LcDREB21位于第4类群,是DREB2类转录因子基因。通过RACE,获得了LcDREB21编码区全长(Genbank登入号:JN860437)。荧光定量PCR结果表明,其表达受干旱和高盐诱导。另外,通过酵母单杂交实验证明了LcDREB21具有转录激活功能;通过表达GFP融合蛋白证明其专一性定位在细胞核中。总之,LcDREB21是一个具有转录激活功能和细胞核专一定位能力的转录因子,可能在植物应对干旱和高盐胁迫的过程中起作用。本实验结果丰富了羊草抗逆基因数据库,为植物改良提供了基因资源。     

鹅p21基因克隆、启动子分析及转录调控研究

【目的】 分析鹅p21基因的结构和启动子活性，探讨p21基因的转录调控机制。【方法】 以泰州鹅为试验对象，通过同源克隆、RACE和生物信息学分析等方法获得鹅p21基因全长序列和5′-侧翼区序列特征；构建6个不同缺失片段的启动子区双荧光素酶报告载体并分析其荧光素酶活性，进而确定p21基因核心启动子区；对核心启动子区转录因子结合位点生肌决定因子(MyoD)(+25~+36 bp)进行定点突变，并构建突变报告基因载体，在C2C12细胞系内初步鉴定鹅p21基因核心转录调控因子。【结果】 鹅p21基因cDNA全长1 943 bp，CDS区大小为453 bp，编码151个氨基酸，蛋白序列包含高度保守的CDI家族结合位点。系统进化树分析表明，鹅p21基因与鸭亲缘关系最近，与鸡和火鸡有较强的进化关系。鹅p21基因5′-侧翼区包含启动子元件，—35~+37 bp是核心启动子区，发挥正向调控作用，结合定点突变技术初步鉴定MyoD是鹅p21基因核心转录调控元件。【结论】 本研究获得了鹅p21基因完整的cDNA序列和启动子区域，MyoD是p21基因核心转录调控因子，为探究p21基因在鹅胚胎期肌肉发育过程中的调控机制提供理论依据。     

基于转录组学比较研究甜高粱幼苗响应干旱和盐胁迫的生理特征北大核心CSCD

探究甜高粱响应干旱与盐胁迫的生理学差异及其分子机制,分析甜高粱在干旱和盐胁迫下的不同基因调控机制和代谢通路,可为饲草作物甜高粱抗逆栽培和育种提供依据。以辽甜1号甜高粱幼苗为材料,使用10%的PEG-6000和0.9%的NaCl模拟中度干旱和盐胁迫,胁迫2和7 d时进行光合气体交换参数、内源激素含量、有机渗透调节物质含量和抗氧化物酶活性测定,同时对幼苗叶片进行转录组测序及生物信息学分析,采用qRT-PCR方法对测序结果进行验证。结果表明:与对照相比,干旱和盐胁迫均显著影响甜高粱幼苗叶片生理指标的变化,但盐胁迫对甜高粱幼苗叶片光合参数和内源激素生长素、细胞分裂素含量抑制程度均高于干旱胁迫,可溶性糖含量在干旱胁迫下显著高于盐胁迫,抗氧化物酶活性和脱落酸含量低于盐胁迫,说明甜高粱对干旱和盐胁迫响应的生理机制不同;利用转录组测序技术鉴定出甜高粱叶片在处理2 d时,干旱和盐胁迫下分别有922和2047个上调基因以及975和1714个下调基因,处理7 d时分别鉴定到157和795个上调基因以及54和722个下调基因;同时鉴定出40个干旱胁迫响应基因和493个盐胁迫响应基因,基于GO富集分析发现甜高粱幼苗在干旱和盐胁迫下响应基因均显著富集于植物逆境响应相关通路,KEGG富集分析发现,干旱胁迫响应基因显著富集于内质网加工和剪切体代谢通路,盐胁迫响应基因显著富集在植物激素信号转导代谢通路;对光合作用、植物激素信号转导、抗氧化物酶和淀粉与蔗糖代谢通路相关差异基因进行分析发现,这些差异基因的表达模式与生理指标变化趋势吻合。因此,甜高粱在转录水平上对盐胁迫的适应稳态落后于干旱胁迫,中度胁迫下甜高粱幼苗对盐胁迫的耐受性要低于干旱胁迫,可溶性糖在甜高粱幼苗抵御干旱胁迫中发挥重要作用,植物激素信号转导和抗氧化物酶活性变化的共同调控是甜高粱幼苗对抗盐胁迫的关键。     

香蕉海藻糖-6-磷酸合成酶基因MaTPS2的克隆及表达特性分析

获得香蕉海藻糖-6-磷酸合成酶基因MaTPS2全长序列并初步探讨MaTPS2基因在植物逆境中的作用。在获得MaTPS2基因全长序列的基础上,采用同源序列比对、聚类分析、实时定量PCR等技术初步探讨香蕉MaTPS2基因的功能。从香蕉A基因组中获得一条3277bp MaTPS2全长序列,包含一个完整的开放阅读框1344bp,编码蛋白含447个氨基酸。其氨基酸理化性质分析表明,MaTPS2蛋白属于不稳定疏水性蛋白,具有一个结构域GT1-TPS。通过与已知植物的TPS氨基酸同源序列比对,同源性达77.95%。系统发育树进化关系分析表明MaTPS2编码的氨基酸序列与银杏TPS氨基酸序列（AAX16014.1）邻近,说明二者的亲缘关系较近。器官特异性分析表明MaTPS2在香蕉植株的根、球茎、假茎、叶、花各器官中均有表达,其中在球茎、假茎和叶中表达量较多,在果实中表达量较低。在ABA、干旱、盐害和枯萎病胁迫处理下MaTPS2表达量升高,在ACC胁迫处理下,MaTPS2表达量显著下降,而低温胁迫条件下,MaTPS2表达不响应。由此推导,MaTPS2可能参与调控香蕉抗旱、抗病机制。激素ACC和ABA胁迫处理说明MaTPS2的表达可能受激素蛋白间接调控。     

紫花苜蓿耐盐性相关NAC转录因子的挖掘及表达分析

NAC(NAM,ATAF1/2,CUC2)是植物特有的一类转录因子家族,可调控植物响应盐胁迫。本研究基于盐胁迫下紫花苜蓿(Medicago sativa L.)耐盐品种GIB和盐敏感品种LS根和叶的转录组数据,通过生物信息学方法筛选出17个差异表达的MsNAC转录因子家族成员,并对其序列特征、系统进化、顺式作用元件及盐胁迫下的表达模式等进行分析。结果显示,17个MsNACs均具有NAC典型的NAM结构域,且蛋白结构相似。系统进化分析将17个MsNACs分成8个亚族,各亚组成员具有高度相似的保守基序。顺式作用元件分析发现,有15个MsNACs具有ABA响应顺式元件ABRE。转录组数据表达模式分析与RT-qPCR验证结果显示MsNAC1,MsNAC2,MsNAC35,MsJUB1和MsNAC40基因在盐处理下GIB叶中上调表达,在LS叶中下调或者未发生变化,推测这些基因参与调控紫花苜蓿的耐盐反应。本研究为进一步研究紫花苜蓿NAC转录因子响应盐胁迫功能提供参考依据。     

家蚕丝腺响应高低温胁迫的转录组分析

丝腺是家蚕合成和分泌丝蛋白的唯一器官，温度在丝腺合成丝蛋白的过程中具有重要调控作用，对蚕丝产量和品质都有重要影响。了解丝腺对饲养温度的响应机制对有效调控熟蚕吐丝有重要借鉴作用，以实用蚕品种贵蚕7号为研究对象，在丝腺发育和丝蛋白合成关键时期(5龄期)分别采用20℃、25℃和30℃不同温度条件进行饲养，并取游走期丝腺进行转录组比较分析。结果发现高温组(30℃)和低温组(20℃)与正常适宜饲养温度组(25℃)相比，共鉴定到109 053个差异表达基因。其中30℃饲养条件下存在67个显著上调表达基因，196个显著下调表达基因；20℃饲养条件下存在101个显著上调表达基因，124个显著下调表达基因。高温组和低温组之间分别存在4个和7个特异差异表达基因。进一步对这些差异基因进行GO和KEGG功能富集分析，发现高温组差异表达基因主要集中于氨基酸代谢、保幼激素代谢和能量代谢等与丝蛋白合成紧密相关的信号通路，低温组差异表达基因主要富集于苯丙素的生物合成、谷胱甘肽代谢等与物质代谢合成相关的信号通路，共同富集的通路主要集中于物质运输和能量代谢。这说明高温和低温胁迫影响丝腺发育或丝蛋白合成分泌的分子途径既存在...     

草地早熟禾NAC基因鉴定及非生物胁迫下表达模式分析

草地早熟禾(Poa pratensis)是世界上广泛应用的草坪草种,其抗逆基因的挖掘对抗逆新品种的选育具有重要意义.NAC基因家族是植物特有的转录因子之一,在植物逆境响应中起到重要作用.以二穗短柄草(Brach ypodium distach yon)的NAC蛋白序列为请求序列,从草地早熟禾转录组数据库序列比对分析得到15个草地早熟禾PpNACs基因的cDNA全长序列,具有完整的开放阅读框.运用生物信息学方法预测早熟禾NAC家族成员的理化性质和保守结构域,发现15个PpNACs转录因子都是酸性蛋白,除PpNAC16外都属于不稳定蛋白,疏水性蛋白,具有相似的结构,通过聚类分析可将15个PpNACs成员分为3类.荧光定量PCR结果显示,在重金属胁迫下显著上调表达的基因有Pp-NAC10、PpNAC18、PpNAC19和PpNAC24,在盐胁迫过程下显著上调表达的基因有PpNAC10和PpNAC24,在干旱胁迫下显著上调表达的基因有PpNAC1、PpNAC6、PpNAC10和PpNAC13,在低温胁迫下显著上调表达的基因有PpNAC1、PpNAC10和PpNAC18,在高温胁迫下显著上调表达的基因有PpNAC20.其中,PpNAC10在重金属、盐、干旱和低温胁迫下均被显著诱导上调表达,Pp-NAC18在重金属、干旱和低温胁迫下均被显著诱导上调表达,PpNAC10和PpNAC18可作为草地早熟禾逆境胁迫响应的候选基因.