过表达ZmIBH1-1提高玉米苗期抗旱性

朱芳芳，董亚辉，任真真，王志勇，苏慧慧，库丽霞，陈彦惠

过表达提高玉米苗期抗旱性

朱芳芳，董亚辉，任真真，王志勇，苏慧慧，库丽霞，陈彦惠

河南农业大学农学院/省部共建小麦玉米作物学国家重点实验室，郑州 450046

【】干旱是严重影响玉米生长发育进程的一个重要因素。挖掘玉米抗旱相关基因，通过转基因功能验证和转录组分析，解析关键基因在响应干旱胁迫过程中的分子调控机制，为抗旱分子育种和遗传改良提供理论依据。以玉米自交系B104（WT）为背景材料，利用农杆菌介导方法构建过表达转基因株系（-OE）；通过对转基因植株进行草铵膦抗性筛选、标记基因和目的基因PCR检测，以及运用实时荧光定量PCR检测目的基因的表达情况，鉴定阳性植株和株系；以WT和-OE转基因株系为材料，通过干旱处理（20% PEG6000），进行表型鉴定和耐旱生理生化指标测定，验证的抗旱功能；通过对干旱胁迫下玉米4叶期转录组的比较分析，鉴定出差异表达的基因（differentially expressed genes，DEGs）；结合DAP-seq（DNA affinity purification sequencing）分析，初步确定ZmIBH1-1蛋白直接调控与抗旱相关的下游靶基因，利用基因组可视化软件IGV（integrative genomics viewer）分析ZmIBH1-1蛋白结合候选靶基因的位置，然后通过Dual-Luciferase试验验证ZmIBH1-1蛋白与靶基因的调控关系。通过玉米遗传转化获得12个转化事件；T3代中，能同时检测到标记基因和目的基因的植株有458个，实时荧光定量PCR检测结果表明，-OE中的表达量显著高于WT，株系3和株系8表达量最高，将其自交获得T4代转基因株系用于后续试验。在干旱胁迫条件下，-OE株系存活率、叶片相对含水量、叶绿素含量、可溶性蛋白含量及其生理生化指标（超氧化物歧化酶、过氧化物酶、过氧化氢酶活性）均显著高于WT，说明玉米中过量表达赋予玉米更高的耐旱性。转录组分析结果表明，WT与-OE株系在干旱胁迫下有1 214个差异表达基因；Gene Ontology（GO）功能富集分析结果表明，差异表达基因主要涉及生物过程、细胞组分和分子功能，如在生物过程中主要涉及到光合作用、应激响应、脱水响应等；KEGG富集分析表明，差异表达基因主要参与植物激素信号传导、新陈代谢等过程。结合转录组显著差异表达基因和DAP-Seq分析所得到ZmIBH1-1蛋白的靶基因，初步确定ZmIBH1-1蛋白直接调控与抗旱相关的11个候选靶基因，包括2个钙信号相关基因、3个半胱氨酸代谢相关基因、1个bHLH转录因子、1个应激响应蛋白、1个谷胱甘肽转移酶、1个氧化还原过程蛋白和2个乙烯响应因子；基因组可视化结果显示ZmIBH1-1蛋白可以结合靶基因启动子区；随后通过Dual-Luciferase试验进一步表明，ZmIBH1-1蛋白可以直接作用于11个候选靶基因，其中，ZmIBH1-1蛋白可以促进、、、、和的表达，抑制、、、和的表达。此外，在干旱胁迫下NAC、WRKY、MYB等转录因子在-OE和WT株系中也存在差异表达。的过表达可以增强玉米苗期的耐旱性；ZmIBH1-1蛋白通过直接调控乙烯信号通路中的和的表达提高玉米的耐旱性；ZmIBH1-1蛋白通过直接调控钙信号相关基因和增强玉米的耐旱性；ZmIBH1-1蛋白可能通过间接调控NAC、WRKY、MYB等转录因子响应干旱胁迫。

玉米；干旱胁迫；；RNA-Seq；转录因子；基因表达

0 引言

