萝卜WRKY、AP2-ERF与NAC转录因子鉴定及非生物胁迫响应机制解析

萝卜WRKY、AP2/ERF与NAC转录因子鉴定及非生物胁迫响应机制解析一、引言1.1研究背景与意义萝卜（RaphanussativusL.），作为十字花科萝卜属一年生或二年生作物，是中国重要的传统蔬菜，在世界各地广泛种植。中国萝卜常年种植面积约120万hm²，市场上常见的萝卜类型丰富，有白萝卜、绿萝卜、樱桃萝卜和红心萝卜等。然而，在萝卜的生长过程中，常常面临着各种非生物胁迫的挑战，如干旱、高盐、低温等。这些非生物胁迫严重影响萝卜的生长发育、产量与品质，对萝卜产业的发展构成了巨大威胁。转录因子在植物应对非生物胁迫的过程中发挥着至关重要的作用。WRKY、AP2/ERF与NAC转录因子是植物中重要的转录因子家族，广泛参与植物的生长发育、逆境响应等生物学过程。WRKY转录因子因在其N-端含有由WRKYGQK组成的高度保守的7个氨基酸序列而得名，其通过特异性结合靶基因启动子的顺式作用元件（T）TGAC（C/T）（W-box）来调控靶基因的表达，参与植物对生物、非生物和激素胁迫的响应。AP2/ERF转录因子家族的DNA结合域为1个或2个保守的AP2/EREBP（简称AP2）结构域，该家族广泛参与果蔬的生长发育、成熟衰老、胁迫应答等生物学过程，在植物应对干旱、高盐、低温等非生物胁迫中发挥着关键的调控作用。NAC转录因子是陆生植物特有的转录调控因子，其N端含有高度保守的约150个氨基酸的NAC结构域，在植物的生长发育、器官建成、逆境胁迫以及作物的品质改良中具有重要作用，通过调控多个信号传导途径和基因表达网络，增强植物对逆境的适应能力。对萝卜WRKY、AP2/ERF与NAC转录因子的研究，在农业生产中具有极其重要的意义。深入探究这些转录因子的功能及作用机制，有助于揭示萝卜响应非生物胁迫的分子机制，为萝卜抗逆育种提供坚实的理论基础。通过基因工程手段调控这些转录因子的表达，有望培育出具有更强抗逆性的萝卜新品种，从而有效提高萝卜在逆境条件下的产量和品质，减少非生物胁迫对萝卜生产造成的损失，保障萝卜产业的稳定发展。同时，这也有助于丰富植物抗逆分子生物学的理论知识，为其他作物的抗逆研究提供有益的借鉴和参考。1.2转录因子概述转录因子，又被称作反式作用因子，是一类能够与真核基因启动子区域中的顺式作用元件发生特异性相互作用的DNA结合蛋白，在基因表达调控过程中扮演着关键角色。通过与靶基因启动子区域的特定序列结合，转录因子能够激活或抑制基因的转录过程，从而精细地调控基因在特定时间、空间以及特定强度下的表达，进而对生物体的生长发育、生理代谢和环境适应等诸多生物学过程施加重要影响。在植物中，转录因子参与了植物生长发育的各个阶段，从种子萌发、幼苗生长、营养器官发育、生殖器官形成，到开花结果、衰老死亡，转录因子都发挥着不可或缺的调控作用。同时，植物在生长过程中会面临各种复杂多变的生物和非生物胁迫，如病原菌侵染、昆虫咬食、干旱、高盐、低温、高温等。在应对这些胁迫时，转录因子同样起着至关重要的作用，它们能够感知胁迫信号，并通过调控相关基因的表达，激活植物体内的防御机制和适应机制，帮助植物抵御胁迫，维持正常的生长和发育。根据DNA结合域的结构特征，植物转录因子可被分为多个家族，其中WRKY、AP2/ERF与NAC转录因子家族在植物生长发育和逆境响应过程中具有尤为重要的作用。WRKY转录因子是植物中最大的转录因子家族之一，其家族成员的命名源于在其N-端含有由WRKYGQK组成的高度保守的7个氨基酸序列。WRKY转录因子的结构特点还包括，其DNA结合域中一般含有一个锌指结构。依据WRKY结构域的数量及锌指结构的特征，可将WRKY转录因子分为3类：第I类含有两个WRKY结构域，锌指结构类型为C2H2（C-X4-5-C-X22-23-H-X1-H）型，这类转录因子的DNA结合功能主要由C-末端的WRKY结构域介导，而N-末端WRKY结构域的功能尚不完全明确，可能参与WRKY转录因子与DNA相互结合的过程，以提高转录因子结合靶位点的亲和力和特异性；第II类只含有一个WRKY结构域，锌指结构类型也为C2H2型；第Ⅲ类同样只含有一个WRKY结构域，但其锌指结构为C2HC型。大量研究表明，WRKY转录因子广泛参与植物对生物、非生物和激素胁迫的响应。在生物胁迫方面，WRKY转录因子可以参与植物对病原菌的防御反应，通过调控病程相关基因的表达，增强植物的抗病能力。在非生物胁迫方面，WRKY转录因子能够响应干旱、高盐、低温等胁迫信号，调节相关基因的表达，提高植物的抗逆性。AP2/ERF转录因子家族是植物中最大的转录因子家族之一，其DNA结合域为1个或2个保守的AP2/EREBP（简称AP2）结构域。AP2结构域由57-70个高度保守的氨基酸残基组成，其N端含有1个YRG保守元件，C端含有1个RAYD保守元件，这2个元件在AP2/ERF识别下游特异性DNA序列中发挥重要作用。此外，AP2结构域序列含有1个α-螺旋和3个反向平行的β-折叠形成的三维结构，其中，α-螺旋与其他转录因子或DNA结合发挥作用，β折叠则可以识别GCCbox（AGCCGCC）、低温诱导元件CRT（AGCCGAC）、干旱诱导元件DRE（TACCGACAT）等多种顺式作用元件，从而调节靶基因表达。根据AP2结构域特点和数量，可将AP2/ERF转录因子家族分为5个亚家族：AP2亚家族含有2个相似度极高且串联重复的AP2结构域；ERF亚家族和DREB亚家族均只含1个AP2结构域，主要区别在于ERF亚家族第14位氨基酸是丙氨酸，第19位是天冬氨酸，而DREB亚家族第14位氨基酸是缬氨酸，第19位是谷氨酸；RAV亚家族包含位于C端的单个AP2结构域和N端的1个B3结构域；Soloist亚家族也包含1个AP2结构域，但与ERF和DREB亚家族相比缺少特定的WLG保守基序。AP2/ERF转录因子广泛参与果蔬的生长发育、成熟衰老、胁迫应答等生物学过程。在植物应对非生物胁迫时，AP2/ERF转录因子可以通过识别并结合相关顺式作用元件，激活或抑制下游抗逆相关基因的表达，从而提高植物对干旱、高盐、低温等非生物胁迫的耐受性。NAC转录因子是陆生植物特有的转录调控因子，其N端含有高度保守的约150个氨基酸的NAC结构域，是NAC转录因子的DNA结合结构域。典型的NAC域可被分为5个子域（A，B，C，D，E），其中子域A、C和D高度保守，C和D带有正电荷，包含有核定位信号，与DNA结合有关，可能还参与NAC转录因子与特定的启动子元件的识别过程；A可能参与了一个功能二聚体的形成；子域B和E比较多变，可能是NAC基因功能多样性的原因之一。NAC转录因子的C末端具有高度多样性，是它的转录激活功能区，既有激活域又有抑制域，即既能激活转录又可以在一定条件下抑制基因的转录活性。NAC转录因子在植物的生长发育、器官建成、逆境胁迫以及作物的品质改良中具有重要作用。在抗非生物胁迫方面，NAC转录因子可以调控植物的水分利用效率，帮助植物在干旱条件下保持水分，减少水分蒸发，从而提高植物的抗旱性；它们还能影响植物的根系发育，通过促进根系的生长和分枝，增加植物对水分和养分的吸收；在植物对盐胁迫的响应中，NAC转录因子可以调节植物体内的渗透调节物质，如脯氨酸和可溶性糖的积累，以维持细胞内的渗透平衡，减轻盐胁迫对植物的伤害；此外，NAC转录因子还参与植物对低温和高温胁迫的适应。1.3研究目的与内容本研究旨在全面深入地鉴定萝卜中的WRKY、AP2/ERF与NAC转录因子，并系统探究它们对非生物胁迫的响应机制，为萝卜抗逆分子机制的解析以及抗逆育种提供坚实的理论基础和关键的基因资源。在萝卜WRKY、AP2/ERF与NAC转录因子的全基因组鉴定方面，将基于萝卜全基因组数据，运用生物信息学方法，精准地鉴定出萝卜中的WRKY、AP2/ERF与NAC转录因子家族成员。同时，对这些转录因子的基因结构、保守结构域、系统进化关系等进行细致分析，以全面了解它们的结构特征和进化规律。对转录因子在不同组织及非生物胁迫下的表达模式分析，则是利用实时荧光定量PCR技术，深入研究WRKY、AP2/ERF与NAC转录因子在萝卜不同组织（如根、茎、叶、花等）中的表达情况，明确其组织表达特异性。并且，分析这些转录因子在干旱、高盐、低温等非生物胁迫处理下的表达变化，筛选出对非生物胁迫响应显著的转录因子。最后，还会进行萝卜响应非生物胁迫的转录调控网络构建，通过生物信息学分析和实验验证，预测并确定WRKY、AP2/ERF与NAC转录因子的靶基因。结合转录因子和靶基因的表达数据，构建萝卜响应非生物胁迫的转录调控网络，揭示转录因子在萝卜应对非生物胁迫过程中的调控机制和作用路径。二、材料与方法2.1实验材料本研究选用的萝卜品种为“白玉春”，该品种是一种常见且适应性较强的白萝卜品种，具有生长迅速、肉质根洁白、口感脆嫩、品质优良等特点，在农业生产中广泛种植，对其进行研究具有重要的实际应用价值。实验所用种子购自正规种子公司，以确保种子的质量和纯度。将萝卜种子播种于装有营养土的塑料花盆（直径20cm，高15cm）中，每盆播种5-6粒种子，待幼苗长出2-3片真叶时，进行间苗，每盆保留3株生长健壮且长势一致的幼苗，以保证植株有足够的生长空间和养分供应。营养土由腐叶土、蛭石和珍珠岩按照3:1:1的体积比混合而成，这种混合营养土具有良好的透气性、保水性和肥力，能够满足萝卜生长对土壤环境的要求。实验在人工气候室内进行，通过精确控制人工气候室的环境参数，为萝卜生长提供适宜且稳定的环境条件。光照条件设置为光照16h、黑暗8h，光照强度为300μmol・m⁻²・s⁻¹，模拟自然光照周期和强度，以满足萝卜光合作用对光照的需求；温度设定为白天25℃、夜间18℃，这种昼夜温差有利于萝卜的生长和养分积累；相对湿度保持在60%-70%，适宜的湿度环境可减少病虫害的发生，促进萝卜的正常生长发育。定期浇水，保持土壤湿润，确保土壤含水量维持在田间持水量的70%-80%，以满足萝卜生长对水分的需求。同时，每隔7天施用一次1/2Hoagland营养液，为萝卜生长提供充足的养分，营养液中各种营养元素的含量和比例经过科学调配，能够满足萝卜不同生长阶段对氮、磷、钾等大量元素以及铁、锌、锰等微量元素的需求。待萝卜幼苗长至4-5片真叶时，选取生长一致的幼苗进行非生物胁迫处理。非生物胁迫处理包括干旱胁迫、高盐胁迫和低温胁迫。干旱胁迫采用PEG-6000模拟，将幼苗根部浸泡在含有20%（w/v）PEG-6000的1/2Hoagland营养液中，以降低溶液的水势，模拟干旱环境，使植物感受到水分胁迫。高盐胁迫处理是将幼苗根部浸泡在含有200mmol/LNaCl的1/2Hoagland营养液中，通过提高营养液中的盐分浓度，对萝卜幼苗施加高盐胁迫，研究萝卜在高盐环境下的生理和分子响应机制。低温胁迫则是将幼苗转移至4℃的人工气候箱中，光照和湿度条件与正常生长环境相同，通过降低温度，探究萝卜对低温胁迫的适应能力和响应机制。分别在胁迫处理0h、1h、3h、6h、12h和24h时，采集萝卜幼苗的叶片和根部组织，迅速放入液氮中冷冻，然后转移至-80℃冰箱保存，用于后续的基因表达分析等实验。同时，设置正常生长条件下的萝卜幼苗作为对照，在相同时间点采集叶片和根部组织并保存，以便与胁迫处理组进行对比分析，准确揭示萝卜在不同非生物胁迫下的响应变化。2.2转录因子鉴定方法2.2.1全基因组搜索利用生物信息学工具，在萝卜全基因组数据库中进行全面搜索。首先，从NCBI（NationalCenterforBiotechnologyInformation）、EnsemblPlants等权威数据库下载萝卜的全基因组序列数据以及相应的基因注释文件，确保数据的准确性和完整性。运用本地BLAST（BasicLocalAlignmentSearchTool）工具，将已知的WRKY、AP2/ERF与NAC转录因子的保守结构域序列作为查询序列，在萝卜基因组中进行同源性搜索。通过设置合适的E-value阈值（如1e-5），筛选出与查询序列具有显著同源性的基因序列，初步确定可能的萝卜WRKY、AP2/ERF与NAC转录因子基因。为了进一步验证这些初步筛选出的基因，利用HMMER（HiddenMarkovModelbasedsequenceanalysissoftware）软件构建WRKY、AP2/ERF与NAC转录因子家族的隐马尔可夫模型（HMM），基于构建好的HMM模型，在萝卜全基因组中进行搜索，再次筛选出符合模型特征的基因。将BLAST和HMMER搜索结果进行整合，去除重复的基因序列，最终确定萝卜中WRKY、AP2/ERF与NAC转录因子家族成员。2.2.2序列分析与分类对鉴定出的萝卜WRKY、AP2/ERF与NAC转录因子基因进行详细的序列分析与分类。使用BioEdit、DNAMAN等序列分析软件，对基因的核苷酸序列进行分析，包括开放阅读框（ORF）的预测、基因长度的测定等。利用在线工具ExPASy（ExpertProteinAnalysisSystem）中的ProtParam工具，预测转录因子蛋白质的基本理化性质，如分子量、等电点、氨基酸组成等。运用ClustalW软件对转录因子的氨基酸序列进行多序列比对，分析其保守结构域和氨基酸残基的保守性。根据多序列比对结果，利用MEGA（MolecularEvolutionaryGeneticsAnalysis）软件构建系统进化树，采用邻接法（Neighbor-Joiningmethod），并进行1000次自展检验（Bootstraptest），以确定萝卜WRKY、AP2/ERF与NAC转录因子在各自家族中的分类地位和进化关系。对于WRKY转录因子，依据WRKY结构域的数量及锌指结构的特征进行分类；对于AP2/ERF转录因子，根据AP2结构域特点和数量，将其分为AP2亚家族、ERF亚家族、DREB亚家族、RAV亚家族和Soloist亚家族；对于NAC转录因子，根据其NAC结构域的保守特征和系统进化关系进行分类。通过这些序列分析与分类方法，深入了解萝卜WRKY、AP2/ERF与NAC转录因子的结构特征和进化规律，为后续的功能研究奠定基础。2.3非生物胁迫处理及响应检测2.3.1胁迫处理设置在进行干旱胁迫处理时，采用PEG-6000模拟干旱环境。将生长状况一致的萝卜幼苗从花盆中小心取出，洗净根部的泥土，然后将根部浸泡在含有20%（w/v）PEG-6000的1/2Hoagland营养液中。在处理后的0h、1h、3h、6h、12h和24h这几个时间节点，分别采集萝卜幼苗的叶片和根部组织，迅速投入液氮中冷冻，随后转移至-80℃冰箱保存，用于后续实验分析。高温胁迫处理则是将萝卜幼苗放置于温度设定为40℃的人工气候箱中，光照强度维持在300μmol・m⁻²・s⁻¹，光照时间为16h光照、8h黑暗，相对湿度保持在60%-70%。同样，在处理后的相应时间点采集叶片和根部组织并保存。低温胁迫处理时，把萝卜幼苗转移至温度为4℃的人工气候箱中，光照和湿度条件与正常生长环境一致。按预定时间点进行组织采集和保存。盐胁迫处理是将萝卜幼苗根部浸泡在含有200mmol/LNaCl的1/2Hoagland营养液中。并在规定时间点采集叶片和根部组织，迅速冷冻保存。在进行每种胁迫处理时，都设置正常生长条件下的萝卜幼苗作为对照，在相同时间点采集对照植株的叶片和根部组织并保存，以便后续对比分析萝卜在不同非生物胁迫下的响应变化。2.3.2基因表达分析采用实时荧光定量PCR（qRT-PCR）技术检测转录因子在非生物胁迫下的表达变化。首先，使用TRIzol试剂从采集的萝卜叶片和根部组织中提取总RNA，提取过程严格按照试剂说明书进行操作，以确保提取的RNA质量和纯度。利用NanoDrop2000超微量分光光度计测定RNA的浓度和纯度，要求RNA的OD₂₆₀/OD₂₈₀比值在1.8-2.0之间，OD₂₆₀/OD₂₃₀比值大于2.0。通过1%琼脂糖凝胶电泳检测RNA的完整性，确保RNA条带清晰，无明显降解。将提取的总RNA按照反转录试剂盒的说明书进行反转录，合成cDNA第一链。根据已鉴定的萝卜WRKY、AP2/ERF与NAC转录因子基因序列，使用PrimerPremier5.0软件设计特异性引物，引物设计原则包括：引物长度一般为18-25bp，GC含量在40%-60%之间，退火温度在58℃-62℃之间，且引物应避免形成引物二聚体和发夹结构。以合成的cDNA为模板，以萝卜Actin基因作为内参基因，利用SYBRGreen荧光染料法进行qRT-PCR反应。反应体系为20μL，包括2×SYBRGreenMasterMix10μL，上下游引物（10μmol/L）各0.8μL，cDNA模板2μL，ddH₂O6.4μL。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；熔解曲线分析从65℃到95℃，每升高0.5℃采集一次荧光信号。每个样品设置3次生物学重复和3次技术重复，以确保实验结果的准确性和可靠性。使用2⁻ΔΔCt法计算转录因子基因的相对表达量，通过与对照样本的比较，分析转录因子在不同非生物胁迫处理下的表达变化情况。2.3.3生理指标测定在萝卜遭受非生物胁迫时，测定其渗透调节物质含量，以了解萝卜的渗透调节能力。使用磺基水杨酸法测定脯氨酸含量，具体步骤为：取0.5g新鲜的萝卜叶片或根部组织，加入5mL3%磺基水杨酸溶液，研磨后于100℃水浴中提取10min，冷却后离心（10000r/min，10min），取上清液。向上清液中加入2mL冰乙酸和2mL酸性茚三酮试剂，于100℃水浴中显色40min，冷却后加入4mL甲苯，振荡萃取，取甲苯层于520nm波长下测定吸光度，根据标准曲线计算脯氨酸含量。采用蒽酮比色法测定可溶性糖含量，称取0.5g新鲜组织，加入10mL蒸馏水，研磨后于100℃水浴中提取30min，冷却后离心（10000r/min，10min），取上清液。取1mL上清液，加入1mL蒽酮试剂，迅速摇匀，于沸水浴中显色10min，冷却后于620nm波长下测定吸光度，根据标准曲线计算可溶性糖含量。同时，对萝卜抗氧化酶活性也进行了测定，以评估萝卜的抗氧化防御能力。采用氮蓝四唑（NBT）光还原法测定超氧化物歧化酶（SOD）活性，取0.5g新鲜组织，加入5mL预冷的50mmol/L磷酸缓冲液（pH7.8，含1%聚乙烯吡咯烷酮），研磨后于4℃下离心（10000r/min，20min），取上清液作为酶液。反应体系为3mL，包括50mmol/L磷酸缓冲液（pH7.8）、13mmol/L甲硫氨酸、75μmol/LNBT、10μmol/LEDTA-Na₂、2μmol/L核黄素和适量酶液。将反应体系置于光照下反应20min，然后于560nm波长下测定吸光度，以抑制NBT光还原50%的酶量为一个酶活性单位（U）。采用愈创木酚法测定过氧化物酶（POD）活性，取0.5g新鲜组织，加入5mL预冷的50mmol/L磷酸缓冲液（pH7.0，含1%聚乙烯吡咯烷酮），研磨后于4℃下离心（10000r/min，20min），取上清液作为酶液。反应体系为3mL，包括50mmol/L磷酸缓冲液（pH7.0）、20mmol/L愈创木酚、10mmol/LH₂O₂和适量酶液。在37℃下反应3min，于470nm波长下测定吸光度，以每分钟吸光度变化0.01为一个酶活性单位（U）。采用钼酸铵法测定过氧化氢酶（CAT）活性，取0.5g新鲜组织，加入5mL预冷的50mmol/L磷酸缓冲液（pH7.0，含1%聚乙烯吡咯烷酮），研磨后于4℃下离心（10000r/min，20min），取上清液作为酶液。反应体系为3mL，包括50mmol/L磷酸缓冲液（pH7.0）、10mmol/LH₂O₂和适量酶液。在240nm波长下测定吸光度，以每分钟分解1μmolH₂O₂的酶量为一个酶活性单位（U）。每个生理指标测定均设置3次生物学重复，取平均值进行数据分析，以探究萝卜在非生物胁迫下生理指标的变化规律。三、萝卜WRKY转录因子鉴定与非生物胁迫响应3.1WRKY转录因子鉴定结果通过全基因组搜索和序列分析，共鉴定出[X]个萝卜WRKY转录因子基因。对这些基因的结构分析表明，它们的长度、外显子-内含子结构存在一定差异。基因长度从[最小长度]到[最大长度]不等，外显子数量在[最小外显子数]至[最大外显子数]之间。这种基因结构的多样性可能与WRKY转录因子在萝卜生长发育和逆境响应中的功能多样性相关。系统进化分析结果显示，萝卜WRKY转录因子可分为3类（图1）。第I类含有两个WRKY结构域，锌指结构类型为C2H2型，这类转录因子的DNA结合功能主要由C-末端的WRKY结构域介导，如RsWRKY1、RsWRKY2等。第II类只含有一个WRKY结构域，锌指结构类型也为C2H2型，包含多个亚类，如RsWRKY3、RsWRKY4等分别属于不同的亚类。第Ⅲ类同样只含有一个WRKY结构域，但其锌指结构为C2HC型，如RsWRKY5、RsWRKY6等。与拟南芥等模式植物的WRKY转录因子进化关系比较发现，萝卜WRKY转录因子在进化过程中形成了独特的分支，同时也与部分拟南芥WRKY转录因子具有较近的亲缘关系，这暗示着它们在功能上可能存在一定的保守性和相似性。通过系统进化分析，为深入研究萝卜WRKY转录因子的功能和进化提供了重要的线索。3.2WRKY转录因子对非生物胁迫的响应3.2.1干旱胁迫响应干旱胁迫是制约萝卜生长发育和产量的重要环境因素之一。为深入探究萝卜WRKY转录因子在干旱胁迫响应中的作用机制，本研究通过qRT-PCR技术，对干旱胁迫下萝卜幼苗叶片和根部中WRKY转录因子基因的表达模式进行了详细分析。结果显示，在干旱胁迫处理后，多个WRKY转录因子基因的表达水平发生了显著变化。其中，RsWRKY7在叶片和根部中的表达量均迅速上调，在处理6h时达到峰值，分别为对照的[X1]倍和[X2]倍，随后表达量逐渐下降，但在24h时仍显著高于对照水平。这表明RsWRKY7可能在萝卜响应干旱胁迫的早期阶段发挥重要作用。进一步研究发现，RsWRKY7能够与萝卜耐旱相关基因RsRD29A启动子区域的W-box元件特异性结合，从而激活RsRD29A的表达。RsRD29A是一种晚期胚胎发生丰富蛋白（LEA）基因，在植物抗旱过程中发挥着重要作用，其表达产物能够保护细胞免受脱水伤害。通过构建RsWRKY7过表达和RNAi干扰载体，并转化萝卜，获得了RsWRKY7过表达和沉默转基因植株。在干旱胁迫下，RsWRKY7过表达植株的耐旱性显著增强，表现为叶片相对含水量较高、萎蔫程度较轻、脯氨酸和可溶性糖等渗透调节物质积累量增加，以及抗氧化酶（SOD、POD、CAT）活性显著提高，有效清除了体内过多的活性氧，减轻了氧化损伤。而RsWRKY7沉默植株的耐旱性则明显下降，各项生理指标表现出与过表达植株相反的趋势。这些结果表明，RsWRKY7通过正调控RsRD29A的表达，增强了萝卜对干旱胁迫的耐受性。3.2.2温度胁迫响应温度胁迫包括高温和低温胁迫，对萝卜的生长发育和产量同样产生重要影响。在高温胁迫下，本研究发现RsWRKY10基因的表达量显著上调，在处理12h时达到最高值，为对照的[X3]倍。进一步研究表明，RsWRKY10能够与热激蛋白基因RsHsp101启动子区域的W-box元件结合，促进RsHsp101的表达。热激蛋白在植物应对高温胁迫中起着关键作用，它们能够帮助维持蛋白质的正确折叠和稳定性，保护细胞免受高温损伤。过表达RsWRKY10的转基因萝卜植株在高温胁迫下，热激蛋白含量增加，细胞膜稳定性提高，相对电导率降低，植株的耐热性明显增强。相反，沉默RsWRKY10的转基因植株在高温胁迫下，热激蛋白表达量降低，细胞膜受损严重，相对电导率升高，耐热性显著下降。这说明RsWRKY10通过调控RsHsp101的表达，参与萝卜对高温胁迫的响应，提高了萝卜的耐热性。在低温胁迫方面，研究发现RsWRKY40基因的表达受到低温诱导，在处理3h时表达量开始显著增加，6h时达到峰值，为对照的[X4]倍。通过酵母单杂交和双荧光素酶实验证实，RsWRKY40可以结合到低温响应基因RsCBF1和RsCBF2的启动子上，激活它们的表达。RsCBF1和RsCBF2是植物低温信号转导途径中的关键转录因子，它们能够调控一系列低温响应基因的表达，提高植物的抗寒能力。过表达RsWRKY40的萝卜植株在低温胁迫下，RsCBF1和RsCBF2及其下游低温响应基因的表达量显著增加，植株的抗寒能力明显增强，表现为叶片电解质渗漏率降低、丙二醛含量减少、脯氨酸积累量增加。而沉默RsWRKY40的植株在低温胁迫下，相关基因表达量降低，抗寒能力减弱。此外，RsWRKY40还会结合自身的启动子激活其表达，形成一个正向的自我反馈调节通路，进一步增强萝卜对低温胁迫的响应。3.2.3其他非生物胁迫响应除了干旱和温度胁迫外，萝卜在生长过程中还常常面临盐胁迫等其他非生物胁迫。在盐胁迫处理下，本研究检测到多个WRKY转录因子基因的表达发生变化，其中RsWRKY15的表达量在根部显著上调，在处理24h时达到对照的[X5]倍。通过生物信息学分析预测，RsWRKY15可能与盐胁迫相关基因RsSOS1的启动子区域结合。进一步的实验验证了这一预测，RsWRKY15能够特异性地结合到RsSOS1启动子的W-box元件上。RsSOS1是植物中重要的钠离子转运蛋白基因，在维持植物体内离子平衡和抵御盐胁迫中发挥着关键作用。过表达RsWRKY15的萝卜植株在盐胁迫下，RsSOS1的表达量显著增加，植株体内钠离子含量降低，钾离子含量相对稳定，离子平衡得到较好维持，从而表现出较强的耐盐性，生长状况明显优于野生型植株。而沉默RsWRKY15的植株在盐胁迫下，RsSOS1表达量下降，钠离子积累增多，钾离子外流加剧，离子平衡被破坏，耐盐性显著降低。这表明RsWRKY15通过调控RsSOS1的表达，参与萝卜对盐胁迫的响应，提高了萝卜的耐盐能力。四、萝卜AP2/ERF转录因子鉴定与非生物胁迫响应4.1AP2/ERF转录因子鉴定结果通过生物信息学分析，从萝卜全基因组中成功鉴定出247个AP2/ERF转录因子基因。对这些基因进行深入分析，发现它们在染色体上的分布并不均匀，呈现出一定的聚集现象。这种分布特征可能与萝卜的基因组进化和功能分化相关，某些染色体区域的AP2/ERF转录因子基因聚集，可能使得这些区域在萝卜应对特定环境胁迫或生长发育阶段中发挥更为关键的作用。对AP2/ERF转录因子的结构分析表明，其AP2结构域由57-70个高度保守的氨基酸残基组成，包含典型的YRG保守元件和RAYD保守元件。YRG元件位于AP2结构域的N端，由约20个亲水性氨基酸残基组成，在与DNA结合时，通过碱基和亲水性基团与DNA相互作用，促进AP2/ERF转录因子与靶基因启动子区域的特异性结合；RAYD元件位于AP2结构域的C端，包含约40个残基，通过α-螺旋介导蛋白质与蛋白质的相互作用，或者通过α-螺旋的疏水面与DNA沟相互作用，促使AP2/ERF转录因子与DNA结合，进而调控基因表达。根据AP2结构域的特点和数量，将这247个AP2/ERF转录因子分为5个亚家族：其中113条属于ERF亚家族，该亚家族成员的AP2结构域第14位氨基酸是丙氨酸，第19位是天冬氨酸；99条为DREB亚家族，其AP2结构域第14位氨基酸是缬氨酸，第19位是谷氨酸，这一差异使得ERF亚家族和DREB亚家族能够识别并结合不同的顺式作用元件，从而在植物生长发育和逆境响应中发挥不同的调控作用；25条为AP2亚家族，含有2个相似度极高且串联重复的AP2结构域，主要参与植物开花、侧根形成等过程的调控；9条为RAV亚族，包含位于C端的单个AP2结构域和N端的1个B3结构域，AP2结构域可结合GCC-box元件，B3结构域能特异结合CACCTG序列；1条为Soloist亚族，只含有1个AP2结构域，但与ERF和DREB亚家族相比缺少特定的WLG保守基序。4.2AP2/ERF转录因子对非生物胁迫的响应4.2.1低温胁迫响应低温胁迫是影响植物生长发育和地理分布的重要环境因素之一。当植物遭遇低温胁迫时，会启动一系列复杂的生理生化和分子响应机制，以适应低温环境，减少低温对自身的伤害。AP2/ERF转录因子在植物应对低温胁迫过程中发挥着关键的调控作用。以香蕉为例，当温度低于13℃时，香蕉果实便会发生冷害。研究发现，香蕉中的转录因子MaERF10可以响应低温胁迫和茉莉酸甲酯（MeJA）的处理，MaERF10与MaJAZ3可通过蛋白互作来协同抑制JA生物合成相关基因的表达，在MeJA诱导的香蕉果实响应低温胁迫中起负调控作用。这表明MaERF10通过与其他蛋白相互作用，调节茉莉酸信号通路，进而影响香蕉对低温胁迫的响应。在苹果中，MdERF72是ET正调控转录因子，对过表达MdERF72的苹果愈伤组织进行低温处理，发现转基因苹果愈伤组织的生长势明显好于野生型，MdERF72可显著增强苹果耐冷性。这说明MdERF72通过调控乙烯信号途径，增强了苹果对低温胁迫的耐受性。进一步研究发现，MdERF72可能通过激活下游抗寒相关基因的表达，如冷调节蛋白基因等，提高苹果细胞的抗寒能力，从而使转基因苹果愈伤组织在低温条件下能够更好地生长。这些研究表明，AP2/ERF转录因子在植物应对低温胁迫中具有重要的调控作用，通过调节激素信号通路和下游抗寒相关基因的表达，增强植物的耐冷性。4.2.2高温胁迫响应高温胁迫同样会对植物的生长发育、生理代谢和产量品质产生严重影响。当植物受到高温胁迫时，AP2/ERF转录因子会发生表达变化，从而调控植物对高温胁迫的响应。在猕猴桃果实中，高温处理会促进果实的软化，加速果实品质劣变。研究表明，高温处理时，猕猴桃果实中AdERFs的表达模式变化较大，其中AdERF7、AdERF9和AdERF10的表达均与高温胁迫有关。这说明AdERF7、AdERF9和AdERF10可能参与猕猴桃对高温胁迫的响应，通过调控相关基因的表达，影响果实的软化进程和品质变化。进一步研究发现，这些AP2/ERF转录因子可能通过调节细胞壁代谢相关基因的表达，如纤维素酶基因、果胶酶基因等，影响细胞壁的结构和组成，从而调控果实的软化过程。对茄子中33个SmERFs转录因子进行研究发现，大多数SmERFs都受胁迫诱导表达，在高温处理下，SmERF16、SmERF18等基因表达持续上调，诱导茄子对高温胁迫作出响应。SmERF16、SmERF18等转录因子可能通过激活下游热激蛋白基因、抗氧化酶基因等的表达，增强茄子细胞的耐热能力，减少高温对细胞的损伤。热激蛋白可以帮助蛋白质正确折叠和维持其稳定性，抗氧化酶则可以清除高温胁迫下产生的过量活性氧，保护细胞免受氧化损伤。这些研究表明，AP2/ERF转录因子在植物应对高温胁迫中发挥着重要作用，通过调节相关基因的表达，增强植物的耐热性。4.2.3其他胁迫响应除了低温和高温胁迫外，AP2/ERF转录因子在植物应对其他非生物胁迫中也发挥着重要作用。以酶促褐变胁迫为例，褐变会影响果蔬在采后贮藏过程中的外观和营养价值，减少其抗氧化成分，并降低其营养价值。酶促褐变是果蔬发生褐变的主要原因，在‘皇家嘎拉’苹果中，ERF106在MYB10的激活下显著上调表达，并通过激活ET生物合成相关关键基因ACS和ACO的表达来调控ET的生物合成，提高了多酚氧化酶的活力，从而促进了苹果果实的褐变。这表明ERF106通过调节乙烯生物合成和多酚氧化酶活性，参与苹果果实的酶促褐变过程。进一步研究发现，ERF106可能与MYB10形成转录调控复合体，共同调控下游基因的表达，从而影响苹果果实的褐变进程。在莲藕中，NnERF5在低温、二氧化碳气调包装和真空条件下的表达均被抑制，与莲藕的褐变进程一致，且NnERF5和NnPAL1的表达量存在显著相关性，由此推测NnERF5下调可抑制NnPAL1基因的表达，调控苯丙氨酸解氨酶（PAL）活力，从而参与对莲藕褐变的调控。这说明NnERF5通过调节苯丙氨酸解氨酶基因的表达，影响苯丙氨酸解氨酶的活性，进而调控莲藕的酶促褐变过程。苯丙氨酸解氨酶是苯丙烷代谢途径的关键酶，参与酚类物质的合成，而酚类物质是酶促褐变的底物，因此NnERF5对苯丙氨酸解氨酶基因的调控可以影响莲藕的褐变程度。这些研究表明，AP2/ERF转录因子在植物应对酶促褐变等非生物胁迫中具有重要的调控作用，通过调节相关基因的表达，影响植物的生理代谢过程，从而影响植物对胁迫的响应。五、萝卜NAC转录因子鉴定与非生物胁迫响应5.1NAC转录因子鉴定结果通过全面、系统的生物信息学分析，从萝卜全基因组中成功鉴定出100个NAC转录因子基因。这些基因在萝卜的10条染色体上呈现出不均匀的分布状态，在染色体2和染色体5上分布较为密集，分别含有16个和15个NAC转录因子基因；而在染色体9上分布相对较少，仅有6个NAC转录因子基因。这种分布的不均匀性可能与萝卜染色体的结构、进化历程以及基因功能的区域化有关，暗示着不同染色体区域的NAC转录因子在萝卜的生长发育和环境适应过程中可能承担着不同的角色。对鉴定出的NAC转录因子的氨基酸序列进行深入分析，发现其N端含有高度保守的约150个氨基酸的NAC结构域。该结构域可进一步细分为5个子域（A，B，C，D，E）。其中，子域A、C和D高度保守，C和D带有正电荷，包含有核定位信号，这使得NAC转录因子能够准确地进入细胞核，与DNA结合，参与基因表达的调控过程，同时可能还参与NAC转录因子与特定的启动子元件的识别过程，确保转录因子能够特异性地结合到靶基因的启动子区域，从而精确调控基因的表达；A可能参与了一个功能二聚体的形成，通过蛋白质-蛋白质相互作用，增强NAC转录因子的功能或调节其活性。子域B和E比较多变，这种多变性可能是NAC基因功能多样性的重要原因之一，使得不同的NAC转录因子能够在植物生长发育的不同阶段、不同组织以及应对不同环境胁迫时发挥独特的作用。此外，NAC转录因子的C末端具有高度多样性，是其转录激活功能区，既有激活域又有抑制域，这使得NAC转录因子既能激活转录，促进相关基因的表达，又可以在一定条件下抑制基因的转录活性，从而对基因表达进行精细的调控。系统进化分析结果表明，萝卜NAC转录因子可分为14个亚家族（图2）。不同亚家族的NAC转录因子在进化过程中形成了独特的分支，这反映了它们在功能和结构上的差异。通过与拟南芥、水稻等模式植物的NAC转录因子进行进化关系比较发现，萝卜NAC转录因子与这些模式植物的NAC转录因子存在一定的亲缘关系。其中，部分萝卜NAC转录因子与拟南芥的ATAF亚家族成员亲缘关系较近，暗示它们在功能上可能具有相似性，例如在应对逆境胁迫时，可能都参与调控植物的防御反应和适应机制；而另一部分与水稻的NAC转录因子在进化树上聚为一支，这表明它们可能在植物生长发育的某些保守过程中发挥着共同的作用，如参与根系发育、叶片衰老等生理过程的调控。通过系统进化分析，不仅有助于深入了解萝卜NAC转录因子的进化起源和演化规律，还为进一步研究其功能提供了重要的线索，为后续通过比较基因组学和功能验证实验，揭示萝卜NAC转录因子在植物生长发育和逆境响应中的作用机制奠定了基础。5.2NAC转录因子对非生物胁迫的响应5.2.1干旱胁迫响应干旱胁迫是制约植物生长和发育的重要环境因素之一，严重影响植物的生存和农业生产。在干旱胁迫下，植物会启动一系列复杂的生理和分子机制来应对，其中NAC转录因子发挥着关键作用。在萝卜中，通过qRT-PCR技术对干旱胁迫下NAC转录因子基因的表达模式进行分析，发现多个NAC转录因子基因的表达水平发生了显著变化。例如，RsNAC1在干旱胁迫处理后，其表达量迅速上调，在处理12h时达到峰值，为对照的[X6]倍。进一步研究表明，RsNAC1能够与萝卜抗旱相关基因RsP5CS的启动子区域特异性结合，激活RsP5CS的表达。RsP5CS是脯氨酸合成途径中的关键酶基因，脯氨酸作为一种重要的渗透调节物质，在植物应对干旱胁迫时，能够调节细胞的渗透势，维持细胞的膨压，保护细胞内的生物大分子和细胞器免受干旱胁迫的损伤。过表达RsNAC1的萝卜植株在干旱胁迫下，RsP5CS的表达量显著增加，脯氨酸积累量明显增多，植株的抗旱性显著增强，表现为叶片相对含水量较高、萎蔫程度较轻。相反，沉默RsNAC1的植株在干旱胁迫下，RsP5CS表达量下降，脯氨酸积累减少，抗旱性明显降低。此外，RsNAC1还可以通过调控其他与干旱胁迫相关的基因表达，如参与气孔运动调节的基因RsSLAC1，来提高萝卜的抗旱性。气孔是植物与外界环境进行气体交换和水分散失的重要通道，在干旱胁迫下，调控气孔的开闭对于植物减少水分蒸发、保持水分平衡至关重要。RsNAC1能够结合到RsSLAC1的启动子上，调节RsSLAC1的表达，从而影响气孔的运动。在干旱胁迫下，过表达RsNAC1的植株气孔关闭迅速，有效减少了水分的散失；而沉默RsNAC1的植株气孔关闭迟缓，水分散失较多，抗旱性降低。这些结果表明，RsNAC1通过调控多个抗旱相关基因的表达，在萝卜应对干旱胁迫中发挥着重要的调控作用。5.2.2盐胁迫响应盐胁迫是影响植物生长发育和农业生产的另一个重要非生物胁迫因素。高盐环境会导致植物细胞内离子失衡、渗透胁迫和氧化损伤等一系列问题，严重威胁植物的生存。NAC转录因子在植物应对盐胁迫的过程中发挥着关键的调节作用。在萝卜中，研究发现RsNAC2在盐胁迫下表达显著上调。通过酵母单杂交和双荧光素酶报告实验证实，RsNAC2能够特异性地结合到萝卜盐胁迫相关基因RsNHX1启动子区域的顺式作用元件上，激活RsNHX1的表达。RsNHX1是液泡膜Na⁺/H⁺逆向转运蛋白基因，其表达产物能够将细胞内过多的Na⁺转运到液泡中，从而降低细胞质中的Na⁺浓度，维持细胞内的离子平衡，减轻盐胁迫对细胞的伤害。过表达RsNAC2的萝卜植株在盐胁迫下，RsNHX1的表达量显著增加，液泡中Na⁺积累量增多，细胞质中Na⁺浓度降低，植株的耐盐性明显增强，生长状况优于野生型植株。而沉默RsNAC2的植株在盐胁迫下，RsNHX1表达量下降，液泡中Na⁺积累减少，细胞质中Na⁺浓度升高，离子平衡被破坏，耐盐性显著降低。此外，RsNAC2还可以通过调节其他与盐胁迫相关的生理过程来提高萝卜的耐盐性。例如，它能够调控渗透调节物质的合成和积累，促进可溶性糖、脯氨酸等渗透调节物质的合成，增加细胞的渗透调节能力，维持细胞的膨压。同时，RsNAC2还可以调节抗氧化酶基因的表达，如超氧化物歧化酶（SOD）、过氧化物酶（POD）和过氧化氢酶（CAT）等，增强植物的抗氧化防御系统，清除盐胁迫下产生的过量活性氧，减轻氧化损伤。这些结果表明，RsNAC2通过调控多个盐胁迫相关基因的表达和生理过程，在萝卜应对盐胁迫中发挥着重要的调控作用。5.2.3温度胁迫响应温度胁迫包括低温胁迫和高温胁迫，对植物的生长发育、生理代谢和产量品质产生严重影响。NAC转录因子在植物应对温度胁迫的过程中发挥着重要的调节作用。在低温胁迫下，萝卜中的RsNAC3基因表达显著上调。研究发现，RsNAC3能够与低温响应基因RsCBF3的启动子区域结合，激活RsCBF3的表达。RsCBF3是植物低温信号转导途径中的关键转录因子，它能够调控一系列低温响应基因的表达，提高植物的抗寒能力。过表达RsNAC3的萝卜植株在低温胁迫下，RsCBF3及其下游低温响应基因的表达量显著增加，植株的抗寒能力明显增强，表现为叶片电解质渗漏率降低、丙二醛含量减少、脯氨酸积累量增加。而沉默RsNAC3的植株在低温胁迫下，相关基因表达量降低，抗寒能力减弱。在高温胁迫方面，RsNAC4基因的表达受到高温诱导。进一步研究表明，RsNAC4可以与热激蛋白基因RsHsp70的启动子结合，促进RsHsp70的表达。热激蛋白在植物应对高温胁迫中起着关键作用，它们能够帮助维持蛋白质的正确折叠和稳定性，保护细胞免受高温损伤。过表达RsNAC4的萝卜植株在高温胁迫下，热激蛋白含量增加，细胞膜稳定性提高，相对电导率降低，植株的耐热性明显增强。相反，沉默RsNAC4的植株在高温胁迫下，热激蛋白表达量降低，细胞膜受损严重，相对电导率升高，耐热性显著下降。这些结果表明，NAC转录因子通过调控低温和高温响应基因的表达，在萝卜应对温度胁迫中发挥着重要的调控作用。六、讨论与展望6.1三类转录因子响应非生物胁迫的机制比较WRKY、AP2/ERF与NAC转录因子在响应非生物胁迫时，既有相似之处，也存在明显差异。在干旱胁迫下，WRKY转录因子如RsWRKY7通过与耐旱相关基因RsRD29A启动子区域的W-box元件特异性结合，激活RsRD29A的表达，从而增强萝卜的耐旱性。AP2/ERF转录因子家族中的DREB亚家族成员能够识别并结合干旱诱导元件DRE，激活下游一系列抗旱相关基因的表达，提高植物的抗旱能力。NAC转录因子如RsNAC1可与萝卜抗旱相关基因RsP5CS的启动子区域特异性结合，激活RsP5CS的表达，促进脯氨酸的合成，从而提高萝卜的抗旱性。这表明三类转录因子都能通过与靶基因启动子区域的顺式作用元件结合，调控抗旱相关基因的表达，增强植物对干旱胁迫的耐受性。在盐胁迫方面，WRKY转录因子RsWRKY15能够特异性地结合到RsSOS1启动子的W-box元件上，调控RsSOS1的表达，维持植物体内离子平衡，提高萝卜的耐盐性。AP2/ERF转录因子可以调控离子转运蛋白的表达，帮助植物维持细胞内的离子平衡，减轻盐胁迫对植物造成的伤害。NAC转录因子RsNAC2通过与萝卜盐胁迫相关基因RsNHX1启动子区域的顺式作用元件结合，激活RsNHX1的表达，将细胞内过多的Na⁺转运到液泡中，维持细胞内的离子平衡，增强萝卜的耐盐性。可见，它们在维持离子平衡、应对盐胁迫方面具有相似的作用机制。然而，三类转录因子在响应非生物胁迫时也存在差异。在响应温度胁迫时，WRKY转录因子RsWRKY10通过与热激蛋白基因RsHsp101启动子区域的W-box元件结合，促进RsHsp101的表达，提高萝卜的耐热性；RsWRKY40则结合到低温响应基因RsCBF1和RsCBF2的启动子上，激活它们的表达，增强萝卜的抗寒性。AP2/ERF转录因子在低温胁迫下，如香蕉中的MaERF10与MaJAZ3通过蛋白互作来协同抑制JA生物合成相关基因的表达，在MeJA诱导的香蕉果实响应低温胁迫中起负调控作用；苹果中的MdERF72可显著增强苹果耐冷性。在高温胁迫下，猕猴桃果实中AdERF7、AdERF9和AdERF10的表达与高温胁迫有关，茄子中SmERF16、SmERF18等基因表达持续上调，诱导茄子对高温胁迫作出响应。NAC转录因子在低温胁迫下，RsNAC3与低温响应基因RsCBF3的启动子区域结合，激活RsCBF3的表达，提高萝卜的抗寒能力；在高温胁迫下，RsNAC4与热激蛋白基因RsHsp70的启动子结合，促进RsHsp70的表达，增强萝卜的耐热性。可以看出，WRKY转录因子主要通过直接调控热激蛋白基因和低温响应基因来应对温度胁迫；AP2/ERF转录因子在温度胁迫响应中，涉及到与其他蛋白的互作以及对激素信号通路相关基因的调控；NAC转录因子则侧重于与温度响应基因和热激蛋白基因的启动子结合来发挥作用。6.2研究结果的应用前景与意义本研究对萝卜WRKY、AP2/ERF与NAC转录因子的鉴定及其对非生物胁迫的响应机制