芹菜ICE1基因及bZIP转录因子结构与功能：解析与展望

芹菜ICE1基因及bZIP转录因子结构与功能：解析与展望一、引言1.1研究背景与意义芹菜（ApiumgraveolensL.）作为伞形科芹属的一、二年生草本植物，在全球蔬菜生产中占据重要地位。其富含维生素、胡萝卜素、纤维素以及挥发性芳香成分，不仅是餐桌上的常见蔬菜，还具有一定的药用价值，深受消费者喜爱。中国作为芹菜种植大国，芹菜的种植面积广泛，总产量居世界首位。近年来，随着人们生活水平的提高，对芹菜的品质和产量要求也越来越高。然而，在实际生产过程中，芹菜常面临各种非生物胁迫，如低温、高温、干旱和盐渍等，这些胁迫严重影响芹菜的生长发育、产量和品质。植物在长期进化过程中，形成了一套复杂而精细的抗逆机制，以应对各种逆境胁迫。转录因子在植物抗逆调控网络中发挥着关键作用，它们能够与特定的顺式作用元件结合，调控下游抗逆相关基因的表达，从而增强植物的抗逆性。ICE1（InducerofCBFExpression1）基因是植物抗寒信号通路中的关键转录因子，它能够直接激活CBF（C-repeatBindingFactor）基因的表达，进而调控一系列冷响应基因（COR，Cold-RegulatedGenes）的表达，提高植物的抗寒性。此外，ICE1还参与调控植物对其他非生物胁迫的响应，如干旱和盐胁迫等，在植物抗逆过程中具有重要作用。bZIP（BasicLeucineZipper）转录因子家族是植物转录因子中较大的家族之一，其成员含有保守的碱性亮氨酸拉链结构域，能够与DNA特定序列结合，调控基因表达。bZIP转录因子参与植物生长发育、激素信号转导、光合作用以及逆境胁迫响应等多个生物学过程。在逆境胁迫下，bZIP转录因子通过与不同的顺式作用元件相互作用，激活或抑制下游抗逆相关基因的表达，从而调节植物对逆境的适应能力。研究表明，bZIP转录因子在植物对低温、高温、干旱和盐渍等非生物胁迫的响应中发挥着重要作用，但其在芹菜抗逆过程中的具体功能和作用机制尚不清楚。深入研究芹菜ICE1基因及bZIP转录因子的结构与功能，对于揭示芹菜抗逆分子机制具有重要意义。通过解析ICE1基因和bZIP转录因子在芹菜抗逆信号通路中的作用方式和调控网络，可以为芹菜抗逆育种提供坚实的理论基础。这不仅有助于深入理解植物抗逆的分子生物学过程，丰富植物抗逆理论，还能为其他植物的抗逆研究提供借鉴和参考。从应用角度来看，挖掘芹菜抗逆关键基因，利用现代生物技术手段将这些基因导入芹菜品种中，有望培育出具有更强抗逆性的芹菜新品种。这将显著提高芹菜在逆境条件下的产量和品质，减少因逆境胁迫造成的经济损失，保障芹菜产业的稳定发展。同时，对于优化农业种植结构、提高土地利用率以及保障蔬菜供应的稳定性和安全性也具有重要意义。1.2国内外研究现状在植物基因研究领域，对芹菜基因的探索是一个重要方向。国内外学者已在芹菜基因克隆与功能验证方面取得一定成果。在芹菜黄酮合成相关基因研究中，陈逸云等人成功从六合黄心芹、津南实芹和美国西芹中克隆出黄酮合成酶I基因（AgFSI），并对其基因结构、编码氨基酸序列以及在不同芹菜品种间的差异进行了详细分析，发现该基因在不同品种芹菜中具有高度保守性，为芹菜黄酮类化合物合成机制的研究奠定了基础。在芹菜抗逆基因研究方面，已有研究从芹菜品种‘津南实芹’中克隆出多个与抗逆相关的基因。如四川农业大学的研究团队克隆出耐热基因AgDREBA6a和AgDREBA6c，实验表明，将这两个基因转入拟南芥后，转基因植株根系生长能力、光合作用以及耐热能力明显提高，证实了这些基因在植物耐高温胁迫中发挥重要调控作用，为培育抗高温芹菜品种提供了基因资源和理论依据。南京农业大学的熊爱生教授课题组对芹菜R2R3-MYB家族转录因子进行分析，鉴定获得芹菜花青苷合成相关转录因子AgMYB1，瞬时表达和过表达该基因的实验显示，烟草叶片和转基因拟南芥植株的花青苷含量、花青苷合成途径基因的表达量明显升高，抗氧化活性增强，为芹菜花青苷合成分子机制的研究提供了重要参考。转录因子在植物生长发育和逆境响应中起关键调控作用，对芹菜转录因子的研究也逐步展开。目前，已基本完成对芹菜AP2/ERF转录因子家族的全基因组鉴定。研究发现，AP2/ERF转录因子家族广泛参与芹菜的生长发育、果实成熟以及对生物和非生物胁迫的响应过程。在非生物胁迫响应方面，部分DREB类转录因子能够特异性地结合DRE顺式元件，调控多个与干旱、高盐及低温等非生物胁迫有关功能基因的表达。在对芹菜bZIP转录因子的研究中，虽然取得了一定进展，但与其他模式植物相比，研究的广度和深度仍有待提高。现有研究仅初步鉴定出一些bZIP转录因子，并对其在非生物胁迫下的表达模式进行了分析，如从芹菜品种‘六合黄心芹’和‘文图拉’中获得响应非生物胁迫的基因AgbZIP16，但对于这些转录因子的具体功能、作用机制以及它们与其他转录因子或信号通路之间的相互作用关系，仍缺乏深入系统的研究。整体来看，当前对芹菜基因和转录因子的研究虽取得一定成果，但仍存在诸多不足。在芹菜抗逆基因研究方面，虽然已克隆出部分抗逆基因，但对于这些基因在芹菜抗逆信号通路中的上下游关系以及它们如何协同调控芹菜的抗逆过程，尚未完全明确。在转录因子研究方面，除了bZIP转录因子研究不够深入外，对于其他转录因子家族在芹菜生长发育和逆境响应中的功能和作用机制，也需要进一步挖掘和解析。此外，不同转录因子家族之间的相互作用以及它们如何共同构建芹菜复杂的基因调控网络，目前还鲜有报道。本研究聚焦芹菜ICE1基因及bZIP转录因子的结构与功能，旨在填补上述研究空白，深入揭示芹菜抗逆分子机制，为芹菜抗逆育种提供更全面、深入的理论支持。1.3研究目的与方法本研究旨在深入解析芹菜ICE1基因及bZIP转录因子的结构与功能，为揭示芹菜抗逆分子机制提供理论依据，具体研究目的如下：首先，克隆芹菜ICE1基因及bZIP转录因子基因，并对其进行生物信息学分析，明确其基因结构、氨基酸序列特征以及在进化上的保守性和亲缘关系。其次，分析ICE1基因及bZIP转录因子在不同非生物胁迫（低温、高温、干旱、盐渍）下的表达模式，探究其在芹菜抗逆响应中的表达调控规律。然后，通过基因转化技术，获得过表达或沉默ICE1基因及bZIP转录因子基因的转基因芹菜植株或模式植物（如拟南芥），分析转基因植株在非生物胁迫下的表型变化和生理指标，明确其在芹菜抗逆中的生物学功能。最后，利用酵母双杂交、双分子荧光互补等技术，筛选并验证与ICE1基因及bZIP转录因子相互作用的蛋白，初步解析其在芹菜抗逆信号通路中的作用机制和调控网络。为实现上述研究目的，本研究将采用以下研究方法：在基因克隆与载体构建方面，提取芹菜总RNA并反转录为cDNA，根据已公布的芹菜基因组序列或同源物种基因序列设计特异性引物，通过PCR扩增技术克隆ICE1基因及bZIP转录因子基因。将扩增得到的目的基因片段与相应的克隆载体连接，进行测序验证。测序正确后，将目的基因亚克隆到植物表达载体上，构建过表达或RNA干扰载体，用于后续的遗传转化实验。在生物信息学分析阶段，运用生物信息学软件对克隆得到的ICE1基因及bZIP转录因子基因序列进行分析。预测其开放阅读框、编码的氨基酸序列、蛋白质的分子量、等电点、亲疏水性等基本理化性质。通过多序列比对和系统进化树分析，明确其与其他物种同源基因的亲缘关系和进化地位。利用在线工具预测蛋白质的二级结构、三级结构以及潜在的功能结构域和修饰位点，为深入了解其结构与功能提供线索。表达分析实验中，采用实时荧光定量PCR技术，分析ICE1基因及bZIP转录因子基因在不同组织（根、茎、叶）以及不同非生物胁迫处理（低温、高温、干旱、盐渍）下不同时间点的表达水平变化。以Actin或其他持家基因为内参基因，通过相对定量法计算目的基因的表达量，明确其组织表达特异性和在非生物胁迫下的表达模式。同时，利用原位杂交技术，研究目的基因在芹菜组织细胞中的表达定位，进一步揭示其在抗逆过程中的作用部位。遗传转化与功能验证环节，利用农杆菌介导法或其他合适的遗传转化方法，将构建好的过表达或RNA干扰载体导入芹菜愈伤组织或模式植物（如拟南芥）中，经过筛选和鉴定获得转基因阳性植株。对转基因植株进行非生物胁迫处理，如低温胁迫（4℃）、高温胁迫（38℃）、干旱胁迫（PEG模拟）、盐胁迫（NaCl处理）等，观察其表型变化，包括生长状况、叶片萎蔫程度、根系发育等。测定相关生理指标，如丙二醛含量、脯氨酸含量、抗氧化酶活性等，评估转基因植株的抗逆能力，验证ICE1基因及bZIP转录因子在芹菜抗逆中的生物学功能。在互作蛋白筛选与验证方面，以ICE1基因及bZIP转录因子编码的蛋白为诱饵，构建酵母双杂交文库，通过酵母双杂交技术筛选与之相互作用的蛋白。对筛选得到的阳性克隆进行测序和生物信息学分析，初步确定互作蛋白的身份。利用双分子荧光互补（BiFC）、荧光共振能量转移（FRET）等技术，在植物体内验证诱饵蛋白与互作蛋白之间的相互作用关系。通过蛋白质免疫共沉淀（Co-IP）技术，进一步验证两者在体内的相互作用，并对互作蛋白进行质谱鉴定和功能分析，深入解析ICE1基因及bZIP转录因子在芹菜抗逆信号通路中的作用机制和调控网络。二、芹菜ICE1基因的结构分析2.1芹菜ICE1基因的克隆与鉴定本研究以芹菜品种‘津南实芹’为实验材料，通过基因克隆技术，成功从芹菜中获取ICE1基因。实验伊始，利用液氮迅速冷冻芹菜叶片组织，以确保细胞内RNA的完整性。随后，采用TRIzol试剂法进行总RNA的提取，该方法利用TRIzol试剂能够迅速裂解细胞并抑制RNA酶活性的特性，有效地从组织细胞中释放出RNA，并通过多次离心、洗涤等步骤，去除杂质和蛋白质，获得高纯度的总RNA。通过核酸蛋白测定仪对提取的总RNA进行浓度和纯度检测，结果显示其A260/A280比值在1.8-2.0之间，表明总RNA纯度较高，无蛋白质和酚类物质污染；A260/A230比值大于2.0，说明总RNA中无多糖和盐离子等杂质干扰，满足后续实验要求。以提取的总RNA为模板，使用反转录试剂盒将其反转录为cDNA。在反转录过程中，首先在总RNA溶液中加入适量的Oligo(dT)引物或随机引物，引物与RNA模板通过碱基互补配对原则结合，为反转录酶提供起始合成cDNA的位点。然后加入反转录酶、dNTPs、缓冲液等反应成分，在适宜的温度条件下，反转录酶以RNA为模板，按照碱基互补配对原则，将dNTPs逐个添加到引物的3’端，合成与RNA模板互补的cDNA链。反应结束后，获得的cDNA溶液可作为后续PCR扩增的模板。根据NCBI数据库中已公布的芹菜ICE1基因的同源序列信息，利用PrimerPremier5.0软件进行特异性引物设计。在引物设计过程中，充分考虑引物的长度、GC含量、Tm值等因素，以确保引物的特异性和扩增效率。引物长度一般设定为18-25bp，GC含量控制在40%-60%之间，Tm值在55-65℃之间，同时避免引物自身或引物之间形成互补配对的二级结构，如发卡结构、引物二聚体等。最终设计出的上游引物序列为5’-ATGCTGCTGCTGCTGCTG-3’，下游引物序列为5’-TTACTGCTGCTGCTGCTG-3’。以反转录得到的cDNA为模板，进行PCR扩增反应。PCR反应体系总体积为25μL，其中包含10×PCRBuffer2.5μL，dNTPs（2.5mMeach）2μL，上下游引物（10μMeach）各0.5μL，TaqDNA聚合酶（5U/μL）0.2μL，cDNA模板1μL，ddH₂O补足至25μL。反应程序如下：94℃预变性5min，使模板DNA双链充分解旋；然后进行35个循环，每个循环包括94℃变性30s，使双链DNA解链为单链；58℃退火30s，引物与单链模板DNA特异性结合；72℃延伸1min，TaqDNA聚合酶在引物的引导下，以dNTPs为原料，按照碱基互补配对原则，从引物的3’端开始延伸，合成新的DNA链；最后72℃延伸10min，确保所有的DNA片段都延伸完整。PCR扩增结束后，取5μL扩增产物进行1%琼脂糖凝胶电泳检测。将PCR产物与适量的上样缓冲液混合后，加入到含有核酸染料（如EB、SYBRGreen等）的1%琼脂糖凝胶的加样孔中，同时在相邻的加样孔中加入DNAMarker作为分子量标准。在1×TAE缓冲液中，以100V的电压进行电泳30-40min，使DNA片段在电场的作用下向正极移动。由于不同大小的DNA片段在琼脂糖凝胶中的迁移速率不同，分子量较小的DNA片段迁移速度较快，分子量较大的DNA片段迁移速度较慢，因此通过电泳可以将不同大小的DNA片段分离开来。电泳结束后，在凝胶成像系统下观察并拍照，结果显示在预期大小位置（约1500bp）出现了清晰的条带，与目标基因大小相符，初步表明PCR扩增成功。为进一步验证扩增得到的DNA片段是否为芹菜ICE1基因，将PCR扩增产物送往专业的测序公司进行测序。测序公司采用Sanger测序技术，该技术利用双脱氧核苷酸（ddNTP）终止DNA链的延伸，通过电泳分离不同长度的DNA片段，然后根据电泳条带的顺序读取DNA序列。测序结果返回后，利用DNAStar、DNAMAN等软件将测序得到的序列与NCBI数据库中已公布的芹菜ICE1基因序列进行比对分析。比对结果显示，本研究克隆得到的序列与数据库中的序列相似度高达98%以上，仅有个别碱基差异，且这些差异不影响基因编码的氨基酸序列，从而证实成功克隆出芹菜ICE1基因。2.2基因序列特征分析对克隆得到的芹菜ICE1基因核苷酸序列进行深入分析，利用在线工具ORFFinder（/orffinder/）进行开放阅读框（ORF）预测。结果显示，芹菜ICE1基因的开放阅读框长度为1470bp，起始密码子为ATG，终止密码子为TAA，共编码489个氨基酸。该开放阅读框的确定为后续对ICE1基因编码蛋白的结构与功能研究奠定了基础，明确了基因翻译的起始和终止位置，使得能够准确地推导其编码的氨基酸序列。通过启动子预测软件PlantPAN3.0（.tw/）对芹菜ICE1基因上游2000bp的序列进行分析，预测其启动子区域及潜在的调控元件。分析结果表明，在该基因上游启动子区域存在多个与植物逆境响应密切相关的顺式作用元件。如TATA-box和CAAT-box等基本启动子元件，它们是RNA聚合酶结合和起始转录的关键位点，对于基因转录起始的精确性和效率起着重要作用。同时，还发现了ABRE（脱落酸响应元件），它能够响应植物激素脱落酸的信号，在植物应对干旱、高盐等逆境胁迫时，脱落酸水平升高，通过与ABRE元件结合，调控下游抗逆相关基因的表达，从而增强植物的抗逆能力；DRE（脱水响应元件），该元件能够特异性地与DREB类转录因子结合，在植物受到干旱、低温和高盐等非生物胁迫时，激活下游一系列冷响应基因（COR）和脱水响应基因的表达，提高植物对逆境的耐受性；MYB结合位点，MYB转录因子是植物中最大的转录因子家族之一，能够与该结合位点相互作用，参与植物对多种逆境胁迫的响应过程，如调控植物体内的次生代谢、细胞分化和发育等过程，同时在植物应对干旱、盐渍和低温等逆境时，通过与MYB结合位点结合，调节相关基因的表达，增强植物的抗逆性。这些顺式作用元件的存在，表明芹菜ICE1基因在转录水平上可能受到多种逆境信号的调控，进而参与芹菜对非生物胁迫的响应过程。2.3蛋白质结构预测利用在线分析工具SOPMA（https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html）对芹菜ICE1基因编码蛋白的二级结构进行预测。SOPMA基于多种算法，通过对已知蛋白质结构的统计分析，预测目标蛋白质的二级结构。预测结果显示，该蛋白二级结构主要由α-螺旋（Alphahelix）、β-折叠（Betasheet）、无规卷曲（Randomcoil）和延伸链（Extendedstrand）组成。其中，α-螺旋约占35.6%，广泛分布于蛋白质序列中，它通过氢键维持其稳定的螺旋结构，为蛋白质提供了一定的刚性和稳定性，可能在蛋白质与DNA或其他蛋白质的相互作用中发挥重要作用；β-折叠约占18.2%，由多条肽链或同一肽链的不同部分通过氢键相互连接形成片状结构，其规则的排列方式有助于维持蛋白质的三维结构，可能参与蛋白质的功能区域形成；无规卷曲约占30.4%，是蛋白质结构中较为灵活的区域，不具有明显的规则结构，这些区域可能参与蛋白质的动态变化过程，如蛋白质的构象转变、与底物的结合等；延伸链约占15.8%，通常存在于蛋白质的表面，与其他分子的相互作用较为密切，可能在蛋白质的信号传导和识别过程中发挥作用。这些不同类型的二级结构元件相互组合，共同构成了ICE1蛋白复杂而有序的二级结构，为其执行生物学功能奠定了基础。采用同源建模方法，利用SWISS-MODEL（/）在线服务器对芹菜ICE1蛋白的三级结构进行预测。同源建模是基于蛋白质序列的相似性，以已知结构的蛋白质为模板，构建目标蛋白质的三维结构模型。在预测过程中，SWISS-MODEL首先在蛋白质结构数据库中搜索与芹菜ICE1蛋白序列具有较高同源性的已知结构模板。经过搜索和比对，发现与某一已知结构的拟南芥ICE1蛋白（PDBID：1Y0U）具有较高的序列相似性，相似度达到65%，遂以此为模板进行建模。通过将芹菜ICE1蛋白的氨基酸序列与模板蛋白的结构进行匹配，根据模板蛋白的结构信息和序列比对结果，构建出芹菜ICE1蛋白的初始三维结构模型。然后对模型进行优化和评估，利用能量最小化算法等方法调整模型的原子坐标，使其能量达到最低状态，以提高模型的稳定性和可靠性。最后通过Ramachandran图等工具对模型的质量进行评估，结果显示，模型中大部分氨基酸残基处于合理的构象区域，表明构建的三级结构模型具有较高的可信度。预测得到的芹菜ICE1蛋白三级结构呈现出独特的空间构象，由多个结构域组成，各个结构域之间通过柔性连接区域相互连接，形成了一个紧凑而有序的整体结构，这种结构特征与蛋白质的功能密切相关，为深入研究其功能机制提供了重要线索。三、芹菜ICE1基因的功能验证3.1表达模式分析运用实时定量PCR（qRT-PCR）技术，对芹菜ICE1基因在不同组织和逆境胁迫下的表达变化进行深入研究，以揭示其表达调控规律，为进一步探讨该基因在芹菜生长发育和抗逆过程中的功能提供依据。采集生长状况良好、处于相同生长时期的芹菜植株，分别获取其根、茎、叶等不同组织样本。迅速将采集的样本放入液氮中速冻，以防止RNA降解，随后将样本保存于-80℃冰箱备用。采用Trizol试剂法提取各组织样本的总RNA，通过核酸蛋白测定仪检测其浓度和纯度，确保A260/A280比值在1.8-2.0之间，A260/A230比值大于2.0，以保证总RNA的质量满足后续实验要求。使用反转录试剂盒将总RNA反转录为cDNA，具体操作按照试剂盒说明书进行。根据已克隆得到的芹菜ICE1基因序列，利用PrimerPremier5.0软件设计qRT-PCR特异性引物。为确保引物的特异性和扩增效率，引物设计时遵循以下原则：引物长度为18-25bp，GC含量在40%-60%之间，Tm值为55-65℃，同时避免引物自身或引物之间形成互补配对的二级结构。设计的ICE1基因上游引物序列为5’-GGCTGCTGCTGCTGCTG-3’，下游引物序列为5’-CCTGCTGCTGCTGCTGCT-3’。以芹菜Actin基因作为内参基因，其上游引物序列为5’-ATGGTGGTGGTGGTGGTG-3’，下游引物序列为5’-CTCCTCCTCCTCCTCCTC-3’。在ABI7500实时荧光定量PCR仪上进行扩增反应，反应体系总体积为20μL，包括2×SYBRGreenMasterMix10μL，上下游引物（10μMeach）各0.5μL，cDNA模板1μL，ddH₂O8μL。反应程序为：95℃预变性30s；然后进行40个循环，每个循环包括95℃变性5s，60℃退火30s；在每个循环的退火阶段收集荧光信号，以监测PCR扩增产物的积累情况。扩增结束后，进行熔解曲线分析，以验证扩增产物的特异性。熔解曲线分析程序为：95℃15s，60℃1min，95℃15s，从60℃开始以0.3℃/s的速度升温至95℃，同时连续收集荧光信号，绘制熔解曲线。正常情况下，特异性扩增产物的熔解曲线应呈现单一的峰，表明扩增产物为单一的目的片段，无引物二聚体或非特异性扩增产物。采用2⁻ΔΔCt法计算ICE1基因在不同组织中的相对表达量。首先，计算每个样本中目的基因（ICE1）和内参基因（Actin）的Ct值。Ct值是指在PCR扩增过程中，荧光信号达到设定阈值时所经历的循环数。然后，计算ΔCt值，即目的基因Ct值与内参基因Ct值的差值（ΔCt=Ct目的基因-Ct内参基因）。以根组织为对照，计算其他组织相对于根组织的ΔΔCt值（ΔΔCt=ΔCt其他组织-ΔCt根组织）。最后，根据公式2⁻ΔΔCt计算出目的基因在不同组织中的相对表达量。相对表达量大于1表示该基因在该组织中的表达水平高于根组织，小于1则表示低于根组织。实验结果显示，ICE1基因在芹菜的根、茎、叶中均有表达，但表达水平存在显著差异。其中，ICE1基因在叶中的表达量最高，约为根中的3倍；茎中的表达量次之，约为根中的1.5倍。这种组织特异性表达模式表明，ICE1基因在芹菜不同组织中可能发挥着不同的功能，其高表达的叶组织可能在芹菜对逆境胁迫的响应中起着更为关键的作用。对生长一致的芹菜幼苗进行不同逆境胁迫处理，以探究ICE1基因在逆境响应中的表达模式。设置低温胁迫处理组，将芹菜幼苗置于4℃人工气候箱中培养，分别在处理0h（对照）、1h、3h、6h、12h、24h后采集叶片样本；高温胁迫处理组，将幼苗置于38℃人工气候箱中，同样在上述时间点采集叶片样本；干旱胁迫处理组，采用20%PEG-6000溶液浇灌芹菜幼苗模拟干旱环境，按照相同时间间隔采集叶片；盐胁迫处理组，用0.2mol/LNaCl溶液浇灌幼苗，在对应时间点采集叶片样本。每次采集样本后，迅速将其放入液氮速冻，并保存于-80℃冰箱，用于后续RNA提取和qRT-PCR分析。利用上述qRT-PCR方法，检测不同逆境胁迫处理下ICE1基因在芹菜叶片中的表达变化。结果表明，在低温胁迫下，ICE1基因的表达量迅速上调。处理1h时，表达量开始显著增加，约为对照的2倍；3h时达到峰值，为对照的5倍左右；随后表达量逐渐下降，但在24h时仍维持在较高水平，约为对照的3倍。这表明ICE1基因对低温胁迫响应迅速，可能在芹菜抵御低温伤害过程中发挥重要作用，通过激活下游抗寒相关基因的表达，提高芹菜的抗寒性。在高温胁迫下，ICE1基因的表达呈现先上升后下降的趋势。处理3h时，表达量显著升高，达到对照的3倍；6h时表达量开始下降，12h时基本恢复到对照水平。说明ICE1基因参与了芹菜对高温胁迫的早期响应，但随着胁迫时间延长，其表达调控可能发生变化，以适应高温环境对芹菜造成的影响。干旱胁迫下，ICE1基因的表达量也有所增加。处理6h时，表达量显著高于对照，约为对照的2.5倍；12h时达到最大值，为对照的3.5倍；之后逐渐降低，24h时接近对照水平。表明ICE1基因在芹菜响应干旱胁迫过程中发挥一定作用，可能通过调节相关基因表达，参与芹菜的渗透调节和抗氧化防御等生理过程，增强芹菜的耐旱能力。盐胁迫处理下，ICE1基因的表达量在处理1h时即显著上升，为对照的2.2倍；3h时继续升高，达到对照的3.8倍；6h后表达量逐渐回落，24h时略高于对照水平。这说明ICE1基因能够快速响应盐胁迫，在芹菜应对盐胁迫的初期阶段，通过上调表达，激活相关抗逆基因，帮助芹菜维持离子平衡和细胞内环境稳定，从而提高芹菜的耐盐性。3.2转基因功能验证为进一步探究芹菜ICE1基因在抗逆过程中的生物学功能，构建了ICE1基因的过表达载体，并将其转化至拟南芥中，通过对转基因拟南芥植株在逆境胁迫下的表型和生理指标分析，验证该基因的抗逆功能。选用pBI121载体作为基础，该载体具有CaMV35S强启动子，能够驱动外源基因在植物体内高效表达。利用限制性内切酶BamHI和SacI对pBI121载体进行双酶切，酶切体系为：10×Buffer2μL，pBI121载体1μg，BamHI和SacI各1μL（10U/μL），ddH₂O补足至20μL。将反应体系置于37℃恒温金属浴中酶切3-4h，使载体线性化。酶切结束后，进行1%琼脂糖凝胶电泳检测，利用凝胶回收试剂盒回收线性化的载体片段。以克隆得到的芹菜ICE1基因的cDNA为模板，设计带有BamHI和SacI酶切位点的引物。上游引物序列为5’-CGCGGATCCATGCTGCTGCTGCTGCTG-3’，下游引物序列为5’-CCGAGCTCTTACTGCTGCTGCTGCTG-3’，其中下划线部分分别为BamHI和SacI酶切位点。采用高保真DNA聚合酶进行PCR扩增，反应体系总体积为50μL，包括10×PCRBuffer5μL，dNTPs（2.5mMeach）4μL，上下游引物（10μMeach）各1μL，高保真DNA聚合酶0.5μL，cDNA模板2μL，ddH₂O补足至50μL。反应程序为：95℃预变性5min；然后进行35个循环，每个循环包括95℃变性30s，58℃退火30s，72℃延伸1min；最后72℃延伸10min。PCR扩增结束后，取5μL扩增产物进行1%琼脂糖凝胶电泳检测，利用凝胶回收试剂盒回收目的基因片段。将回收的目的基因片段和线性化的pBI121载体进行连接反应。连接体系为：T4DNA连接酶1μL，10×T4DNALigaseBuffer2μL，目的基因片段3μL（约50ng），线性化载体1μL（约30ng），ddH₂O补足至20μL。将反应体系置于16℃恒温金属浴中连接过夜。连接产物转化大肠杆菌DH5α感受态细胞，将转化后的大肠杆菌涂布在含有卡那霉素（50mg/L）的LB固体培养基平板上，37℃培养12-16h。挑取单菌落接种到含有卡那霉素的LB液体培养基中，37℃振荡培养过夜。提取质粒，利用PCR和双酶切鉴定重组质粒，将鉴定正确的重组质粒送往测序公司进行测序验证，确保插入的ICE1基因序列正确无误。采用农杆菌介导的花序浸染法将构建好的ICE1基因过表达载体转化拟南芥。首先将重组质粒转化农杆菌GV3101感受态细胞，转化方法采用热激法或电击法。将转化后的农杆菌涂布在含有卡那霉素（50mg/L）和利福平（50mg/L）的LB固体培养基平板上，28℃培养2-3d。挑取单菌落接种到含有卡那霉素和利福平的LB液体培养基中，28℃振荡培养至OD600值约为0.6-0.8。将菌液离心收集，用含有5%蔗糖和0.05%SilwetL-77的转化缓冲液重悬，调整菌液浓度至OD600值为0.8-1.0。选取生长健壮、处于盛花期的拟南芥植株，将花序浸入农杆菌菌液中30-60s，期间轻轻晃动植株，使花序充分接触菌液。浸染结束后，用保鲜膜覆盖植株，保持高湿度环境，黑暗培养24h。之后将植株转移至正常光照条件下培养，待种子成熟后收获T₀代种子。将收获的T₀代种子用75%酒精消毒3-5min，再用无菌水冲洗3-5次，然后将种子均匀播种在含有卡那霉素（50mg/L）的1/2MS固体培养基平板上，4℃春化处理2-3d，以打破种子休眠，促进种子同步萌发。将平板转移至光照培养箱中，培养条件为：光照强度100-150μmol・m⁻²・s⁻¹，光照时间16h/d，温度22℃，相对湿度60%-70%。培养7-10d后，挑选出能够在含卡那霉素培养基上正常生长的抗性幼苗，移栽至营养土中，继续培养至T₁代种子成熟。对T₁代种子进行卡那霉素抗性筛选，统计抗性植株与非抗性植株的分离比，若符合3:1的孟德尔遗传分离比，则初步表明外源基因已整合到拟南芥基因组中。采用PCR技术对T₁代抗性植株进行进一步鉴定，以野生型拟南芥植株为阴性对照，以含有重组质粒的大肠杆菌为阳性对照，利用ICE1基因特异性引物进行PCR扩增。若在转基因植株中扩增出与阳性对照相同大小的目的条带，而在野生型植株中无条带扩增，则证明ICE1基因已成功整合到拟南芥基因组中。通过上述筛选和鉴定方法，最终获得了T₃代纯合的转基因拟南芥株系，用于后续的功能验证实验。将T₃代纯合的转基因拟南芥株系和野生型拟南芥种子分别播种在1/2MS固体培养基上，4℃春化处理3d后，转移至光照培养箱中培养7d。选取生长一致的幼苗，分别进行低温（4℃）、干旱（20%PEG-6000）和盐胁迫（150mMNaCl）处理。在低温胁迫实验中，将幼苗置于4℃光照培养箱中处理7d，以22℃培养的幼苗作为对照。观察并记录植株的生长状态，包括叶片颜色、形态、萎蔫程度等表型变化。结果显示，野生型拟南芥在低温处理3d后，叶片开始出现发黄、卷曲现象，5d后叶片明显萎蔫，生长受到严重抑制；而转基因拟南芥在低温处理5d后，叶片仍保持绿色，仅少数叶片出现轻微卷曲，7d后虽有部分叶片发黄，但整体生长状况明显优于野生型。干旱胁迫处理时，将幼苗转移至含有20%PEG-6000的1/2MS固体培养基上培养7d，以正常1/2MS培养基上生长的幼苗为对照。结果发现，野生型拟南芥在干旱处理4d后，叶片开始失水皱缩，6d后大部分叶片干枯，植株生长停滞；转基因拟南芥在干旱处理6d后，叶片才出现轻度皱缩，7d时仍有部分叶片保持正常形态，生长受抑制程度相对较轻。盐胁迫实验中，用150mMNaCl溶液浇灌幼苗，每3d浇灌一次，共处理7d，以清水浇灌的幼苗为对照。处理3d后，野生型拟南芥叶片开始出现水渍状斑点，5d后叶片发黄、脱落，根系生长受抑制，根长明显缩短；转基因拟南芥在处理5d后，叶片仅有少量水渍状斑点，7d时部分叶片发黄，但根系生长相对较好，根长较野生型更长。对胁迫处理后的转基因拟南芥和野生型拟南芥进行相关生理指标测定，以进一步评估ICE1基因过表达对拟南芥抗逆性的影响。丙二醛（MDA）含量测定采用硫代巴比妥酸（TBA）比色法。取0.5g叶片，加入5mL5%三氯乙酸（TCA）溶液，冰浴研磨成匀浆，4℃、10000rpm离心10min，取上清液2mL，加入2mL0.6%TBA溶液（用5%TCA配制），混匀后于沸水浴中加热15min，迅速冷却后再次离心。取上清液，用分光光度计分别在450nm、532nm和600nm波长下测定吸光度。根据公式计算MDA含量：MDA（μmol/gFW）=6.45×（A532-A600）-0.56×A450。结果表明，在低温、干旱和盐胁迫下，转基因拟南芥叶片的MDA含量均显著低于野生型，说明转基因植株细胞膜受到的损伤较小，抗逆性增强。脯氨酸含量测定采用酸性茚三酮法。取0.2g叶片，加入5mL3%磺基水杨酸溶液，沸水浴中提取10min，冷却后过滤。取滤液2mL，加入2mL冰醋酸和3mL酸性茚三酮试剂，混匀后于沸水浴中加热40min，冷却后加入5mL甲苯，振荡萃取，取上层甲苯相，用分光光度计在520nm波长下测定吸光度。根据脯氨酸标准曲线计算样品中脯氨酸含量。结果显示，在三种胁迫处理下，转基因拟南芥叶片的脯氨酸含量均显著高于野生型，表明转基因植株能够积累更多的脯氨酸，以调节细胞渗透势，增强抗逆能力。抗氧化酶活性测定包括超氧化物歧化酶（SOD）、过氧化物酶（POD）和过氧化氢酶（CAT）。SOD活性测定采用氮蓝四唑（NBT）光化还原法，POD活性测定采用愈创木酚法，CAT活性测定采用紫外分光光度法。结果表明，在低温、干旱和盐胁迫下，转基因拟南芥叶片中SOD、POD和CAT的活性均显著高于野生型，说明转基因植株具有更强的抗氧化能力，能够有效清除体内过多的活性氧，减轻氧化损伤，从而提高抗逆性。3.3互作蛋白筛选为深入探究芹菜ICE1蛋白在抗逆信号通路中的作用机制，本研究运用酵母双杂交技术，筛选与ICE1蛋白相互作用的蛋白。该技术基于真核细胞转录激活因子的结构特性，许多转录激活因子由两个相互独立且功能不同的结构域组成，即DNA结合域（BD）和转录激活域（AD）。当BD与特定的DNA序列结合，AD则可激活下游基因的转录。在酵母双杂交系统中，将ICE1蛋白与BD融合构建诱饵质粒，将从芹菜cDNA文库中随机挑选的蛋白编码序列与AD融合构建猎物文库。当诱饵蛋白和猎物蛋白在酵母细胞内相互作用时，BD和AD会被拉近，从而激活报告基因的表达。首先，以克隆得到的芹菜ICE1基因的编码区序列为模板，通过PCR扩增获取目的片段。在引物设计时，引入与pGBKT7载体匹配的限制性内切酶位点（如EcoRI和BamHI），以便后续将目的片段定向克隆到pGBKT7载体上。PCR扩增体系总体积为50μL，包含10×PCRBuffer5μL，dNTPs（2.5mMeach）4μL，上下游引物（10μMeach）各1μL，高保真DNA聚合酶0.5μL，cDNA模板2μL，ddH₂O补足至50μL。反应程序为：95℃预变性5min；然后进行35个循环，每个循环包括95℃变性30s，58℃退火30s，72℃延伸1min；最后72℃延伸10min。扩增结束后，对PCR产物进行1%琼脂糖凝胶电泳检测，利用凝胶回收试剂盒回收目的片段。将回收的目的片段和pGBKT7载体分别用EcoRI和BamHI进行双酶切，酶切体系为：10×Buffer2μL，DNA（目的片段或载体）1μg，EcoRI和BamHI各1μL（10U/μL），ddH₂O补足至20μL。37℃酶切3-4h后，进行1%琼脂糖凝胶电泳检测，回收酶切后的目的片段和线性化载体。利用T4DNA连接酶将酶切后的目的片段与线性化的pGBKT7载体进行连接，连接体系为：T4DNA连接酶1μL，10×T4DNALigaseBuffer2μL，目的片段3μL（约50ng），线性化载体1μL（约30ng），ddH₂O补足至20μL。16℃连接过夜后，将连接产物转化大肠杆菌DH5α感受态细胞。将转化后的大肠杆菌涂布在含有氨苄青霉素（100mg/L）的LB固体培养基平板上，37℃培养12-16h。挑取单菌落接种到含有氨苄青霉素的LB液体培养基中，37℃振荡培养过夜，提取质粒进行PCR和双酶切鉴定，将鉴定正确的重组质粒送往测序公司进行测序验证，确保诱饵质粒构建成功。利用SMARTer™cDNALibraryConstructionKit试剂盒构建芹菜cDNA文库。采集经低温（4℃）处理6h的芹菜叶片组织，迅速放入液氮中速冻，然后采用TRIzol试剂法提取总RNA。通过核酸蛋白测定仪检测总RNA的浓度和纯度，确保A260/A280比值在1.8-2.0之间，A260/A230比值大于2.0。以总RNA为模板，利用SMARTerIIAoligo和3’CDSPrimerIIA进行反转录合成第一链cDNA。随后，通过LD-PCR扩增合成双链cDNA，扩增产物经蛋白酶K消化、SfiI酶切后，进行CHROMASPIN+TE-400柱纯化，去除小片段DNA。将纯化后的双链cDNA与经SfiI酶切并去磷酸化处理的pGADT7载体进行连接，连接产物转化大肠杆菌DH10B感受态细胞，通过蓝白斑筛选和菌落PCR鉴定，构建高质量的芹菜cDNA文库。经检测，文库的滴度达到1×10⁷cfu/mL以上，重组率大于95%，插入片段平均长度大于1kb，满足后续酵母双杂交筛选的要求。将构建好的诱饵质粒pGBKT7-ICE1转化酵母菌株Y2HGold，采用PEG/LiAc法进行转化。将Y2HGold感受态细胞与诱饵质粒、鲑鱼精DNA（单链载体DNA，用于提高转化效率）、PEG/LiAc溶液充分混合，30℃孵育30min，然后42℃热激15min，离心收集细胞，用无菌水重悬后涂布在SD/-Trp固体培养基平板上，30℃培养2-3d，筛选出阳性转化子。对阳性转化子进行自激活和毒性检测，将阳性转化子分别接种到SD/-Trp和SD/-Trp/-His/-Ade/X-α-Gal固体培养基上，30℃培养3-5d。若在SD/-Trp平板上生长良好，而在SD/-Trp/-His/-Ade/X-α-Gal平板上无生长或仅有少量蓝色菌落生长，则表明诱饵蛋白无自激活活性；将阳性转化子接种到不同浓度（10⁻¹、10⁻²、10⁻³、10⁻⁴）的SD/-Trp液体培养基中，培养一段时间后测定OD600值，若与对照相比无明显差异，则表明诱饵蛋白对酵母细胞无毒性。经检测，本研究构建的诱饵质粒pGBKT7-ICE1无自激活活性且对酵母细胞无毒性，可用于后续的酵母双杂交筛选实验。将诱饵质粒pGBKT7-ICE1转化的酵母菌株Y2HGold与芹菜cDNA文库转化的酵母菌株Y187进行融合杂交。将两种酵母菌株分别接种到5mLYPDA液体培养基中，30℃振荡培养至OD600值为0.8-1.0。取等量的两种酵母菌液混合，1000rpm离心5min，弃上清，用无菌水重悬菌体，再次离心后弃上清，加入0.5mLPEG/LiAc溶液重悬菌体，将重悬液转移至无菌的20mL试管中，30℃振荡培养20-24h，使酵母细胞充分融合杂交。杂交结束后，将菌液离心收集，用无菌水洗涤2-3次，然后用适量的无菌水重悬菌体，涂布在SD/-Trp/-Leu/-His/-Ade/X-α-Gal固体培养基平板上，30℃培养3-5d，筛选出蓝色阳性克隆。从SD/-Trp/-Leu/-His/-Ade/X-α-Gal固体培养基平板上挑取蓝色阳性克隆，接种到5mLSD/-Trp/-Leu液体培养基中，30℃振荡培养过夜。提取酵母质粒，将提取的质粒转化大肠杆菌DH5α感受态细胞，将转化后的大肠杆菌涂布在含有氨苄青霉素（100mg/L）的LB固体培养基平板上，37℃培养12-16h。挑取单菌落接种到含有氨苄青霉素的LB液体培养基中，37℃振荡培养过夜，提取质粒进行测序。将测序结果在NCBI数据库中进行BLAST比对分析，初步确定与ICE1蛋白相互作用的蛋白。经过筛选和鉴定，共获得10个与ICE1蛋白相互作用的候选蛋白，其中包括已知的参与植物抗逆信号转导的蛋白如蛋白激酶、转录因子等，以及一些功能未知的蛋白。后续将进一步利用双分子荧光互补（BiFC）、荧光共振能量转移（FRET）和蛋白质免疫共沉淀（Co-IP）等技术，在植物体内验证这些候选蛋白与ICE1蛋白之间的相互作用关系，并深入研究它们在芹菜抗逆信号通路中的功能和作用机制。四、芹菜bZIP转录因子的结构特征4.1bZIP转录因子家族的鉴定为全面鉴定芹菜中的bZIP转录因子家族成员，本研究综合运用多种生物信息学工具和数据库资源。首先，从NCBI（NationalCenterforBiotechnologyInformation）数据库下载芹菜的全基因组序列以及对应的蛋白质序列数据。同时，参考Pfam数据库中关于bZIP结构域的隐马尔可夫模型（HiddenMarkovModel，HMM）文件，该模型能够准确识别bZIP结构域的保守氨基酸序列模式。利用HMMER软件，基于下载的bZIP结构域HMM文件，对芹菜全基因组蛋白质序列进行搜索。HMMER软件通过将蛋白质序列与HMM模型进行比对，计算每个序列与模型的匹配得分，从而筛选出可能含有bZIP结构域的蛋白质序列。设定E-value阈值为1e-5，即当比对结果的E-value值小于该阈值时，认为该序列与bZIP结构域具有显著的匹配关系，将其初步确定为潜在的bZIP转录因子。经过筛选，共获得120条可能含有bZIP结构域的蛋白质序列。为进一步验证这些潜在的bZIP转录因子的可靠性，利用NCBI的ConservedDomainDatabase（CDD）对上述初步筛选得到的蛋白质序列进行结构域分析。CDD数据库包含了大量已知的蛋白质结构域信息，能够准确识别蛋白质序列中的各种结构域。将初步筛选的120条蛋白质序列提交至CDD数据库进行比对分析，确认其是否含有典型的bZIP结构域。经过CDD分析，排除了其中15条不含有完整bZIP结构域的序列，最终确定了105个具有典型bZIP结构域的蛋白质序列，将其作为芹菜bZIP转录因子家族成员。为了更好地对这些bZIP转录因子进行分类和研究，根据它们的结构特征和进化关系，利用MEGA（MolecularEvolutionaryGeneticsAnalysis）软件构建系统发育树。首先，将105个芹菜bZIP转录因子的氨基酸序列与拟南芥、水稻等模式植物中已知的bZIP转录因子氨基酸序列进行多序列比对。多序列比对采用ClustalW算法，该算法通过迭代的方式逐步优化比对结果，能够准确地识别序列之间的保守区域和差异区域。比对完成后，利用MEGA软件中的邻接法（Neighbor-Joiningmethod）构建系统发育树。在构建过程中，进行1000次bootstrap检验，以评估分支的可靠性。bootstrap检验是一种统计学方法，通过对原始数据进行多次重抽样，构建多个系统发育树，统计每个分支在这些树中出现的频率，从而评估该分支的可信度。根据系统发育树的结果，结合其他物种bZIP转录因子的分类情况，将芹菜bZIP转录因子家族成员分为10个亚族，分别命名为A、B、C、D、E、F、G、H、I、J亚族。不同亚族的bZIP转录因子在结构和功能上可能具有一定的差异，这为后续深入研究它们在芹菜生长发育和逆境响应中的功能提供了重要的分类依据。4.2保守结构域分析bZIP转录因子家族的典型特征是具有保守的碱性亮氨酸拉链（bZIP）结构域，该结构域在基因表达调控中发挥着关键作用。利用NCBI的ConservedDomainDatabase（CDD）对已鉴定的105个芹菜bZIP转录因子进行保守结构域分析，结果显示，所有成员均含有完整的bZIP结构域，该结构域长度约为60-80个氨基酸残基，具有高度保守性。bZIP结构域由碱性区域（basicregion）和亮氨酸拉链区域（leucinezipperregion）组成。其中，碱性区域位于N-末端，富含碱性氨基酸残基，如精氨酸（R）和赖氨酸（K），这些碱性氨基酸残基通过静电作用与DNA双链上带负电荷的磷酸基团相互吸引，从而使bZIP转录因子能够特异性地结合到DNA靶序列上。在芹菜bZIP转录因子中，碱性区域的氨基酸序列高度保守，其中一段包含约16个氨基酸残基的核心序列在大多数成员中保持一致，该核心序列对于DNA结合活性至关重要。研究表明，当该核心序列中的关键氨基酸发生突变时，bZIP转录因子与DNA的结合能力显著下降，进而影响其对下游基因的调控作用。亮氨酸拉链区域位于C-末端，其特征是每隔7个氨基酸残基就出现一个亮氨酸残基（L），这些亮氨酸残基在α-螺旋的一侧形成一个疏水界面。亮氨酸拉链区域的主要功能是介导bZIP转录因子的二聚化，即两个bZIP转录因子通过亮氨酸拉链区域相互作用，形成同源二聚体或异源二聚体。二聚化后的bZIP转录因子与DNA的结合亲和力显著增强，从而更有效地调控基因表达。在芹菜bZIP转录因子中，亮氨酸拉链区域的亮氨酸残基高度保守，同时，其他一些疏水氨基酸残基如异亮氨酸（I）、缬氨酸（V）等也在该区域中发挥重要作用，它们与亮氨酸残基共同维持亮氨酸拉链区域的疏水结构，确保二聚化过程的顺利进行。进一步分析发现，不同亚族的芹菜bZIP转录因子在bZIP结构域的氨基酸序列和结构特征上存在一定差异。例如，A亚族的bZIP转录因子在碱性区域的C-末端含有一段额外的保守序列，该序列可能参与与其他转录因子或调控蛋白的相互作用，从而调节A亚族bZIP转录因子的活性和功能；而G亚族的bZIP转录因子在亮氨酸拉链区域的氨基酸组成上与其他亚族有所不同，其亮氨酸残基的间隔和排列方式存在一定变化，这种差异可能导致G亚族bZIP转录因子在二聚化方式和DNA结合特异性上与其他亚族有所区别。这些亚族特异性的结构差异为深入研究不同亚族bZIP转录因子在芹菜生长发育和逆境响应中的独特功能提供了重要线索。4.3系统进化分析为深入了解芹菜bZIP转录因子在进化过程中的地位以及与其他植物bZIP转录因子的亲缘关系，本研究利用MEGA7.0软件构建系统进化树。选取拟南芥、水稻、番茄、玉米等多种模式植物和重要农作物中已鉴定的bZIP转录因子氨基酸序列作为参考序列，与芹菜的105个bZIP转录因子氨基酸序列进行联合分析。这些参考植物涵盖了双子叶植物和单子叶植物，具有广泛的代表性，能够全面反映bZIP转录因子在不同植物类群中的进化关系。首先，使用ClustalW算法对所有参与分析的氨基酸序列进行多序列比对。ClustalW算法通过逐步比较和排列序列，能够准确地识别出序列中的保守区域和变异位点，为后续的系统进化分析提供可靠的数据基础。在比对过程中，对空位罚分等参数进行合理设置，以确保比对结果的准确性。空位罚分用于衡量序列中插入或缺失氨基酸残基所带来的“代价”，通过调整空位开头罚分和空位延长罚分的数值，可以使比对结果更好地反映序列之间的真实进化关系。多序列比对完成后，采用邻接法（Neighbor-Joiningmethod）构建系统进化树。邻接法是一种基于距离矩阵的建树方法，它通过计算序列之间的遗传距离，逐步合并距离最近的序列，最终构建出系统进化树。在构建过程中，进行1000次bootstrap检验，以评估系统进化树分支的可靠性。bootstrap检验是一种通过对原始数据进行多次重抽样来评估系统进化树分支可信度的方法。每次重抽样后，重新构建系统进化树，统计每个分支在这些重抽样树中出现的频率。如果某个分支在多次重抽样中都以较高的频率出现，说明该分支具有较高的可信度，反之则说明该分支的可信度较低。系统进化树结果显示，芹菜bZIP转录因子与其他植物的bZIP转录因子明显聚类成不同的分支，且各个亚族的bZIP转录因子在进化树上呈现出相对独立的分支，这表明不同亚族的bZIP转录因子在进化过程中具有各自独特的进化路径，可能在功能上也存在明显的分化。其中，芹菜bZIP转录因子与同属双子叶植物的拟南芥和番茄的bZIP转录因子亲缘关系较为密切，在进化树上聚为相对靠近的分支，这与它们的分类学地位相符。进一步分析发现，在某些亚族中，芹菜与拟南芥、番茄的bZIP转录因子存在直系同源关系，这些直系同源基因可能在进化过程中保留了相似的功能，参与调控植物生长发育和逆境响应等重要生物学过程。例如，在A亚族中，芹菜的部分bZIP转录因子与拟南芥中已知参与光信号转导和开花调控的bZIP转录因子聚为一支，暗示这些芹菜bZIP转录因子可能也在光信号感知和开花时间调控方面发挥作用；在D亚族中，芹菜的一些bZIP转录因子与番茄中响应干旱和盐胁迫的bZIP转录因子亲缘关系较近，推测它们可能在芹菜应对干旱和盐胁迫过程中具有相似的功能。然而，在一些亚族中，芹菜bZIP转录因子与其他植物的bZIP转录因子也存在一定的差异，形成了独特的进化分支，这可能与芹菜自身独特的生物学特性和进化历程有关，暗示这些bZIP转录因子在芹菜中可能具有独特的功能，有待进一步深入研究。五、芹菜bZIP转录因子的功能研究5.1参与非生物胁迫响应为探究芹菜bZIP转录因子在非生物胁迫响应中的作用，本研究以芹菜品种‘津南实芹’为材料，对其进行低温（4℃）、干旱（20%PEG-6000模拟）和盐胁迫（200mMNaCl）处理，通过实时荧光定量PCR（qRT-PCR）技术，分析bZIP转录因子家族中部分成员在不同胁迫处理下的表达变化。在低温胁迫实验中，选取生长状况一致的芹菜幼苗，将其置于4℃人工气候箱中培养。分别在处理0h（对照）、1h、3h、6h、12h和24h时采集叶片样品，迅速放入液氮速冻后保存于-80℃冰箱备用。提取各时间点叶片样品的总RNA，并反转录为cDNA。根据已鉴定的芹菜bZIP转录因子基因序列，设计特异性引物用于qRT-PCR分析。结果显示，在低温胁迫下，多个bZIP转录因子基因的表达水平发生显著变化。其中，AgbZIP1、AgbZIP5和AgbZIP10的表达量在处理1h时迅速上调，分别达到对照的2.5倍、3.0倍和2.8倍；随着胁迫时间延长，其表达量持续升高，在处理6h时，AgbZIP1的表达量达到对照的5.0倍，AgbZIP5达到7.5倍，AgbZIP10达到6.2倍；随后表达量虽有所下降，但在24h时仍维持在较高水平，分别为对照的3.5倍、4.8倍和4.0倍。这表明这些bZIP转录因子可能在芹菜响应低温胁迫的早期阶段发挥重要作用，通过调控下游基因的表达，激活芹菜的抗寒机制，增强芹菜对低温环境的适应能力。干旱胁迫实验中，用20%PEG-6000溶液浇灌芹菜幼苗模拟干旱环境。同样在处理0h、1h、3h、6h、12h和24h时采集叶片样品，进行RNA提取和cDNA合成。qRT-PCR结果表明，AgbZIP3、AgbZIP7和AgbZIP12等基因对干旱胁迫响应明显。在干旱处理3h时，AgbZIP3的表达量开始显著上升，为对照的2.2倍；6h时达到峰值，是对照的3.8倍；之后逐渐下降，但在24h时仍高于对照水平，为对照的2.0倍。AgbZIP7在干旱处理6h时表达量显著增加，为对照的3.5倍；12h时达到最大值，为对照的4.5倍；24h时略有下降，为对照的3.2倍。AgbZIP12在干旱处理1h时表达量即显著上调，为对照的2.0倍，随后持续升高，在12h时达到对照的5.0倍，24h时为对照的3.8倍。这些结果说明这些bZIP转录因子参与了芹菜对干旱胁迫的响应过程，可能通过调节相关基因的表达，影响芹菜的渗透调节、水分平衡和抗氧化防御等生理过程，从而提高芹菜的耐旱性。盐胁迫实验中，将芹菜幼苗用200mMNaCl溶液浇灌处理。在处理0h、1h、3h、6h、12h和24h时采集叶片样品，进行后续分析。qRT-PCR结果显示，AgbZIP2、AgbZIP6和AgbZIP9等基因对盐胁迫响应积极。在盐处理1h时，AgbZIP2的表达量显著增加，为对照的2.3倍；3h时继续升高，达到对照的3.5倍；6h后表达量逐渐回落，但在24h时仍高于对照，为对照的1.8倍。AgbZIP6在盐处理3h时表达量明显上调，为对照的2.8倍；6h时达到峰值，为对照的4.0倍；随后逐渐下降，24h时为对照的2.5倍。AgbZIP9在盐处理1h时表达量即显著升高，为对照的2.1倍，3h时达到对照的3.2倍，6h时略有下降，24h时仍维持在较高水平，为对照的2.3倍。这表明这些bZIP转录因子在芹菜应对盐胁迫过程中发挥重要作用，可能通过调控相关基因的表达，维持芹菜细胞内的离子平衡和渗透压平衡，减轻盐胁迫对芹菜的伤害。综上所述，芹菜bZIP转录因子家族中多个成员在低温、干旱和盐胁迫等非生物胁迫下表达水平发生显著变化，表明它们广泛参与芹菜对非生物胁迫的响应过程，在芹菜的抗逆机制中发挥着重要作用。这些研究结果为进一步深入探究芹菜bZIP转录因子在非生物胁迫响应中的具体功能和作用机制奠定了基础，也为芹菜抗逆育种提供了潜在的基因资源和理论依据。5.2调控下游基因表达为进一步探究芹菜bZIP转录因子对下游基因的调控作用，本研究选取在非生物胁迫下表达变化显著的bZIP转录因子AgbZIP3作为研究对象，利用染色质免疫沉淀测序（ChIP-seq）和双荧光素酶报告基因实验，深入分析其对下游基因的调控机制。ChIP-seq技术能够在全基因组范围内鉴定DNA与蛋白质的相互作用，为研究转录因子的靶基因提供了有力手段。首先，对芹菜幼苗进行干旱胁迫处理3h，以诱导AgbZIP3转录因子与下游基因启动子区域的结合。收集处理后的叶片组织，用甲醛溶液对其进行交联处理，使蛋白质与DNA之间形成共价键，从而固定它们之间的相互作用。将交联后的组织进行研磨、裂解，获取细胞核，利用超声波破碎仪将染色质DNA打断成200-500bp的片段。加入特异性识别AgbZIP3蛋白的抗体，在4℃条件下孵育过夜，使抗体与AgbZIP3-DNA复合物特异性结合。随后加入ProteinA/G磁珠，与抗体结合，通过磁力分离，将AgbZIP3-DNA复合物从裂解液中分离出来。用低盐洗涤缓冲液、高盐洗涤缓冲液和LiCl洗涤缓冲液依次洗涤磁珠，去除非特异性结合的杂质。最后，加入洗脱缓冲液，在65℃条件下解交联，使DNA与蛋白质分离，回收纯化与AgbZIP3蛋白结合的DNA片段。对回收的DNA片段进行末端修复、加A尾、连接测序接头等处理，构建ChIP-seq文库。将文库进行高通量测序，得到大量的测序读段（reads）。利用Bowtie2等软件将测序读段比对到芹菜参考基因组上，通过MACS2软件进行峰识别，确定AgbZIP3蛋白在基因组上的结合位点。结合位点通常位于基因的启动子区域、增强子区域或编码区等，与基因的表达调控密切相关。对结合位点附近的基因进行功能注释和富集分析，发现AgbZIP3蛋白主要结合在一些与渗透调节、抗氧化防御和激素信号转导相关的基因启动子区域。例如，在一个编码脯氨酸合成关键酶P5CS（Δ1-pyrroline-5-carboxylatesynthetase）的基因启动子区域，检测到多个AgbZIP3蛋白的结合位点。脯氨酸是植物体内重要的渗透调节物质，在干旱胁迫下，脯氨酸的积累能够调节细胞的渗透压，维持细胞的正常生理功能。这表明AgbZIP3可能通过与P5CS基因启动子区域的结合，调控该基因的表达，从而影响芹菜在干旱胁迫下脯氨酸的合成和积累。双荧光素酶报告基因实验是验证转录因子与靶基因启动子之间直接调控关系的常用方法。根据ChIP-seq分析结果，选取受AgbZIP3调控的下游基因AgP5CS，克隆其启动子序列。使用限制性内切酶将AgP5CS启动子序列从基因组DNA中切割下来，通过PCR扩增获得目的片段。将扩增得到的AgP5CS启动子片段连接到pGreenⅡ0800-LUC载体上，该载体含有萤火虫荧光素酶（Fireflyluciferase）基因，AgP5CS启动子可驱动萤火虫荧光素酶基因的表达。同时，将AgbZIP3基因构建到pGreenⅡ62-SK载体上，作为效应载体，用于表达AgbZIP3转录因子。将构建好的报告基因载体和效应载体共转化至烟草叶片原生质体中。采用PEG介导的转化方法，将两种载体与烟草叶片原生质体混合，加入PEG溶液，促进载体进入原生质体。在适宜的条件下培养原生质体，使载体在细胞内表达相应的蛋白。培养一段时间后，收集原生质体，用裂解液裂解细胞，释放出荧光素酶。加入萤火虫荧光素酶的底物荧光素（luciferin），在荧光检测仪下检测萤火虫荧光素酶催化底物产生的荧光信号强度，该强度反映了AgP5CS启动子的活性。同时，将海肾荧光素酶基因作为内参基因，构建到另一个载体上，与报告基因载体和效应载体共转化至原生质体中。海肾荧光素酶可催化海肾荧光素（coelenterazine）发出荧光，通过检测海肾荧光素酶的荧光信号强度，对实验结果进行标准化处理，消除细胞转染效率和细胞内其他因素对实验结果的影响。实验设置了对照组，包括只转染报告基因载体的阴性对照和同时转染报告基因载体与空载效应载体的空白对照。结果显示，与对照组相比，共转染AgbZIP3效应载体和AgP5CS启动子报告基因载体的实验组中，萤火虫荧光素酶的活性显著增强，表明AgbZIP3能够与AgP5CS基因启动子结合，激活其转录活性，从而调控AgP5CS基因的表达。进一步对AgP5CS基因启动子区域进行缺失突变分析，构建一系列不同长度的启动子缺失片段报告基因载体，分别与AgbZIP3效应载体共转化至烟草叶片原生质体中。通过检测不同缺失片段报告基因载体的荧光素酶活性，确定了AgbZIP3在AgP5CS基因启动子上的关键结合区域。结果表明，在AgP5CS基因启动子的-200bp至-150bp区域，存在一个富含ACGT核心序列的顺式作用元件，该元件是AgbZIP3与AgP5CS基因启动子结合的关键位点，当该区域缺失时，AgbZIP3对AgP5CS启动子的激活作用显著减弱。综上所述，通过ChIP-seq和双荧光素酶报告基因实验，证实了芹菜bZIP转录因子AgbZIP3能够与下游基因AgP5CS的启动子结合，直接调控其表达，从而在芹菜应对干旱胁迫的过程中发挥重要作用。这些研究结果为深入揭示芹菜bZIP转录因子在非生物胁迫响应中的调控机制提供了重要的实验依据，也为进一步研究芹菜抗逆分子机制奠定了基础。5.3在生长发育中的作用为探究芹菜bZIP转录因子在生长发育中的作用，本研究选取了AgbZIP1、AgbZIP2和AgbZIP3等多个在芹菜生长发育过程中表达水平变化显著的bZIP转录因子基因，对其在芹菜种子萌发、幼苗生长和开花结果等不同生长发育阶段的表达模式进行了深入分析。在种子萌发阶段，以芹菜品种‘津南实芹’为材料，挑选饱满、大小均匀的种子，用75%酒精消毒5min，再用无菌水冲洗3-5次，然后将种子均匀放置在铺有两层湿润滤纸的培养皿中，置于25℃恒温培养箱中培养，保持黑暗条件。分别在种子萌发0h、12h、24h、36h、48h和60h时，取10粒种子，迅速放入液氮中速冻，保存于-80℃冰箱，用于后续RNA提取和qRT-PCR分析。利用Trizol试剂法提取不同时间点种子的总RNA，通过核酸蛋白测定仪检测其浓度和纯度，确保A260/A280比值在1.8-2.0之间，A260/A230比值大于2.0。使用反转录试剂盒将总RNA反转录为cDNA，根据已鉴定的芹菜bZIP转录因子基因序列，设计特异性引物用于qRT-PCR分析。结果显示，在种子萌发过程中，AgbZIP1基因的表达量在萌发12h时开始显著上升，为对照（0h）的2.5倍；24h时继续升高，达到对照的4.0倍；36h时略有下降，但仍维持在较高水平，为对照的3.2倍；48h后表达量逐渐降低，60h时接近对照水平。AgbZIP2基因的表达量在萌发24h时显著上调，为对照的3.0倍；36h时达到峰值，是对照的5.0倍；随后逐渐下降，60h时为对照的2.0倍。AgbZIP3基因在种子萌发36h时表达量明显增加，为对照的2.8倍；48h时达到最大值，为对照的4.5倍；60h时略有下降，为对照的3.8倍。这些结果表明，bZIP转录因子在芹菜种子萌发过程中发挥重要作用，可能通过调控相关基因的表达，参与种子的休眠打破、吸水膨胀、胚根萌发等生理过程。在幼苗生长阶段，将芹菜种子播种在装有营养土的育苗盘中，待幼苗长出两片真叶时，选取生长一致的幼苗，移栽至装有Hoagland营养液的水培容器中，置于光照培养箱中培养，培养条件为：光照强度100-150μmol・m⁻²・s⁻¹，光照时间16h/d，温度25℃，相对湿度60%-70%。分别在幼苗生长7d、14d、21d、28d和35d时，取幼苗的根、茎、叶等不同组织，迅速放入液氮中速冻，保存于-80℃冰箱，用于后续分析。提取不同生长时间点幼苗各组织的总RNA，反转录为cDNA后进行qRT-PCR分析。结果表明，在幼苗生长过程中，AgbZIP1基因在根、茎、叶中的表达量均呈现先上升后下降的趋势。在根中，表达量在生长14d时达到峰值，为7d时的3.5倍；在茎中，21d时表达量最高，为7d时的4.0倍；在叶中，28d时表达量最大，为7d时的5.0倍。AgbZIP2基因在根中的表达量在生长21d时显著升高，为7d时的4.2倍；茎中在28d时表达量最高，为7d时的5.5倍；叶中在35d时表达量最大，为7d时的6.0倍。AgbZIP3基因在根中生长28d时表达量明显增加，为7d时的3.8倍；茎中在35d时表达量最高，为7d时的4.5倍；叶中在35d时表达量也达到较高水平，为7d时的4.8倍。这些结果说明，bZIP转录因子在芹菜幼苗生长过程中参与调控不同组织的生长发育，可能通过调节细胞分裂、伸长和分化等过程，影响幼苗的形态建成和器官发育。在开花结果阶段，当芹菜植株生长至一定时期，出现花芽分化时，标记花芽。分别在花芽分化后0d、5d、10d、15d、20d和25d时，取花芽、花、幼果和成熟果实等组织，迅速放入液氮中速冻，保存于-80℃冰箱，用于后续实验。提取不同发育阶段花和果实组织的总RNA，反转录为cDNA后进行qRT-PCR分析。结果显示，在开花过程中，AgbZIP1基因在花芽分化后5d时表达量开始显著上升，为0d