【研究意义】玉米是世界第一大粮食作物，其产量占全球总谷物产量的37.2%[1]。玉米产量的提高对保障国家粮食安全至关重要。干旱、盐碱、矿物质缺乏等非生物胁迫对玉米正常的生长发育会造成极大的影响，进而导致玉米产量的损失。其中，干旱是制约玉米产量提高的重要因素之一。因此，挖掘玉米抗旱耐旱的关键基因，从分子水平上揭示玉米抗旱机制，为培育抗旱耐旱玉米新品种提供理论依据，对保障粮食产量稳定具有十分重要的意义。【前人研究进展】为了适应复杂多变的自然环境，在漫长的进化过程中，植物自身形成了一种特异的适应机制，植物在应对干旱等逆境时，会启动相应的功能基因，调节其生理生化水平，以适应不良环境[2-3]。已有的研究表明，植物中bHLH、NAC、WRKY、bZIP、NF-Y、MYB等转录因子家族的基因通过不同的调控路径，在植物抵御逆境胁迫反应中发挥着重要的作用[4-10]。在干旱等胁迫条件下，多个玉米转录因子家族基因的过量表达可以提高植物的抗旱耐旱特性。例如，玉米苗期过表达可以提高玉米水分利用率并诱导干旱响应基因的表达，进而提高玉米干旱耐受性[4]；的表达受干旱、高盐等胁迫条件诱导，过表达可以增强拟南芥干旱耐受性[5]。作为干旱胁迫和高温胁迫的正调控因子，通过调控ABA信号路径相关基因，过表达提高了拟南芥中超氧化物歧化酶（superoxide dismutase，SOD）、过氧化物酶（peroxidase，POD）和过氧化氢酶（catalase，CAT）活性，降低了活性氧（reactive oxygen species，ROS）含量，转基因拟南芥表现出耐旱性和耐热性[6]。过表达提高干旱胁迫下转基因拟南芥中POD和CAT的活性，降低了ROS积累，通过调控胁迫相关基因提高转基因拟南芥的抗旱性[7]。植物bZIP转录因子调节多种功能，包括植物发育、胁迫反应和信号传导等过程，过表达的拟南芥比野生型的电解质渗漏低，水分流失慢，干旱和盐耐受性增强[8]。干旱处理24 h后，地上部的表达量升高5.49倍，推测该基因可能参与了玉米干旱胁迫响应[9]。干旱和盐胁迫诱导的表达，拟南芥过表达该基因呈现生长性能增强，存活率高，CAT、POD和SOD活性升高，增强了转基因拟南芥对干旱和盐胁迫的耐受性[10]。除去转录因子，还有与干旱胁迫响应相关的基因的报道。例如，干旱处理后，过表达转基因玉米成活率显著提高，而敲除突变体对干旱表现出更加敏感，说明增强了玉米干旱耐受性[11]。ZHANG等[12]发现过表达的转基因玉米根毛长度增加，对水分亏缺的耐性提升。DING等[13]发现的遗传变异与玉米苗期抗旱性显著关联。玉米突变体表皮蜡质减少、角质层渗透性增加，降低了幼苗抗旱性[14]。通过诱导胁迫响应基因和提高ROS清除酶活性，可能在干旱和盐胁迫耐受中起作用[15]。WANG等[16]鉴定到一个玉米Ⅰ类SUMO结合酶基因（），该基因在干旱胁迫下表达上调，拟南芥中过表达该基因则增强其抗旱性。在干旱条件下，过表达的转基因株系通过提高SOD和CAT的活性，增加气孔关闭，减少ROS的积累正向调控植物耐旱性[17]。bHLH家族是仅次于MYB家族的第二大家族[18]，在植物抗旱中发挥着重要作用。与编码典型的bHLH转录因子，可以促进转基因拟南芥体内类黄酮生物合成、脱落酸（abscisic acid，ABA）信号通路、脯氨酸生物合成及ROS清除酶等基因表达上调，增强拟南芥耐旱性[19-20]。胡杨基因被干旱诱导后表达上调，主要通过调节气孔密度和大小、光合作用及生长发育等从而提高其抗旱性[21]。茶树在干旱胁迫下有39个bHLH转录因子表达上调[22]。拟南芥是耐旱性、耐盐性等多个逆境胁迫信号的调节因子[23]。在水稻中，高量表达显著提高了水稻的抗旱性[24]。ZmPTF1是一种磷酸饥饿诱导的bHLH转录因子，过表达株系改善玉米根系发育，增加ABA含量，激活ABA-、CBF4-、ATAF2-和NAC30介导的胁迫响应，提高玉米的耐旱性[25]。【本研究切入点】编码一个bHLH型转录因子，在叶夹角形成发育中负向调控玉米叶夹角的大小[26]。然而，关于是否可以提高植物的耐旱性，目前尚未报道。【拟解决的关键问题】本研究以B104（WT）及B104为背景的过表达转基因（-OE）株系为材料，对玉米幼苗进行干旱胁迫处理，结合RNA-Seq和DAP-Seq数据分析，解析响应干旱胁迫的分子机制，为挖掘玉米抗旱基因及选育抗旱型玉米新品种奠定理论基础。

1 材料与方法

试验于2018年6月—2020年10月在河南农业大学农学院/省部共建小麦玉米作物学国家重点实验室完成。

1.1 试验材料

玉米遗传转化所用农杆菌（）菌株EHA105购自河南三瑞生物科技有限公司；pFGC5941植物表达载体购自BioVector NTCC典型培养物保藏中心；过表达重组质粒由河南农业大学省部共建小麦玉米作物学国家重点实验室陈彦惠教授课题组构建，T-DNA区结构如图1所示；玉米遗传转化由北京博美兴奥科技有限公司完成；T1—T4自交纯合及阳性转基因植株鉴定由河南农业大学省部共建小麦玉米作物学国家重点实验室陈彦惠教授课题组完成。过表达转基因株系命名为-OE。

LB：T-DNA的左边重复序列；Tnos：终止子；Bar：除草剂筛选基因；PMAS：MAS启动子；P35S：35S启动子；ETMV：TMV增强子；RB：T-DNA的右边重复序列

1.2 转基因植株的检测

采用除草剂抗性筛选、标记和目的基因PCR鉴定和目的基因的荧光定量PCR表达分析3种方法对转基因植株进行检测鉴定。除草剂抗性筛选：利用WT和-OE为材料，从植株3片叶开始，每隔3—5 d向叶片喷施草铵膦（300 mg·L-1），共喷施3次，处理一周后观察抗性表型，初步筛选转基因阳性植株。标记和目的基因PCR鉴定：以提取的植株叶片DNA为模板，利用标记基因特异引物（Bar-F和Bar-R）和特异引物（ZmIBH-F和ZmIBH-R）进行扩增检测，根据扩增条带鉴定出转基因阳性植株，所用引物序列见表1。目的基因的荧光定量PCR表达分析：用Trizol法提取玉米第4片完全展开叶的RNA，按PrimeScripTMRT reagent Kit（TaKaRa）的操作步骤进行RNA反转录，用定量引物（qIBH-F和qIBH-R）和内参引物（Tublin-F和Tublin-R）进行荧光定量PCR检测，采用2-ΔΔCt法进行定量分析[27]检测阳性植株。每个样品3次生物学重复，3次技术重复，引物序列见表1。将转基因不同世代经过PCR检测到的阳性单株移栽到大田，在开花期进行人工套袋自交获得种子。

表1 试验所用引物

1.3 干旱胁迫下WT和ZmIBH1-1-OE的表型鉴定

