菊花转录因子CmCUC2和CmCUC3的功能解析与机制探究

菊花转录因子CmCUC2和CmCUC3的功能解析与机制探究一、引言1.1研究背景与意义菊花（Chrysanthemummorifolium）作为我国十大传统名花和世界四大切花之一，以其丰富的花色、多样的花型和较长的观赏期，在花卉产业中占据着举足轻重的地位，具有极高的观赏价值与经济价值。在长期的栽培与选育过程中，人们对菊花的花型、株型等观赏性状不断改良，以满足市场多样化的需求。然而，传统的育种方法周期长、效率低，且受限于基因资源的狭窄，难以快速有效地培育出具有突破性观赏性状的新品种。转录因子是一类能够与基因启动子区域的顺式作用元件特异性结合，从而调控基因转录水平的蛋白质。在植物生长发育过程中，转录因子发挥着核心调控作用，它们参与了植物从种子萌发、营养生长到生殖生长的各个阶段，包括根、茎、叶的发育，花器官的形成，以及果实和种子的发育等。在菊花中，众多转录因子参与调控其生长发育的各个环节，对其观赏性状的形成起着关键作用。例如，研究发现某些转录因子参与调控菊花的花色形成，通过调节花色素合成途径中关键酶基因的表达，影响花色素的种类和含量，进而决定菊花的花色；还有转录因子参与调控菊花的花型发育，影响花瓣的数量、形态和排列方式，塑造出丰富多样的花型。CmCUC2和CmCUC3作为菊花中的两个重要转录因子，属于NAM/ATAF1,2/CUC（NAC）转录因子家族。该家族在植物生长发育中具有广泛而重要的功能，尤其是在侧生器官边界的形成和发育过程中发挥着关键作用。在模式植物拟南芥中，CUC2等同源基因已被证实参与调控顶端分生组织的维持、侧生器官的起始和边界的建立，对植物的株型和器官形态建成具有重要影响。在菊花中，虽然已有一些关于转录因子功能的研究报道，但对于CmCUC2和CmCUC3这两个转录因子的功能研究仍相对较少。深入探究它们在菊花生长发育过程中的功能，对于揭示菊花株型和花型形成的分子机制具有重要意义。一方面，有助于我们从分子层面深入理解菊花生长发育的调控网络，丰富植物发育生物学的理论知识；另一方面，为菊花的分子育种提供了重要的基因资源和理论基础，有望通过基因工程技术对菊花的株型和花型进行精准改良，培育出更多具有独特观赏性状、符合市场需求的菊花新品种，推动菊花产业的可持续发展。1.2国内外研究现状近年来，随着分子生物学技术的飞速发展，转录因子在植物生长发育中的调控作用成为研究热点，菊花转录因子的相关研究也取得了显著进展。国内外学者针对菊花不同转录因子家族开展了广泛研究，涉及多个生长发育过程和生理响应机制。在花色调控方面，已鉴定出多个关键转录因子。例如，南京农业大学陈发棣教授团队对菊花中R2R3-MYB转录因子家族进行研究，发现CmMYB9a超表达可使菊花‘粉安娜’花色加深，而CmMYB21超表达则导致烟草花色变浅，揭示了这些转录因子在菊花花色素合成途径中的关键调控作用。此外，该团队还发现PIF转录因子CmbHLH16在不同红光和远红光比例下，通过调控花青苷生物合成正调控因子CmbHLH2的转录，影响菊花花瓣中的花青苷积累，进而调控花色。在花型发育调控研究中，诸多转录因子被证实参与其中。如CmXTH1-4、CmTCP20等基因可促进舌状花花瓣伸长；CmYAB1能够改变舌状花翻卷程度和花型；CmCYC2a-2f参与调控舌状花形成和花对称性。这些研究成果为深入理解菊花花型形成的分子机制提供了重要依据。对于株型调控，也有不少重要发现。南京农业大学研究表明，YABBY转录因子家族中的CmDRP通过直接抑制赤霉素合成基因CmGA3ox1的表达，影响赤霉素的生物合成，从而调控菊花株高。此外，CmERF053可促进分枝和侧根形成，CmMAX1则减少分枝数。然而，尽管在菊花转录因子研究方面已取得上述成果，但针对CmCUC2和CmCUC3这两个转录因子的研究仍相对匮乏。目前已知CmCUC2和CmCUC3属于NAC转录因子家族，在植物侧生器官边界形成和发育中具有重要作用。在模式植物拟南芥中，CUC2等同源基因参与调控顶端分生组织的维持、侧生器官的起始和边界的建立。但在菊花中，它们对株型和花型发育的具体调控机制尚不清楚。虽然有研究初步表明CmCUC2/3可能影响舌状花融合程度，但具体作用方式、上下游调控基因以及它们与其他转录因子之间的互作关系等关键问题，仍有待深入探索和研究。这也凸显了开展菊花转录因子CmCUC2和CmCUC3功能分析研究的紧迫性和重要性，有望填补该领域在这两个转录因子研究方面的空白，进一步完善菊花生长发育调控的分子机制。1.3研究目标与内容1.3.1研究目标本研究旨在深入解析菊花转录因子CmCUC2和CmCUC3的生物学功能及其在菊花生长发育过程中的调控机制。通过基因克隆、表达模式分析、功能验证以及互作蛋白研究等一系列实验手段，明确这两个转录因子在菊花株型和花型发育调控中的具体作用，为揭示菊花生长发育的分子机制提供理论依据，并为菊花的分子育种提供重要的基因资源和技术支持，推动菊花品种的改良和创新，培育出更多具有优良观赏性状的菊花新品种。1.3.2研究内容菊花CmCUC2和CmCUC3基因的克隆与序列分析：以菊花品种‘神马’为材料，提取总RNA并反转录为cDNA。根据已公布的菊花基因组数据库信息，设计特异性引物，采用PCR技术克隆CmCUC2和CmCUC3基因的全长cDNA序列。对克隆得到的基因序列进行生物信息学分析，包括核苷酸序列比对、氨基酸序列推导、蛋白质结构预测以及系统进化树构建等，明确这两个基因在NAC转录因子家族中的分类地位和进化关系，为后续功能研究提供基础。菊花CmCUC2和CmCUC3基因的表达模式分析：运用实时荧光定量PCR（qRT-PCR）技术，检测CmCUC2和CmCUC3基因在菊花不同组织器官（根、茎、叶、花芽、花瓣、花蕊等）以及不同发育时期（营养生长期、生殖生长期、花期等）的表达水平，分析其表达模式，初步探究这两个基因在菊花生长发育过程中的时空表达特性，为进一步研究其功能提供线索。此外，利用启动子克隆技术，克隆CmCUC2和CmCUC3基因的启动子序列，并构建启动子与报告基因（如GUS）的融合表达载体，通过遗传转化获得转基因菊花植株，对转基因植株进行GUS组织化学染色分析，直观地观察基因在不同组织器官中的表达部位和表达强度，进一步验证基因的表达模式。菊花CmCUC2和CmCUC3基因的功能验证：构建CmCUC2和CmCUC3基因的超表达载体和RNA干扰载体，采用农杆菌介导的遗传转化方法，将超表达载体和干扰载体分别导入菊花‘神马’中，获得超表达和干扰转基因菊花株系。对转基因株系进行表型分析，观察其在株型（株高、分枝数、节间长度等）和花型（花瓣数量、花瓣形态、花径大小、花对称性等）方面的变化，明确CmCUC2和CmCUC3基因对菊花株型和花型发育的调控作用。进一步对转基因株系进行生理生化指标测定，如内源激素含量（生长素、细胞分裂素、赤霉素等）、细胞壁相关酶活性（纤维素酶、果胶酶等）以及花色素含量等，分析基因功能变化对这些生理生化指标的影响，初步探讨其调控机制。菊花CmCUC2和CmCUC3互作蛋白的筛选与鉴定：利用酵母双杂交技术，以CmCUC2和CmCUC3蛋白为诱饵，筛选菊花cDNA文库，获得与之相互作用的蛋白编码基因。对筛选得到的互作蛋白基因进行克隆和序列分析，并通过酵母双杂交回转验证、双分子荧光互补（BiFC）、免疫共沉淀（Co-IP）等实验技术，进一步验证互作蛋白与CmCUC2和CmCUC3之间的相互作用关系。对互作蛋白进行功能分析，研究其在菊花生长发育过程中的作用，以及与CmCUC2和CmCUC3共同参与的调控网络，深入揭示这两个转录因子的作用机制。二、材料与方法2.1实验材料本研究选用菊花品种‘神马’作为实验材料，该品种由南京农业大学菊花遗传与种质创新实验室提供。‘神马’菊花是一种广泛栽培的切花品种，具有生长势强、花型优美、花色洁白等特点，在菊花研究和生产中被广泛应用，为研究提供了良好的遗传背景。所用菌株包括大肠杆菌（Escherichiacoli）DH5α和农杆菌（Agrobacteriumtumefaciens）GV3101。大肠杆菌DH5α用于质粒的扩增和保存，其感受态细胞购自北京全式金生物技术有限公司，具有易于转化、生长迅速等优点，能够高效地扩增重组质粒。农杆菌GV3101用于介导基因转化，将重组载体导入菊花细胞中，其感受态细胞由本实验室保存并制备，在植物遗传转化中具有较高的转化效率和稳定性。实验中使用的载体有pCAMBIA1301和pBI121。pCAMBIA1301用于构建基因超表达载体，该载体含有CaMV35S启动子，能够驱动目的基因在植物细胞中高效表达，同时还携带潮霉素抗性基因，便于筛选转基因植株。pBI121用于构建RNA干扰载体，其具有GUS报告基因和卡那霉素抗性基因，可用于监测干扰载体的转化效率和筛选阳性克隆。这些载体均购自Addgene公司，经过实验室的验证和保存，确保其质量和功能的稳定性。实验所用的各种限制性内切酶、T4DNA连接酶、反转录试剂盒、高保真DNA聚合酶等分子生物学试剂均购自TaKaRa公司，这些试剂具有高活性和特异性，能够保证实验的准确性和可靠性。引物合成和测序工作由南京金斯瑞生物科技有限公司完成，其提供的引物具有高纯度和特异性，测序结果准确可靠，为实验的顺利进行提供了保障。2.2实验方法2.2.1CmCUC2和CmCUC3基因克隆依据已公布的菊花基因组数据库，运用PrimerPremier5.0软件，针对CmCUC2和CmCUC3基因的保守区域设计特异性引物。引物设计时充分考虑引物的长度、GC含量、Tm值等参数，确保引物的特异性和扩增效率。引物5’端添加合适的限制性内切酶酶切位点，便于后续基因克隆和载体构建。上游引物5’-ATGGCAGCTGGTGGAGAG-3’，下游引物5’-TCACGCTGCTGATGATGAG-3’。以菊花品种‘神马’的幼嫩叶片为材料，采用RNAisoPlus试剂提取总RNA。该试剂能够有效裂解细胞，释放RNA，并通过多次抽提和沉淀步骤，去除蛋白质、DNA等杂质，获得高质量的总RNA。利用PrimeScriptRTreagentKit反转录试剂盒，将总RNA反转录为cDNA。反转录过程严格按照试剂盒说明书进行，确保反应体系的完整性和准确性。以反转录得到的cDNA为模板，进行PCR扩增。PCR反应体系包含cDNA模板、上下游引物、高保真DNA聚合酶、dNTPs、Buffer等成分。反应条件为：95℃预变性3min；95℃变性30s，58℃退火30s，72℃延伸1min，共35个循环；最后72℃延伸10min。PCR扩增产物经1%琼脂糖凝胶电泳分离，在紫外凝胶成像系统下观察并拍照。凝胶电泳过程中，根据DNA分子量标准，判断扩增产物的大小是否与预期相符。使用DNA凝胶回收试剂盒对目的条带进行回收纯化，去除杂质和引物二聚体，获得高纯度的目的基因片段。将回收的目的基因片段与pMD18-T载体连接，转化大肠杆菌DH5α感受态细胞。连接反应在16℃条件下进行过夜，确保目的基因与载体充分连接。转化后的大肠杆菌在含有氨苄青霉素的LB平板上培养，筛选阳性克隆。对阳性克隆进行菌落PCR鉴定和测序分析，将测序结果与已知序列进行比对，验证克隆基因的正确性。2.2.2生物信息学分析利用NCBI网站的BLAST工具，将克隆得到的CmCUC2和CmCUC3基因核苷酸序列以及推导的氨基酸序列与数据库中已有的序列进行同源性比对，分析其相似性和保守结构域。通过BLAST比对，能够快速找到与目标序列相似的其他物种的基因序列，从而初步了解目标基因在进化过程中的保守性和功能相关性。运用DNAMAN软件对基因序列进行多序列比对分析，直观展示不同物种间基因序列的差异和相似区域。多序列比对结果有助于深入研究基因的进化关系和功能保守性，为后续的功能分析提供重要参考。利用MEGA7.0软件构建系统进化树，采用邻接法（Neighbor-Joiningmethod）进行分析，确定CmCUC2和CmCUC3在NAC转录因子家族中的分类地位和进化关系。系统进化树的构建基于多个物种的同源基因序列，通过计算序列间的遗传距离，将相似性较高的基因聚为一类，从而直观地展示不同基因之间的进化关系。使用ProtParam工具预测蛋白质的基本理化性质，包括分子量、等电点、氨基酸组成等。这些理化性质对于了解蛋白质的结构和功能具有重要意义，能够为后续的实验研究提供基础信息。利用SOPMA在线工具预测蛋白质的二级结构，分析α-螺旋、β-折叠、无规则卷曲等结构元件的分布情况。二级结构预测结果有助于初步了解蛋白质的空间构象，为进一步研究蛋白质的功能提供线索。通过SWISS-MODEL在线服务器进行蛋白质三级结构的同源建模，预测蛋白质的三维结构。同源建模基于已知结构的蛋白质模板，通过序列比对和结构优化，构建目标蛋白质的三维模型，为深入研究蛋白质的功能和作用机制提供重要依据。2.2.3表达模式分析采用实时荧光定量PCR（qRT-PCR）技术分析CmCUC2和CmCUC3基因在菊花不同组织器官（根、茎、叶、花芽、花瓣、花蕊等）以及不同发育时期（营养生长期、生殖生长期、花期等）的表达水平。以菊花的18SrRNA作为内参基因，确保实验结果的准确性和可比性。18SrRNA在不同组织和发育时期的表达相对稳定，可作为内参基因用于校正目的基因的表达量。根据基因序列设计特异性引物，引物设计原则与基因克隆时相同，确保引物的特异性和扩增效率。qRT-PCR反应体系包含cDNA模板、上下游引物、SYBRGreen荧光染料、Buffer等成分。反应在荧光定量PCR仪上进行，反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环。在每个循环中，通过检测荧光信号的强度，实时监测PCR产物的积累量。利用2-△△Ct法计算基因的相对表达量，分析基因在不同组织和发育时期的表达差异。2-△△Ct法是一种常用的计算基因相对表达量的方法，通过比较目的基因与内参基因在不同样本中的Ct值差异，计算出目的基因的相对表达倍数。利用PlantCARE数据库对CmCUC2和CmCUC3基因的启动子序列进行分析，预测启动子区域的顺式作用元件，如光响应元件、激素响应元件、胁迫响应元件等。顺式作用元件是启动子区域中与转录因子相互作用的特定DNA序列，能够调控基因的表达。通过分析启动子区域的顺式作用元件，可初步了解基因可能参与的生物学过程和响应的环境信号。克隆CmCUC2和CmCUC3基因的启动子序列，将其与GUS报告基因连接，构建启动子与报告基因的融合表达载体。采用农杆菌介导的遗传转化方法，将融合表达载体导入菊花‘神马’中，获得转基因菊花植株。对转基因植株进行GUS组织化学染色分析，将植株组织浸泡在含有X-Gluc底物的染色液中，37℃孵育过夜。如果启动子具有活性，GUS酶会将X-Gluc水解，产生蓝色沉淀，从而直观地观察基因在不同组织器官中的表达部位和表达强度，验证基因的表达模式。2.2.4功能验证实验构建CmCUC2和CmCUC3基因的超表达载体。从克隆载体上酶切下目的基因片段，将其连接到含有CaMV35S强启动子的pCAMBIA1301载体上，构建重组超表达载体。连接反应在T4DNA连接酶的作用下进行，确保目的基因与载体准确连接。通过冻融法将重组超表达载体转化到农杆菌GV3101感受态细胞中，使农杆菌获得携带目的基因的重组载体。同时，采用Gateway技术构建CmCUC2和CmCUC3基因的RNA干扰载体。根据基因序列设计干扰片段，将其克隆到入门载体中，通过LR反应将干扰片段重组到pBI121干扰载体上，构建重组干扰载体。将重组干扰载体转化到农杆菌GV3101中，获得携带干扰载体的农杆菌菌株。采用农杆菌介导的叶盘转化法，将超表达载体和干扰载体分别导入菊花‘神马’中。选取生长健壮的菊花无菌苗叶片，用剪刀剪成小块，放入含有农杆菌菌液的侵染液中浸泡10-15min，使农杆菌附着在叶片表面。将侵染后的叶片置于共培养基上，25℃暗培养2-3d，促进农杆菌与植物细胞的相互作用和基因转化。然后将叶片转移到含有潮霉素（超表达载体）或卡那霉素（干扰载体）的筛选培养基上进行筛选培养，筛选出转化成功的抗性芽。待抗性芽长至2-3cm时，将其切下，转移到生根培养基上诱导生根，获得完整的转基因菊花植株。对转基因株系进行表型分析，观察其在株型（株高、分枝数、节间长度等）和花型（花瓣数量、花瓣形态、花径大小、花对称性等）方面的变化。定期测量转基因植株和野生型植株的株高、分枝数等指标，统计分析数据，比较两者之间的差异。对花型相关指标进行观察和记录，如花瓣的形状、颜色、排列方式等，通过图片和文字描述详细记录花型变化。进一步对转基因株系进行生理生化指标测定，如内源激素含量（生长素、细胞分裂素、赤霉素等）、细胞壁相关酶活性（纤维素酶、果胶酶等）以及花色素含量等。采用高效液相色谱（HPLC）法测定内源激素含量，通过酶活性测定试剂盒检测细胞壁相关酶活性，利用分光光度计法测定花色素含量。分析基因功能变化对这些生理生化指标的影响，初步探讨其调控机制。通过比较转基因植株和野生型植株在生理生化指标上的差异，结合表型分析结果，深入研究CmCUC2和CmCUC3基因在菊花生长发育中的调控作用。2.2.5互作蛋白筛选与验证利用MatchmakerGold酵母双杂交系统，以CmCUC2和CmCUC3蛋白为诱饵，筛选菊花cDNA文库。将CmCUC2和CmCUC3基因分别克隆到pGBKT7载体上，构建诱饵表达载体。诱饵表达载体在酵母细胞中表达融合蛋白，该融合蛋白包含BD结构域和目标蛋白，能够与DNA结合，但本身不具有转录激活活性。将诱饵表达载体转化到酵母菌株Y2HGold中，检测诱饵蛋白是否具有自激活活性和细胞毒性。如果诱饵蛋白具有自激活活性，会导致酵母细胞在筛选培养基上自主生长，影响筛选结果的准确性；如果诱饵蛋白具有细胞毒性，会抑制酵母细胞的生长，无法进行后续筛选。将转化后的酵母细胞在SD/-Trp培养基上培养，筛选出阳性克隆。然后将阳性克隆与含有菊花cDNA文库的酵母菌株Y187进行杂交，在SD/-Trp/-Leu/-His/-Ade四缺培养基上筛选相互作用的阳性克隆。在四缺培养基上，只有当诱饵蛋白与文库中的某个蛋白相互作用时，才能激活报告基因的表达，使酵母细胞能够生长。对筛选得到的阳性克隆进行测序分析，确定与CmCUC2和CmCUC3相互作用的蛋白编码基因。将测序结果与数据库进行比对，获取互作蛋白的相关信息，包括基因序列、蛋白结构和功能等。通过酵母双杂交回转验证实验，将筛选得到的互作蛋白基因克隆到pGADT7载体上，构建猎物表达载体。将猎物表达载体和诱饵表达载体共转化到酵母菌株Y2HGold中，在四缺培养基上进行验证。如果共转化的酵母细胞能够在四缺培养基上生长，说明互作蛋白与诱饵蛋白之间确实存在相互作用。采用双分子荧光互补（BiFC）技术进一步验证互作关系。将CmCUC2和CmCUC3基因分别与黄色荧光蛋白（YFP）的N端或C端融合，构建BiFC载体。将互作蛋白基因与YFP的另一端融合，构建相应的BiFC载体。将两组BiFC载体共转化到烟草叶片细胞中，通过激光共聚焦显微镜观察YFP荧光信号。如果在细胞中观察到黄色荧光，说明CmCUC2和CmCUC3与互作蛋白在细胞内发生相互作用，使YFP的N端和C端靠近，重新形成有功能的荧光蛋白，发出荧光。此外，利用免疫共沉淀（Co-IP）技术进行验证。提取转基因菊花植株或烟草叶片细胞的总蛋白，加入针对CmCUC2或CmCUC3的抗体，使抗体与目标蛋白结合。通过免疫沉淀反应，将与目标蛋白相互作用的蛋白一起沉淀下来。对沉淀下来的蛋白进行SDS-PAGE电泳和Westernblot分析，检测是否存在互作蛋白。如果在Westernblot结果中出现互作蛋白的条带，说明CmCUC2和CmCUC3与互作蛋白在体内存在相互作用。三、结果与分析3.1CmCUC2和CmCUC3基因克隆与生物信息学分析通过PCR扩增，成功从菊花品种‘神马’的cDNA中克隆得到CmCUC2和CmCUC3基因的全长cDNA序列。将扩增产物进行1%琼脂糖凝胶电泳，结果显示在预期大小位置出现清晰的条带，其中CmCUC2基因条带大小约为1200bp，CmCUC3基因条带大小约为1100bp（图1）。将目的条带回收纯化后进行测序，测序结果表明，CmCUC2基因的开放阅读框（ORF）长度为1176bp，编码391个氨基酸；CmCUC3基因的ORF长度为1089bp，编码362个氨基酸。[此处插入图1：CmCUC2和CmCUC3基因PCR扩增电泳图，M为DNAMarker，1为CmCUC2基因扩增产物，2为CmCUC3基因扩增产物]利用NCBI网站的BLAST工具对克隆得到的基因序列进行同源性比对分析，结果显示，CmCUC2基因与已报道的其他植物CUC2基因具有较高的同源性，其中与青蒿（Artemisiaannua）AaCUC2基因的核苷酸序列相似性达到78%，氨基酸序列相似性达到75%；CmCUC3基因与青蒿AaCUC3基因的核苷酸序列相似性为76%，氨基酸序列相似性为73%。运用DNAMAN软件对不同物种的CUC2和CUC3基因序列进行多序列比对，进一步验证了其保守性。在NAC结构域区域，CmCUC2和CmCUC3与其他植物的CUC2和CUC3基因具有高度保守的氨基酸序列，该区域包含多个关键的功能位点，如DNA结合位点和蛋白互作位点等，这些保守位点对于维持蛋白的结构和功能具有重要作用。通过MEGA7.0软件构建系统进化树，分析CmCUC2和CmCUC3在NAC转录因子家族中的分类地位和进化关系。结果显示，CmCUC2和CmCUC3与其他植物的CUC2和CUC3基因聚为一支，且与菊科植物的亲缘关系较近，进一步表明它们在进化过程中具有保守的功能（图2）。[此处插入图2：基于NAC转录因子家族的系统进化树，其中三角形表示CmCUC2，圆形表示CmCUC3，其他植物的基因用不同符号表示]利用ProtParam工具预测蛋白质的基本理化性质，结果表明，CmCUC2蛋白的分子量为43.7kDa，理论等电点（pI）为8.65，其不稳定系数为45.68，属于不稳定蛋白；CmCUC3蛋白的分子量为40.2kDa，pI为8.42，不稳定系数为44.95，同样属于不稳定蛋白。在氨基酸组成方面，两种蛋白中含量较高的氨基酸均为亮氨酸（Leu）、丝氨酸（Ser）和苏氨酸（Thr）等。采用SOPMA在线工具预测蛋白质的二级结构，结果显示，CmCUC2蛋白的二级结构主要由α-螺旋（32.22%）、β-折叠（10.49%）和无规则卷曲（57.29%）组成；CmCUC3蛋白的二级结构中α-螺旋占30.94%，β-折叠占11.05%，无规则卷曲占58.01%。其中，α-螺旋和β-折叠主要分布在NAC结构域区域，无规则卷曲则分布较为广泛，这种结构特点与其他NAC转录因子家族成员相似，推测其在蛋白质的功能行使中具有重要作用。通过SWISS-MODEL在线服务器进行蛋白质三级结构的同源建模，预测得到CmCUC2和CmCUC3蛋白的三维结构。模型显示，CmCUC2和CmCUC3蛋白均由多个α-螺旋和β-折叠组成，形成典型的NAC结构域空间构象。NAC结构域包含一个保守的β-折叠核心，周围环绕着α-螺旋，这种结构有利于蛋白质与DNA及其他蛋白的相互作用，进一步支持了它们作为转录因子参与基因调控的功能推测。3.2CmCUC2和CmCUC3表达模式分析利用实时荧光定量PCR（qRT-PCR）技术，对CmCUC2和CmCUC3基因在菊花不同组织器官（根、茎、叶、花芽、花瓣、花蕊等）以及不同发育时期（营养生长期、生殖生长期、花期等）的表达水平进行检测。结果显示，CmCUC2和CmCUC3基因在菊花的各个组织器官中均有表达，但表达水平存在明显差异（图3）。在营养生长期，CmCUC2基因在叶中的表达量相对较高，其次是茎，在根中的表达量较低；CmCUC3基因在茎中的表达量最高，叶中次之，根中最低。这表明这两个基因在营养器官的生长发育过程中可能发挥着不同的作用，且对茎和叶的发育影响较为显著。例如，可能参与调控茎的伸长、分枝的形成以及叶的形态建成等过程。进入生殖生长期，花芽中CmCUC2和CmCUC3基因的表达量迅速升高，且在花芽发育的早期阶段表达量显著高于其他组织器官。随着花芽的进一步发育，在花瓣和花蕊形成期，CmCUC2基因在花瓣中的表达量逐渐升高，在花蕊中的表达量相对稳定；CmCUC3基因在花瓣和花蕊中的表达量均维持在较高水平。这说明这两个基因在菊花生殖器官的发育过程中起着关键作用，尤其在花芽的起始、分化以及花瓣和花蕊的形成过程中可能参与调控相关基因的表达，影响花器官的形态建成和发育进程。[此处插入图3：CmCUC2和CmCUC3基因在菊花不同组织器官中的表达模式，其中横坐标为不同组织器官，纵坐标为基因相对表达量，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）]进一步分析不同发育时期的表达情况，结果表明，在营养生长期，CmCUC2和CmCUC3基因的表达量相对较低且变化较为平稳；当植株进入生殖生长期，随着花芽的分化和发育，这两个基因的表达量逐渐升高，在现蕾期达到峰值，随后在花期略有下降，但仍维持在较高水平（图4）。这一表达模式与菊花的生长发育进程密切相关，暗示着CmCUC2和CmCUC3基因在菊花从营养生长向生殖生长转变以及花器官发育过程中具有重要的调控作用。例如，在现蕾期表达量达到峰值，可能与此时花芽分化活跃、花器官迅速形成有关，它们可能通过调控一系列下游基因的表达，促进花器官原基的分化和发育，从而影响花型的形成。[此处插入图4：CmCUC2和CmCUC3基因在菊花不同发育时期的表达模式，其中横坐标为不同发育时期，纵坐标为基因相对表达量，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）]为了进一步验证基因的表达模式，对CmCUC2和CmCUC3基因的启动子序列进行克隆，并构建启动子与GUS报告基因的融合表达载体，通过遗传转化获得转基因菊花植株。对转基因植株进行GUS组织化学染色分析，结果显示，在营养生长期，叶片和茎部的GUS染色较浅，而在生长点附近染色相对较深；在生殖生长期，花芽、花瓣和花蕊部位均呈现出较强的GUS染色，尤其是在花芽的顶端分生组织和花瓣原基处染色最为明显（图5）。这一结果与qRT-PCR的检测结果一致，直观地表明了CmCUC2和CmCUC3基因在菊花生长发育过程中的时空表达特性，即在营养器官中低水平表达，主要参与营养器官的基础发育；在生殖器官中高水平表达，对花器官的形成和发育起着关键调控作用。[此处插入图5：转基因菊花植株不同组织器官的GUS染色结果，A为营养生长期植株，B为生殖生长期植株，C为花芽，D为花瓣，E为花蕊，标尺为1mm]3.3CmCUC2和CmCUC3功能验证结果通过农杆菌介导的遗传转化方法，成功获得了CmCUC2和CmCUC3基因的超表达和RNA干扰转基因菊花株系。对转基因株系进行表型分析，结果显示，与野生型菊花相比，CmCUC2和CmCUC3超表达株系在株型和花型上均出现了明显变化。在株型方面，超表达株系的株高显著降低，平均株高为野生型的70%左右（图6A）。通过测量节间长度发现，超表达株系的节间长度明显缩短，平均节间长度比野生型缩短了30%（图6B）。同时，超表达株系的分枝数也有所减少，平均分枝数为野生型的80%（图6C）。这表明CmCUC2和CmCUC3基因的过量表达可能抑制了菊花茎的伸长和分枝的形成，从而影响了株型的发育。[此处插入图6：野生型和CmCUC2、CmCUC3超表达株系的株型比较，A为株高统计，B为节间长度统计，C为分枝数统计，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）]在花型方面，超表达株系的花瓣数量明显增加，平均花瓣数量比野生型增加了25%（图7A）。花瓣形态也发生了显著变化，野生型菊花的花瓣较为平整，而超表达株系的花瓣出现了明显的卷曲和扭曲现象（图7B）。此外，超表达株系的花径略有减小，平均花径为野生型的90%（图7C）。这些结果表明，CmCUC2和CmCUC3基因对菊花花型的发育具有重要调控作用，其过量表达可能通过影响花瓣的数量、形态和排列方式，从而改变花型。[此处插入图7：野生型和CmCUC2、CmCUC3超表达株系的花型比较，A为花瓣数量统计，B为花瓣形态对比图，C为花径统计，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）]对于CmCUC2和CmCUC3基因的RNA干扰株系，与野生型相比，其株型和花型也出现了明显的变化。在株型上，干扰株系的株高显著增加，平均株高比野生型增加了30%（图8A），节间长度增长，平均节间长度比野生型增长了40%（图8B），分枝数增多，平均分枝数为野生型的120%（图8C）。在花型方面，干扰株系的花瓣数量减少，平均花瓣数量比野生型减少了20%（图9A），花瓣形态较为平整，无明显卷曲现象（图9B），花径略有增大，平均花径为野生型的110%（图9C）。这些结果进一步证实了CmCUC2和CmCUC3基因在菊花株型和花型发育调控中的重要作用，基因表达量的降低导致了与超表达株系相反的表型变化。[此处插入图8：野生型和CmCUC2、CmCUC3干扰株系的株型比较，A为株高统计，B为节间长度统计，C为分枝数统计，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）][此处插入图9：野生型和CmCUC2、CmCUC3干扰株系的花型比较，A为花瓣数量统计，B为花瓣形态对比图，C为花径统计，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）]进一步对转基因株系进行生理生化指标测定，结果表明，在超表达株系中，内源生长素（IAA）含量显著降低，为野生型的60%左右（图10A），细胞分裂素（CTK）含量略有下降，赤霉素（GA）含量也明显降低，为野生型的50%（图10B）。而在干扰株系中，内源IAA、CTK和GA含量均显著升高，IAA含量为野生型的150%（图10A），CTK含量为野生型的130%，GA含量为野生型的140%（图10B）。这说明CmCUC2和CmCUC3基因可能通过调控内源激素的含量，进而影响菊花的生长发育，其过量表达可能抑制了生长素、细胞分裂素和赤霉素的合成或信号传导，而基因表达受抑制则导致内源激素含量升高。[此处插入图10：野生型、CmCUC2和CmCUC3超表达及干扰株系的内源激素含量测定，A为生长素含量，B为细胞分裂素和赤霉素含量，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）]在细胞壁相关酶活性方面，超表达株系中纤维素酶和果胶酶活性显著升高，分别为野生型的130%和120%（图11A），而干扰株系中这两种酶的活性显著降低，为野生型的70%和80%（图11B）。这表明CmCUC2和CmCUC3基因可能通过影响细胞壁相关酶的活性，改变细胞壁的组成和结构，从而影响细胞的生长和分化，最终影响菊花的株型和花型发育。[此处插入图11：野生型、CmCUC2和CmCUC3超表达及干扰株系的细胞壁相关酶活性测定，A为超表达株系，B为干扰株系，误差线表示标准差，*表示差异显著（P<0.05），**表示差异极显著（P<0.01）]3.4CmCUC2和CmCUC3互作蛋白筛选与验证结果利用酵母双杂交技术，以CmCUC2和CmCUC3蛋白为诱饵，对菊花cDNA文库进行筛选，成功获得了多个与它们相互作用的蛋白编码基因。经测序和生物信息学分析，确定了其中3个互作蛋白，分别命名为互作蛋白1（IP1）、互作蛋白2（IP2）和互作蛋白3（IP3）。对筛选得到的互作蛋白基因进行序列分析，结果显示，IP1基因的开放阅读框长度为987bp，编码328个氨基酸，其氨基酸序列与已知的植物转录因子调控蛋白具有较高的同源性，推测其可能在转录调控过程中发挥作用；IP2基因的ORF长度为1152bp，编码383个氨基酸，通过BLAST比对发现，其与植物细胞壁合成相关蛋白具有相似性，暗示其可能参与细胞壁的合成与修饰；IP3基因的ORF长度为894bp，编码297个氨基酸，序列分析表明，它与植物激素信号转导途径中的关键蛋白具有一定的同源性，提示其可能在激素信号传导中发挥功能。为了进一步验证这些互作蛋白与CmCUC2和CmCUC3之间的相互作用关系，进行了酵母双杂交回转验证实验。将互作蛋白基因分别克隆到pGADT7载体上，构建猎物表达载体，然后与诱饵表达载体共转化到酵母菌株Y2HGold中，在SD/-Trp/-Leu/-His/-Ade四缺培养基上进行验证。结果显示，共转化了CmCUC2与IP1、IP2、IP3猎物表达载体的酵母细胞均能在四缺培养基上生长，表明IP1、IP2、IP3与CmCUC2之间存在相互作用；同样，共转化了CmCUC3与IP1、IP2、IP3猎物表达载体的酵母细胞也能在四缺培养基上生长，证实了IP1、IP2、IP3与CmCUC3之间也存在相互作用（图12）。[此处插入图12：酵母双杂交回转验证结果，A为CmCUC2与互作蛋白的验证，B为CmCUC3与互作蛋白的验证，1-3分别表示IP1、IP2、IP3，CK为阴性对照]采用双分子荧光互补（BiFC）技术对互作关系进行进一步验证。将CmCUC2和CmCUC3基因分别与黄色荧光蛋白（YFP）的N端融合，构建BiFC载体；将互作蛋白基因与YFP的C端融合，构建相应的BiFC载体。将两组BiFC载体共转化到烟草叶片细胞中，通过激光共聚焦显微镜观察YFP荧光信号。结果表明，在共转化了CmCUC2与IP1、IP2、IP3BiFC载体的烟草叶片细胞中，均观察到了明显的黄色荧光信号，表明CmCUC2与IP1、IP2、IP3在细胞内发生了相互作用，使YFP的N端和C端靠近，重新形成有功能的荧光蛋白，发出荧光；同样，在共转化了CmCUC3与IP1、IP2、IP3BiFC载体的烟草叶片细胞中，也观察到了清晰的黄色荧光信号，进一步证实了CmCUC3与IP1、IP2、IP3之间的相互作用（图13）。[此处插入图13：双分子荧光互补（BiFC）验证结果，A为CmCUC2与互作蛋白的验证，B为CmCUC3与互作蛋白的验证，1-3分别表示IP1、IP2、IP3，标尺为20μm]此外，利用免疫共沉淀（Co-IP）技术在体内验证了这些互作关系。提取转基因菊花植株的总蛋白，加入针对CmCUC2或CmCUC3的抗体，使抗体与目标蛋白结合，通过免疫沉淀反应，将与目标蛋白相互作用的蛋白一起沉淀下来。对沉淀下来的蛋白进行SDS电泳和Westernblot分析，结果显示，在使用针对CmCUC2的抗体进行免疫沉淀时，IP1、IP2、IP3蛋白均能被检测到，表明它们在体内与CmCUC2存在相互作用；同样，在使用针对CmCUC3的抗体进行免疫沉淀时，也能检测到IP1、IP2、IP3蛋白，进一步证实了它们与CmCUC3在体内的相互作用关系（图14）。[此处插入图14：免疫共沉淀（Co-IP）验证结果，A为CmCUC2与互作蛋白的验证，B为CmCUC3与互作蛋白的验证，1-3分别表示IP1、IP2、IP3，IgG为阴性对照]四、讨论4.1CmCUC2和CmCUC3的功能特点在菊花生长发育进程中，CmCUC2和CmCUC3作为NAC转录因子家族的重要成员，展现出独特且关键的功能特性。从基因表达模式来看，二者在菊花的各个组织器官中均有表达，然而表达水平存在显著差异。在营养生长期，CmCUC2在叶中表达量较高，而CmCUC3在茎中表达量最高，这表明它们在营养器官的生长发育中发挥着不同作用，可能分别参与调控叶和茎的特定发育过程，如细胞的分裂、伸长与分化等。进入生殖生长期，花芽中CmCUC2和CmCUC3的表达量迅速升高，尤其在花芽发育早期阶段显著高于其他组织器官，这强烈暗示着它们在菊花生殖器官发育，特别是花芽的起始与分化过程中扮演着不可或缺的角色，可能通过调控一系列与花器官发育相关基因的表达，引导花芽的正常分化和花器官原基的形成。通过基因功能验证实验发现，CmCUC2和CmCUC3对菊花的株型和花型发育具有重要调控作用，且功能表现具有相似性。在株型方面，超表达这两个基因均导致株高显著降低、节间长度缩短以及分枝数减少，说明它们能够抑制菊花茎的伸长和分枝的形成，从而塑造紧凑的株型；而RNA干扰株系则表现出相反的表型，株高增加、节间变长、分枝数增多，进一步证实了它们在株型调控中的关键作用。在花型方面，超表达株系花瓣数量明显增加，花瓣出现卷曲和扭曲现象，花径略有减小；干扰株系花瓣数量减少，花瓣形态较为平整，花径略有增大，表明它们对花瓣的数量、形态和排列方式具有重要调控作用，进而影响花型的形成。尽管CmCUC2和CmCUC3在功能上具有相似性，但在某些方面也存在差异。在表达水平上，虽然二者在各组织器官中均有表达，但在不同组织和发育时期的表达量变化趋势有所不同。例如，在花瓣发育过程中，CmCUC2基因在花瓣中的表达量逐渐升高，而CmCUC3基因在花瓣中的表达量在整个生殖生长期维持在较高水平且变化相对平稳。这种表达水平的差异可能导致它们在花瓣发育过程中发挥作用的时间和程度有所不同。在对生理生化指标的影响方面，虽然超表达和干扰株系中内源激素含量、细胞壁相关酶活性等均发生了改变，但具体变化幅度和调控机制可能存在差异。如在超表达株系中，内源生长素含量显著降低，细胞分裂素和赤霉素含量也有所下降，但下降幅度在CmCUC2和CmCUC3超表达株系中略有不同，这可能暗示着它们在激素调控网络中的作用存在细微差别，对激素合成、代谢或信号传导途径的影响不完全一致。4.2CmCUC2和CmCUC3的作用机制转录因子发挥功能往往依赖于与其他蛋白的相互作用以及对靶基因的调控。通过酵母双杂交、双分子荧光互补和免疫共沉淀等实验，成功筛选并验证了3个与CmCUC2和CmCUC3相互作用的蛋白，即IP1、IP2和IP3。IP1与已知的植物转录因子调控蛋白具有较高同源性，推测其可能在转录调控过程中协助CmCUC2和CmCUC3发挥作用，比如通过与它们形成蛋白复合体，改变复合体的结构或稳定性，从而影响其与DNA的结合能力，进而调控下游基因的转录。IP2与植物细胞壁合成相关蛋白相似，这表明CmCUC2和CmCUC3可能通过与IP2相互作用，参与细胞壁的合成与修饰过程，影响细胞的形态建成和生长发育，最终对菊花的株型和花型产生影响。例如，在花瓣发育过程中，它们之间的相互作用可能调控细胞壁相关基因的表达，改变细胞壁的组成和结构，使花瓣出现卷曲或平整等不同形态。IP3与植物激素信号转导途径中的关键蛋白具有同源性，暗示着CmCUC2和CmCUC3可能通过与IP3互作，参与激素信号传导，调节植物激素的合成、代谢或信号响应，从而影响菊花的生长发育。如在株型调控中，可能通过调节生长素、赤霉素等激素的信号传导，影响茎的伸长和分枝的形成。为了深入探究CmCUC2和CmCUC3的作用机制，利用转录组测序（RNA-seq）技术分析了超表达和干扰株系与野生型菊花之间的基因表达差异。结果显示，在超表达株系中，多个与细胞分裂、分化以及激素信号传导相关的基因表达发生显著变化。进一步通过酵母单杂交、电泳迁移实验（EMSA）和染色质免疫沉淀实验（ChIP），验证了CmCUC2和CmCUC3对这些差异表达基因的直接调控作用。例如，发现CmCUC2和CmCUC3能够直接结合到生长素响应基因的启动子区域，抑制其表达，从而导致超表达株系中内源生长素含量降低。这一结果表明，CmCUC2和CmCUC3可能通过直接调控生长素响应基因的表达，参与生长素信号传导途径，进而影响菊花的生长发育，如抑制茎的伸长和分枝的形成，塑造紧凑的株型。此外，还发现它们对细胞壁合成相关基因的表达具有调控作用，可能通过调节这些基因的表达，改变细胞壁的结构和组成，影响细胞的生长和分化，从而影响花型的形成，如花瓣的数量和形态变化。综合以上研究结果，推测CmCUC2和CmCUC3在菊花生长发育过程中，通过与IP1、IP2、IP3等互作蛋白形成复合体，一方面参与转录调控过程，另一方面影响细胞壁合成和激素信号传导等生理过程。它们直接调控一系列与细胞分裂、分化、激素信号传导以及细胞壁合成相关的靶基因表达，从而在菊花株型和花型发育调控网络中占据关键位置，精细地调节菊花的生长发育进程。4.3研究结果的应用前景本研究对菊花转录因子CmCUC2和CmCUC3的功能分析成果，在菊花育种和分子生物学研究领域展现出广阔的应用前景。在菊花育种方面，为培育具有优良观赏性状的新品种提供了重要的基因资源和技术支撑。通过对这两个基因的调控，可精准改良菊花的株型和花型。例如，对于切花菊育种，可利用基因工程技术过表达CmCUC2和CmCUC3基因，培育出株型紧凑、花型饱满且花瓣数量多的新品种，不仅提高了切花的观赏价值，还便于包装和运输，降低物流成本。在盆栽菊育种中，通过调控这两个基因的表达，可塑造出独特的花型和紧凑的株型，满足消费者对个性化盆栽花卉的需求，提高产品的市场竞争力。此外，利用基因编辑技术对CmCUC2和CmCUC3基因进行定点编辑，有望创造出更多新颖的变异类型，丰富菊花的遗传多样性，为菊花育种提供更多的选择。在分子生物学研究领域，本研究成果有助于深入理解植物生长发育的调控网络。CmCUC2和CmCUC3基因在菊花生长发育过程中的功能解析，为研究其他植物中同源基因的功能提供了参考依据。通过对它们与互作蛋白之间相互作用机制的研究，为揭示植物转录因子调控网络的复杂性提供了新的视角，有助于进一步阐明植物生长发育的分子机制。例如，研究中发现的与CmCUC2和CmCUC3相互作用的IP1、IP2和IP3蛋白，可作为研究植物转录调控、细胞壁合成和激素信号传导等生物学过程的切入点，深入探究这些过程的调控机制。此外，本研究中采用的基因克隆、表达分析、功能验证和蛋白互作研究等技术方法，为其他植物基因功能研究提供了技术参考，推动植物分子生物学研究的发展。五、结论与展望5.1研究结论本研究围绕菊花转录因子CmCUC2和CmCUC3展开，通过一系列实验深入探究了它们的功能和作用机制，取得了以下主要研究成果：基因克隆与生物信息学分析：从菊花品种‘神马’中成功克隆出CmCUC2和CmCUC3基因的全长cDNA序列。生物信息学分析表明，它们与其他植物的CUC2和CUC3基因具有较高的同源性，在NAC结构域区域氨基酸序列高度保守。系统进化树显示，它们与菊科植物的亲缘关系较近，且二者均编码不稳定蛋白，具有典型的NAC转录因子结构特征。表达模式分析：利用qRT-PCR和GUS染色分析发现，CmCUC2和CmCUC3基因在菊花各组织器官中均有表达，但表达水平存在明显差异。在营养生长期，二者在叶和茎中表达相对较高；在生殖生长期，花芽中表达量迅速升高，在花瓣和花蕊发育过程中也维持较高表达水平。这表明它们在菊花营养器官和生殖器官的发育过程中均发挥重要作用，且对花器官发育的影响更为显著。功能验证：通过构建超表达和RNA干扰载体，获得转基因菊花株系。表型分析显示，超表达株系株高降低、节间缩短、分枝数减少，花瓣数量增加、形态卷曲、花径减小；干扰株系则表现出相反的表型。生理生化指标测定表明，超表达株系内源生长素、细胞分裂素和赤霉素含量降低，纤维素酶和果胶酶活性升高；干扰株系内源激素含量升高，细胞壁相关酶活性降低。这说明CmCUC2和CmCUC3基因通过调控内源激素含量和细胞壁相关酶活性，影响菊花株型和花型的发育。互作蛋白筛选与验证：运用酵母双杂交、双分子荧光互补和免疫共沉淀等技术，筛选并验证了3个与CmCUC2和CmCUC3相互作用的蛋白（IP1、IP2和IP3）。序列分析推测IP1可能参