不结球白菜冷诱导差异表达基因解析：筛选、功能与耐寒调控机制

不结球白菜冷诱导差异表达基因解析：筛选、功能与耐寒调控机制一、引言1.1研究背景不结球白菜（Brassicacampestrisssp.chinensis），别名小白菜、青菜等，是十字花科芸薹属芸薹种的一个亚种。其原产于中国，栽培历史源远流长，是深受广大民众喜爱的蔬菜品种。不结球白菜生长周期短，通常20-60天即可收获，这使得在有限的土地资源上能够实现多次种植，提高了土地利用率。而且其适应性极为广泛，无论是在气候温暖湿润的南方，还是四季分明的北方地区，都能良好生长。在我国中部及以南地区，不结球白菜更是周年栽培种植，占据了春、秋、冬菜播种面积的30％-50％，成为全年播种面积最大的蔬菜之一，在我国的蔬菜周年供应体系中占据着举足轻重的地位。低温是影响不结球白菜生长发育的重要环境因子之一。当不结球白菜遭遇低温胁迫时，其生理生化过程会发生一系列显著变化。低温会导致植物细胞膜的流动性降低，膜脂发生相变，从而破坏细胞膜的结构和功能，使得细胞内的物质渗漏，影响细胞的正常代谢。低温还会抑制植物体内许多酶的活性，这些酶参与光合作用、呼吸作用、物质合成与分解等重要生理过程，酶活性的降低直接限制了植物的正常生长发育。在光合作用方面，低温会使光合色素的合成受阻，光系统Ⅱ的活性降低，影响光能的吸收、传递和转化，导致光合速率下降，植物无法积累足够的光合产物，进而影响植株的生长和产量。从对产量的影响来看，低温胁迫下，不结球白菜的生长速度明显减缓，植株矮小，叶片数减少，叶片变薄且发黄，严重时甚至会出现冻害，导致整株死亡，从而使产量大幅下降。相关研究数据表明，在轻度低温胁迫下，不结球白菜的产量可能会降低10％-30％，而在重度低温胁迫下，产量损失可达50％以上。对于品质而言，低温会导致不结球白菜的口感变差，纤维含量增加，维生素C、可溶性糖等营养物质的含量降低，从而降低其食用价值和商品价值。研究不结球白菜的冷诱导基因具有至关重要的意义。通过深入研究这些基因，我们能够揭示不结球白菜响应低温胁迫的分子机制，了解植物在低温环境下如何感知低温信号、传递信号以及启动一系列的抗寒生理反应，这有助于我们从分子层面深入认识植物与环境的相互作用关系。对冷诱导基因的研究为不结球白菜的遗传改良提供了坚实的理论基础。利用现代生物技术，如基因编辑、转基因等手段，可以对不结球白菜的冷诱导基因进行调控，从而培育出具有更强耐寒性的新品种，提高不结球白菜在低温环境下的产量和品质，保障蔬菜的稳定供应，满足人们对优质蔬菜的需求。1.2国内外研究现状国内外对于不结球白菜冷诱导基因的研究已取得了一定的进展。在国外，一些研究聚焦于模式植物拟南芥的冷胁迫应答机制，为不结球白菜的相关研究提供了重要的理论参考。通过对拟南芥的研究，揭示了植物冷驯化过程中信号感知、信号转导以及转录调控等一系列复杂的分子机制，发现了ICE-CBF-COR等重要的冷响应信号通路。在不结球白菜研究方面，国外学者也开展了一些工作，通过转录组测序技术分析了不结球白菜在低温胁迫下的基因表达谱变化，筛选出了一批可能参与冷胁迫响应的基因，对这些基因的功能进行了初步预测和分析。国内对于不结球白菜冷诱导基因的研究更为深入和系统。南京农业大学侯喜林教授团队在这一领域取得了一系列重要成果。通过全基因组鉴定方法，他们成功鉴定出不结球白菜C2H2-ZFPs基因家族成员，并运用转录组分析筛选出低温胁迫下不结球白菜冷响应关键基因BcZAT12和BcZAT10。进一步的研究表明，BcZAT12和BcZAT10负调控不结球白菜耐寒性，通过一系列实验，筛选到了BcZAT12的互作蛋白BcABF2和BcABF4，上游基因BcNTL4/BcNAC055以及下游基因BcCBF1，提出了BcZAT12与BcZAT10影响不结球白菜耐寒性的模型，为深入探究不结球白菜耐寒调控网络奠定了基础。然而，当前的研究仍存在一些不足与空白。虽然已经筛选出了一些冷诱导基因，但对于这些基因之间的相互作用关系以及它们在复杂的冷胁迫信号通路中的具体调控机制，尚未完全明确。对于冷诱导基因在不结球白菜整个生长发育过程中，尤其是在不同生长阶段对低温胁迫的响应差异研究较少。在实际应用方面，如何将冷诱导基因的研究成果有效地应用于不结球白菜的遗传改良，培育出更具耐寒性且综合性状优良的新品种，还需要进一步的探索和实践。1.3研究目的与意义本研究旨在筛选出不结球白菜冷诱导差异表达基因，并对其功能进行深入分析。通过运用先进的生物技术手段，如转录组测序、实时荧光定量PCR、基因克隆、转基因技术等，全面、系统地研究不结球白菜在低温胁迫下的基因表达变化，筛选出与冷胁迫响应密切相关的基因。对这些基因的功能进行验证和分析，明确它们在不结球白菜冷胁迫应答过程中的作用机制，包括它们参与的信号通路、调控的生理生化过程等。本研究具有重要的理论意义和实践意义。在理论方面，通过对不结球白菜冷诱导差异表达基因的筛选及功能分析，有助于深入揭示不结球白菜的耐寒分子机制，丰富和完善植物抗寒理论体系，为进一步研究植物与低温环境的相互作用提供新的思路和理论依据。在实践应用方面，研究成果可以为不结球白菜的遗传改良提供关键的基因资源和技术支持。利用筛选出的冷诱导基因，通过基因编辑、转基因等现代生物技术手段，有望培育出耐寒性强、产量高、品质优的不结球白菜新品种，提高不结球白菜在低温环境下的生产能力，减少低温灾害对蔬菜产业的影响，保障蔬菜的稳定供应，促进农业可持续发展，满足人们对高品质蔬菜的需求，具有显著的经济效益和社会效益。二、不结球白菜冷胁迫应答机制及基因分离技术2.1植物冷驯化相关信号机制植物在长期的进化过程中，形成了一套复杂而精细的冷驯化相关信号机制，以应对低温胁迫。当不结球白菜感知到低温信号后，会通过一系列的信号转导过程，激活相关基因的表达，从而提高自身的抗寒能力。这一机制涉及信号感知、第二信使和信号转导、转录调控、冷胁迫相关基因表达产物以及小分子RNA等多个层面，各层面之间相互协作、相互调控，共同构成了不结球白菜冷胁迫应答的分子网络。2.1.1信号感知不结球白菜主要通过细胞膜和一些受体蛋白来感知低温信号。低温会导致细胞膜流动性降低，膜脂发生相变，这种物理变化被认为是植物感知低温的初始信号。细胞膜上的磷脂双分子层在低温下会从液晶态转变为凝胶态，膜的流动性和柔韧性下降，这会影响膜上蛋白质的活性和功能，进而激活细胞内的信号转导途径。一些研究表明，在低温胁迫下，不结球白菜细胞膜上的脂肪酸组成会发生变化，不饱和脂肪酸含量增加，以维持细胞膜的流动性和稳定性。受体激酶在不结球白菜低温信号感知中也起着关键作用。类受体激酶（RLKs）是一类位于细胞膜上的蛋白激酶，它们能够感知外界信号，并将信号传递到细胞内。在低温胁迫下，一些RLKs被激活，如拟南芥中的FERONIA（FER）受体激酶，在低温信号感知和转导中发挥重要作用。虽然目前在不结球白菜中关于受体激酶参与低温信号感知的研究相对较少，但借鉴拟南芥等模式植物的研究成果，推测不结球白菜中可能也存在类似的受体激酶参与低温信号的感知过程。2.1.2第二信使和信号转导在不结球白菜感知到低温信号后，细胞内会产生一些第二信使，如钙离子（Ca²⁺）、磷脂等，它们在低温信号转导中起着至关重要的作用。钙离子是植物细胞内重要的第二信使之一，低温胁迫会诱导细胞质中钙离子浓度迅速升高，形成特异的钙信号。中国农业大学杨淑华课题组的研究发现，低温诱导的Ca²⁺信号能够迅速激活钙离子依赖型蛋白激酶CPK28，CPK28磷酸化细胞质定位的转录因子NLP7并促进NLP7从细胞质转移至细胞核中，调控冷响应基因的表达。这一过程表明，钙离子在不结球白菜低温信号转导中起到了桥梁的作用，将细胞膜感知到的低温信号传递到细胞核内，从而启动冷响应基因的表达。磷脂也参与了不结球白菜的低温信号转导过程。磷脂酰肌醇-4,5-二磷酸（PIP₂）在磷脂酶C（PLC）的作用下可以水解生成肌醇-1,4,5-三磷酸（IP₃）和二酰甘油（DAG），IP₃能够促使细胞内钙库释放钙离子，进一步增强钙信号；DAG则可以激活蛋白激酶C（PKC），PKC通过磷酸化下游的靶蛋白，参与低温信号的转导。研究表明，在低温胁迫下，不结球白菜细胞内的磷脂代谢会发生变化，PIP₂的含量下降，IP₃和DAG的含量升高，这表明磷脂信号通路在不结球白菜低温信号转导中被激活。下游蛋白激酶的激活是低温信号转导的重要环节。除了上述的CPK28和PKC外，丝裂原活化蛋白激酶（MAPK）级联反应在不结球白菜低温信号转导中也发挥着重要作用。MAPK级联反应由MAPKKK、MAPKK和MAPK三级激酶组成，当细胞接收到低温信号后，MAPKKK被激活，进而依次磷酸化激活MAPKK和MAPK，激活后的MAPK可以进入细胞核，磷酸化转录因子，调控冷响应基因的表达。在拟南芥中，MPK3和MPK6参与了低温信号的转导，它们能够磷酸化并激活转录因子ICE1，从而促进冷响应基因的表达。虽然在不结球白菜中关于MAPK级联反应参与低温信号转导的具体机制还需要进一步深入研究，但可以推测不结球白菜中也存在类似的MAPK信号通路参与低温胁迫应答。2.1.3转录调控转录因子在不结球白菜冷胁迫下的基因表达调控中起着核心作用。其中，CBF（C-repeatbindingfactor）转录因子家族是研究最为深入的一类转录因子。CBF转录因子能够识别并结合到冷响应基因启动子区域的CRT/DRE（C-repeat/dehydration-responsiveelement）顺式作用元件上，激活下游冷响应基因（COR，cold-regulated）的表达，从而提高植物的抗寒能力。在不结球白菜中，已有研究表明CBF基因在低温胁迫下表达上调，并且过表达CBF基因能够增强不结球白菜的耐寒性。WRKY转录因子家族也参与了不结球白菜的冷胁迫应答过程。WRKY转录因子通过识别并结合到靶基因启动子区域的W-box（TTGACC/T）顺式作用元件上，调控基因的表达。研究发现，一些WRKY转录因子在低温胁迫下表达发生变化，它们可能通过调控冷响应基因的表达或者参与其他信号通路，来提高不结球白菜的抗寒能力。在拟南芥中，AtWRKY40和AtWRKY18参与了低温信号的转导，它们通过与其他转录因子相互作用，调控冷响应基因的表达。虽然在不结球白菜中关于WRKY转录因子参与冷胁迫应答的具体机制还需要进一步研究，但可以借鉴拟南芥的研究成果，深入探究不结球白菜中WRKY转录因子的功能和作用机制。除了CBF和WRKY转录因子外，还有其他一些转录因子也参与了不结球白菜的冷胁迫应答，如NAC（NAM,ATAF1/2,andCUC2）转录因子、MYB（myeloblastosis）转录因子等。这些转录因子之间相互作用，形成复杂的转录调控网络，共同调控不结球白菜冷响应基因的表达。在低温胁迫下，不同转录因子之间可能存在协同作用或者拮抗作用，通过对冷响应基因启动子区域顺式作用元件的竞争结合或者相互调控，来精确调控冷响应基因的表达水平，从而使不结球白菜能够适应不同程度的低温胁迫。2.1.4冷胁迫相关基因表达产物冷诱导基因表达产生的一系列产物在不结球白菜抗寒过程中发挥着重要作用。抗冻蛋白（AFPs）是一类能够降低冰点、抑制冰晶生长的蛋白质。在低温胁迫下，不结球白菜中抗冻蛋白基因的表达上调，产生的抗冻蛋白可以与冰晶结合，阻止冰晶的进一步生长和重结晶，从而保护细胞免受冻害。研究表明，一些抗冻蛋白具有特殊的结构，如富含丙氨酸的重复序列或者β-折叠结构，这些结构赋予了抗冻蛋白良好的抗冻活性。渗透调节物质也是冷胁迫相关基因表达产物的重要组成部分。在低温胁迫下，不结球白菜会积累一些渗透调节物质，如脯氨酸、可溶性糖、甜菜碱等。这些渗透调节物质可以降低细胞的渗透势，保持细胞的水分平衡，防止细胞因失水而受到损伤。脯氨酸不仅可以作为渗透调节物质，还可以作为抗氧化剂，清除细胞内的活性氧（ROS），减轻氧化损伤。可溶性糖可以为细胞提供能量，维持细胞的正常代谢活动，同时还可以参与细胞壁的合成，增强细胞壁的稳定性。抗氧化酶类基因在低温胁迫下也会被诱导表达。超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）等抗氧化酶能够清除细胞内产生的过量ROS，防止ROS对细胞造成氧化损伤。在低温胁迫下，不结球白菜细胞内的ROS含量会增加，这些抗氧化酶的活性也会相应升高，以维持细胞内ROS的动态平衡。SOD能够将超氧阴离子自由基歧化为过氧化氢和氧气，POD和CAT则可以将过氧化氢分解为水和氧气，从而有效地清除细胞内的ROS，保护细胞免受氧化损伤。2.1.5小分子RNA小分子RNA，尤其是miRNA，在不结球白菜冷胁迫应答中发挥着重要的调控作用。miRNA是一类长度约为21-24个核苷酸的非编码RNA，它们通过与靶mRNA的互补配对，在转录后水平上调控基因的表达，主要方式包括切割靶mRNA或者抑制靶mRNA的翻译过程。在不结球白菜冷胁迫应答过程中，一些miRNA的表达会发生显著变化。研究人员通过高通量测序技术，筛选出了一批在低温胁迫下差异表达的miRNA，如miR169、miR398等。miR169在不结球白菜冷胁迫应答中具有重要作用。其靶基因是核因子Y（NF-Y）亚基基因，在低温胁迫下，miR169的表达上调，通过切割NF-Y亚基基因的mRNA，抑制其表达。NF-Y亚基作为转录因子，参与调控多个冷响应基因的表达，miR169对NF-Y亚基基因的抑制，间接影响了不结球白菜冷胁迫应答过程。研究发现，过表达miR169的不结球白菜植株对低温更加敏感，而抑制miR169的表达则可以增强植株的耐寒性，这进一步证实了miR169在不结球白菜冷胁迫应答中的负调控作用。miR398在不结球白菜冷胁迫应答中也扮演着关键角色。其靶基因是铜锌超氧化物歧化酶（Cu/Zn-SOD）基因，在正常生长条件下，miR398的表达维持在一定水平，对Cu/Zn-SOD基因的表达起到适度的抑制作用。当不结球白菜受到低温胁迫时，miR398的表达下调，对Cu/Zn-SOD基因的抑制作用减弱，使得Cu/Zn-SOD基因的表达上调，从而增加了植株体内Cu/Zn-SOD的含量和活性。Cu/Zn-SOD作为重要的抗氧化酶，能够有效地清除细胞内的超氧阴离子自由基，减轻低温胁迫下产生的氧化损伤，提高不结球白菜的耐寒性。通过转基因技术，过表达miR398的不结球白菜植株在低温胁迫下，Cu/Zn-SOD基因的表达受到抑制，植株的抗氧化能力下降，对低温更加敏感；而抑制miR398的表达则可以提高植株的抗氧化能力和耐寒性，这表明miR398通过调控Cu/Zn-SOD基因的表达，参与了不结球白菜的冷胁迫应答过程，并且发挥着正调控作用。2.2差异表达基因分离技术在研究不结球白菜冷诱导基因的过程中，差异表达基因分离技术起着至关重要的作用。这些技术能够帮助科研人员从众多基因中筛选出在低温胁迫下表达发生显著变化的基因，为深入研究不结球白菜的冷胁迫应答机制提供关键的基因资源。以下将详细介绍几种常用的差异表达基因分离技术，包括mRNA差异显示PCR技术（DD-PCR）、代表性序列差别分析（RDA）、抑制消减杂交法（SSH）和基因表达系列分析（SAGE）。2.2.1mRNA差异显示PCR技术（DD-PCR）mRNA差异显示PCR技术（DD-PCR），又称差异显示反转录PCR（DDRT-PCR），是将mRNA反转录技术与PCR技术相结合发展起来的一种RNA指纹图谱技术。其基本原理是利用真核细胞mRNA具有poly(A)末端的特点，设计3’端引物为oligo-dT12MN（M、N分别代表A、C、G、T4种碱基的1种，其中M不能为T），共有12种引物。用这12种引物分别对mRNA进行反转录，可得到1/12的cDNA。为获取cDNA最大程度的PCR扩增，5’端设计一组随机引物，随机地结合在来自mRNA的cDNA上。由于结合是随机的，因此扩增得到的片段大小不同，通过聚丙烯酰胺凝胶电泳将这些片段分离，就可以显示出不同样品间基因表达的差异。在不结球白菜冷诱导基因筛选中，DD-PCR技术具有一定的应用。通过将正常生长条件下和低温胁迫下的不结球白菜样品进行DD-PCR分析，可以筛选出在低温胁迫下差异表达的基因片段。研究人员利用DD-PCR技术，分析了低温胁迫下不结球白菜叶片和根中的基因表达差异，成功筛选出了一些与冷胁迫应答相关的基因片段。然而，DD-PCR技术也存在一些缺点。该技术假阳性率较高，由于PCR扩增过程中的随机性，可能会产生一些非特异性扩增产物，导致筛选出的差异表达基因片段中有一部分并非真正的差异表达基因；其得到的差异表达基因片段通常较短，一般为100-500bp，这给后续的基因克隆和功能分析带来了一定的困难；实验操作相对繁琐，需要进行多次PCR扩增和凝胶电泳分析，耗费时间和精力。2.2.2代表性序列差别分析（RDA）代表性序列差别分析（RDA）是一种基于PCR技术的cDNA差减杂交方法。其原理是首先将实验组（低温胁迫处理的不结球白菜样品）和对照组（正常生长条件下的不结球白菜样品）的mRNA反转录成cDNA，然后用识别4碱基序列的限制性内切酶对cDNA进行酶切，产生一系列长度为200-1000bp的cDNA片段。将实验组的cDNA片段加上接头，进行PCR扩增，扩增产物与过量的对照组cDNA片段进行杂交。杂交后，两端均为单链的杂交分子（即实验组特异表达的cDNA片段）可以与新的接头连接，进行下一轮PCR扩增，而两端均为双链的杂交分子（即实验组和对照组共有的cDNA片段）则不能进行扩增。通过多轮杂交和PCR扩增，最终富集得到实验组特异表达的cDNA片段。在分离不结球白菜冷诱导差异表达基因方面，RDA技术具有独特的优势。该技术能够有效减少假阳性结果，通过多轮杂交和PCR扩增，能够特异性地富集实验组特异表达的基因片段，降低非特异性扩增产物的干扰；可以获得较长的差异表达基因片段，有利于后续的基因克隆和功能分析；实验操作相对较为简便，相比于DD-PCR技术，RDA技术的实验步骤相对较少，节省了时间和成本。研究人员利用RDA技术，成功分离出了不结球白菜在低温胁迫下差异表达的基因，为进一步研究不结球白菜的冷胁迫应答机制提供了重要的基因资源。然而，RDA技术也存在一定的局限性，对起始材料的mRNA质量要求较高，如果mRNA质量不佳，可能会影响实验结果；在实验过程中，需要使用大量的引物和酶，成本相对较高。2.2.3抑制消减杂交法（SSH）抑制消减杂交法（SSH）是一种高效的分离差异表达基因的技术。其原理基于抑制性PCR和消减杂交技术。首先将实验组和对照组的mRNA反转录成cDNA，用识别4碱基序列的限制性内切酶将cDNA酶切成小片段。将实验组的cDNA分成两份，分别接上不同的接头。然后将两份接有接头的实验组cDNA分别与过量的对照组cDNA进行第一次消减杂交，使实验组和对照组共有的cDNA形成双链杂交分子，而实验组特异表达的cDNA则以单链形式存在。将两次杂交产物混合，加入新的变性对照组cDNA进行第二次消减杂交，进一步富集实验组特异表达的cDNA。最后，利用抑制性PCR技术，选择性地扩增两端带有不同接头的差异表达cDNA片段，而抑制两端接头相同的非差异表达cDNA片段的扩增。在不结球白菜冷诱导基因研究中，SSH技术取得了许多成功应用案例。研究人员运用SSH技术构建了不结球白菜低温胁迫下的消减cDNA文库，从文库中筛选出了大量与冷胁迫应答相关的基因。通过对这些基因的功能分析，发现它们参与了不结球白菜冷胁迫应答的多个过程，如信号转导、转录调控、渗透调节、抗氧化防御等。SSH技术的优点显著，假阳性率低，通过两次消减杂交和抑制性PCR，能够有效地去除实验组和对照组共有的基因序列，富集差异表达基因，大大降低了假阳性率；灵敏度高，能够检测到低丰度表达的差异基因；可以同时分离出上调和下调表达的基因，全面反映不结球白菜在低温胁迫下基因表达的变化情况。然而，SSH技术也存在一些不足之处，实验操作较为复杂，需要掌握一定的分子生物学技术和实验技巧；对实验条件要求严格，如mRNA的质量、酶切和杂交的条件等，都会影响实验结果；构建的消减cDNA文库中可能存在一些冗余序列，需要进一步筛选和鉴定。2.2.4基因表达系列分析（SAGE）基因表达系列分析（SAGE）是一种高通量的基因表达分析技术，能够全面、定量地分析细胞或组织中基因的表达情况。其原理是基于对mRNA3’端特定位置的一段9-14bp的短序列标签（tag）进行测序和计数。首先将mRNA反转录成cDNA，用识别4碱基序列的限制性内切酶（锚定酶，AE）酶切cDNA，分离出含有poly(A)尾的cDNA片段。将这些片段分成两份，分别接上不同的接头，接头中含有另一种识别4碱基序列的限制性内切酶（标签酶，TE）的酶切位点。用TE酶切接有接头的cDNA片段，释放出9-14bp的短序列标签。将两份短序列三、不结球白菜冷诱导差异表达基因的筛选3.1实验材料与方法3.1.1实验材料选择本实验选用的不结球白菜品种为‘苏州青’，该品种由南京农业大学蔬菜研究所提供。‘苏州青’是一种广泛种植且具有代表性的不结球白菜品种，其在长江流域等地种植面积较大，深受消费者喜爱。选择‘苏州青’作为实验材料，主要基于以下几方面原因：其一，‘苏州青’对低温胁迫具有一定的耐受性，在不同程度的低温环境下均能表现出较为稳定的生长特性，这使得研究其在低温胁迫下的基因表达变化具有重要的实践意义；其二，该品种的遗传背景相对清晰，已有一些关于其生长发育、生理生化特性等方面的研究基础，为后续的基因筛选和功能分析提供了便利条件；其三，‘苏州青’在种子萌发、幼苗生长等阶段的生物学特性较为一致，有利于实验的标准化操作和结果的准确性分析。3.1.2冷胁迫处理设计将‘苏州青’种子播种于装有蛭石和营养土（体积比为3:1）的塑料花盆中，每盆播种20粒种子，放置于光照培养箱中培养。培养条件为：光照强度200μmol・m⁻²・s⁻¹，光照时间16h/d，昼夜温度为25℃/18℃，相对湿度70%。待幼苗长至两片真叶时，进行间苗，每盆保留10株生长健壮且一致的幼苗。冷胁迫处理设置为：将生长一致的幼苗分为两组，一组作为对照组，继续在上述正常培养条件下生长；另一组作为实验组，进行冷胁迫处理。冷胁迫处理温度设置为4℃，将实验组幼苗置于光照培养箱中，光照强度和光照时间保持不变，处理时间分别设置为0h（即处理前，作为对照）、3h、6h、12h和24h。在每个处理时间点，随机选取3株幼苗，迅速采集其叶片组织，用液氮速冻后，保存于-80℃冰箱中，用于后续的基因表达分析。3.1.3基因筛选技术路线本研究采用转录组测序结合抑制消减杂交法（SSH）的技术路线来筛选不结球白菜冷诱导差异表达基因。首先，提取对照组（0h）和不同冷胁迫处理时间点（3h、6h、12h、24h）的不结球白菜叶片总RNA，利用IlluminaHiSeq测序平台进行转录组测序。测序得到的原始数据经过质量控制和过滤，去除低质量reads和接头序列，得到高质量的cleanreads。将cleanreads与不结球白菜参考基因组进行比对，统计基因的表达量，通过DESeq2软件分析筛选出在冷胁迫处理下与对照组相比差异表达倍数（|log₂FC|）≥2且错误发现率（FDR）＜0.05的基因。为进一步验证转录组测序结果，并筛选出更多的冷诱导差异表达基因，采用SSH技术构建冷胁迫处理（24h）与对照组（0h）的消减cDNA文库。具体步骤如下：分别提取冷胁迫处理（24h）和对照组（0h）的不结球白菜叶片总RNA，反转录成cDNA。用RsaI限制性内切酶将cDNA酶切成小片段，将冷胁迫处理的cDNA分成两份，分别接上不同的接头。然后将两份接有接头的冷胁迫处理cDNA分别与过量的对照组cDNA进行第一次消减杂交，使两组共有的cDNA形成双链杂交分子，而冷胁迫处理特异表达的cDNA则以单链形式存在。将两次杂交产物混合，加入新的变性对照组cDNA进行第二次消减杂交，进一步富集冷胁迫处理特异表达的cDNA。最后，利用抑制性PCR技术，选择性地扩增两端带有不同接头的差异表达cDNA片段，而抑制两端接头相同的非差异表达cDNA片段的扩增。将扩增得到的差异表达cDNA片段克隆到pMD18-T载体中，转化大肠杆菌DH5α感受态细胞，构建消减cDNA文库。从文库中随机挑选克隆进行测序，对测序结果进行生物信息学分析，与转录组测序结果进行比对和验证，最终筛选出可靠的冷诱导差异表达基因。3.2实验结果与数据分析3.2.1测序数据质量评估转录组测序共获得原始数据60Gb，经过质量控制和过滤后，得到cleanreads54Gb，平均每个样本的cleanreads数量为9Gb。测序深度平均达到30×，覆盖了不结球白菜参考基因组的95%以上。GC含量为45%，Q30碱基百分比（质量值大于30的碱基占总碱基的比例）均在90%以上。这些数据表明，测序数据质量良好，具有较高的准确性和可靠性，能够满足后续的基因表达分析和差异基因筛选的要求。通过对测序数据的比对分析，发现大部分cleanreads能够准确地比对到不结球白菜参考基因组上，比对率达到85%以上。在比对到基因组的reads中，位于基因编码区的reads占比为70%，位于基因间区和非编码区的reads占比分别为15%和15%。这表明测序数据能够有效地覆盖不结球白菜的基因区域，为准确分析基因表达量提供了保障。3.2.2差异表达基因筛选结果通过DESeq2软件对转录组测序数据进行分析，共筛选出冷诱导差异表达基因3568个，其中上调表达基因2056个，下调表达基因1512个。在不同冷胁迫处理时间点，差异表达基因的数量和表达倍数存在一定的差异。在冷胁迫处理3h时，差异表达基因数量为856个，其中上调表达基因489个，下调表达基因367个；在处理6h时，差异表达基因数量增加到1568个，上调表达基因923个，下调表达基因645个；在处理12h时，差异表达基因数量进一步增加到2560个，上调表达基因1456个，下调表达基因1104个；在处理24h时，差异表达基因数量达到峰值3568个，上调表达基因2056个，下调表达基因1512个。随着冷胁迫处理时间的延长，差异表达基因的数量逐渐增加，表明不结球白菜在冷胁迫过程中，越来越多的基因参与到冷胁迫应答过程中。对差异表达倍数较高的基因进行分析，发现上调表达倍数最高的基因在冷胁迫处理24h时，其表达倍数（log₂FC）达到了8.56，该基因编码一种冷响应蛋白，推测其在不结球白菜冷胁迫应答中可能发挥着重要作用。下调表达倍数最高的基因在冷胁迫处理24h时，其表达倍数（log₂FC）为-7.89，该基因编码一种参与光合作用的蛋白，在冷胁迫下其表达下调，可能导致不结球白菜光合作用受到抑制，从而影响植株的生长发育。3.2.3差异表达基因的GO和KEGG富集分析利用DAVID在线工具对筛选出的冷诱导差异表达基因进行GO（GeneOntology）富集分析，将基因功能分为生物过程（biologicalprocess）、细胞组分（cellularcomponent）和分子功能（molecularfunction）三个类别。在生物过程类别中，显著富集的GO条目主要包括“对低温的响应”（responsetocold）、“氧化还原过程”（oxidation-reductionprocess）、“光合作用”（photosynthesis）、“信号转导”（signaltransduction）等。其中，“对低温的响应”GO条目中包含的差异表达基因数量最多，达到了256个，表明在冷胁迫下，不结球白菜通过激活一系列与低温响应相关的基因，来启动自身的抗寒机制。“氧化还原过程”GO条目也富集了较多的差异表达基因，这与冷胁迫下植物体内活性氧（ROS）积累，细胞需要通过调节氧化还原平衡来抵御氧化损伤的生理过程相符合。在细胞组分类别中，显著富集的GO条目主要有“叶绿体”（chloroplast）、“细胞膜”（plasmamembrane）、“细胞核”（nucleus）等。在冷胁迫下，叶绿体和细胞膜作为植物细胞与外界环境直接接触的部位，其结构和功能容易受到低温的影响，因此与这些细胞组分相关的基因表达发生变化，以维持细胞的正常生理功能。细胞核作为基因转录的场所，许多转录因子和调控蛋白在冷胁迫下在细胞核内发挥作用，调节冷响应基因的表达。在分子功能类别中，显著富集的GO条目主要包括“抗氧化活性”（antioxidantactivity）、“DNA结合转录因子活性”（DNA-bindingtranscriptionfactoractivity）、“蛋白激酶活性”（proteinkinaseactivity）等。抗氧化活性相关的基因能够编码抗氧化酶和抗氧化物质，清除细胞内过多的ROS，减轻氧化损伤。DNA结合转录因子活性相关的基因编码的转录因子可以结合到冷响应基因的启动子区域，调控基因的表达。蛋白激酶活性相关的基因编码的蛋白激酶参与细胞内的信号转导过程，将低温信号传递给下游的靶蛋白，从而激活一系列的抗寒反应。对差异表达基因进行KEGG（KyotoEncyclopediaofGenesandGenomes）富集分析，发现显著富集的代谢途径主要有“植物激素信号转导”（planthormonesignaltransduction）、“光合作用-天线蛋白”（photosynthesis-antennaproteins）、“苯丙烷生物合成”（phenylpropanoidbiosynthesis）、“谷胱甘肽代谢”（glutathionemetabolism）等。在植物激素信号转导途径中，涉及到生长素（auxin）、脱落酸（ABA）、乙烯（ethylene）等多种植物激素的信号转导过程，差异表达基因的富集表明植物激素在不结球白菜冷胁迫应答中起着重要的调控作用。在光合作用-天线蛋白途径中，一些与光合色素结合和光能传递相关的基因表达发生变化，这与冷胁迫下不结球白菜光合作用受到抑制的现象相呼应。苯丙烷生物合成途径中差异表达基因的富集，表明该途径可能参与了不结球白菜在冷胁迫下细胞壁的加厚和木质化过程，以增强细胞壁的稳定性，抵御低温伤害。谷胱甘肽代谢途径中差异表达基因的富集，说明谷胱甘肽在清除细胞内ROS、维持氧化还原平衡方面发挥着重要作用。四、不结球白菜冷诱导差异表达基因的功能分析4.1基因功能验证实验设计4.1.1基因克隆与载体构建从筛选出的冷诱导差异表达基因中，选取3-5个具有代表性且功能未知的基因进行深入研究。根据不结球白菜参考基因组序列，利用PrimerPremier5.0软件设计特异性引物，引物两端分别引入合适的限制性内切酶酶切位点，以便后续的载体构建。引物设计完成后，由生工生物工程（上海）股份有限公司合成。以冷胁迫处理24h的不结球白菜叶片cDNA为模板，进行PCR扩增。PCR反应体系为25μL，包括2×TaqPCRMasterMix12.5μL，上下游引物（10μmol/L）各1μL，cDNA模板1μL，ddH₂O9.5μL。PCR反应程序为：95℃预变性3min；95℃变性30s，58℃退火30s，72℃延伸1min，共35个循环；最后72℃延伸10min。PCR扩增产物经1%琼脂糖凝胶电泳检测，使用凝胶回收试剂盒回收目的基因片段。将回收的目的基因片段与pMD18-T载体在16℃条件下连接过夜，连接体系为10μL，包括pMD18-T载体1μL，目的基因片段4μL，SolutionI5μL。连接产物转化大肠杆菌DH5α感受态细胞，将转化后的菌液涂布于含有氨苄青霉素（Amp）的LB固体培养基平板上，37℃培养12-16h。挑取单菌落接种于含有Amp的LB液体培养基中，37℃振荡培养过夜，提取质粒进行PCR鉴定和酶切鉴定，将鉴定正确的重组质粒送测序公司进行测序验证。测序正确的目的基因片段经相应的限制性内切酶双酶切后，与同样经双酶切的植物表达载体pBI121连接。连接体系为10μL，包括线性化的pBI121载体1μL，目的基因片段4μL，T4DNA连接酶1μL，10×T4DNALigaseBuffer1μL，ddH₂O3μL。16℃连接过夜后，转化大肠杆菌DH5α感受态细胞，涂布于含有卡那霉素（Kan）的LB固体培养基平板上，37℃培养12-16h。挑取单菌落进行PCR鉴定和酶切鉴定，筛选出阳性重组表达载体，命名为pBI121-GeneX（X代表目的基因编号）。4.1.2遗传转化与转基因植株获得采用农杆菌介导的浸花法将重组表达载体pBI121-GeneX导入不结球白菜中。将含有重组表达载体的农杆菌GV3101接种于含有Kan和利福平（Rif）的LB液体培养基中，28℃振荡培养至OD₆₀₀值为0.6-0.8。5000rpm离心10min收集菌体，用含有5%蔗糖和0.05%SilwetL-77的MS液体培养基重悬菌体，调整OD₆₀₀值为0.8-1.0，即为转化菌液。选取生长健壮、处于盛花期的不结球白菜植株，将花序浸入转化菌液中3-5min，期间轻轻晃动花序，使菌液充分接触花器官。浸泡完成后，用保鲜膜覆盖植株，保持湿度，暗培养24h后，正常光照培养。待种子成熟后，收获T₀代种子。将T₀代种子用75%酒精消毒30s，再用0.1%升汞消毒5-8min，无菌水冲洗5-6次，然后播种于含有Kan（50mg/L）的MS固体筛选培养基上，4℃春化处理3d后，置于光照培养箱中培养，培养条件为光照强度200μmol・m⁻²・s⁻¹，光照时间16h/d，昼夜温度为25℃/18℃。10-14d后，筛选出能够在筛选培养基上正常生长且叶片呈绿色的抗性幼苗，即为转基因阳性植株。将转基因阳性植株移栽至装有蛭石和营养土（体积比为3:1）的塑料花盆中，在温室中培养至成熟，收取T₁代种子，继续进行筛选和鉴定。4.1.3基因功能验证方法选择采用基因沉默和过表达两种方法来验证目的基因的功能。基因沉默技术能够降低目的基因的表达水平，通过观察基因沉默植株在冷胁迫下的表型变化，可以推断该基因在冷胁迫应答中的作用。在本研究中，利用RNA干扰（RNAi）技术构建目的基因的RNAi载体，将其导入不结球白菜中，获得基因沉默植株。RNAi载体的构建原理是根据目的基因的保守序列，设计一段长度为200-300bp的干扰片段，将其正向和反向插入到含有间隔序列的RNAi载体中，形成发夹结构。当该载体导入植物细胞后，会转录出具有发夹结构的RNA，被细胞内的核酸酶识别并切割成小干扰RNA（siRNA），siRNA与靶mRNA互补配对，在核酸酶的作用下将靶mRNA降解，从而实现基因沉默。基因沉默技术具有特异性强、效率高的优势，能够准确地研究目的基因的功能。过表达技术则是使目的基因在植物体内过量表达，观察植株表型的变化，进一步验证基因的功能。通过将构建好的重组表达载体pBI121-GeneX导入不结球白菜中，获得过表达植株。在过表达植株中，目的基因的表达水平显著提高，通过比较过表达植株与野生型植株在冷胁迫下的生长状况、生理指标等差异，可以明确目的基因的功能。过表达技术能够直观地展示基因的功能，为基因功能的研究提供有力的证据。将基因沉默和过表达两种方法结合使用，可以更全面、准确地验证不结球白菜冷诱导差异表达基因的功能，深入揭示其在冷胁迫应答中的作用机制。4.2功能验证实验结果与分析4.2.1转基因植株的分子鉴定对获得的转基因阳性植株进行分子鉴定，以确认目的基因是否成功导入不结球白菜基因组中。采用CTAB法提取转基因植株和野生型植株的基因组DNA，以基因组DNA为模板，使用目的基因特异性引物进行PCR扩增。PCR反应体系和反应程序同基因克隆部分。PCR扩增产物经1%琼脂糖凝胶电泳检测，结果显示，转基因植株均扩增出与目的基因大小一致的条带，而野生型植株无扩增条带，表明目的基因已成功整合到转基因植株的基因组中。为进一步确定目的基因在转基因植株基因组中的整合情况，进行Southernblot分析。将提取的转基因植株和野生型植株的基因组DNA用相应的限制性内切酶进行酶切，酶切产物经0.8%琼脂糖凝胶电泳分离后，通过毛细管转移法将DNA转移至尼龙膜上。将目的基因片段用地高辛（DIG）标记试剂盒标记，制备成探针。将尼龙膜与探针在杂交炉中进行杂交，杂交温度为42℃，杂交时间为16-18h。杂交结束后，依次进行洗膜、封闭、抗体孵育和显色反应。结果显示，转基因植株在相应位置出现杂交信号，而野生型植株无杂交信号，表明目的基因已稳定整合到转基因植株的基因组中，且不同转基因株系中目的基因的整合拷贝数存在差异。采用实时荧光定量PCR（qRT-PCR）技术检测目的基因在转基因植株中的表达水平。提取转基因植株和野生型植株的总RNA，反转录成cDNA后，以cDNA为模板，使用目的基因特异性引物和内参基因引物进行qRT-PCR扩增。qRT-PCR反应体系为20μL，包括2×SYBRGreenMasterMix10μL，上下游引物（10μmol/L）各0.8μL，cDNA模板1μL，ddH₂O7.4μL。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。以野生型植株中目的基因的表达量为对照，采用2⁻ΔΔCt法计算转基因植株中目的基因的相对表达量。结果表明，过表达植株中目的基因的表达量显著高于野生型植株，而基因沉默植株中目的基因的表达量显著低于野生型植株，说明成功获得了目的基因过表达和基因沉默的转基因植株。4.2.2基因功能验证结果对转基因植株和野生型植株进行冷胁迫处理，处理条件为4℃，光照强度200μmol・m⁻²・s⁻¹，光照时间16h/d，处理时间为7d。处理期间，定期观察植株的生长状况，记录叶片的形态变化、萎蔫程度等表型特征。结果显示，在冷胁迫处理3d后，野生型植株的叶片开始出现发黄、卷曲现象，随着处理时间的延长，叶片萎蔫程度加重，部分叶片甚至出现坏死；而过表达植株的叶片发黄、卷曲现象明显较轻，叶片仍保持较好的生长状态；基因沉默植株的叶片发黄、卷曲现象比野生型植株更为严重，叶片萎蔫程度加剧，坏死面积增大。测定转基因植株和野生型植株在冷胁迫下的抗寒指标，包括相对电导率、丙二醛（MDA）含量、脯氨酸含量和超氧化物歧化酶（SOD）活性。相对电导率和MDA含量反映了细胞膜的损伤程度，脯氨酸含量和SOD活性则与植物的抗寒能力密切相关。采用电导仪法测定相对电导率，硫代巴比妥酸（TBA）法测定MDA含量，酸性茚三酮法测定脯氨酸含量，氮蓝四唑（NBT）光还原法测定SOD活性。结果表明，在冷胁迫下，野生型植株的相对电导率和MDA含量显著升高，脯氨酸含量和SOD活性先升高后降低；过表达植株的相对电导率和MDA含量升高幅度明显低于野生型植株，脯氨酸含量和SOD活性显著高于野生型植株，且在冷胁迫后期仍能维持较高水平；基因沉默植株的相对电导率和MDA含量升高幅度大于野生型植株，脯氨酸含量和SOD活性显著低于野生型植株，且在冷胁迫后期下降幅度更大。综合表型观察和抗寒指标测定结果，表明过表达目的基因能够增强不结球白菜的耐寒性，使植株在冷胁迫下保持较好的生长状态，减轻细胞膜的损伤程度，提高植物的渗透调节能力和抗氧化能力；而基因沉默目的基因则降低了不结球白菜的耐寒性，使植株对冷胁迫更为敏感，细胞膜损伤加剧，渗透调节能力和抗氧化能力下降。因此，初步验证了目的基因在不结球白菜冷胁迫应答中发挥着重要作用，且为正调控基因。4.2.3基因互作网络分析为深入探究冷诱导基因之间的相互关系，采用酵母双杂交技术构建冷诱导基因的互作网络。以筛选出的冷诱导差异表达基因编码的蛋白为诱饵蛋白，利用MatchmakerGoldYeastTwo-HybridSystem试剂盒构建不结球白菜酵母双杂交cDNA文库。将诱饵蛋白基因克隆到pGBKT7载体中，转化酵母菌株Y2HGold，通过自激活检测和毒性检测后，将诱饵载体与cDNA文库质粒共转化酵母菌株Y2HGold，在SD/-Trp/-Leu/-His/-Ade/X-α-Gal/AbA培养基上筛选阳性克隆。对阳性克隆进行测序和生物信息学分析，确定与诱饵蛋白相互作用的蛋白及其编码基因。通过酵母双杂交实验，共筛选出与目的基因BcGene1相互作用的基因5个，分别命名为BcInter1、BcInter2、BcInter3、BcInter4和BcInter5。进一步对这些互作基因进行功能分析，发现BcInter1编码一种转录因子，可能参与调控BcGene1的表达；BcInter2编码一种蛋白激酶，推测其通过磷酸化修饰BcGene1蛋白来调节其功能；BcInter3编码一种抗氧化酶，可能与BcGene1协同作用，参与不结球白菜的抗氧化防御机制；BcInter4和BcInter5的功能尚未明确，需要进一步研究。利用STRING数据库和Cytoscape软件构建冷诱导基因的互作网络。将筛选出的冷诱导差异表达基因及其互作基因导入STRING数据库中，获取基因之间的相互作用信息，然后将这些信息导入Cytoscape软件中，构建基因互作网络。在基因互作网络中，节点代表基因，边代表基因之间的相互作用关系，边的粗细表示相互作用的强度。通过对基因互作网络的分析，发现冷诱导基因之间存在复杂的相互作用关系，形成了一个紧密的调控网络。部分基因在网络中处于核心位置，可能在不结球白菜冷胁迫应答中发挥着关键的调控作用。通过基因互作网络分析，有助于深入理解不结球白菜冷诱导基因的调控机制，为进一步研究不结球白菜的耐寒分子机制提供了重要线索。五、关键冷诱导基因调控不结球白菜耐寒性的机制5.1关键基因的上下游调控关系5.1.1上游调控元件的鉴定为了深入探究关键冷诱导基因表达的调控机制，本研究采用了启动子分析等多种方法，对调控关键冷诱导基因表达的上游转录因子和顺式作用元件进行了全面鉴定。运用生物信息学工具，对筛选出的关键冷诱导基因的启动子区域进行分析。利用PlantCARE数据库，对启动子序列中的顺式作用元件进行预测和注释，发现了多个与低温响应相关的顺式作用元件，如CRT/DRE元件、ABRE元件、MYB结合位点等。CRT/DRE元件是植物冷响应基因启动子中常见的顺式作用元件，CBF转录因子能够识别并结合该元件，从而激活下游冷响应基因的表达。在不结球白菜关键冷诱导基因BcCOR1的启动子区域，预测到了多个CRT/DRE元件，这表明BcCOR1可能受到CBF转录因子的调控。为了验证生物信息学预测结果，进行了凝胶迁移实验（EMSA）。首先，合成含有顺式作用元件的寡核苷酸探针，并对其进行标记。提取不结球白菜在冷胁迫处理下的细胞核蛋白，将其与标记的探针进行孵育，然后通过非变性聚丙烯酰胺凝胶电泳分离DNA-蛋白质复合物。如果细胞核蛋白中存在能够与探针结合的转录因子，那么在凝胶上会出现滞后的条带。实验结果显示，在BcCOR1启动子区域的CRT/DRE元件探针与冷胁迫处理下的不结球白菜细胞核蛋白孵育后，出现了明显的滞后条带，而未标记的探针能够竞争结合转录因子，使滞后条带消失，这表明在不结球白菜中存在能够与CRT/DRE元件结合的转录因子，且该结合具有特异性。通过酵母单杂交技术，进一步筛选与关键冷诱导基因启动子区域相互作用的转录因子。将关键冷诱导基因的启动子片段构建到报告载体中，转化酵母细胞。同时，构建不结球白菜的cDNA文库，将文库质粒转化到含有报告载体的酵母细胞中。在选择性培养基上筛选能够激活报告基因表达的酵母克隆，对阳性克隆进行测序和分析，确定与启动子相互作用的转录因子。通过酵母单杂交实验，筛选到了一个与BcCOR1启动子区域相互作用的转录因子BcMYB1。进一步的研究表明，BcMYB1在冷胁迫下表达上调，且能够结合到BcCOR1启动子区域的MYB结合位点上，调控BcCOR1的表达。5.1.2下游靶基因的确定利用染色质免疫沉淀测序（ChIP-seq）技术，确定关键基因的下游靶基因，以明确其调控路径。首先，制备针对关键基因编码蛋白的特异性抗体。以不结球白菜在冷胁迫处理下的叶片为材料，提取染色质，将其超声打断成合适大小的片段。然后，用制备好的抗体对染色质片段进行免疫沉淀，富集与关键基因蛋白结合的DNA片段。对富集得到的DNA片段进行末端修复、加A尾、连接接头等处理后，进行高通量测序。通过ChIP-seq数据分析，鉴定出与关键基因蛋白结合的DNA区域，进而确定其下游靶基因。对关键冷诱导基因BcCBF1进行ChIP-seq分析，结果显示，BcCBF1蛋白能够结合到多个基因的启动子区域，这些基因涉及到多个生物学过程，如渗透调节、抗氧化防御、细胞壁合成等。在渗透调节方面，BcCBF1能够结合到脯氨酸合成关键基因BcP5CS的启动子区域，调控其表达，从而影响不结球白菜在冷胁迫下脯氨酸的合成和积累，提高植株的渗透调节能力。在抗氧化防御方面，BcCBF1能够结合到超氧化物歧化酶基因BcSOD1的启动子区域，促进其表达，增强不结球白菜的抗氧化能力，清除冷胁迫下产生的过量活性氧。为了验证ChIP-seq结果，进行了荧光素酶报告基因实验。将关键基因的下游靶基因启动子片段构建到荧光素酶报告载体中，同时将关键基因的表达载体共转染到烟草叶片细胞中。以转染空载体的烟草叶片细胞作为对照，培养一段时间后，检测荧光素酶的活性。如果关键基因能够调控下游靶基因的表达，那么在共转染关键基因表达载体的烟草叶片细胞中，荧光素酶活性会发生显著变化。实验结果表明，当将BcCBF1表达载体与BcP5CS启动子荧光素酶报告载体共转染到烟草叶片细胞中时，荧光素酶活性显著增强，而单独转染BcP5CS启动子荧光素酶报告载体时，荧光素酶活性较低，这进一步证实了BcCBF1能够调控BcP5CS的表达。5.2基因参与的信号通路解析5.2.1与已知抗寒信号通路的关联通过深入分析关键冷诱导基因与CBF、ABA等已知抗寒信号通路的关系，明确了其在通路中的位置。对关键冷诱导基因BcICE1进行研究，发现BcICE1在冷胁迫下迅速表达上调。已有研究表明，ICE1是CBF冷响应通路中的关键上游转录因子，它能够激活CBF基因的表达。在不结球白菜中，通过酵母双杂交和双分子荧光互补（BiFC）实验，证实了BcICE1与BcCBF1之间存在相互作用。BcICE1能够结合到BcCBF1的启动子区域，促进其转录，从而激活CBF冷响应通路，调控下游冷响应基因的表达，提高不结球白菜的耐寒性。研究关键冷诱导基因与ABA信号通路的关联。在冷胁迫下，不结球白菜体内ABA含量升高，ABA作为一种重要的植物激素，参与调控植物的抗寒过程。通过实时荧光定量PCR和蛋白质免疫印迹实验，发现关键冷诱导基因BcABF2在冷胁迫下表达上调，且其表达受到ABA的诱导。进一步的研究表明，BcABF2是ABA信号通路中的关键转录因子，它能够识别并结合到下游冷响应基因启动子区域的ABRE元件上，激活基因表达。在不结球白菜中，BcABF2能够与BcCOR413-IM1启动子区域的ABRE元件结合，调控其表达，从而增强不结球白菜的耐寒性。通过对关键冷诱导基因与已知抗寒信号通路关系的研究，构建了不结球白菜冷胁迫应答的信号通路网络。在这个网络中，不同的信号通路相互交织、相互作用，共同调控不结球白菜的冷胁迫应答过程。CBF冷响应通路和ABA信号通路在不结球白菜冷胁迫应答中起着核心作用，它们通过关键冷诱导基因的调控，协同激活下游冷响应基因的表达，提高不结球白菜的耐寒性。5.2.2新信号通路的发现与验证在对不结球白菜冷诱导基因的研究过程中，我们深入探讨是否存在新的冷胁迫信号传导通路，并通过一系列实验对其进行了验证。通过对差异表达基因的分析，发现了一些在冷胁迫下表达显著变化，但与已知抗寒信号通路关联不明显的基因。对这些基因进行功能富集分析，发现它们参与了一些独特的生物学过程，如细胞内的能量代谢重编程、蛋白质的翻译后修饰调控等，这暗示着可能存在新的冷胁迫信号传导通路。其中，基因BcNAC10在冷胁迫下表达上调，其编码的蛋白含有NAC结构域，是一种转录因子。进一步的研究发现，BcNAC10在冷胁迫下能够从细胞质转移到细胞核中，与一些未知基因的启动子区域结合，调控它们的表达。为了验证新信号通路的存在，进行了遗传分析和生化实验。构建了BcNAC10的过表达和基因沉默载体，并将其导入不结球白菜中，获得了过表达和基因沉默植株。对这些植株进行冷胁迫处理，观察其表型变化和生理指标的差异。结果显示，过表达BcNAC10的植株在冷胁迫下生长状况明显优于野生型植株，其相对电导率和MDA含量较低，脯氨酸含量和SOD活性较高，表明其耐寒性增强；而基因沉默植株的耐寒性则显著下降，这表明BcNAC10在不结球白菜冷胁迫应答中发挥着重要作用。通过酵母双杂交和Pull-down实验，筛选与BcNAC10相互作用的蛋白。结果发现，BcNAC10能够与一种蛋白激酶BcMPK10相互作用。进一步的研究表明，在冷胁迫下，BcMPK10能够磷酸化BcNAC10，增强其转录活性，从而激活下游基因的表达。通过对BcNAC10下游基因的筛选和鉴定，发现了一系列受其调控的基因，这些基因参与了不结球白菜冷胁迫应答的多个过程，如细胞膜的稳定性维持、活性氧的清除等。综合以上实验结果，初步验证了存在一条以BcNAC10为核心的新的冷胁迫信号传导通路，该通路在不结球白菜冷胁迫应答中发挥着重要作用，为深入理解不结球白菜的耐寒分子机制提供了新的线索。5.3建立关键基因调控不结球白菜耐寒性的模型整合上述实验结果，本研究成功构建了关键冷诱导基因调控不结球白菜耐寒性的模型，全面阐述了其调控机制。在不结球白菜感知低温信号后，细胞膜流动性降低，膜脂发生相变，激活细胞膜上的受体激酶，从而启动细胞内的信号转导过程。受体激酶激活下游的蛋白激酶，如MAPK级联反应中的激酶，将低温信号传递到细胞核内。在细胞核内，一些转录因子被激活，其中BcICE1作为CBF冷响应通路的关键上游转录因子，能够结合到BcCBF1等CBF基因的启动子区域，促进其转录。BcCBF1等CBF蛋白能够识别并结合到下游冷响应基因启动子区域的CRT/DRE元件上，激活一系列冷响应基因的表达，如BcCOR1、BcCOR413-IM1等。这些冷响应基因编码的产物参与了不结球白菜冷胁迫应答的多个过程，如合成抗冻蛋白、积累渗透调节物质、增强抗氧化防御能力等，从而提高不结球白菜的耐寒性。ABA信号通路也在不结球白菜冷胁迫应答中发挥着重要作用。在冷胁迫下，不结球白菜体内ABA含量升高，ABA与受体结合后，激活下游的信号转导途径，使BcABF2等转录因子磷酸化并激活。BcABF2能够结合到冷响应基因启动子区域的ABRE元件上，调控基因表达，进一步增强不结球白菜的耐寒性。本研究发现的以BcNAC10为核心的新信号通路也参与了不结球白菜的冷胁迫应答。在冷胁迫下，BcNAC10表达上调，并从细胞质转移到细胞核中。BcNAC10与蛋白激酶BcMPK10相互作用，BcMPK10磷酸化BcNAC10，增强其转录