干旱胁迫处理：将-OE和WT玉米籽粒置于10%的H2O2中浸泡洗涤20 min，无菌水洗涤3次，每次5 min。消毒后的种子放在无菌湿润的发芽纸中，28℃黑暗培养2 d，挑选发芽均一的籽粒种植蛭石中，在（30±2）℃恒温培养室中培养（12 h黑暗/12 h光照）。玉米幼苗生长至2叶1心转入Hoagland营养液，4叶1心时采用添加20% PEG-6000的Hoagland营养液进行干旱处理，同时设置对照（Hoagland营养液），试验重复3次。

耐旱表型鉴定：选取对照和干旱处理的WT和-OE植株进行存活率统计和相对含水量测定。

相对含水量（relative water content，RWC）的测定：干旱处理24 h后取新鲜叶片擦干净后称重（fresh weight，FM），将叶片放入水中5—6 h，使叶片吸水达到饱和状态后，取出叶片并擦干叶片表面水分后再称重（saturation weight，TM），再将叶片放入烘箱，105℃杀青30 min，然后在80℃环境下烘至恒重，称重（dry weight，DM）。叶片相对含水量RWC（%）=（FM-DM）/（TM-DM）×100。

1.4 生理生化指标的测定

干旱处理24 h后取处理组和对照组的植株叶片，用于生理生化指标的测定，每个指标3次生物学重复。用韩赞平[28]方法进行SOD、POD、CAT活性测定；用叶绿素用丙酮-乙醇混合液的萃取方法[29]进行叶绿素含量（chlorophyll contents，Cht）和类胡萝卜素含量（carotenoids，Car）测定，称取0.2 g新鲜玉米叶片，加入80%预冷的丙酮，研磨成匀浆，于6 000 r/min离心15 min，抽取上清液测定其在665、649和470 nm波长下的吸光值，试验在黑暗环境中进行。计算公式为：

叶绿素a浓度（mg·L-1）Ca＝13.95A665-6.8A649

叶绿素b浓度（mg·L-1）Cb=24.96 A649-7.32A665

叶绿素浓度（mg·L-1）Cht=Ca+Cb

类胡萝卜素浓度（mg·L-1）Car＝（1000A470-2.05Ca- 114.8Cb）/248

采用考马斯亮蓝染色法[30]，对可溶性蛋白（soluble protein，SP）含量进行测定，称取0.1 g新鲜玉米叶片，采用Tris-HCL缓冲液（pH6.8）进行可溶性蛋白研磨提取，通过测定595 nm处吸光值测定可溶性蛋白含量。

1.5 数据分析

应用Microsoft Excel 2010进行数据分析和作图。

1.6 转录组测序及分析

以正常、干旱处理的WT和-OE植株的叶片（2个生物学重复）为材料，委托武汉希望组生物科技有限公司完成RNA提取和文库构建，然后通过MGI-T7平台对文库进行PE150测序，共构建8个cDNA文库。使用在线软件Trim Galore（www.bioinformatics.babraham.ac.uk/projects/trim_galore/）处理测序产出的原始数据，去除含有接头、poly-N及低质量reads后获得高质量clean reads，并对Q20、Q30、GC含量以及重复序列进行统计；将上述高质量的clean reads与玉米基因组B73_V4（AGPV4版本）进行比对，通过FPKM（Fragments per Kilobase per Millon Mapped Fragments）对基因表达量进行标准化，以FPKM＞1为表达标准，随后根据正常与干旱处理的WT和-OE的比较，用DEseq（|log2FC(fold change)|≥1和-value/FDR＜0.05）确定DEGs。利用Uniprot、Swissprot、COG、NR、GO和KEGG等数据库对DEGs进行功能注释，利用在线工具WEGO（http://wego.genomics.org.cn）和Gene Ontology（http://www.geneontology.org/）对差异基因进行GO功能分类[31]。通过KEGG分析DEGs主要参与的代谢途径和信号通路（https://www.kegg.jp/）。

1.7 ZmIBH1-1对其直接作用的下游靶基因的调控分析

为了明确在玉米发育过程中的调控网络，Cao等[26]采用DAP-seq技术证明了ZmIBH1-1蛋白通过绑定4个Cis-elements（NNCAAGTNG、CANGTN、CTTCGNN和GGNGGAGA）直接作用于启动子区域的靶基因有1 188个，结合RNA-Seq数据明确了调控玉米叶夹角发育的靶基因。为了进一步明确响应干旱胁迫所涉及的调控路径，首先用VLOOKUP对干旱胁迫下的RNA-Seq分析获得的差异表达基因和ZmIBH1-1蛋白DAP-seq获得的1 188个靶基因[26]进行交集分析，初步确定干旱胁迫响应的候选靶基因。然后用基因组可视化软件IGV（integrative genomics viewer）分析ZmIBH1-1蛋白结合候选靶基因的位置。

采用Dual-Luciferase试验进一步验证ZmIBH1-1与靶基因的调控关系。将靶基因的启动子（含有ZmIBH1-1的结合位点）克隆到pGreenII0800-luc载体上（含有报告基因）；将基因克隆到pCAMBIA1300载体上。将上述重组载体分别转入农杆菌GV3101中，采用3种组合（融合有靶基因启动子的双荧光素酶报告载体、过表达转录因子载体以及融合有靶基因启动子的双荧光素酶报告载体、双荧光素酶报告载体（空载））侵染烟草（）叶片。培养（14 h光照/10 h黑暗）2 d后，提取烟草叶片蛋白（Cat#E1910，Promega），使用Glomax®20/20生物/化学发光检测仪（Cat#E5311，Promega）测定萤火虫荧光素酶（firefly luciferase，LUC）和海肾荧光素酶（renilla luciferase，REN），每个样品3个生物学重复。

