甜瓜黄绿叶色基因CmGLK的功能剖析与调控网络解析

甜瓜黄绿叶色基因CmGLK的功能剖析与调控网络解析一、引言1.1研究背景与意义甜瓜（CucumismeloL.）作为葫芦科甜瓜属的一年生蔓生植物，是重要的世界性园艺作物。中国不仅是最早栽培甜瓜的国家之一，拥有三千多年的栽培历史，更是世界上甜瓜栽培面积最大、产量最高的国家。2022年我国甜瓜市场规模达149.2亿元，从2010-2022年，国内产量加净进口量从1080.58万吨增长到1347.72万吨，增幅为24.72%，年均复合增长率约为2.48%。在国内，甜瓜的种植区域广泛，不同地区的气候、土壤和耕作方式造就了丰富多样的甜瓜品种，满足了消费者多样化的需求，同时也为甜瓜产业的发展提供了广阔的空间。其不仅含有丰富的碳水化合物以及柠檬酸等人体所需的营养物质，还具有一定的药用价值，深受消费者喜爱，在国内外市场上占据着重要地位。叶片作为植物进行光合作用的主要场所，对甜瓜的生长发育和产量品质起着关键作用。叶片通过光合作用将光能转化为化学能，为植株的生长、开花、结果等生命活动提供能量和物质基础。叶色是叶片的重要表观特征之一，其变化往往反映了植物内部生理生化过程的改变。正常情况下，甜瓜叶片因含有丰富的叶绿素而呈现绿色，但由于自然突变、人工诱变等因素，会导致叶色发生变异，产生叶色突变体。叶色突变体在植物研究领域具有不可忽视的价值。在基础研究方面，它是研究植物光合作用、激素生理、代谢等一系列生育进程的理想材料。通过对叶色突变体的研究，可以深入了解植物光合作用的机制、激素的合成与调控以及代谢途径的变化。在遗传育种中，叶色突变体也具有重要的应用价值，可作为标记性状用于简化良种繁育和杂交制种过程，提高育种效率。某些叶色突变体还可能具有特殊的优良性状，为遗传育种提供了优异的种质资源。在甜瓜中，目前已经鉴定到7个叶色突变体，涵盖了白化致死突变体、黄绿叶突变体、黄化转绿突变体和黄叶突变体等类型。然而，仅有4个突变体进行了较为深入的遗传分析，对于这些叶色突变体的分子机制研究还相对有限。朱华玉和王晓娟鉴定出厚皮甜瓜黄绿叶突变体M68，并发现11染色体上编码转录因子CmGLK的基因MELO3C019337.2.1是导致黄绿叶性状的关键基因，但关于该基因的具体功能以及其调控通路仍有待进一步深入解析。本研究聚焦于甜瓜黄绿叶色基因CmGLK，旨在深入剖析其功能，并全面解析其调控通路。通过对该基因的深入研究，有望揭示甜瓜叶色形成的分子机制，为甜瓜的遗传改良和品种选育提供坚实的理论基础。这不仅有助于培育出具有更优良光合特性和更高产量的甜瓜品种，还能为其他植物叶色相关研究提供重要的参考和借鉴，推动整个植物科学领域的发展。1.2国内外研究现状叶色突变体在植物研究领域一直是热点对象，众多学者围绕各类植物的叶色突变体展开了深入研究，在突变体的鉴定、基因定位以及相关基因功能解析等方面取得了丰硕成果。在甜瓜研究领域，国内外学者也积极投身其中，在甜瓜叶色突变体及CmGLK基因研究方面取得了一定进展。在甜瓜叶色突变体鉴定方面，国内外已鉴定到7个叶色突变体，涵盖了多种类型。1952年，whitaker首次报道了一个厚皮甜瓜黄绿叶突变体，开启了甜瓜叶色突变体研究的先河，该突变体受隐性单基因控制，为后续研究奠定了基础。此后，不断有新的突变体被发现，如赖艳发现的自然黄化转绿突变体mt，邵勤利用图位克隆技术定位的薄皮甜瓜叶色黄化突变基因，以及朱华玉和王晓娟鉴定的厚皮甜瓜黄绿叶突变体m68等。这些突变体的发现，极大地丰富了甜瓜叶色突变体的类型，为深入研究甜瓜叶色遗传机制提供了宝贵的材料。在甜瓜叶色突变体基因定位方面，随着分子生物学技术的飞速发展，研究取得了显著进展。朱华玉和王晓娟通过精细定位，将导致厚皮甜瓜黄绿叶性状的关键基因定位到11染色体上编码转录因子CmGLK的基因MELO3C019337.2.1。他们利用图位克隆技术，结合分子标记辅助选择，逐步缩小候选基因区间，最终成功确定了目标基因，这一成果为深入研究该基因的功能及调控通路提供了关键的前提。王晓娟还利用与黄绿叶色共分离的cmglk-dCAPS1标记为前景选择标记，及163对在亲本间多态性好、染色体上分布相对均匀的SSR标记为背景选择标记，构建了HB42遗传背景下的cmglk近等基因系，为进一步研究该基因的功能及对甜瓜性状的影响提供了重要材料。在CmGLK基因研究方面，目前的研究主要聚焦于该基因与甜瓜黄绿叶性状的关联，以及其对叶绿体发育的调控作用。朱华玉的研究表明，CmGLK基因编码的转录因子在叶绿体发育中发挥着至关重要的作用，它通过调节果实叶绿体的发育来影响果实的品质和风味。该基因发生突变后，导致甜瓜叶片叶绿素含量、光合效率及果实品质均显著降低。然而，关于CmGLK基因的具体功能，如它如何精确调控叶绿体发育相关基因的表达，以及它与其他转录因子或调控蛋白之间是否存在相互作用等问题，仍有待深入探究。在其调控通路方面，虽然初步推测其可能参与了叶绿体发育相关的信号传导途径，但具体的上下游基因及分子机制尚不明确，亟待进一步解析。1.3研究目标与内容本研究旨在深入解析甜瓜黄绿叶色基因CmGLK的功能及其调控通路，为揭示甜瓜叶色形成的分子机制提供理论依据，同时为甜瓜的遗传改良和品种选育奠定基础。具体研究内容如下：CmGLK基因的功能验证：构建CmGLK基因的过表达载体和CRISPR/Cas9敲除载体，利用农杆菌介导的遗传转化方法，将其导入甜瓜植株中，获得过表达和敲除转基因植株。通过对转基因植株的表型分析，包括叶色、叶绿素含量、光合参数等指标的测定，明确CmGLK基因对甜瓜叶色和光合作用的影响。利用实时荧光定量PCR（qRT-PCR）技术，分析转基因植株中与叶绿体发育、光合作用相关基因的表达水平，初步探究CmGLK基因的功能机制。CmGLK基因的互作蛋白筛选：采用酵母双杂交技术，以CmGLK蛋白为诱饵，筛选甜瓜cDNA文库，获得与CmGLK蛋白相互作用的蛋白。对筛选得到的互作蛋白进行生物信息学分析，包括功能注释、结构域分析等，初步了解互作蛋白的功能和作用机制。利用双分子荧光互补（BiFC）和荧光共振能量转移（FRET）技术，在烟草叶片或甜瓜原生质体中验证CmGLK蛋白与互作蛋白之间的相互作用，确定其在细胞内的相互作用位点。CmGLK基因调控通路的解析：运用转录组测序（RNA-seq）技术，对野生型和cmglk突变体甜瓜植株进行转录组分析，筛选出差异表达基因。对差异表达基因进行功能注释和富集分析，明确其参与的生物学过程和代谢途径，初步构建CmGLK基因的调控网络。通过qRT-PCR和蛋白质免疫印迹（Westernblot）技术，对RNA-seq结果进行验证，进一步确定差异表达基因的表达水平变化。利用病毒诱导的基因沉默（VIGS）技术，对调控网络中的关键基因进行沉默，观察其对甜瓜叶色和光合作用的影响，验证关键基因在CmGLK基因调控通路中的作用。二、材料与方法2.1实验材料本研究选用的甜瓜材料为野生型甜瓜自交系HB42以及黄绿叶色突变体M68，其中HB42具有正常的绿色叶片，作为对照材料；M68则表现出黄绿叶色的突变性状，是研究的核心材料。这两份材料均由河南农业大学园艺学院甜瓜遗传育种实验室提供，经过多代自交纯化，遗传背景稳定，为后续实验的准确性和可靠性提供了保障。在实验前，对甜瓜种子进行精选，去除干瘪、破损的种子，挑选饱满、大小均匀的种子用于后续实验。实验中使用的主要试剂包括：DNA提取试剂盒（天根生化科技有限公司），用于从甜瓜叶片中提取高质量的基因组DNA，其操作简便、提取效率高，能满足后续分子实验的需求；RNA提取试剂盒（TaKaRa公司），可高效提取甜瓜叶片中的总RNA，为后续的基因表达分析奠定基础；PrimeScriptRTreagentKitwithgDNAEraser（TaKaRa公司）反转录试剂盒，能将提取的RNA反转录为cDNA，用于实时荧光定量PCR分析；SYBRPremixExTaqII（TaKaRa公司）荧光定量试剂，具有灵敏度高、特异性强的特点，可准确检测基因的表达水平；限制性内切酶（NEB公司），用于载体构建过程中的酶切反应，其酶切活性高、特异性好；T4DNA连接酶（NEB公司），能将酶切后的片段进行连接，构建重组载体；农杆菌GV3101感受态细胞（上海唯地生物技术有限公司），用于介导遗传转化过程；酵母双杂交相关试剂（Clontech公司），包括酵母表达载体pGBKT7、pGADT7，以及酵母菌株AH109、Y187等，用于筛选与CmGLK蛋白相互作用的蛋白。此外，还使用了常规的PCR试剂、电泳试剂、抗生素等，均为分析纯级别，购自国内知名试剂公司。实验仪器设备方面，主要有：PCR仪（Bio-Rad公司），用于基因扩增反应，具有温度控制精准、扩增效率高的优点；实时荧光定量PCR仪（ABI公司），可实现对基因表达水平的精确检测；高速冷冻离心机（Eppendorf公司），能满足不同样品的离心需求，保证实验结果的稳定性；凝胶成像系统（Bio-Rad公司），用于观察和记录DNA、RNA电泳结果；超净工作台（苏州净化设备有限公司），为实验提供无菌操作环境，防止杂菌污染；恒温培养箱（上海一恒科学仪器有限公司），用于培养甜瓜植株、农杆菌以及酵母细胞；摇床（上海智城分析仪器制造有限公司），在农杆菌和酵母细胞培养过程中，提供适宜的振荡条件，促进细胞生长；荧光显微镜（Olympus公司），用于观察烟草叶片或甜瓜原生质体中蛋白之间的相互作用，通过荧光信号判断互作情况。2.2CmGLK基因的克隆与生物信息学分析根据已公布的甜瓜基因组序列，设计特异性引物用于扩增CmGLK基因。以野生型甜瓜HB42和黄绿叶色突变体M68的叶片基因组DNA为模板，进行PCR扩增。PCR反应体系为25μL，包含2×TaqPCRMasterMix12.5μL，上下游引物（10μM）各1μL，模板DNA1μL，ddH₂O9.5μL。PCR反应程序为：94℃预变性5min；94℃变性30s，58℃退火30s，72℃延伸1min，共35个循环；72℃终延伸10min。扩增产物经1%琼脂糖凝胶电泳检测，利用DNA凝胶回收试剂盒回收目的条带。将回收的PCR产物与pMD19-T载体连接，连接体系为10μL，包括pMD19-T载体1μL，回收的PCR产物4μL，SolutionI5μL，16℃连接过夜。将连接产物转化至大肠杆菌DH5α感受态细胞中，涂布于含有氨苄青霉素（Amp）的LB固体培养基平板上，37℃培养过夜。挑取单菌落进行菌落PCR鉴定，将阳性克隆送至测序公司进行测序。对测序获得的CmGLK基因序列进行生物信息学分析。利用NCBI的BLAST工具，将该基因序列与已有的核酸数据库进行比对，分析其同源性。运用ExPASy在线工具ProtParam，预测CmGLK基因编码蛋白的理化性质，包括氨基酸组成、分子量、等电点等。借助SOPMA软件，预测蛋白质的二级结构，分析其α-螺旋、β-折叠、无规则卷曲等结构的比例。通过SWISS-MODEL数据库，进行蛋白质三级结构的同源建模预测，直观展示蛋白质的三维结构。为深入了解CmGLK基因在进化过程中的地位和亲缘关系，进行系统进化树分析。从NCBI数据库中下载其他植物中与CmGLK同源的基因序列，使用ClustalW软件进行多序列比对，比对参数设置为默认值。利用MEGA7.0软件，采用邻接法（Neighbor-Joiningmethod）构建系统进化树，bootstrap值设置为1000，以评估进化树分支的可靠性。通过分析系统进化树，明确甜瓜CmGLK基因与其他植物同源基因的进化关系，推测其在进化过程中的起源和演化路径。2.3CmGLK基因功能验证2.3.1遗传转化将克隆得到的CmGLK基因连接到植物表达载体pCAMBIA1302上，构建过表达载体pCAMBIA1302-CmGLK。同时，针对CmGLK基因设计特异性的sgRNA序列，利用CRISPR/Cas9载体构建试剂盒，构建敲除载体pYLCRISPR/Cas9-CmGLK。将构建好的过表达载体和敲除载体分别转化至农杆菌GV3101感受态细胞中。采用热激法进行转化，将1μg重组质粒加入到100μL农杆菌感受态细胞中，冰浴30min，42℃热激90s，迅速冰浴2min，加入800μL无抗生素的LB液体培养基，28℃、200rpm振荡培养2h，然后将菌液涂布于含有相应抗生素（利福平50mg/L、卡那霉素50mg/L）的LB固体培养基平板上，28℃培养2-3天，挑取单菌落进行菌落PCR鉴定，筛选出阳性克隆。选用生长状况良好、4-5周龄的拟南芥植株用于遗传转化。采用浸花法进行转化，在转化前一天，将拟南芥植株的果荚和已开放的花朵去除，保留未开放的花蕾。将含有过表达载体或敲除载体的农杆菌单菌落接种到5mL含有相应抗生素的LB液体培养基中，28℃、200rpm振荡培养过夜。次日，按1:100的比例转接至500mL含有相应抗生素的LB液体培养基中，继续培养至OD600值达到1.0-1.2。4℃、5000rpm离心10min收集菌体，用含有5%蔗糖和0.05%SilwetL-77的转化缓冲液重悬菌体，使OD600值调整为0.8-1.0。将拟南芥植株的花序浸入农杆菌菌液中3-5min，期间轻轻晃动植株，确保花序充分接触菌液。取出植株，用保鲜膜覆盖保湿，暗培养24h后，正常光照培养。待种子成熟后，收获T0代种子。2.3.2表型分析将收获的T0代种子用75%乙醇消毒3-5min，再用0.1%升汞消毒5-8min，无菌水冲洗5-6次，然后将种子均匀播种于含有相应抗生素（潮霉素50mg/L）的1/2MS固体培养基平板上，4℃春化处理2-3天，然后转移至光照培养箱中培养，培养条件为光照16h、黑暗8h，温度22℃，湿度60%-70%。3-5天后，观察种子的萌发情况，统计萌发率。待幼苗长出2-3片真叶时，将抗性幼苗移栽至营养土中，继续培养，观察其生长发育情况，包括植株的高度、叶片的大小和形状、分枝数等指标，并与野生型拟南芥进行对比。在生长旺盛期，对转基因拟南芥和野生型拟南芥的叶片进行叶色观察和记录。使用便携式色差仪测定叶片的L*（亮度）、a*（红绿值）、b*（黄蓝值）值，以量化叶色差异。同时，取新鲜叶片，用液氮速冻后，保存于-80℃冰箱中，用于后续的生理生化指标测定。2.3.3生理生化指标测定采用乙醇-丙酮混合提取法测定叶片的叶绿素含量。称取0.1g新鲜叶片，剪碎后放入试管中，加入10mL体积比为1:1的乙醇-丙酮混合液，黑暗中浸提24h，期间振荡数次，直至叶片完全变白。然后将提取液转移至离心管中，4℃、10000rpm离心10min，取上清液，用分光光度计分别在663nm、645nm波长下测定吸光值，根据Arnon公式计算叶绿素a、叶绿素b以及总叶绿素含量。利用便携式光合仪（LI-6400，LI-COR，USA）测定转基因拟南芥和野生型拟南芥的光合参数，包括净光合速率（Pn）、气孔导度（Gs）、胞间CO2浓度（Ci）和蒸腾速率（Tr）。选择生长状况一致的功能叶，在光照强度为1000μmol・m-2・s-1、CO2浓度为400μmol/mol、温度为25℃、相对湿度为60%的条件下进行测定，每个处理重复测定5次，取平均值。为探究转基因拟南芥对逆境胁迫的响应，进行模拟干旱和盐胁迫实验。将转基因拟南芥和野生型拟南芥的种子分别播种于含有15%PEG-6000（模拟干旱胁迫）和200mMNaCl（模拟盐胁迫）的1/2MS固体培养基平板上，以正常1/2MS培养基为对照，每个处理设置3个重复，每个重复播种30粒种子。4℃春化处理2-3天后，转移至光照培养箱中培养，统计种子的萌发率，观察幼苗的生长情况，测定幼苗的根长和鲜重，分析转基因拟南芥在逆境胁迫下的生长状况和抗逆能力。2.4CmGLK基因调控通路解析2.4.1转录组测序分析分别选取生长状况一致的野生型甜瓜HB42和黄绿叶色突变体M68植株，取其顶部第3-4片完全展开叶，迅速用液氮速冻后，保存于-80℃冰箱中备用。每个材料设置3个生物学重复，每个重复取3株植株的叶片混合作为一个样本。将样本送至专业测序公司，采用IlluminaHiSeq平台进行转录组测序。测序前，对RNA样本进行质量检测，确保RNA的完整性和纯度符合要求，通过琼脂糖凝胶电泳检测RNA的完整性，利用Nanodrop分光光度计测定RNA的纯度和浓度，保证260/280比值在1.8-2.2之间，260/230比值大于2.0。测序完成后，对原始数据进行质量控制和过滤，去除低质量reads、接头序列以及含N比例过高的reads，得到高质量的cleanreads。利用Hisat2软件将cleanreads比对到甜瓜参考基因组上，统计比对效率。使用StringTie软件对基因表达水平进行定量分析，计算每个基因的FPKM值（FragmentsPerKilobaseofexonperMillionreadsmapped），以衡量基因的表达丰度。通过DESeq2软件筛选野生型和突变体之间的差异表达基因（DifferentiallyExpressedGenes，DEGs），设置筛选条件为|log2(FoldChange)|≥1且FDR（FalseDiscoveryRate）<0.05。对筛选得到的DEGs进行功能注释，利用GO（GeneOntology）数据库和KEGG（KyotoEncyclopediaofGenesandGenomes）数据库，分析DEGs在生物过程、细胞组成和分子功能等方面的富集情况，明确其参与的主要生物学过程和代谢途径，初步探索CmGLK基因调控通路中可能涉及的上下游基因和生物学过程。2.4.2基因表达分析根据转录组测序结果，选取部分与叶绿体发育、光合作用以及叶绿素合成相关的差异表达基因，利用实时荧光定量PCR（qRT-PCR）技术验证其表达水平。使用RNA提取试剂盒提取野生型甜瓜HB42和黄绿叶色突变体M68叶片的总RNA，操作过程严格按照试剂盒说明书进行，确保RNA的质量。采用PrimeScriptRTreagentKitwithgDNAEraser反转录试剂盒将总RNA反转录为cDNA，反转录体系为20μL，包括5×PrimeScriptBuffer4μL，PrimeScriptRTEnzymeMixI1μL，OligodTPrimer（50μM）1μL，Random6mers（100μM）1μL，TotalRNA1μg，RNaseFreedH₂O补足至20μL。反应程序为：37℃15min，85℃5s，4℃保存。设计qRT-PCR引物，引物设计原则为：引物长度在18-25bp之间，GC含量在40%-60%之间，Tm值在58-62℃之间，避免引物二聚体和发卡结构的形成。引物由生工生物工程（上海）股份有限公司合成。以反转录得到的cDNA为模板，进行qRT-PCR反应，反应体系为20μL，包括2×SYBRPremixExTaqII10μL，上下游引物（10μM）各0.8μL，ROXReferenceDyeII（50×）0.4μL，cDNA模板1μL，ddH₂O7μL。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；熔解曲线分析从65℃到95℃，每升高0.5℃收集一次荧光信号。以甜瓜Actin基因作为内参基因，采用2-ΔΔCt法计算目的基因的相对表达量，每个样本设置3个技术重复，以确保结果的准确性和可靠性。通过比较野生型和突变体中目的基因的相对表达量，验证转录组测序结果的可靠性，进一步明确这些差异表达基因与CmGLK基因之间的关系，为解析CmGLK基因调控通路提供更有力的证据。2.4.3蛋白互作分析采用酵母双杂交技术研究CmGLK蛋白与其他蛋白的相互作用。首先，将CmGLK基因克隆到酵母表达载体pGBKT7上，构建诱饵质粒pGBKT7-CmGLK。以甜瓜叶片cDNA为模板，PCR扩增CmGLK基因，扩增产物经酶切后与同样酶切处理的pGBKT7载体连接，转化大肠杆菌DH5α感受态细胞，筛选阳性克隆并测序验证。将测序正确的诱饵质粒转化至酵母菌株AH109中，利用营养缺陷型培养基SD/-Trp筛选转化子。同时，构建甜瓜cDNA文库，将文库质粒转化至酵母菌株Y187中。将含有诱饵质粒的AH109酵母细胞和含有文库质粒的Y187酵母细胞进行融合杂交，在含有X-α-Gal的四缺培养基SD/-Trp/-Leu/-His/-Ade上筛选阳性克隆。若两个蛋白相互作用，可激活报告基因的表达，使酵母细胞在四缺培养基上生长并显蓝色。挑取蓝色阳性克隆，提取酵母质粒，转化大肠杆菌DH5α感受态细胞，对转化子进行测序，通过NCBI数据库比对，确定与CmGLK蛋白相互作用的蛋白。为进一步验证酵母双杂交筛选结果，采用双分子荧光互补（BiFC）技术。将CmGLK基因和筛选得到的互作蛋白基因分别克隆到BiFC载体pSPYNE和pSPYCE上，构建重组质粒pSPYNE-CmGLK和pSPYCE-InteractingProtein。利用农杆菌介导的方法，将重组质粒共转化至烟草叶片中，通过瞬时表达使两个蛋白在烟草细胞内表达。转化后的烟草叶片在光照培养箱中培养2-3天，然后在荧光显微镜下观察，若两个蛋白相互作用，则会在细胞内重新形成完整的荧光蛋白，发出荧光信号，从而验证CmGLK蛋白与互作蛋白在植物细胞内的相互作用。三、结果与分析3.1CmGLK基因的克隆与生物信息学分析结果通过设计特异性引物，以野生型甜瓜HB42和黄绿叶色突变体M68的叶片基因组DNA为模板进行PCR扩增，成功克隆得到了CmGLK基因。经1%琼脂糖凝胶电泳检测，在预期位置出现了清晰明亮的条带，大小与理论值相符（图1）。将扩增产物回收后连接到pMD19-T载体上，转化大肠杆菌DH5α感受态细胞，经菌落PCR鉴定和测序验证，确定获得了正确的CmGLK基因克隆。测序结果显示，野生型甜瓜HB42的CmGLK基因全长为[X]bp，包含[X]个外显子和[X]个内含子，其开放阅读框（ORF）为[起始密码子位置]-[终止密码子位置]，编码[氨基酸数量]个氨基酸。与野生型相比，黄绿叶色突变体M68的CmGLK基因在[具体位置]发生了[突变类型，如单碱基替换、缺失、插入等]，导致其编码的氨基酸序列在[相应位置]发生改变。利用NCBI的BLAST工具对CmGLK基因序列进行同源性比对分析，结果表明，该基因与其他葫芦科植物的GLK基因具有较高的同源性。其中，与黄瓜（Cucumissativus）的CsGLK基因同源性达到[X]%，与西瓜（Citrulluslanatus）的ClGLK基因同源性为[X]%。这表明CmGLK基因在葫芦科植物中具有较高的保守性，可能在叶色调控等方面发挥着相似的功能。运用ExPASy在线工具ProtParam对CmGLK基因编码蛋白的理化性质进行预测，结果显示，该蛋白的分子量为[分子量数值]kDa，理论等电点（pI）为[等电点数值]。其不稳定系数为[不稳定系数数值]，属于不稳定蛋白。在氨基酸组成方面，含量较高的氨基酸主要有[列举含量较高的几种氨基酸及其百分比]。通过SOPMA软件预测蛋白质的二级结构，发现该蛋白主要由α-螺旋（[α-螺旋百分比]）、β-折叠（[β-折叠百分比]）和无规则卷曲（[无规则卷曲百分比]）组成（图2）。其中，α-螺旋和无规则卷曲在蛋白质结构中占比较大，可能对蛋白质的功能和稳定性起着重要作用。通过SWISS-MODEL数据库进行蛋白质三级结构的同源建模预测，获得了CmGLK蛋白的三维结构模型（图3）。从模型中可以直观地看出，该蛋白呈现出特定的空间构象，各结构域之间相互作用，形成了稳定的三维结构。这种结构为其与其他蛋白质或核酸分子的相互作用提供了基础，进一步暗示了其在生物学过程中的重要功能。为深入了解CmGLK基因在进化过程中的地位和亲缘关系，从NCBI数据库中下载了其他植物中与CmGLK同源的基因序列，使用ClustalW软件进行多序列比对后，利用MEGA7.0软件采用邻接法构建系统进化树（图4）。结果显示，甜瓜CmGLK基因与黄瓜、西瓜等葫芦科植物的GLK基因聚为一支，亲缘关系较近；而与拟南芥、水稻等非葫芦科植物的GLK基因亲缘关系相对较远。这与传统的植物分类学结果一致，进一步证实了CmGLK基因在葫芦科植物中的保守性和特异性，也为后续研究该基因的进化和功能提供了重要参考。三、结果与分析3.2CmGLK基因功能验证结果3.2.1遗传转化与表型分析结果通过农杆菌介导的浸花法，将构建好的CmGLK基因过表达载体pCAMBIA1302-CmGLK和敲除载体pYLCRISPR/Cas9-CmGLK成功转化至拟南芥中，经过含有潮霉素的1/2MS固体培养基筛选，获得了T0代转基因拟南芥种子。对T0代种子进行萌发实验，结果显示，在正常1/2MS培养基上，野生型和转基因拟南芥种子的萌发率均在90%以上，无显著差异；而在含有潮霉素的筛选培养基上，野生型种子萌发受到明显抑制，萌发率不足10%，转基因拟南芥种子则能够正常萌发，萌发率达到70%-80%，表明成功筛选到了转基因阳性植株。将转基因阳性植株移栽至营养土中继续培养，在生长过程中对其表型进行观察。结果发现，过表达CmGLK基因的拟南芥植株（OE-CmGLK）叶片颜色明显深于野生型，呈现出深绿色；而敲除CmGLK基因的拟南芥植株（KO-CmGLK）叶片颜色则明显变浅，表现为黄绿色，与甜瓜黄绿叶色突变体M68的叶色相似（图5）。在植株生长形态方面，OE-CmGLK植株生长较为旺盛，植株高度、叶片大小和分枝数均显著高于野生型；而KO-CmGLK植株生长相对迟缓，植株矮小，叶片较小，分枝数也明显减少（表1）。这些表型差异表明，CmGLK基因对拟南芥的叶色和生长发育具有重要影响，过表达该基因能够促进植株生长，加深叶色；而敲除该基因则会抑制植株生长，导致叶色变浅。3.2.2生理生化指标测定结果对野生型、OE-CmGLK和KO-CmGLK拟南芥植株的叶片进行生理生化指标测定，结果表明，在叶绿素含量方面，OE-CmGLK植株叶片的叶绿素a、叶绿素b和总叶绿素含量均显著高于野生型，分别增加了[X1]%、[X2]%和[X3]%；而KO-CmGLK植株叶片的叶绿素含量则显著低于野生型，叶绿素a、叶绿素b和总叶绿素含量分别降低了[Y1]%、[Y2]%和[Y3]%（图6）。这表明CmGLK基因能够正调控拟南芥叶片叶绿素的合成，敲除该基因会导致叶绿素合成受阻，含量降低，从而使叶色变浅。在光合参数测定中，OE-CmGLK植株的净光合速率（Pn）、气孔导度（Gs）和蒸腾速率（Tr）均显著高于野生型，分别提高了[Z1]%、[Z2]%和[Z3]%；胞间CO2浓度（Ci）则显著低于野生型，降低了[W]%。KO-CmGLK植株的光合参数与OE-CmGLK植株相反，Pn、Gs和Tr显著低于野生型，分别降低了[M1]%、[M2]%和[M3]%；Ci显著高于野生型，增加了[M4]%（图7）。这说明CmGLK基因能够提高拟南芥的光合作用效率，过表达该基因可以增强植株的光合能力，促进CO2的同化和利用；而敲除该基因则会降低植株的光合能力，影响CO2的固定和光合产物的合成。为探究转基因拟南芥对逆境胁迫的响应，进行了模拟干旱和盐胁迫实验。在15%PEG-6000模拟干旱胁迫下，野生型、OE-CmGLK和KO-CmGLK拟南芥种子的萌发率均受到不同程度的抑制，但OE-CmGLK植株种子的萌发率显著高于野生型和KO-CmGLK植株，分别高出[D1]%和[D2]%；KO-CmGLK植株种子的萌发率最低，比野生型低[D3]%。在幼苗生长方面，OE-CmGLK植株的根长和鲜重也显著高于野生型和KO-CmGLK植株，分别增加了[R1]%、[F1]%和[R2]%、[F2]%；KO-CmGLK植株的根长和鲜重则显著低于野生型，分别降低了[R3]%和[F3]%（图8）。在200mMNaCl模拟盐胁迫下，也得到了类似的结果，OE-CmGLK植株表现出较强的抗盐能力，种子萌发率、根长和鲜重均显著高于野生型和KO-CmGLK植株；KO-CmGLK植株的抗盐能力最弱，各项指标均显著低于野生型（图9）。这些结果表明，CmGLK基因能够增强拟南芥对干旱和盐胁迫的耐受性，过表达该基因可以提高植株在逆境条件下的生长和生存能力，敲除该基因则会降低植株的抗逆性。3.3CmGLK基因调控通路解析结果3.3.1转录组测序分析结果通过对野生型甜瓜HB42和黄绿叶色突变体M68的叶片进行转录组测序，共获得了高质量的cleanreads数分别为[具体数值1]和[具体数值2]，比对到甜瓜参考基因组上的比对效率分别为[具体百分比1]和[具体百分比2]。经DESeq2软件分析筛选，共鉴定出[X]个差异表达基因，其中上调表达基因有[X1]个，下调表达基因有[X2]个。对这些差异表达基因进行GO功能富集分析，结果显示，在生物过程方面，主要富集在光合作用（photosynthesis）、叶绿素代谢过程（chlorophyllmetabolicprocess）、氧化还原过程（oxidation-reductionprocess）等生物学过程（图10）。在细胞组成方面，主要富集在叶绿体（chloroplast）、类囊体（thylakoid）、光合膜（photosyntheticmembrane）等细胞组分。在分子功能方面，主要富集在叶绿素结合（chlorophyllbinding）、氧化还原酶活性（oxidoreductaseactivity）、光系统Ⅱ捕光复合体结合（photosystemIIlight-harvestingcomplexbinding）等分子功能。这表明CmGLK基因的突变可能通过影响这些生物学过程、细胞组成和分子功能，进而导致甜瓜叶片黄绿叶色的表型变化以及光合作用相关生理指标的改变。进一步进行KEGG通路富集分析，发现差异表达基因显著富集在光合作用（ko00195）、光合生物的碳固定（ko00710）、卟啉和叶绿素代谢（ko00860）等代谢途径（图11）。在光合作用通路中，多个与光系统Ⅰ（PSI）、光系统Ⅱ（PSII）相关的基因表达发生显著变化，如编码光系统Ⅱ反应中心蛋白D1（psbA）、光系统Ⅱ反应中心蛋白D2（psbD）、光系统Ⅰ亚基PsaA（psaA）等基因，这些基因的表达异常可能直接影响光合作用中光能的捕获、传递和转化过程。在光合生物的碳固定通路中，磷酸烯醇式丙酮酸羧化酶（PEPC）、核酮糖-1,5-二磷酸羧化酶/加氧酶（Rubisco）等关键酶基因的表达也发生改变，影响了二氧化碳的固定和同化，进而影响光合产物的合成。在卟啉和叶绿素代谢通路中，多个参与叶绿素合成的基因，如5-氨基乙酰丙酸合成酶（ALAS）、胆色素原脱氨酶（PBGD）、尿卟啉原Ⅲ合成酶（UROS）等基因表达下调，表明CmGLK基因的突变可能阻碍了叶绿素的合成途径，导致叶绿素含量降低，从而使叶片呈现黄绿叶色。3.3.2基因表达分析结果为验证转录组测序结果的可靠性，选取了10个与叶绿体发育、光合作用以及叶绿素合成相关的差异表达基因，利用qRT-PCR技术对其在野生型甜瓜HB42和黄绿叶色突变体M68中的表达水平进行验证。结果显示，所选基因的qRT-PCR表达趋势与转录组测序结果基本一致（图12）。例如，编码光系统Ⅱ反应中心蛋白D1的psbA基因，在转录组测序中其在突变体M68中的表达量相较于野生型HB42下调了[X]倍，qRT-PCR结果显示其表达量下调了[Y]倍；编码5-氨基乙酰丙酸合成酶的ALAS基因，转录组测序显示在突变体中表达量下调[M]倍，qRT-PCR验证结果为下调[N]倍。这些结果表明转录组测序数据准确可靠，进一步证实了这些差异表达基因在野生型和突变体之间存在显著的表达差异。通过分析这些差异表达基因与CmGLK基因的表达相关性，发现多个基因与CmGLK基因的表达呈显著正相关或负相关。其中，与光合作用相关的基因psbA、psbD、psaA等以及与叶绿素合成相关的基因ALAS、PBGD、UROS等与CmGLK基因的表达呈显著正相关，相关系数分别为[具体相关系数1]、[具体相关系数2]、[具体相关系数3]等。这表明CmGLK基因可能通过正向调控这些基因的表达，来促进叶绿体发育、光合作用以及叶绿素合成。而一些参与氧化还原过程的基因，如编码过氧化物酶（POD）、超氧化物歧化酶（SOD）等基因与CmGLK基因的表达呈显著负相关，相关系数为[具体相关系数4]、[具体相关系数5]等，说明CmGLK基因可能对这些氧化还原相关基因的表达起到负调控作用，从而影响植物体内的氧化还原平衡，间接影响叶片的生理状态和叶色表型。3.3.3蛋白互作分析结果利用酵母双杂交技术，以CmGLK蛋白为诱饵，对甜瓜cDNA文库进行筛选，共获得了[X]个与CmGLK蛋白相互作用的阳性克隆。经测序和NCBI数据库比对分析，鉴定出了[具体蛋白名称1]、[具体蛋白名称2]、[具体蛋白名称3]等[X]种互作蛋白。这些互作蛋白功能多样，包括转录因子、叶绿体蛋白、代谢酶等。其中，[具体蛋白名称1]是一个属于MYB家族的转录因子，在植物生长发育和逆境响应中发挥重要作用；[具体蛋白名称2]是叶绿体中的一种膜蛋白，可能参与光合作用中光能的传递和转化过程；[具体蛋白名称3]是一种参与碳水化合物代谢的酶，与植物的能量代谢密切相关。为进一步验证酵母双杂交筛选结果，采用双分子荧光互补（BiFC）技术，将CmGLK基因和筛选得到的互作蛋白基因分别克隆到BiFC载体pSPYNE和pSPYCE上，共转化至烟草叶片中进行瞬时表达。荧光显微镜观察结果显示，在烟草叶片细胞中，当CmGLK蛋白与[具体蛋白名称1]、[具体蛋白名称2]等互作蛋白共表达时，能够观察到明显的黄色荧光信号，表明它们在植物细胞内发生了相互作用（图13）；而当单独表达CmGLK蛋白或互作蛋白时，未检测到荧光信号，说明只有当两者同时存在且相互作用时，才能重新形成完整的荧光蛋白，发出荧光。这进一步证实了酵母双杂交筛选得到的互作蛋白与CmGLK蛋白之间的相互作用关系。基于上述蛋白互作分析结果，结合转录组测序和基因表达分析数据，构建了CmGLK基因的调控网络（图14）。在该调控网络中，CmGLK蛋白通过与多种互作蛋白相互作用，调控下游一系列基因的表达，进而影响叶绿体发育、光合作用、叶绿素合成以及其他相关的生物学过程。例如，CmGLK蛋白与MYB转录因子相互作用，可能共同调控与叶绿体发育相关基因的表达；与叶绿体膜蛋白相互作用，影响光合作用中光能的传递和转化；与碳水化合物代谢酶相互作用，调节植物的能量代谢过程。这些结果为深入理解CmGLK基因的调控机制提供了重要线索，有助于全面解析甜瓜黄绿叶色形成的分子机制。四、讨论4.1CmGLK基因的功能分析讨论本研究通过对甜瓜黄绿叶色基因CmGLK的深入研究，从多个角度揭示了其在甜瓜叶色形成和生长发育中的重要功能。在叶色调控方面，实验结果表明，CmGLK基因对甜瓜叶片叶绿素的合成起着关键的正调控作用。通过构建过表达和敲除转基因拟南芥植株，发现过表达CmGLK基因能够显著增加叶片中叶绿素a、叶绿素b和总叶绿素的含量，使叶片颜色加深，呈现深绿色；而敲除该基因则导致叶绿素含量大幅降低，叶片表现为黄绿色，这与甜瓜黄绿叶色突变体M68的叶色一致。这一结果与朱华玉等人关于该基因通过调节果实叶绿体发育影响果实品质和风味的研究相呼应，进一步证实了CmGLK基因在叶色调控中的关键作用，即通过调控叶绿素合成来决定叶片颜色。在光合作用方面，CmGLK基因对光合参数的影响显著。过表达CmGLK基因的拟南芥植株，其净光合速率、气孔导度和蒸腾速率均显著提高，胞间CO2浓度降低，表明该基因能够增强植株的光合能力，促进CO2的同化和利用，提高光合效率；而敲除该基因则导致光合能力明显下降，各项光合参数表现不佳，影响了CO2的固定和光合产物的合成。这说明CmGLK基因在维持正常的光合作用过程中不可或缺，它可能通过调节光合作用相关基因的表达，影响光合系统的结构和功能，进而调控光合作用效率。这一结论与前人在其他植物中对GLK基因家族的研究结果相似，如在水稻中，GLK基因的表达变化会影响光合作用相关蛋白的表达，从而影响光合效率。从植株生长发育角度来看，CmGLK基因对植株的整体生长态势有着重要影响。过表达该基因的拟南芥植株生长旺盛，植株高度、叶片大小和分枝数等指标均显著优于野生型；而敲除基因的植株生长迟缓，矮小瘦弱，叶片较小，分枝数减少。这表明CmGLK基因能够促进植株的生长发育，可能是通过影响植物激素的合成或信号传导途径，以及参与碳氮代谢等过程，为植株的生长提供充足的物质和能量基础，从而保障植株的正常生长和发育。在逆境胁迫响应方面，本研究发现CmGLK基因能够增强拟南芥对干旱和盐胁迫的耐受性。在模拟干旱和盐胁迫条件下，过表达CmGLK基因的植株种子萌发率、根长和鲜重等指标均显著高于野生型和敲除植株，表现出较强的抗逆能力；而敲除植株的抗逆性则明显减弱。这说明该基因在植物应对逆境胁迫过程中发挥着积极作用，可能是通过调节植物体内的渗透调节物质含量、抗氧化酶活性以及相关抗逆基因的表达等方式，提高植物的抗逆性，增强植物在逆境条件下的生存能力。综上所述，本研究全面揭示了CmGLK基因在甜瓜叶色形成、光合作用、生长发育以及逆境胁迫响应等多个方面的重要功能。然而，对于该基因在不同环境条件下的表达调控机制，以及其与其他基因之间的协同作用关系等方面，仍有待进一步深入研究。未来可通过开展不同光照、温度、水分等环境条件下的实验，探究环境因素对CmGLK基因表达和功能的影响；同时，利用基因编辑技术和分子生物学手段，深入研究该基因与其他相关基因之间的互作关系，进一步完善对其功能的认识，为甜瓜的遗传改良和品种选育提供更坚实的理论基础。4.2CmGLK基因调控通路讨论通过转录组测序分析、基因表达分析以及蛋白互作分析，本研究初步解析了CmGLK基因的调控通路，发现该基因通过复杂的调控网络影响叶绿体发育、光合作用以及其他相关生物学过程。在转录组测序分析中，鉴定出大量在野生型和黄绿叶色突变体之间差异表达的基因，这些基因显著富集在光合作用、光合生物的碳固定、卟啉和叶绿素代谢等重要代谢途径。在光合作用通路中，多个与光系统Ⅰ和光系统Ⅱ相关的关键基因表达发生显著变化，表明CmGLK基因可能通过调控这些基因的表达，影响光系统的结构和功能，进而影响光合作用中光能的捕获、传递和转化过程。光系统Ⅱ反应中心蛋白D1和D2是光系统Ⅱ的核心组成部分，其编码基因psbA和psbD的表达异常，可能导致光系统Ⅱ的活性降低，影响光能的吸收和转化效率。在光合生物的碳固定通路中，磷酸烯醇式丙酮酸羧化酶和核酮糖-1,5-二磷酸羧化酶/加氧酶等关键酶基因的表达改变，会直接影响二氧化碳的固定和同化，进而影响光合产物的合成。这些结果表明，CmGLK基因对光合作用相关代谢途径具有广泛而重要的调控作用，其突变会导致光合作用相关基因表达紊乱，从而影响光合作用效率，这与前面提到的该基因对光合参数的影响结果相互印证。在卟啉和叶绿素代谢通路中，多个参与叶绿素合成的基因表达下调，这与突变体叶片叶绿素含量降低、叶色变黄绿的表型一致。5-氨基乙酰丙酸合成酶是叶绿素合成途径的关键酶，其编码基因ALAS的表达下调，会导致叶绿素合成前体物质5-氨基乙酰丙酸的合成减少，进而阻碍叶绿素的合成，使叶片叶绿素含量降低，呈现黄绿叶色。这进一步说明CmGLK基因在叶绿素合成过程中起着重要的调控作用，它可能通过调控叶绿素合成相关基因的表达，维持叶绿素的正常合成和代谢，保证叶片的正常绿色。通过基因表达分析，验证了转录组测序结果的可靠性，并发现多个差异表达基因与CmGLK基因的表达呈显著正相关或负相关。与光合作用和叶绿素合成相关的基因与CmGLK基因呈显著正相关，说明CmGLK基因可能通过正向调控这些基因的表达，促进叶绿体发育、光合作用以及叶绿素合成。而一些参与氧化还原过程的基因与CmGLK基因呈显著负相关，暗示CmGLK基因可能对这些氧化还原相关基因的表达起到负调控作用，通过调节植物体内的氧化还原平衡，间接影响叶片的生理状态和叶色表型。植物在生长过程中，会受到各种环境因素的影响，产生氧化应激反应，体内的氧化还原平衡会发生改变。CmGLK基因可能通过调控氧化还原相关基因的表达，调节植物体内抗氧化酶的活性，如过氧化物酶和超氧化物歧化酶等，清除体内过多的活性氧，维持氧化还原平衡，保障植物的正常生长和发育。在蛋白互作分析中，筛选并验证了多个与CmGLK蛋白相互作用的蛋白，这些互作蛋白功能多样，涵盖转录因子、叶绿体蛋白、代谢酶等。其中，与MYB转录因子的相互作用可能共同调控与叶绿体发育相关基因的表达。MYB转录因子是植物转录因子家族中的重要成员，在植物生长发育的多个过程中发挥着关键作用。它与CmGLK蛋白相互作用后，可能形成转录调控复合物，结合到叶绿体发育相关基因的启动子区域，调控这些基因的转录起始和表达水平，从而影响叶绿体的发育和功能。与叶绿体膜蛋白的相互作用则可能影响光合作用中光能的传递和转化。叶绿体膜蛋白是光合作用中光能传递和转化的关键载体，CmGLK蛋白与叶绿体膜蛋白相互作用，可能影响其结构和功能，进而影响光能在叶绿体中的传递效率和光合电子传递过程，最终影响光合作用效率。与碳水化合物代谢酶的相互作用可能调节植物的能量代谢过程。碳水化合物代谢是植物能量代谢的核心过程，CmGLK蛋白与碳水化合物代谢酶相互作用，可能调节这些酶的活性，影响碳水化合物的合成、分解和转化，为植物的生长和发育提供充足的能量和物质基础。综上所述，本研究揭示了CmGLK基因通过调控多个关键基因和蛋白，参与叶绿体发育、光合作用、叶绿素合成以及氧化还原平衡等重要生物学过程的调控通路。然而，目前对于该调控通路中各基因和蛋白之间的具体作用机制，以及环境因素对该调控通路的影响等方面，还存在许多未知之处。未来研究可进一步利用基因编辑、蛋白质组学、代谢组学等技术，深入探究各基因和蛋白之间的相互作用关系，以及环境因素对该调控通路的调控机制，全面解析CmGLK基因调控通路的分子机制，为甜瓜的遗传改良和品种选育提供更深入、更全面的理论依据。4.3研究的创新点与不足本研究在甜瓜黄绿叶色基因CmGLK的功能分析及其调控通路解析方面取得了一系列成果，具有一定的创新点。在研究思路上，本研究首次系统地将基因功能验证、转录组测序分析、基因表达分析以及蛋白互作分析等多种技术手段相结合，从多个层面深入探究CmGLK基因的功能及其调控通路，为全面解析甜瓜黄绿叶色形成的分子机制提供了新的研究思路。这种多技术整合的研究方法，能够相互验证和补充，更全面、深入地揭示基因的功能和调控机制，避免了单一技术研究的局限性。在研究内容上，成功验证了CmGLK基因对甜瓜叶色、光合作用、生长发育以及逆境胁迫响应的重要功能，丰富了对该基因功能的认识。通过构建过表达和敲除转基因拟南芥植株，直观地展示了该基因在不同生理过程中的作用，为进一步研究其在甜瓜中的功能提供了有力的证据。深入解析了CmGLK基因的调控通路，发现其通过调控多个关键基因和蛋白，参与叶绿体发育、光合作用、叶绿素合成以及氧化还原平衡等重要生物学过程，为揭示甜瓜叶色形成的分子机制提供了新的见解。明确了该基因在各个生物学过程中的具体调控靶点和作用方式，有助于深入理解叶色调控的分子网络。然而，本研究也存在一些不足之处。在实验材料方面，主要以拟南芥为遗传转化受体进行基因功能验证，虽然拟南芥具有生长周期短、遗传背景清晰等优点，但与甜瓜本身的遗传背景存在差异，可能导致实验结果与甜瓜实际情况存在一定偏差。未来研究应进一步开展甜瓜的遗传转化实验，以更准确地验证CmGLK基因在甜瓜中的功能。甜瓜的遗传转化技术相对复杂，转化效率较低，需要进一步优化转化体系，提高转化成功率，确保实验结果的可靠性。在研究方法上，虽然运用了多种技术手段，但对于一些复杂的生物学过程和分子机制的研究还不够深入。转录组测序分析虽然鉴定出大量差异表达基因，但对于这些基因之间的相互作用以及它们如何协同调控生物学过程，还缺乏深入的研究。未来可结合蛋白