2 结果

2.1 ZmIBH1-1-OE转基因玉米阳性植株的鉴定

通过玉米遗传转化共获得12个转化事件。T3代时，每个转化事件种植5个株系，每个株系种植10个单株，共600个单株，进行鉴定（图2）。结果表明，T3代株系中约80%具有筛选标记基因和目的基因。标记基因和目的基因同时被检测到的植株有458个，其中，来自于2个（3、8）转化事件的3个株系，共30个单株全部为阳性。

对T3代12个独立事件的转基因株系在4叶期进行qRT-PCR检测，结果表明，-OE株系中的表达量显著高于WT，且在3、8转化事件中该基因的表达量最高（图3）。结合抗除草剂鉴定和PCR检测结果，说明这两个转化事件已经纯合。转化事件3和8套袋自交获得的T4代转基因株系用于后续试验研究。

A：喷洒或涂抹除草剂草胺膦的表型鉴定；B：Bar引物PCR检测；C：IBH1-1引物PCR检测。WT：野生型单株；T：转基因单株；1：阳性对照（质粒DNA）；2：空白对照；M：DL5000 DNA Marker；其余泳道代表转基因植株

WT：野生型；1—12：12个转化事件的ZmIBH1-1-OE株系（n=3，±SD，**P＜0.01）

2.2 过表达ZmIBH1-1增强玉米耐旱性

-OE和WT植株模拟干旱处理后发现，-OE株系比WT植株抗旱性强（图4-A）。具体表现为：干旱胁迫下，-OE株系的存活率显著高于WT（图4-B），-OE株系的叶片相对含水量显著高于WT（图4-C），-OE株系的SOD、CAT及POD酶活性均显著高于WT（图4-D），类胡萝卜素、总叶绿素、可溶性蛋白含量也均显著高于WT（图4-E）。以上结果说明过表达增强玉米的耐旱性。

2.3 干旱胁迫下ZmIBH1-1-OE和WT叶片转录组分析

为明确-OE植株在干旱胁迫下转录水平变化，对正常和干旱胁迫下WT和-OE株系的幼苗进行转录组测序。去除低质量reads后，8个样品均有超过91%的clean reads。将clean reads与玉米B73基因组参考序列V4进行比对，8个样品均有超过88%的unique reads（表2）；PCA分析结果显示不同样品间存在差异，同一样品不同生物学重复聚集（图5-A）；相关性分析表明，生物学重复间相关系数高（图5-B）。以上结果表明转录组数据稳定可靠，可用于后续分析。

A：B104和ZmIBH1-1-OE株系在PEG6000胁迫下的表型；B：B104和ZmIBH1-1-OE株系在PEG6000胁迫下的存活率（n=3，±SD，**P＜0.01，n.s不显著）；C：B104和ZmIBH1-1-OE株系在PEG6000胁迫下的叶片平均含水量（n=3，±SD，**P＜0.01，n.s不显著）；D：B104和ZmIBH1-1-OE株系在PEG6000胁迫下的POD、SOD、CAT酶活性测定（n=3，±SD，**P＜0.01，n.s不显著）；E：B104和ZmIBH1-1-OE株系在PEG6000胁迫下的总叶绿素（Cht）、类胡萝卜素（Car）和可溶性蛋白（SP）含量的变化测定（n=3，±SD，**P＜0.01，n.s不显著）

表2 RNA-seq数据reads数总结

转录组分析结果表明，干旱胁迫下的表达量显著高于正常条件，且在-OE株系中表达量显著高于WT（图5-C）；WT与-OE株系在干旱胁迫下有1 214个基因表达差异显著。将1 214个差异表达基因进行GO注释分析（电子附图1-A），这些基因主要涉及到生物过程、细胞组分和分子功能。KEGG代谢通路分析结果表明，差异表达基因主要参与植物激素信号传导、新陈代谢等路径（电子附图1-C）。1 214个差异表达基因中还包括Ca+通道蛋白及NAC、WRKY、MYB等类型转录因子（电子附图1-B和电子附表1）。

A：转录组数据PCA分析。p代表PEG处理，0 h代表未处理，1和2代表2个生物学重复；B：转录组相关性分析热图；C：正常和干旱条件下，B104和ZmIBH1-1-OE株系中ZmIBH1-1的表达量（n=2，±SD，**p＜0.01）

2.4 ZmIBH1-1蛋白直接作用的干旱胁迫响应的靶基因分析

为进一步解析ZmIBH1-1蛋白直接作用的干旱响应靶基因，对RNA-Seq差异表达基因和DAP-Seq分析得到的ZmIBH1-1蛋白的靶基因[26]进行分析，初步确定11个ZmIBH1-1可能直接调控的与抗旱相关的候选靶基因（电子附表2），包括2个钙信号相关基因（：Zm00001d025340；：Zm00001d051676）、3个半胱氨酸代谢相关基因（：Zm00001d034736；：Zm00001d049110；：Zm00001d038173）、1个bHLH转录因子（：Zm00001d011847）、1个应激响应蛋白（：Zm00001d047302）、1个谷胱甘肽转移酶（：Zm00001d036951）、1个氧化还原过程蛋白（：Zm00001d029519）和2个乙烯响应因子（：Zm00001d017466；：Zm00001d028974）。ZmIBH1-1的DAP-Seq结果显示ZmIBH1-1在靶基因启动子区域存在peaks富集（图6），说明ZmIBH1-1可以直接结合靶基因启动子。

A—D分别代表ZmIBH1-1分别结合ZmEIN3、ZmERF-107、ZmCLPB3和ZmAGD12的启动子区域。红框代表ZmIBH1-1结合峰所在位置

随后采用Dual-Luciferase试验进一步验证ZmIBH1-1对候选靶基因的调控作用（图7-A）。结果表明，ZmIBH1-1分别显著增强由、、、、和的启动子启动的报告基因的表达（图7-B），而显著降低由、、、和的启动子启动的报告基因的表达（图7-C）；RNA-Seq结果表明，在干旱胁迫下，、、、、和在-OE株系中的表达量显著高于WT，而、、、和在-OE株系中的表达量显著低于WT。以上结果表明ZmIBH1-1可以直接调控靶基因的表达。

A：35s启动子启动海肾荧光素酶（REN）作为内参，基因启动子启动报告基因萤火虫荧光素酶（LUC）基因，ZmIBH1-1作为效应因子，空载体为对照，在本氏烟草中瞬时表达。B—C：LUC/REN的比值代表启动子的活性，每个试验重复3次（n=3，±SD，**p＜0.01）。NC：阴性对照

3 讨论

本研究结果表明，过表达可以增强玉米的耐旱性，ZmIBH1-1蛋白直接结合11个靶基因的启动子调控靶基因的表达，其中包括2个乙烯响应基因和，因此推测可能通过乙烯信号通路来响应干旱胁迫。植物在受到干旱胁迫时植物体内的内源乙烯含量会有很大程度的提升[32]，过量的乙烯会诱发植物细胞的衰老和程序性死亡[33-34]。在玫瑰中促进表达，而可以显著抑制乙烯受体蛋白基因的表达[35]，ETRs结合乙烯，减少或抑制乙烯信号转导诱发的细胞死亡[36]。ZmEIN3作为EIN3-like转录因子，是乙烯信号通路中的重要转录因子[37]。本研究结果表明，在干旱胁迫下，的表达被显著抑制。说明可能通过直接抑制的表达，间接促进玉米中乙烯受体蛋白基因的表达，从而减少或抑制乙烯诱发的细胞程序性死亡，提高玉米的抗旱性。在逆境胁迫下，植物通过转录因子调控相关基因的表达提高植物体抵御逆境胁迫的能力[4-10]。ERF转录因子在响应逆境胁迫中起到重要作用[38-39]。ERF家族成员通过特异结合基因启动子GCC-box，从而参与乙烯应答及非生物胁迫[38]。水稻中，过表达ERF转录因子基因（）和可以调节内源乙烯等激素合成的功能，以抵抗淹水胁迫[39]。本研究中ZmIBH1-1蛋白通过抑制乙烯响应因子基因的表达参与玉米对干旱胁迫响应的调节。

钙离子作为第二信使，在植物应对非生物胁迫时（盐、干旱、低温等）发挥着重要作用[40-41]。干旱缺水引起细胞质内钙离子浓度变化从而激活钙依赖性蛋白激酶（calmodulin-dependent protein kinases，CPKs）信号导致ABA释放，ABA浓度的积累导致气孔关闭减少水分损失[42-43]。本研究中，在干旱胁迫下，-OE植株的RWC显著高于WT。RNA-Seq分析发现多个与钙信号相关的基因在WT和-OE株系中存在差异表达，其中，ZmIBH1-1直接作用于和的启动子调控其表达，间接调控钙依赖性蛋白激酶Zm00001d014773以及钙调素蛋白Zm00001d040323和Zm00001d028948。说明ZmIBH1-1可能通过调控钙信号相关基因的表达使气孔关闭减少水分蒸腾，从而提高玉米- OE植株的抗旱性。

已有的研究表明，植物中bHLH、NAC、WRKY、bZIP、NF-Y、MYB等转录因子家族的基因通过不同的调控路径，在植物抵御逆境胁迫反应中发挥着重要的作用[4-10]。MAO等[4]克隆了基因，该基因过表达明显提高了玉米苗期的抗旱性和水分利用效率，并诱导应答抗旱基因在干旱胁迫下表达上调。Ren等[44]对NAC型转录因子ZmNST3展开了研究，敲除导致玉米的耐旱耐盐能力降低。MYB蛋白是植物抵御各种逆境的重要转录因子。如、A、、、、等都能通过不同的方式和路径提高植物的抗旱性[45-49]。干旱和盐胁迫可以诱导的表达，拟南芥过表达该基因，CAT、POD和SOD活性升高，增强了转基因拟南芥对干旱和盐胁迫的耐受性[10]。植物特有转录因子WRKY在植物非生物胁迫响应过程中发挥重要作用[50]。在拟南芥中过表达、，通过提高SOD、POD和CAT活性降低ROS含量，进而提高转基因拟南芥的耐旱性和耐热性[6-7]。本研究中，转录组分析表明，在干旱胁迫下，WT与-OE株系的差异表达基因中包括8个NAC、8个WRKY、11个MYB类型转录因子（电子附图1-B和电子附表1）。说明ZmIBH1-1可能通过间接调控NAC、MYB和WRKY转录因子参与干旱胁迫响应。

4 结论

过表达可以增强玉米的耐旱性；ZmIBH1-1可以直接结合11个靶基因的启动子区域，调控11个靶基因的表达；ZmIBH1-1通过直接调控乙烯信号通路中的基因和的表达提高玉米的耐旱性；ZmIBH1-1通过直接调控钙信号相关基因和增强玉米的耐旱性；ZmIBH1-1可能通过间接调控NAC、WRKY、MYB等转录因子响应干旱胁迫。

[1] QU D Y. Food and Agriculture Organization of the United Nations. Viale delle Terme di Caracalla, Rome, Italy. 2016 [June 2021], http://www.fao.org/statistics/databases/en/.

[2] 代宇佳, 罗晓峰, 周文冠, 陈锋, 帅海威, 杨文钰, 舒凯. 生物和非生物逆境胁迫下的植物系统信号. 植物学报, 2019, 54(2): 102-111.

Dai Y J, Luo X F, Zhou W G, Chen F, Shuai H W, Yang W Y, Shu K. Plant systemic signaling under biotic and abiotic stresses conditions. Chinese Bulletin of Botany, 2019, 54(2): 102-111. (in Chinese)

[3] Bartels D, Sunkar R. Drought and salt tolerance in plants. Critical Reviews in Plant Sciences, 2005, 24(1): 23-58.

[4] Mao H, Wang H, Liu S, Li Z, Yang X, Yan J, Li J, Tran L, Qin F. A transposable element in a NAC gene is associated with drought tolerance in maize seedlings. Nature Communications, 2015, 6: 8326.

[5] Mao H, Yu L, Han R, Li Z, Liu H. ZmNAC55, a maize stress-responsive NAC transcription factor, confers drought resistance in transgenic. Plant Physiology & Biochemistry, 2016, 105: 55-66.

[6] Wang C T, Ru J N, Liu Y W, Li M, Zhao D, Yang J F, Fu J D, Xu Z S. Maize WRKY transcription factor ZmWRKY106 confers drought and heat tolerance in transgenic plants. International Journal of Molecular Sciences, 2018, 19(10): 3046.

[7] Wang C T, Ru J N, Liu Y W, Yang J F, Meng L, Xu Z S, Fu J D. The maize WRKY transcription factor ZmWRKY40 confers drought resistance in transgenic. International Journal of Molecular Sciences, 2018, 19(9): 2580.

[8] Ying S, Zhang D F, Jing F, Shi Y S, Song Y C, Wang T Y, Yu L. Cloning and characterization of a maize bZIP transcription factor, ZmbZIP72, confers drought and salt tolerance in transgenic. Planta, 2012, 235(2): 253-266.

[9] 刘彦丹, 英生, 张登峰, 石云素, 宋燕春, 白志川, 王天宇, 黎裕. 玉米逆境胁迫响应基因ZmbZIP71的克隆与表达分析. 植物遗传资源学报, 2011, 12(5): 775-781.

Liu Y D, Ying S, Zhang D F, Shi Y S, Song Y C, Bai Z C, Wang T Y, Li Y. Isolation and expression analysis of a stress-responsive gene ZmbZIP71 in maize (L.). Journal of Plant Genetic Resources, 2011, 12(5): 775-781. (in Chinese)

[10] Wu J, Jiang Y, Liang Y, Chen L, Chen W, Cheng B. Expression of the maize MYB transcription factor ZmMYB3R enhances drought and salt stress tolerance in transgenic plants. Plant Physiology and Biochemistry, 2019, 137: 179-188.

[11] Zhang H, Xiang Y, He N, Liu X, Dai M. Enhanced vitamin C production mediated by an ABA-induced PTP-like nucleotidase improves drought tolerance ofand maize. Molecular Plant, 2020, 13(5): 760-776.

[12] Zhang X, Mi Y, Mao H, Liu S, Qin F. Genetic variation in ZmTIP1 contributes to root hair elongation and drought tolerance in maize. Plant Biotechnology Journal, 2020, 18(5): 1271-1283.

[13] Ding S, He F, Tang W, Du H, Wang H. Identification of maize CC-type glutaredoxins that are associated with response to drought stress. Genes, 2019, 10(8): 610.

[14] Li L, Du Y, He C, Dietrich C R, Zheng J. Maize glossy6 is involved in cuticular wax deposition and drought tolerance. Journal of Experimental Botany, 2019, 70(12): 3089-3099.

[15] Zhou L, Zhou J, Xiong Y, Liu C, Wang J, Wang G, Cai Y, Wu K. Overexpression of a maize plasma membrane intrinsic protein ZmPIP1;1 confers drought and salt tolerance in. Plos One, 2018, 13(6): e198639.

[16] Wang H, Wang M, Xia Z. The maize class-I SUMO conjugating enzyme ZmSCE1d is involved in drought stress response. International Journal of Molecular Sciences, 2019, 21(1): 29.

[17] Liang Y, Jiang Y, Du M, Li B, Wu J. ZmASR3 from the maize ASR gene family positively regulates drought tolerance in transgenic. International Journal of Molecular Sciences, 2019, 20(9): 2278.

[18] Feller A, Machemer K, Braun E L, Grotewold E. Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant Journal for Cell & Molecular Biology, 2011, 66(1): 94-116.

[19] Wang F, Zhu H, Kong W, Peng R, Liu Q, Yao Q. The antirrhinum AmDEL gene enhances flavonoids accumulation and salt and drought tolerance in transgenic. Planta, 2016, 244(1): 59-73.

[20] Wang F, Zhu H, Chen D, Li Z, Peng R, Yao Q. A grape bHLH transcription factor gene, VvbHLH1, increases the accumulation of flavonoids and enhances salt and drought tolerance in transgenic. Plant Cell Tissue & Organ Culture, 2016, 125(2): 387-398.

[21] DONG Y, WANG C, HAN X, TANG S, LIU S, XIA X, YIN W. A novel bHLH transcription factor PebHLH35 fromconfers drought tolerance through regulating stomatal development, photosynthesis and growth in. Biochemical and Biophysical Research Communications, 2014, 450(1): 453-458.

[22] Cui X, Wang Y X, Liu Z W, Wang W L, Li H, Zhuang J. Transcriptome-wide identification and expression profile analysis of the bHLH family genes in. Functional and Integrative Genomics, 2018, 18(5): 489-503.

[23] Liu W, Tai H, Li S, Gao W, Zhao M, Xie C, Li W X. bHLH122is important for drought and osmotic stress resistance inand in the repression of ABA catabolism. New Phytologist, 2014, 201(4): 1192-1204.

[24] Seo J S, Joo J, Kim M J, Kim Y K, Nahm B H, Sang I S, Cheong J J, Lee J S, Kim J K, Yang D C. OsbHLH148, a basic helix-loop-helix protein, interacts with OsJAZ proteins in a jasmonate signaling pathway leading to drought tolerance in rice. Plant Journal for Cell & Molecular Biology, 2011, 65(6): 907-921.

[25] Li Z, Liu C, Zhang Y, Wang B, Ran Q, Zhang J. The bHLH family member ZmPTF1 regulates drought tolerance in maize by promoting root development and ABA synthesis. Journal of Experimental Botany, 2019, 70(19): 5471-5486.

[26] Cao Y, Zeng H, Ku L X, Ren Z, Han Y, Su H, Dou D, Liu H, Dong Y, Zhu F. ZmIBH1-1 regulates plant architecture in maize. Journal of Experimental Botany, 2020, 71(10): 2943-2955.

[27] Kenneth J L, Thomas D S. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCTmethod. Methods, 2002, 25: 402-408.

[28] 韩赞平. 玉米种子活力相关性状QTL定位及相关基因的克隆[D]. 郑州: 河南农业大学, 2014.

Han Z P. QTL mapping of seed vigor related traits and related gene cloning in maize [D]. Zhengzhou: Henan Agricultural University, 2014. (in Chinese)

[29] 姚晓云, 蓝海军, 邓伟, 陈红萍, 罗晨曦, 况震, 罗宗铭, 王记林, 陈大洲. 水稻淡白叶突变体的叶绿素含量测定及农艺性状比较分析. 江西农业学报, 2020, 32(12): 12-15.

Yao X Y, Lan H J, Deng W, Chen H P, Luo C X, Kuang Z, Luo Z M, Wang L J, Chen D Z. Determination of chlorophyll content and comparative analysis of agronomic traits of pale-white- leaf mutant in rice. Acta Agriculturae Jiangxi, 2020, 32(12): 12-15. (in Chinese)

[30] 焦洁. 考马斯亮蓝G-250染色法测定苜蓿中可溶性蛋白含量. 农业工程技术, 2016, 36(17): 33-34.

Jiao J. Determination of soluble protein content in Alfalfa by Coomassie brilliant blue G-250 staining. Agricultural Engineering Technology, 2016, 36(17): 33-34. (in Chinese)

[31] Ye J. WEGO: a web tool for plotting GO annotations. Nucleic Acids Research, 2006, 34(Web Server issue): W293-W297.

[32] Pierik R, Sasidharan R, Voesenek L. Growth control by ethylene: Adjusting phenotypes to the environment. Journal of Plant Growth Regulation, 2007, 26(2): 188-200.

[33] Lanahan M B. The never ripe mutation blocks ethylene perception in tomato. The Plant Cell, 1994, 6(4): 521-530.

[34] Grbi V, Bleecker A B. Ethylene regulates the timing of leaf senescence in. The Plant Journal, 1995, 8(4): 595-602.

[35] Luo J, Ma N, Pei H, Chen J, Li J, Gao J. A DELLA gene, RhGAI1, is a direct target of EIN3 and mediates ethylene-regulated rose petal cell expansion via repressing the expression of RhCesA2. Journal of Experimental Botany, 2013, 64(16): 5075-5084.

[36] Hua J, Meyerowitz E M. Ethylene responses are negatively regulated by a receptor gene family in. Cell (Cambridge), 1998, 94(2): 261-271.

[37] Guo H W. Paradigms and paradox in the ethylene signaling pathway and interaction network. Molecular Plant, 2011, 4(4): 626-634.

[38] OHME-TAKAGI M, SHINSHI H. Ethylene-inducible DNA binding proteins that interact with an ethylene-responsive element. The Plant Cell, 1995, 7(2): 173-182.

[39] HATTORI Y, NAGAI K, FURUKAWA S, SONG XJ, KAWANO R, SAKAKIBARA H, WU J, MATSUMOTO T, YOSHIMURA A, KITANO H, MATSUOKA M, MORI H, ASHIKARI M. The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature, 2009, 460(7258): 1026-1030.

[40] Knight H. Calcium signaling during abiotic stress in plants. International Review of Cytology-a Survey of Cell Biology, 1999, 195: 269-324.

[41] Reddy A, Ali G S, Celesnik H, Day I S. Coping with stresses: Roles of calcium- and calcium/calmodulin-regulated gene expression. The Plant Cell, 2011, 23(6): 2010-2032.

[42] Dubrovina A S, Kiselev K V, Khristenko V S, Aleynova O A. VaCPK20, a calcium-dependent protein kinase gene of wild grapevineRupr., mediates cold and drought stress tolerance. Journal of Plant Physiology, 2015, 185: 1-12.

[43] Zou J J, Li X D, Ratnasekera D, Wang C, Liu W X, Song L F, Zhang W Z, Wu W H.CALCIUM-DEPENDENT PROTEIN KINASE8 and CATALASE3 function in abscisic acid-mediated signaling and H2O2homeostasis in stomatal guard cells under drought stress. The Plant Cell, 2015, 27(5): 1445-1460.

[44] Ren Z, Zhang D, Cao L, Zhang W, Ku L. Functions and regulatory framework of ZmNST3 in maize under lodging and drought stress. Plant Cell and Environment, 2020, 43(9): 2272-2286.

[45] Oh JE, Kwon Y, Kim JH, Noh H, Hong SW, Lee H. A dual role for MYB60 in stomatal regulation and root growth ofunder drought stress. Plant Molecular Biology, 2011, 77(1/2): 91-103.

[46] Nakabayashi R, Yonekura-Sakakibara K, Urano K, Suzuki M, Yamada Y, Nishizawa T, Matsuda F, Kojima M, Sakakibara H, Shinozaki K. Enhancement of oxidative and drought tolerance inby overaccumulation of antioxidant flavonoids. The Plant Journal, 2014, 77(3): 367-379.

[47] Dai X Y, Xu Y Y, Ma Q B, Xu W Y, Wang T, Xue Y B, Chong K. Overexpression of an R1R2R3 MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic. Plant Physiology, 2007, 143(4): 1739-1751.

[48] Hong S, Chen S, Jiang J, Chen F, Chen Y, Gu C, Li P, Song A, Zhu X, Gao H. Heterologous expression of the chrysanthemum R2R3-MYB Transcription Factor CmMYB2 enhances drought and salinity tolerance, increases hypersensitivity to ABA and delays flowering in. Molecular Biotechnology, 2012, 51(2): 160-173.

[49] Qin Y, Wang M, Tian Y, He W, Lu H, Xia G. Over-expression of TaMYB33 encoding a novel wheat MYB transcription factor increases salt and drought tolerance in. Molecular Biology Reports, 2012, 39(6): 7183-7192.

[50] Finatto T, Viana V E, Woyann L G, Busanello C, Oliveira A. Can WRKY transcription factors help plants to overcome environmental challenges? Genetics & Molecular Biology, 2018, 41(3): 533-544.

附表1 正常和干旱条件下，B104和ZmIBH1-1过表达株系中胁迫和光合作用相关的差异表达基因Supplementary table 1 Differentially expressed genes related to stressand photosynthesisin B104and ZmIBH1-1overexpression linesundernormaland droughtstress

?

?

?

?

?

附表2 干旱条件下，ZmIBH1-1直接调控的差异表达基因Supplementary table 2 DEGsdirectly regulated by ZmIBH1-1underdroughtstress

Over-expression ofto improve drought resistance in maize seedlings

ZHU FangFang, DONG YaHui, REN ZhenZhen, WANG ZhiYong, SU HuiHui, KU LiXia, CHEN YanHui

College of Agronomy, Henan Agricultural University/National Key Laboratory of Wheat and Maize Crop Science, Zhengzhou 450046

【】Drought is an important factor that affects the growth and development of maize seriously. Through the mining of genes related to drought resistance in maize, transgene function verification and transcriptome analysis, analyzing the molecular regulation mechanism of key genes in response to drought stress, this paper provides the theoretical basis for drought resistance molecular breeding and genetic improvement. 【】In this study, the maize inbred line B104 (wild-type, WT) was used as the background to construct theoverexpression (-OE) transgenic line by-mediated method. The transgenic plants and lines were identified via screening transgenic plants for glufosinate-ammonium resistance, PCR detection of marker gene and target gene and expression analysis of target gene by qRT-PCR. We used the-OE and WT transgenic lines as materials. Through drought treatment (20% PEG6000), phenotype identification and drought tolerance physiological and biochemical index determination were carried out to verify the drought resistance function of; RNA-Seq was used to identify differentially expressed genes (DEGs) under drought stress at the 4-leaf stage; Combined with DAP-seq (DNA affinity purification sequencing) analysis, it is preliminarily determined that ZmIBH1-1 protein directly regulates downstream target genes related to drought resistance, and IGV (Integrative Genomics Viewer) was used to analyze the position of the ZmIBH1-1 protein binding candidate target gene, and then the Dual-Luciferase assay was used to verify the regulatory relationship between ZmIBH1-1 protein and target genes. 【】12 transformation events were obtained by genetic transformation of maize. In the T3generation, there were 458 plants in which the marker geneand the target genewere simultaneously detected. The results of qRT-PCR showed that the expression level of-OE lines was significantly higher than that of WT and the expression levels of transformation events 3 and 8 were the highest, which were self-crossed to obtain T4generation for subsequent experiments. Under drought stress, the survival rate, the relative water content, the chlorophyll content, soluble protein content and the physiological and biochemical indicators (superoxide dismutase, peroxidase, catalase activity) of-OE were higher than those of WT significantly, which indicating that the overexpression of-OE lines. Gene Ontology (GO) analysis showed that DEGs were mainly involved in biological processes, cell components and molecular functions, such as photosynthesis, stress response, dehydration response, etc. in biological processes; KEGG enrichment analysis showed that DEGs were mainly involved in the signal transduction of plant hormones, the metabolism and other processes. Combining the significantly DEGs of RNA-Seq and the target genes of ZmIBH1-1 obtained from DAP-seq analysis, it is preliminarily identified 11 candidate target genes related to drought resistance that may be directly regulated by ZmIBH1-1, including 2 calcium signal related genes, 3 cysteine metabolism related genes, 1 bHLH transcription factor, 1 stress response protein, 1 glutathione transferase, 1 redox process protein and 2 ethylene response factor; Integrative genomics viewer showed that ZmIBH1-1 protein could bind to the promoters of the target genes; Subsequent Dual-Luciferase assay further showed that ZmIBH1-1 protein can directly act on 11 candidate target genes, of which, ZmIBH1-1 directly binds to the promoters of,,,,andto promote their expression, and directly binds to the promoters of,,,andto repress their expression. In addition, transcription factors such as NAC, WRKY and MYB also differentially expressed between WT and-OE under drought stress. 【】The overexpression ofcan enhance the drought tolerance of maize; ZmIBH1-1 improves the drought tolerance of maize by directly regulating the expression of genesandand; ZmIBH1-1 may indirectly regulate NAC, WRKY, MYB and other transcription factors in response to drought stress.

maize; drought stress;; RNA-Seq; transcription factor; gene expression

10.3864/j.issn.0578-1752.2021.21.002

2021-04-25；

2021-06-16

国家自然科学基金（31871639）

朱芳芳，E-mail：1261546452@qq.com。董亚辉，E-mail：1018661461@qq.com。朱芳芳和董亚辉为同等贡献作者。通信作者陈彦惠，E-mail：chy9890@163.com

（责任编辑 李莉）