甜瓜品种响应低温胁迫的转录组解析：差异与机制洞察

甜瓜品种响应低温胁迫的转录组解析：差异与机制洞察一、引言1.1研究背景甜瓜（CucumismeloL.）作为葫芦科中重要的园艺作物，在全球范围内广泛种植。中国作为世界上最大的甜瓜生产国，甜瓜的种植面积和产量均居世界首位，其香甜的口感和丰富的营养价值深受消费者喜爱。甜瓜起源于热带及亚热带地区，属于典型的喜温作物，对温度要求较高，其生长发育的适宜温度一般在25-35℃之间。在低温环境下，甜瓜的生长和发育会受到显著影响，导致产量和品质下降。据相关研究表明，在低温胁迫下，甜瓜的产量可降低10%-30%，果实的可溶性固形物含量也会明显降低，严重影响了甜瓜产业的经济效益。在我国，尤其是北方地区，早春季节的低温冷害以及设施栽培中夜间温度过低等问题，经常给甜瓜生产带来严重的威胁。华北地区早春茬甜瓜栽培时，倒春寒等低温天气会严重抑制植株的生长，使植株矮小、叶片发黄、光合作用减弱，导致果实发育不良，品质下降，同时也会推迟果品的上市期，给瓜农造成较大的经济损失。而在设施栽培中，由于设施的保温性能有限，夜间温度常常难以满足甜瓜生长的需求，使得甜瓜容易遭受低温胁迫，影响其生长发育和产量。此外，低温还会使甜瓜更容易受到病虫害的侵袭，进一步降低产量和品质。不同甜瓜品种对低温胁迫的响应存在显著差异，这种差异不仅体现在形态和生理生化指标上，还体现在基因表达水平上。研究不同甜瓜品种在低温胁迫下的转录组，有助于揭示甜瓜耐低温的分子机制，挖掘耐低温相关基因，为甜瓜耐低温品种的选育提供理论依据和基因资源。通过对不同甜瓜品种转录组的分析，可以了解在低温胁迫下哪些基因被激活或抑制，这些基因参与了哪些生物学过程和代谢途径，从而深入理解甜瓜对低温胁迫的适应机制。这对于培育具有强耐冷性的甜瓜新品种，提高甜瓜在低温环境下的产量和品质，推动甜瓜产业的可持续发展具有重要的现实意义。1.2国内外研究现状1.2.1低温胁迫对植物的影响研究低温胁迫作为一种重要的非生物胁迫，对植物的生长发育、生理生化等方面均产生显著影响。在生长发育方面，低温会抑制植物细胞的分裂和伸长，导致植株矮小、根系发育不良、叶片生长受阻等。研究发现，在低温条件下，黄瓜幼苗的株高、茎粗和叶面积的增长速率明显降低，根系的生长也受到抑制，根长和根表面积减小。低温还会影响植物的生殖生长，导致花芽分化异常、开花延迟、授粉受精不良等问题，从而降低结实率和产量。从生理生化角度来看，低温胁迫会破坏植物细胞膜的结构和功能，使细胞膜的透性增加，细胞内物质外渗，导致植物代谢紊乱。低温会影响植物的光合作用，使光合色素含量降低，光合酶活性下降，进而抑制光合电子传递和碳同化过程，减少光合产物的积累。低温还会使植物体内的活性氧（ROS）积累，引发氧化胁迫，对植物细胞造成损伤。为了应对氧化胁迫，植物会启动抗氧化防御系统，提高超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）等抗氧化酶的活性，以清除过量的ROS。低温胁迫还会影响植物体内的激素平衡，如脱落酸（ABA）、生长素（IAA）、细胞分裂素（CTK）等激素的含量和分布发生变化，进而调控植物的生长发育和抗逆反应。在甜瓜研究中，低温胁迫同样对其生长和发育造成诸多不利影响。低温会使甜瓜植株生长缓慢，叶片发黄、卷曲，光合作用减弱，导致果实发育不良，品质下降。有研究表明，低温胁迫下，甜瓜果实的可溶性固形物含量、维生素C含量和可溶性糖含量均显著降低，果实硬度增加，口感变差。不同甜瓜品种对低温胁迫的耐受性存在差异，一些耐低温品种能够通过调节自身的生理生化过程，如提高抗氧化酶活性、积累渗透调节物质等，来减轻低温对植株的伤害，维持相对正常的生长和发育。1.2.2转录组测序技术在植物研究中的应用转录组测序技术（RNA-seq）是一种利用高通量测序技术对细胞或组织中全部或部分mRNA进行测序分析的技术，能够全面、快速地获取基因表达信息，在植物研究中具有广泛的应用。在解析植物基因表达调控方面，转录组测序技术可以帮助研究人员了解在不同生长发育阶段、不同环境条件下植物基因的表达模式和调控机制。通过对不同组织或器官的转录组分析，可以确定哪些基因在特定组织或器官中特异性表达，从而揭示这些基因在植物生长发育中的功能。对拟南芥根、茎、叶等不同组织的转录组测序分析，发现了大量在根中特异性表达的基因，这些基因参与了根的生长、发育和对养分的吸收等过程。在逆境胁迫研究中，转录组测序技术可以鉴定出受胁迫诱导或抑制表达的基因，进一步分析这些基因参与的生物学过程和信号转导途径，有助于揭示植物对逆境胁迫的响应机制。在干旱胁迫下，通过对小麦转录组的分析，发现了一系列与干旱响应相关的基因，这些基因参与了渗透调节、抗氧化防御、激素信号转导等过程。在挖掘关键基因方面，转录组测序技术为筛选和鉴定植物重要性状相关基因提供了有力工具。通过比较不同品种或基因型植物在相同环境条件下的转录组差异，或者同一品种在不同环境条件下的转录组变化，可以发现与目标性状相关的差异表达基因。对水稻耐盐品种和盐敏感品种进行转录组测序分析，筛选出了多个与耐盐性相关的基因，为水稻耐盐分子育种提供了重要的基因资源。在甜瓜耐低温研究中，转录组测序技术也发挥了重要作用。通过对耐低温甜瓜品种和低温敏感品种在低温胁迫下的转录组分析，能够挖掘出与耐低温相关的关键基因，解析甜瓜耐低温的分子机制。已有研究利用转录组测序技术，在甜瓜中鉴定出了一些受低温诱导表达的基因，这些基因涉及ABA信号传导、逆境蛋白合成、转录因子调控等多个方面，为深入研究甜瓜耐低温机制提供了重要线索。1.3研究目的与意义1.3.1研究目的本研究旨在运用转录组测序技术，深入剖析不同甜瓜品种在低温胁迫下的基因表达谱，筛选出耐低温相关的关键基因，揭示甜瓜响应低温胁迫的分子机制，为甜瓜耐低温品种的培育提供坚实的理论依据和丰富的基因资源。具体研究目的如下：分析不同甜瓜品种在低温胁迫下的生理生化特性：通过对耐低温甜瓜品种和低温敏感品种在低温胁迫下的形态指标、生理生化指标进行测定和分析，明确不同品种对低温胁迫的耐受性差异，以及低温胁迫对甜瓜生长发育、细胞膜稳定性、抗氧化系统、渗透调节等生理过程的影响。解析不同甜瓜品种在低温胁迫下的转录组特征：利用转录组测序技术，全面获取不同甜瓜品种在低温胁迫前后的基因表达信息，鉴定差异表达基因，分析差异表达基因的功能注释、富集的生物学过程和代谢途径，从而深入了解甜瓜在低温胁迫下的基因表达调控网络。挖掘甜瓜耐低温相关基因并验证其功能：基于转录组数据分析结果，筛选出与甜瓜耐低温密切相关的候选基因，采用实时荧光定量PCR、基因克隆、转基因等技术对候选基因的表达模式和功能进行验证，为进一步揭示甜瓜耐低温的分子机制奠定基础。1.3.2研究意义理论意义丰富植物抗逆分子生物学理论：深入研究甜瓜在低温胁迫下的转录组变化，有助于揭示植物响应低温胁迫的分子调控机制，进一步丰富植物抗逆分子生物学的理论体系。通过分析甜瓜在低温胁迫下基因表达的变化规律，以及这些基因所参与的生物学过程和代谢途径，能够为理解植物如何感知低温信号、传递信号以及启动防御机制提供新的视角和证据，为深入研究植物抗逆性的分子基础提供参考。完善甜瓜遗传学研究：不同甜瓜品种对低温胁迫的响应存在差异，研究这些差异背后的遗传基础，有助于完善甜瓜遗传学研究。通过鉴定与甜瓜耐低温相关的基因和分子标记，可以为甜瓜遗传图谱的构建、基因定位和克隆等研究提供重要的数据支持，推动甜瓜遗传学研究的深入发展，加深对甜瓜遗传多样性和进化的认识。实践意义指导甜瓜耐低温品种的选育：筛选出的耐低温相关基因和分子标记，为甜瓜耐低温品种的选育提供了重要的理论依据和技术手段。在甜瓜育种过程中，可以利用这些基因和标记进行分子标记辅助选择，提高育种效率，加快耐低温甜瓜品种的选育进程，培育出适应低温环境的优质甜瓜品种，满足市场对不同类型甜瓜的需求，推动甜瓜产业的可持续发展。提高甜瓜生产的经济效益：低温冷害是影响甜瓜产量和品质的重要因素之一。通过培育耐低温甜瓜品种，可以有效减轻低温对甜瓜生产的影响，提高甜瓜在低温环境下的产量和品质，降低生产成本，增加瓜农的经济收入。耐低温甜瓜品种的推广应用，还可以扩大甜瓜的种植区域和种植季节，提高土地资源的利用率，促进农业产业结构的调整和优化，为农村经济发展做出贡献。二、材料与方法2.1试验材料2.1.1甜瓜品种选择选用耐低温品种‘伊丽莎白’和低温敏感品种‘羊角蜜’作为试验材料。‘伊丽莎白’具有早熟、高产、优质、适应性广、抗性强等特点，特别是耐弱光能力强，在低温环境下仍能保持较好的生长态势，果实品质稳定。‘羊角蜜’果实呈长锥形，一端大一端稍细而尖，似羊角状，其口感清甜脆嫩，深受市场欢迎，但对低温较为敏感，在低温胁迫下生长发育容易受到抑制，产量和品质会受到较大影响。选择这两个品种进行研究，能够形成鲜明对比，有利于深入分析不同甜瓜品种对低温胁迫响应的差异，挖掘耐低温相关基因和分子机制。2.1.2材料培养与处理将‘伊丽莎白’和‘羊角蜜’甜瓜种子用55℃温水浸种10分钟，待水温降低至30℃左右时，继续浸种3-4小时，以打破种子休眠并进行消毒。浸种后，用湿纱布包好种子，置于28-30℃恒温培养箱中催芽，待种子露白后，播于装有基质的50孔穴盘中，每穴播1粒种子，播种深度约1-1.5cm。播种后，浇透水，保持基质湿润，并将穴盘置于光照培养箱中培养，光照强度为3000-5000lx，光照时间为16h/d，昼夜温度分别为28℃和22℃。待甜瓜幼苗长至三叶一心时，选取生长健壮、长势一致的幼苗进行低温胁迫处理。将一部分幼苗移入人工气候箱中，设置温度为10℃，光照强度为3000lx，光照时间为12h/d，作为低温处理组；另一部分幼苗继续在原光照培养箱中生长，温度为28℃，光照强度和光照时间不变，作为对照组。分别在低温处理0h、6h、12h、24h和48h时，采集幼苗的叶片，迅速放入液氮中速冻，然后转移至-80℃冰箱中保存，用于后续的生理生化指标测定和转录组测序分析。2.2生理指标测定2.2.1测定指标及方法丙二醛（MDA）含量：采用硫代巴比妥酸（TBA）比色法测定。称取0.5g甜瓜叶片，加入5mL5%三氯乙酸（TCA）和少量石英砂，在冰浴条件下研磨成匀浆，然后在4℃、10000rpm条件下离心10min，取上清液备用。取2mL上清液，加入2mL0.6%TBA溶液（用5%TCA配制），混合均匀后，在沸水浴中加热15min，迅速冷却后再次离心。取上清液，用分光光度计分别测定其在450nm、532nm和600nm波长下的吸光值。根据公式计算MDA含量，以nmol/gFW表示。公式如下：MDA含量(nmol/gFW)=\frac{6.45\times(A_{532}-A_{600})-0.56\timesA_{450}\timesV_{t}}{V_{s}\timesFW}其中，A_{450}、A_{532}和A_{600}分别为450nm、532nm和600nm波长下的吸光值；V_{t}为提取液总体积（mL）；V_{s}为测定用提取液体积（mL）；FW为样品鲜重（g）。游离脯氨酸含量：采用酸性茚三酮法测定。称取0.2g甜瓜叶片，加入5mL3%磺基水杨酸溶液，研磨成匀浆后，在100℃水浴中提取10min，然后冷却并在4℃、8000rpm条件下离心10min，取上清液备用。取2mL上清液，加入2mL冰醋酸和2mL2.5%酸性茚三酮溶液，混合均匀后，在100℃水浴中加热30min，迅速冷却后加入4mL甲苯，振荡萃取5min，静置分层后，取上层甲苯相，用分光光度计测定其在520nm波长下的吸光值。根据标准曲线计算游离脯氨酸含量，以μg/gFW表示。超氧化物歧化酶（SOD）活性：采用氮蓝四唑（NBT）光还原法测定。称取0.5g甜瓜叶片，加入5mL预冷的50mmol/L磷酸缓冲液（pH7.8，含1%聚乙烯吡咯烷酮），在冰浴条件下研磨成匀浆，然后在4℃、12000rpm条件下离心20min，取上清液作为酶液。取3mL反应混合液（含50mmol/L磷酸缓冲液（pH7.8）、13mmol/L甲硫氨酸、75μmol/LNBT、10μmol/LEDTA-Na₂和2μmol/L核黄素），加入50μL酶液，以不加酶液的反应混合液作为对照，在光照强度为4000lx的条件下反应20min，然后用分光光度计测定其在560nm波长下的吸光值。SOD活性以抑制NBT光还原50%为一个酶活性单位（U），用U/gFW表示。过氧化物酶（POD）活性：采用愈创木酚法测定。称取0.5g甜瓜叶片，加入5mL50mmol/L磷酸缓冲液（pH7.0），在冰浴条件下研磨成匀浆，然后在4℃、10000rpm条件下离心15min，取上清液作为酶液。取3mL反应混合液（含50mmol/L磷酸缓冲液（pH7.0）、20mmol/L愈创木酚、10mmol/LH₂O₂），加入50μL酶液，以不加酶液的反应混合液作为对照，在37℃条件下反应5min，然后加入1mL2mol/L硫酸终止反应，用分光光度计测定其在470nm波长下的吸光值。POD活性以每分钟吸光值变化0.01为一个酶活性单位（U），用U/gFW表示。2.2.2数据处理与分析使用Excel2019软件对生理指标测定数据进行整理和初步计算，计算每个处理的平均值和标准差。利用SPSS26.0统计分析软件进行单因素方差分析（One-WayANOVA），比较不同甜瓜品种在不同处理时间下生理指标的差异显著性。当P<0.05时，认为差异显著；当P<0.01时，认为差异极显著。采用Duncan氏新复极差法进行多重比较，进一步明确各处理间的差异情况。通过Origin2021软件绘制柱状图、折线图等，直观展示不同甜瓜品种在低温胁迫下生理指标的变化趋势，以便更清晰地分析低温胁迫对不同甜瓜品种生理指标的影响。2.3转录组测序及分析2.3.1RNA提取与质量检测采用Trizol试剂法提取不同处理时间下‘伊丽莎白’和‘羊角蜜’甜瓜叶片的总RNA。具体步骤如下：取0.1g左右的甜瓜叶片，迅速放入液氮中研磨成粉末状，将研磨好的样品转移至无RNase的1.5mL离心管中，加入1mLTrizol试剂，剧烈振荡混匀，室温静置5min，使样品充分裂解。加入0.2mL氯仿，盖紧管盖，手动剧烈振荡15s，室温孵育2-3min，然后在4℃、12000rpm条件下离心15min。离心后，溶液分为三层，上层为无色水相，中层为白色蛋白层，下层为红色酚氯仿相，小心吸取上层水相转移至新的无RNase离心管中。加入等体积的异丙醇，轻轻颠倒混匀，室温静置10min，使RNA沉淀，4℃、12000rpm条件下离心10min，此时RNA沉淀在管底形成胶状沉淀。弃去上清液，加入1mL75%乙醇（用DEPC水配制），轻轻颠倒洗涤RNA沉淀，4℃、7500rpm条件下离心5min，弃去上清液，重复洗涤一次。将离心管置于超净工作台中，敞口干燥5-10min，使RNA沉淀中的乙醇挥发干净。加入适量无RNase的水，用枪头反复吹打，使RNA沉淀充分溶解，获得的RNA溶液保存于-80℃冰箱中备用。使用超微量紫外分光光度计测定RNA溶液在260nm和280nm波长下的吸光值，计算A260/A280的比值，以评估RNA的纯度。一般来说，高质量的RNA样品其A260/A280比值应在1.8-2.1之间。同时，通过测定260nm处的吸光值，根据公式计算RNA的浓度。公式如下：RNA浓度(μg/μL)=A_{260}\times稀释倍数\times40\div1000取适量RNA样品进行1%琼脂糖凝胶电泳检测，在紫外透射光下观察RNA条带的完整性。完整的RNA样品应呈现出清晰的28S和18S核糖体RNA条带，且28S条带的亮度约为18S条带的2倍，同时无明显的弥散现象。只有RNA质量符合要求，即纯度高、完整性好的样品，才用于后续的转录组测序分析。2.3.2转录组测序将质量合格的RNA样品送至上**因公司进行转录组测序。测序平台采用IlluminaHiSeqXTen，该平台具有高通量、高准确性和高稳定性的特点，能够快速、准确地获取大量的测序数据。测序策略为双端测序（Paired-End，PE），测序读长为150bp。在测序过程中，首先将RNA样品进行片段化处理，然后利用随机引物合成cDNA第一链，再通过DNA聚合酶合成cDNA第二链。对合成的双链cDNA进行末端修复、加A尾和连接测序接头等操作，构建成测序文库。将测序文库进行定量和质量检测后，在IlluminaHiSeqXTen测序平台上进行测序，得到原始测序数据（RawReads）。通过严格的质量控制流程，去除低质量reads、接头序列和污染序列等，得到高质量的测序数据（CleanReads）。最终，每个样品获得的CleanReads数据量均达到6Gb以上，Q30碱基百分比（测序错误率低于0.1%的碱基所占比例）均在90%以上，为后续的数据分析提供了充足且可靠的数据基础。2.3.3数据分析流程数据过滤与质量评估：对测序得到的原始数据进行质量控制，使用FastQC软件对原始数据进行质量评估，查看数据的碱基质量分布、GC含量分布、测序接头污染等情况。利用Trimmomatic软件去除低质量碱基（质量值小于20）、接头序列以及含有N碱基比例过高的reads，得到高质量的CleanReads，确保后续分析的数据准确性。参考基因组比对：将CleanReads与甜瓜参考基因组（Cucumismelov3.6.1）进行比对，使用HISAT2软件进行比对分析。HISAT2是一种高效的比对工具，能够快速准确地将测序reads定位到参考基因组上。通过比对，确定每个reads在基因组上的位置，计算基因的表达量，常用的表达量衡量指标为每千碱基转录本百万映射reads数（FPKM）。差异表达基因分析：利用DESeq2软件对不同处理组（低温处理组和对照组）之间的基因表达数据进行差异分析。以|log2(FoldChange)|≥1且FDR（FalseDiscoveryRate）<0.05为筛选标准，鉴定出在低温胁迫下差异表达的基因。FoldChange表示两组基因表达量的比值，反映了基因表达的变化倍数；FDR是对多重假设检验进行校正后的P值，用于控制假阳性率。基因功能注释：对差异表达基因进行功能注释，将差异表达基因的序列与多个公共数据库进行比对，包括NR（NCBI非冗余蛋白质数据库）、Swiss-Prot（高质量的蛋白质序列数据库）、GO（GeneOntology，基因本体数据库）、KEGG（KyotoEncyclopediaofGenesandGenomes，京都基因与基因组百科全书）等。通过比对，获得差异表达基因的功能信息，了解这些基因参与的生物学过程、分子功能和细胞组成等。代谢通路富集分析：基于KEGG数据库，对差异表达基因进行代谢通路富集分析，使用clusterProfiler软件进行分析。通过富集分析，确定哪些代谢通路在低温胁迫下发生了显著变化，筛选出与低温胁迫响应密切相关的代谢通路，如植物激素信号转导、抗氧化防御系统、碳水化合物代谢等通路，进一步揭示甜瓜在低温胁迫下的生理代谢变化机制。2.3.4荧光定量PCR验证为了验证转录组测序结果的可靠性，选取部分在转录组分析中筛选出的差异表达基因进行荧光定量PCR（qRT-PCR）验证。根据转录组测序得到的基因序列，利用PrimerPremier5.0软件设计特异性引物，引物设计原则包括引物长度一般为18-25bp，GC含量在40%-60%之间，引物Tm值在58-62℃之间，避免引物二聚体和发夹结构的形成等。引物合成由**工生物工程（上海）股份有限公司完成。以提取的总RNA为模板，使用反转录试剂盒（TaKaRaPrimeScriptRTreagentKitwithgDNAEraser）将RNA反转录成cDNA。反转录反应体系和条件按照试剂盒说明书进行操作。qRT-PCR反应采用SYBRGreen荧光染料法，在实时荧光定量PCR仪（ABI7500）上进行。反应体系为20μL，包括2×SYBRGreenMasterMix10μL，上下游引物（10μmol/L）各0.8μL，cDNA模板1μL，ddH₂O7.4μL。反应程序为：95℃预变性30s；95℃变性5s，60℃退火延伸34s，共40个循环；熔解曲线分析从60℃到95℃，每升高0.5℃采集一次荧光信号。每个样品设置3次生物学重复和3次技术重复。以甜瓜的β-actin基因作为内参基因，用于校正和标准化目的基因的表达量。采用2^(-ΔΔCt)法计算目的基因的相对表达量。将qRT-PCR结果与转录组测序结果进行相关性分析，若两者结果趋势一致，则表明转录组测序结果可靠，为后续深入研究甜瓜耐低温的分子机制提供有力的数据支持。三、结果与分析3.1低温胁迫下甜瓜表型与生理生化指标变化3.1.1形态变化观察在正常生长条件下，耐低温品种‘伊丽莎白’和低温敏感品种‘羊角蜜’的甜瓜幼苗均生长健壮，叶片翠绿、舒展，无明显差异。随着低温胁迫时间的延长，两个品种的形态均发生了明显变化。在低温处理6h后，‘羊角蜜’叶片开始出现轻微的萎蔫现象，叶片边缘逐渐失水卷曲，颜色也由深绿变为淡绿。而‘伊丽莎白’叶片仍保持相对挺立，仅叶片边缘稍有卷曲，颜色变化不明显。当低温处理12h时，‘羊角蜜’叶片萎蔫加剧，部分叶片出现发黄现象，生长点生长缓慢。‘伊丽莎白’叶片虽然也出现了一定程度的萎蔫，但发黄现象不明显，生长点仍能保持相对正常的生长。在低温处理24h后，‘羊角蜜’叶片严重萎蔫，大部分叶片发黄，甚至出现坏死斑，植株生长明显受到抑制。相比之下，‘伊丽莎白’叶片虽然也受到一定程度的损伤，但仍有部分叶片保持绿色，植株的生长虽有所减缓，但未受到严重抑制。48h时，‘羊角蜜’植株几乎停止生长，叶片大量枯黄、脱落，表现出对低温胁迫的高度敏感性。‘伊丽莎白’植株虽也遭受了一定程度的伤害，但仍有部分叶片维持一定的生理功能，表现出较强的耐低温能力。通过对两个品种在低温胁迫下形态变化的观察可以直观地看出，‘伊丽莎白’在低温环境下能够更好地维持植株的正常形态和生长，对低温胁迫具有较强的耐受性；而‘羊角蜜’在低温胁迫下形态变化明显，生长受到严重抑制，对低温较为敏感。这些形态变化为进一步研究不同甜瓜品种在低温胁迫下的生理生化响应和分子机制提供了重要的表型依据。3.1.2MDA含量变化丙二醛（MDA）是膜脂过氧化的最终产物，其含量的高低可以反映植物细胞膜脂过氧化的程度和植物受到逆境胁迫伤害的程度。在正常生长条件下，‘伊丽莎白’和‘羊角蜜’甜瓜幼苗叶片中的MDA含量较低，且两者之间无显著差异。随着低温胁迫时间的延长，两个品种叶片中的MDA含量均呈上升趋势（图1）。在低温处理6h时，‘羊角蜜’叶片中的MDA含量显著增加，较对照增加了50.3%，而‘伊丽莎白’叶片中的MDA含量虽也有所增加，但增加幅度相对较小，较对照增加了25.6%。处理12h时，‘羊角蜜’叶片中的MDA含量继续快速上升，达到对照的2.8倍，而‘伊丽莎白’叶片中的MDA含量为对照的1.8倍。当低温处理24h后，‘羊角蜜’叶片中的MDA含量急剧上升，是对照的4.2倍，‘伊丽莎白’叶片中的MDA含量则为对照的2.5倍。在低温处理48h时，‘羊角蜜’叶片中的MDA含量高达对照的5.5倍，而‘伊丽莎白’叶片中的MDA含量为对照的3.2倍。经统计学分析，在低温胁迫的各个时间点，‘羊角蜜’叶片中的MDA含量均显著高于‘伊丽莎白’（P<0.05）。这表明在低温胁迫下，‘羊角蜜’细胞膜脂过氧化程度更为严重，细胞膜受到的损伤更大，而‘伊丽莎白’能够较好地维持细胞膜的稳定性，对低温胁迫具有更强的耐受性。由此可见，MDA含量的变化与甜瓜品种的耐低温能力密切相关，耐低温品种在低温胁迫下能够有效抑制膜脂过氧化，降低MDA的积累，从而减轻低温对植株的伤害。3.1.3游离脯氨酸含量变化游离脯氨酸是植物体内重要的渗透调节物质之一，在逆境胁迫下，植物会积累大量的游离脯氨酸来调节细胞的渗透势，维持细胞的正常生理功能，提高植物的抗逆性。在正常生长条件下，‘伊丽莎白’和‘羊角蜜’甜瓜幼苗叶片中的游离脯氨酸含量较低，且品种间差异不显著。随着低温胁迫时间的延长，两个品种叶片中的游离脯氨酸含量均呈现出先上升后下降的趋势（图2）。在低温处理6h时，‘伊丽莎白’叶片中的游离脯氨酸含量迅速上升，较对照增加了1.2倍，而‘羊角蜜’叶片中的游离脯氨酸含量增加幅度相对较小，较对照增加了0.8倍。处理12h时，‘伊丽莎白’叶片中的游离脯氨酸含量继续上升，达到对照的3.5倍，‘羊角蜜’叶片中的游离脯氨酸含量也显著增加，为对照的2.5倍。当低温处理24h后，‘伊丽莎白’叶片中的游离脯氨酸含量达到峰值，是对照的5.2倍，‘羊角蜜’叶片中的游离脯氨酸含量也在此时达到较高水平，为对照的3.8倍。此后，随着低温胁迫时间的进一步延长，两个品种叶片中的游离脯氨酸含量均开始下降。在低温处理48h时，‘伊丽莎白’叶片中的游离脯氨酸含量仍显著高于对照，为对照的3.0倍，‘羊角蜜’叶片中的游离脯氨酸含量为对照的2.0倍。通过对不同处理时间下两个品种游离脯氨酸含量的比较分析可知，在低温胁迫前期，‘伊丽莎白’叶片中游离脯氨酸含量的上升速度和积累量均显著高于‘羊角蜜’（P<0.05）。这表明‘伊丽莎白’在低温胁迫下能够更快地积累游离脯氨酸，从而更有效地调节细胞渗透势，维持细胞的正常生理功能，增强对低温胁迫的抵抗能力。而‘羊角蜜’在游离脯氨酸的积累速度和积累量上相对较弱，对低温胁迫的适应能力较差。因此，游离脯氨酸含量的变化在甜瓜耐低温过程中发挥着重要作用，耐低温品种通过快速积累游离脯氨酸来提高自身的耐低温能力。3.1.4SOD和POD活性变化超氧化物歧化酶（SOD）和过氧化物酶（POD）是植物抗氧化防御系统中的关键酶，它们能够催化活性氧的歧化反应，将超氧阴离子自由基转化为过氧化氢和氧气，以及将过氧化氢分解为水和氧气，从而清除植物体内过量的活性氧，保护植物细胞免受氧化损伤。在正常生长条件下，‘伊丽莎白’和‘羊角蜜’甜瓜幼苗叶片中的SOD和POD活性相对稳定，且两个品种之间无显著差异。在低温胁迫下，两个品种叶片中的SOD和POD活性均发生了明显变化。随着低温胁迫时间的延长，‘伊丽莎白’叶片中的SOD活性呈现出先上升后略有下降的趋势（图3）。在低温处理6h时，‘伊丽莎白’叶片中的SOD活性显著升高，较对照增加了35.6%，而‘羊角蜜’叶片中的SOD活性虽也有所升高，但增加幅度较小，较对照增加了18.5%。处理12h时，‘伊丽莎白’叶片中的SOD活性继续升高，达到对照的1.8倍，‘羊角蜜’叶片中的SOD活性为对照的1.3倍。当低温处理24h后，‘伊丽莎白’叶片中的SOD活性达到峰值，是对照的2.2倍，此时‘羊角蜜’叶片中的SOD活性也有所增加，为对照的1.6倍。此后，随着低温胁迫时间的进一步延长，‘伊丽莎白’叶片中的SOD活性略有下降，但仍显著高于对照，而‘羊角蜜’叶片中的SOD活性下降较为明显。在低温处理48h时，‘伊丽莎白’叶片中的SOD活性为对照的1.8倍，‘羊角蜜’叶片中的SOD活性仅为对照的1.2倍。对于POD活性，在低温胁迫下，‘伊丽莎白’叶片中的POD活性同样呈现出先上升后下降的趋势（图4）。在低温处理6h时，‘伊丽莎白’叶片中的POD活性显著升高，较对照增加了42.8%，‘羊角蜜’叶片中的POD活性增加幅度相对较小，较对照增加了25.6%。处理12h时，‘伊丽莎白’叶片中的POD活性继续上升，达到对照的2.5倍，‘羊角蜜’叶片中的POD活性为对照的1.8倍。当低温处理24h后，‘伊丽莎白’叶片中的POD活性达到峰值，是对照的3.2倍，‘羊角蜜’叶片中的POD活性也在此时达到较高水平，为对照的2.2倍。随后，随着低温胁迫时间的延长，两个品种叶片中的POD活性均开始下降。在低温处理48h时，‘伊丽莎白’叶片中的POD活性仍显著高于对照，为对照的2.0倍，‘羊角蜜’叶片中的POD活性为对照的1.5倍。经统计学分析，在低温胁迫的大部分时间点，‘伊丽莎白’叶片中的SOD和POD活性均显著高于‘羊角蜜’（P<0.05）。这表明在低温胁迫下，‘伊丽莎白’能够更有效地激活抗氧化防御系统，提高SOD和POD的活性，及时清除体内过量的活性氧，从而减轻氧化损伤，维持细胞的正常生理功能，表现出较强的耐低温能力。而‘羊角蜜’在低温胁迫下抗氧化酶活性的提升幅度相对较小，对活性氧的清除能力较弱，导致细胞受到的氧化损伤较大，生长发育受到严重抑制。因此，SOD和POD活性的变化与甜瓜品种的耐低温能力密切相关，抗氧化酶活性的提高是甜瓜耐低温的重要生理机制之一。3.2转录组测序结果3.2.1测序数据质量评估对‘伊丽莎白’和‘羊角蜜’甜瓜在低温胁迫不同时间点（0h、6h、12h、24h、48h）的叶片样本进行转录组测序，共获得了30个测序文库，每个文库的原始数据量均在6Gb以上。经过严格的数据过滤和质量控制，去除低质量reads、接头序列以及含有N碱基比例过高的reads后，得到高质量的CleanReads，各文库的CleanReads数据量均达到5.8Gb以上，Q30碱基百分比均在90%以上，表明测序数据质量良好，能够满足后续的分析需求。对CleanReads的碱基质量分布进行分析，结果显示，各样本的碱基质量值主要集中在30-40之间，表明测序过程中碱基识别的准确性较高。GC含量分布分析结果表明，各样本的GC含量在45%-50%之间，处于合理范围内，说明测序数据不存在明显的GC偏好性。此外，通过对测序数据的重复率进行评估，发现各样本的重复率均低于5%，表明测序数据的重复性较好，能够真实反映基因的表达情况。3.2.2与参考基因组比对情况将获得的CleanReads与甜瓜参考基因组（Cucumismelov3.6.1）进行比对，结果显示，各样本的比对率均在85%以上，其中‘伊丽莎白’和‘羊角蜜’在不同处理时间下的平均比对率分别为87.6%和86.8%。唯一比对率（即reads唯一比对到基因组上的比例）在75%-80%之间，表明大部分测序数据能够准确地定位到参考基因组上，为后续的基因表达分析提供了可靠的基础。对不同样本在参考基因组上的覆盖度进行分析，结果表明，各样本在参考基因组上的覆盖度较为均匀，大部分基因区域均能被有效覆盖。在不同处理时间下，‘伊丽莎白’和‘羊角蜜’的基因覆盖度差异不显著，说明两个品种在基因覆盖方面不存在明显的差异，测序数据能够全面地反映两个品种在低温胁迫下的基因表达情况。3.2.3差异表达基因筛选与分析以|log2(FoldChange)|≥1且FDR<0.05为筛选标准，对‘伊丽莎白’和‘羊角蜜’在低温胁迫不同时间点与对照（0h）之间的基因表达数据进行差异分析，共筛选出了大量的差异表达基因。随着低温胁迫时间的延长，两个品种的差异表达基因数量均呈现出先增加后减少的趋势。在低温处理6h时，‘伊丽莎白’和‘羊角蜜’分别筛选出1568个和2135个差异表达基因；处理12h时，差异表达基因数量达到峰值，分别为2567个和3089个；随后，随着低温胁迫时间的进一步延长，差异表达基因数量逐渐减少。在低温处理48h时，‘伊丽莎白’和‘羊角蜜’的差异表达基因数量分别为1892个和2256个。对差异表达基因的上调和下调情况进行分析，结果显示，在低温胁迫下，‘伊丽莎白’和‘羊角蜜’的上调基因和下调基因数量在不同处理时间点存在一定差异。在大部分处理时间点，‘羊角蜜’的上调基因数量均高于‘伊丽莎白’，而‘伊丽莎白’的下调基因数量相对较多。例如，在低温处理12h时，‘羊角蜜’的上调基因数量为1856个，下调基因数量为1233个；‘伊丽莎白’的上调基因数量为1489个，下调基因数量为1078个。这表明在低温胁迫下，‘羊角蜜’可能通过上调更多的基因来响应低温胁迫，而‘伊丽莎白’则可能通过下调某些基因来维持自身的生理平衡。3.2.4差异表达基因聚类分析对筛选出的差异表达基因进行聚类分析，以展示不同基因在不同处理时间和品种中的表达模式。通过聚类分析，将差异表达基因分为多个聚类簇，每个聚类簇中的基因具有相似的表达模式。结果显示，在低温胁迫下，不同聚类簇中的基因表达模式存在明显差异。一些聚类簇中的基因在低温处理早期（6h、12h）表达量显著上调，随后逐渐下降，表明这些基因可能在低温胁迫早期发挥重要作用，参与了甜瓜对低温信号的感知和早期响应过程。例如，在聚类簇1中，包含了一些与植物激素信号转导、转录因子调控相关的基因，这些基因在低温处理6h和12h时表达量显著上调，推测它们可能通过调节植物激素信号通路和基因转录，来启动甜瓜对低温胁迫的防御反应。另一些聚类簇中的基因在低温处理后期（24h、48h）表达量持续升高，说明这些基因可能在甜瓜对低温胁迫的长期适应过程中发挥关键作用。如聚类簇2中包含了一些与抗氧化防御、渗透调节相关的基因，这些基因在低温处理24h和48h时表达量逐渐增加，表明它们可能参与了甜瓜在低温胁迫后期对氧化损伤的修复和渗透调节，以维持细胞的正常生理功能。此外，还发现一些基因在‘伊丽莎白’和‘羊角蜜’中的表达模式存在显著差异，这些基因可能是导致两个品种耐低温能力不同的关键因素。例如，在聚类簇3中，某些基因在‘伊丽莎白’中表达量上调，而在‘羊角蜜’中表达量下调，进一步验证了两个品种在低温胁迫下基因表达调控的差异。3.2.5差异表达基因GO功能富集分析为了明确差异表达基因的功能类别，对筛选出的差异表达基因进行GO（GeneOntology）功能富集分析。GO功能注释将基因的功能分为细胞组分（CellularComponent）、分子功能（MolecularFunction）和生物过程（BiologicalProcess）三个类别。在细胞组分方面，差异表达基因主要富集在细胞膜、叶绿体、线粒体等细胞结构相关的条目上。例如，在低温胁迫下，与细胞膜相关的基因表达发生显著变化，这可能与细胞膜的稳定性和物质运输功能在低温环境下的调整有关。细胞膜是细胞与外界环境的屏障，低温胁迫可能导致细胞膜的流动性降低和透性改变，通过调节与细胞膜相关的基因表达，甜瓜可以维持细胞膜的正常功能，减少低温对细胞的伤害。在分子功能方面，差异表达基因主要富集在氧化还原酶活性、转录因子活性、离子结合等条目上。其中，氧化还原酶活性相关基因的富集表明在低温胁迫下，甜瓜通过调节氧化还原酶的表达来应对活性氧的积累，维持细胞内的氧化还原平衡。转录因子活性相关基因的富集则说明转录因子在调控低温胁迫相关基因表达过程中发挥着重要作用。转录因子可以与靶基因的启动子区域结合，调控基因的转录起始，从而影响基因的表达水平。在低温胁迫下，转录因子通过激活或抑制下游靶基因的表达，参与了甜瓜对低温胁迫的响应和适应过程。在生物过程方面，差异表达基因主要富集在植物激素信号转导、氧化应激响应、碳水化合物代谢、光合作用等生物过程相关的条目上。植物激素信号转导相关基因的富集表明植物激素在甜瓜响应低温胁迫过程中起到重要的信号传递作用。例如，脱落酸（ABA）作为一种重要的植物激素，在低温胁迫下，ABA信号通路相关基因的表达发生变化，ABA通过与受体结合，激活下游的信号转导途径，调控相关基因的表达，从而影响植物的生长发育和抗逆性。氧化应激响应相关基因的富集进一步证实了低温胁迫会导致甜瓜体内活性氧积累，引发氧化应激，而甜瓜通过激活氧化应激响应相关基因的表达，增强抗氧化防御能力，减轻氧化损伤。碳水化合物代谢相关基因的变化可能与甜瓜在低温胁迫下能量代谢的调整有关。在低温环境下，甜瓜可能通过调节碳水化合物的合成、分解和转运，为细胞提供足够的能量，以维持正常的生理功能。光合作用相关基因的富集表明低温胁迫对甜瓜的光合作用产生了显著影响。光合作用是植物生长发育的基础，低温可能抑制光合酶的活性、影响光合色素的合成和光合电子传递过程，通过调节光合作用相关基因的表达，甜瓜试图维持一定的光合能力，以满足自身生长和抗逆的需求。3.2.6差异基因的KEGGPathway富集分析基于KEGG（KyotoEncyclopediaofGenesandGenomes）数据库，对差异表达基因进行代谢通路富集分析，以确定显著富集的代谢通路，并探讨这些通路在低温胁迫响应中的作用。结果显示，差异表达基因显著富集在植物激素信号转导、苯丙氨酸代谢、谷胱甘肽代谢、光合作用-天线蛋白、碳代谢等多个代谢通路上。在植物激素信号转导通路中，多个与生长素（IAA）、脱落酸（ABA）、乙烯（ETH）等激素信号转导相关的基因表达发生显著变化。如在ABA信号通路中，PYR/PYL受体基因、PP2C蛋白磷酸酶基因和SnRK2蛋白激酶基因的表达量均发生了改变。在低温胁迫下，ABA含量升高，ABA与PYR/PYL受体结合，抑制PP2C的活性，从而激活SnRK2蛋白激酶，SnRK2进一步磷酸化下游的转录因子和离子通道蛋白等，调控相关基因的表达，使植物产生一系列生理生化变化，以适应低温环境。这表明植物激素信号转导通路在甜瓜响应低温胁迫过程中起着关键的调控作用，通过调节激素信号的传递和响应，甜瓜能够协调自身的生长发育和抗逆反应。苯丙氨酸代谢通路的富集表明该通路在甜瓜应对低温胁迫中具有重要意义。苯丙氨酸是植物体内多种次生代谢产物的前体，在低温胁迫下，苯丙氨酸代谢途径中的关键酶基因表达上调，促进苯丙氨酸向木质素、黄酮类等次生代谢产物的转化。木质素的合成有助于增强细胞壁的强度和稳定性，提高植物对低温胁迫的抵抗能力；黄酮类物质则具有抗氧化活性，能够清除体内过量的活性氧，减轻低温胁迫引起的氧化损伤。因此，苯丙氨酸代谢通路的激活可能是甜瓜适应低温胁迫的一种重要机制。谷胱甘肽代谢通路的显著富集说明谷胱甘肽在甜瓜抵御低温胁迫中发挥着重要作用。谷胱甘肽是一种重要的抗氧化剂，在谷胱甘肽代谢通路中，谷胱甘肽合成酶基因和谷胱甘肽还原酶基因的表达上调，使得谷胱甘肽的合成和再生能力增强。在低温胁迫下，甜瓜体内活性氧积累，谷胱甘肽通过与活性氧反应，将其还原为无害的物质，从而保护细胞免受氧化损伤。此外，谷胱甘肽还参与了植物体内的信号转导过程，调节相关基因的表达，进一步增强植物的抗逆性。光合作用-天线蛋白通路的富集表明低温胁迫对甜瓜的光合作用产生了显著影响。在该通路中，多个与光合天线蛋白相关的基因表达下调，如LHCb1、LHCb2等基因。光合天线蛋白是光合作用中捕获光能的重要结构，其表达下调可能导致光能捕获效率降低，进而影响光合作用的正常进行。这可能是甜瓜在低温胁迫下为了减少能量消耗和避免光氧化损伤而做出的适应性调节。同时，也反映出光合作用在甜瓜响应低温胁迫过程中是一个重要的生理过程，受到低温的显著影响。碳代谢通路的富集说明碳代谢在甜瓜应对低温胁迫中起着重要作用。在低温胁迫下，碳代谢通路中的多个基因表达发生变化，如参与糖酵解、三羧酸循环和磷酸戊糖途径的关键酶基因。这些基因表达的改变可能导致碳水化合物的代谢途径发生调整，使甜瓜能够更有效地利用碳水化合物来提供能量和维持细胞的正常生理功能。例如，在低温胁迫下，糖酵解途径的关键酶基因表达上调，促进葡萄糖的分解，为细胞提供更多的能量；而三羧酸循环中的一些酶基因表达下调，可能是为了减少能量的过度消耗，以维持细胞的能量平衡。3.2.7转录因子分析转录因子在植物应对逆境胁迫的基因表达调控过程中发挥着关键作用。通过对差异表达基因的分析，共鉴定出了多个参与低温胁迫响应的转录因子家族，包括AP2/ERF、bZIP、MYB、NAC等。这些转录因子家族在不同的生物学过程中发挥着重要作用，它们通过与靶基因的启动子区域结合，调控基因的转录起始，从而影响基因的表达水平，参与了甜瓜对低温胁迫的响应和适应过程。在AP2/ERF转录因子家族中，多个成员在低温胁迫下表达显著上调。AP2/ERF转录因子家族是植物特有的一类转录因子，其成员含有一个或两个AP2结构域，在植物生长发育、激素信号转导和逆境胁迫响应等过程中发挥着重要作用。在低温胁迫下，AP2/ERF转录因子可能通过与下游靶基因启动子区域的顺式作用元件结合，激活或抑制相关基因的表达，从而调控甜瓜对低温胁迫的响应。例如，一些AP2/ERF转录因子可能参与了植物激素信号转导途径，通过调节ABA、ETH等激素信号通路相关基因的表达，来影响甜瓜的抗低温能力。bZIP转录因子家族在低温胁迫下也有多个成员的表达发生显著变化。bZIP转录因子含有一个高度保守的碱性亮氨酸拉链结构域，能够与DNA上的特定序列结合，调控基因的转录。在低温胁迫下，bZIP转录因子可能通过与其他转录因子或蛋白相互作用，形成转录调控复合物，共同调节下游靶基因的表达。一些bZIP转录因子可能参与了氧化应激响应和碳水化合物代谢等过程，通过调控相关基因的表达，增强甜瓜在低温胁迫下的抗氧化能力和能量代谢水平。MYB转录因子家族在植物生长发育和逆境胁迫响应中也具有重要作用。在低温胁迫下，MYB转录因子家族中的一些成员表达上调，这些MYB转录因子可能通过与靶基因启动子区域的MYB结合位点结合，激活或抑制相关基因的表达。例如，一些MYB转录因子可能参与了苯丙氨酸代谢途径的调控，通过调节苯丙氨酸解氨酶等关键酶基因的表达，促进木质素和黄酮类等次生代谢产物的合成，增强甜瓜对低温胁迫的抵抗能力。NAC转录因子家族在植物逆境胁迫响应中同样发挥着重要作用。在低温胁迫下，多个NAC转录因子的表达发生显著变化。NAC转录因子含有一个保守的NAC结构域，能够与DNA结合，调控基因的转录。一些NAC转录因子可能参与了植物激素信号转导、氧化应激响应和细胞凋亡等过程，通过调节相关基因的表达，提高甜瓜在低温胁迫下的抗逆性。例如，某些NAC转录因子可能通过激活抗氧化酶基因的表达，增强甜瓜的抗氧化防御能力，减轻低温胁迫引起的氧化损伤。3.3荧光定量PCR验证结果为了验证转录组测序结果的可靠性，选取了10个在转录组分析中差异表达较为显著的基因进行荧光定量PCR（qRT-PCR）验证，这些基因涉及植物激素信号转导、抗氧化防御、渗透调节等多个与低温胁迫响应密切相关的生物学过程。qRT-PCR结果显示，10个基因在‘伊丽莎白’和‘羊角蜜’两个品种中的表达趋势与转录组测序结果基本一致（图5）。以基因CM01为例，在转录组测序中，该基因在‘伊丽莎白’低温处理6h时表达量上调，为对照的2.5倍；在qRT-PCR验证中，其表达量同样上调，为对照的2.3倍。基因CM05在转录组测序中，在‘羊角蜜’低温处理12h时表达量下调，为对照的0.4倍；qRT-PCR结果显示其表达量为对照的0.45倍，表达趋势与转录组测序结果相符。通过对qRT-PCR和转录组测序结果进行相关性分析，得到两者的相关系数R²=0.86（P<0.01），表明qRT-PCR结果与转录组测序结果具有高度的相关性。这说明转录组测序结果准确可靠，能够真实地反映不同甜瓜品种在低温胁迫下的基因表达变化情况，为后续深入研究甜瓜耐低温的分子机制提供了坚实的数据基础。四、讨论4.1低温胁迫下甜瓜品种苗期生理变化与耐低温机制在本研究中，通过对耐低温品种‘伊丽莎白’和低温敏感品种‘羊角蜜’在低温胁迫下的生理指标测定和转录组分析，发现不同甜瓜品种在低温胁迫下苗期生理变化存在显著差异，这些差异与甜瓜的耐低温机制密切相关。从生理指标变化来看，MDA含量是衡量植物细胞膜脂过氧化程度和膜损伤程度的重要指标。在低温胁迫下，‘羊角蜜’叶片中的MDA含量显著高于‘伊丽莎白’，且随着胁迫时间的延长，‘羊角蜜’MDA含量的上升幅度更大。这表明‘羊角蜜’在低温胁迫下细胞膜受到的损伤更为严重，膜脂过氧化程度更高，而‘伊丽莎白’能够较好地维持细胞膜的稳定性，减轻膜脂过氧化程度，从而表现出较强的耐低温能力。这可能是因为‘伊丽莎白’具有更有效的细胞膜保护机制，如含有更多的不饱和脂肪酸，使细胞膜在低温下仍能保持较好的流动性和稳定性，或者具有更强的抗氧化防御系统，能够及时清除膜脂过氧化产生的自由基，减少对细胞膜的损伤。游离脯氨酸作为一种重要的渗透调节物质，在植物应对逆境胁迫中发挥着关键作用。在低温胁迫下，‘伊丽莎白’叶片中游离脯氨酸含量的上升速度和积累量均显著高于‘羊角蜜’。这使得‘伊丽莎白’能够更有效地调节细胞渗透势，维持细胞的膨压和正常生理功能，增强对低温胁迫的抵抗能力。‘伊丽莎白’可能具有更高效的脯氨酸合成途径，或者在低温胁迫下能够更迅速地启动脯氨酸合成相关基因的表达，从而促进脯氨酸的积累。而‘羊角蜜’在游离脯氨酸积累方面相对较弱，可能导致其在低温胁迫下细胞渗透调节能力不足，生长发育受到抑制。SOD和POD是植物抗氧化防御系统中的关键酶，能够清除体内过量的活性氧，保护细胞免受氧化损伤。在低温胁迫下，‘伊丽莎白’叶片中的SOD和POD活性显著高于‘羊角蜜’，且在胁迫前期活性上升更为迅速。这表明‘伊丽莎白’能够更有效地激活抗氧化防御系统，及时清除低温胁迫下产生的过量活性氧，维持细胞内的氧化还原平衡，从而减轻氧化损伤，保证细胞的正常生理功能。‘伊丽莎白’可能具有更高效的抗氧化酶基因表达调控机制，或者其抗氧化酶本身具有更高的活性和稳定性，能够在低温环境下更好地发挥作用。而‘羊角蜜’在低温胁迫下抗氧化酶活性提升幅度较小，对活性氧的清除能力较弱，导致细胞受到的氧化损伤较大，生长发育受到严重抑制。结合转录组数据进一步分析，在低温胁迫下，两个品种的差异表达基因数量和表达模式存在明显差异。GO功能富集分析和KEGGPathway富集分析结果表明，这些差异表达基因主要参与了植物激素信号转导、氧化应激响应、碳水化合物代谢、光合作用等多个与耐低温密切相关的生物学过程和代谢通路。在植物激素信号转导通路中，ABA信号通路在甜瓜响应低温胁迫中起着重要作用。ABA作为一种重要的逆境激素，能够调节植物对低温胁迫的响应。在低温胁迫下，‘伊丽莎白’中与ABA信号转导相关的基因表达变化更为显著，可能通过激活ABA信号通路，促进下游抗逆基因的表达，从而增强对低温胁迫的耐受性。如ABA信号通路中的PYR/PYL受体基因、PP2C蛋白磷酸酶基因和SnRK2蛋白激酶基因的表达量在‘伊丽莎白’中发生了明显改变，这些基因的协同作用可能有助于‘伊丽莎白’感知低温信号，并启动一系列生理生化反应来适应低温环境。氧化应激响应相关基因的表达变化也与甜瓜的耐低温能力密切相关。在低温胁迫下，‘伊丽莎白’中参与氧化应激响应的基因表达上调更为明显，这与‘伊丽莎白’较高的SOD和POD活性相呼应。这些基因可能通过编码抗氧化酶、抗氧化剂合成相关蛋白等，增强‘伊丽莎白’的抗氧化防御能力，减轻低温胁迫引起的氧化损伤。例如，谷胱甘肽代谢通路中的相关基因在‘伊丽莎白’中表达上调，促进了谷胱甘肽的合成和再生，增强了其清除活性氧的能力。碳水化合物代谢通路在甜瓜应对低温胁迫中也具有重要意义。在低温胁迫下，‘伊丽莎白’中参与碳水化合物代谢的基因表达发生了适应性变化，可能通过调节碳水化合物的合成、分解和转运，为细胞提供足够的能量，以维持正常的生理功能。如糖酵解途径的关键酶基因表达上调，促进了葡萄糖的分解，为细胞提供更多的能量；而三羧酸循环中的一些酶基因表达下调，可能是为了减少能量的过度消耗，以维持细胞的能量平衡。光合作用相关基因的表达变化反映了低温胁迫对甜瓜光合作用的影响。在低温胁迫下，两个品种中与光合作用-天线蛋白通路相关的基因表达均下调，但‘伊丽莎白’下调的幅度相对较小。这表明‘伊丽莎白’在低温胁迫下能够更好地维持光合作用，为植物生长和抗逆提供必要的物质和能量基础。‘伊丽莎白’可能通过调节光合作用相关基因的表达，优化光合机构的功能，提高光能利用效率，从而在一定程度上减轻低温对光合作用的抑制。综上所述，不同甜瓜品种在低温胁迫下苗期生理变化的差异是其耐低温能力不同的重要体现。耐低温品种‘伊丽莎白’在低温胁迫下能够通过维持细胞膜稳定性、快速积累游离脯氨酸、高效激活抗氧化防御系统等生理机制，以及调节植物激素信号转导、氧化应激响应、碳水化合物代谢、光合作用等相关基因的表达，来适应低温环境，表现出较强的耐低温能力。而低温敏感品种‘羊角蜜’在这些方面相对较弱，导致其在低温胁迫下生长发育受到严重抑制。本研究结果为深入理解甜瓜耐低温机制提供了重要的理论依据，也为甜瓜耐低温品种的选育提供了有益的参考。4.2差异表达基因与低温胁迫响应通路在低温胁迫下，植物会通过一系列复杂的信号转导途径和基因表达调控来应对逆境，维持自身的生长和发育。对不同甜瓜品种在低温胁迫下的转录组数据分析发现，多个差异表达基因参与了低温胁迫响应通路，这些通路在调节甜瓜耐低温性中发挥着关键作用。4.2.1脱落酸信号转导通路脱落酸（ABA）作为一种重要的植物激素，在植物应对低温胁迫过程中起着核心的信号传递作用。在本研究中，转录组数据分析显示，多个与ABA信号转导相关的基因在低温胁迫下表达发生显著变化。ABA信号转导通路主要包括ABA的感知、信号传递和下游基因的表达调控等环节。在正常生长条件下，ABA含量较低，PYR/PYL受体处于非活性状态，PP2C蛋白磷酸酶能够抑制SnRK2蛋白激酶的活性。当植物受到低温胁迫时，体内ABA含量迅速升高，ABA与PYR/PYL受体结合，形成ABA-PYR/PYL复合物，该复合物与PP2C蛋白磷酸酶结合，抑制其活性，从而解除对SnRK2蛋白激酶的抑制。激活的SnRK2蛋白激酶进一步磷酸化下游的转录因子和离子通道蛋白等，调控相关基因的表达，使植物产生一系列生理生化变化，以适应低温环境。在本研究中，耐低温品种‘伊丽莎白’在低温胁迫下，PYR/PYL受体基因、PP2C蛋白磷酸酶基因和SnRK2蛋白激酶基因的表达量均发生了明显改变。其中，PYR/PYL受体基因的表达上调，可能增强了对ABA的感知能力，使得植物能够更迅速地响应低温信号。PP2C蛋白磷酸酶基因的表达下调，减少了对SnRK2蛋白激酶的抑制作用，从而促进了ABA信号的传递。SnRK2蛋白激酶基因的表达上调，进一步增强了下游基因的磷酸化水平，激活了一系列与耐低温相关的基因表达。这些基因的协同作用可能有助于‘伊丽莎白’感知低温信号，并启动一系列生理生化反应来适应低温环境。相比之下，低温敏感品种‘羊角蜜’在这些基因的表达变化上相对较弱，可能导致其ABA信号转导通路的激活不够迅速和有效，从而影响了其对低温胁迫的耐受性。4.2.2逆境蛋白合成通路逆境蛋白是植物在逆境胁迫下产生的一类蛋白质，它们在保护植物细胞、维持细胞正常生理功能和提高植物抗逆性方面发挥着重要作用。在低温胁迫下，甜瓜中多个与逆境蛋白合成相关的基因表达发生显著变化。其中，脱水素（DHN）基因在低温胁迫下表达上调。脱水素是一种富含亲水性氨基酸的蛋白质，能够与细胞内的水分子结合，形成水合层，从而保护细胞免受脱水伤害。在低温环境下，植物细胞容易失水，脱水素的积累可以提高细胞的保水能力，维持细胞的膨压和正常生理功能。研究表明，过表达脱水素基因的植物在低温胁迫下能够更好地保持细胞膜的完整性，减少细胞内物质的外渗，从而提高植物的耐低温能力。富含甘氨酸的RNA结合蛋白（GRP）基因在低温胁迫下也表达上调。GRP能够与RNA结合，参与mRNA的加工、转运和稳定性调节等过程。在低温胁迫下，GRP可能通过调节与耐低温相关基因的mRNA稳定性和翻译效率，来调控这些基因的表达水平，从而增强植物的耐低温能力。热激蛋白（HSP）基因在低温胁迫下同样表达上调。HSP具有分子伴侣的功能，能够帮助其他蛋白质正确折叠、组装和转运，维持蛋白质的结构和功能稳定性。在低温环境下，蛋白质的结构和功能容易受到影响，HSP的表达上调可以保护细胞内的蛋白质免受低温损伤，确保细胞内各种生理生化反应的正常进行。4.2.3其他低温胁迫响应通路除了脱落酸信号转导通路和逆境蛋白合成通路外，还有许多其他通路也参与了甜瓜对低温胁迫的响应。在苯丙氨酸代谢通路中，多个关键酶基因的表达上调，促进了苯丙氨酸向木质素、黄酮类等次生代谢产物的转化。木质素的合成有助于增强细胞壁的强度和稳定性，提高植物对低温胁迫的抵抗能力。黄酮类物质则具有抗氧化活性，能够清除体内过量的活性氧，减轻低温胁迫引起的氧化损伤。在谷胱甘肽代谢通路中，谷胱甘肽合成酶基因和谷胱甘肽还原酶基因的表达上调，使得谷胱甘肽的合成和再生能力增强。谷胱甘肽是一种重要的抗氧化剂，在低温胁迫下，它能够与活性氧反应，将其还原为无害的物质，从而保护细胞免受氧化损伤。谷胱甘肽还参与了植物体内的信号转导过程，调节相关基因的表达，进一步增强植物的抗逆性。在光合作用相关通路中，多个与光合天线蛋白、光合电子传递和碳同化相关的基因表达下调。这可能是甜瓜在低温胁迫下为了减少能量消耗和避免光氧化损伤而做出的适应性调节。在低温环境下，光合作用受到抑制，光能利用效率降低，通过下调这些基因的表达，甜瓜可以减少光合机构的活性，降低能量消耗，同时避免过多的光能转化为活性氧，对细胞造成氧化损伤。然而，耐低温品种‘伊丽莎白’在这些基因表达下调的程度上相对较小，说明其能够更好地维持光合作用，为植物生长和抗逆提供必要的物质和能量基础。4.3转录因子在低温胁迫响应中的调控作用转录因子在植物响应低温胁迫的基因表达调控网络中扮演着核心角色，它们能够特异性地结合到靶基因的启动子区域，通过激活或抑制基因转录，协调植物体内一系列生理生化反应，从而增强植物对低温胁迫的耐受性。在本研究中，通过转录组数据分析，鉴定出多个在低温胁迫下差异表达的转录因子家族，这些转录因子家族在甜瓜应对低温胁迫过程中发挥着关键作用。AP2/ERF转录因子家族是植物特有的一类转录因子，在低温胁迫响应中具有重要功能。该家族成员含有保守的AP2结构域，能够识别并结合靶基因启动子区域的顺式作用元件，调控基因表达。在低温胁迫下，本研究中多个AP2/ERF转录因子基因在耐低温品种‘伊丽莎白’中表达显著上调，可能通过激活下游与低温耐受相关基因的表达，增强了‘伊丽莎白’对低温胁迫的抵抗能力。研究表明，AP2/ERF转录因子可以调控植物激素信号转导途径，如通过调节ABA、ETH等激素信号通路相关基因的表达，影响植物对低温胁迫的响应。AP2/ERF转录因子还可能参与调控逆境蛋白合成、抗氧化防御系统等相关基因的表达，从而提高植物的耐低温能力。bZIP转录因子家族在植物生长发育和逆境响应中也起着重要作用。其成员通过碱性亮氨酸拉链结构域与DNA上的特定序列结合，形成同源或异源二聚体，调控基因转录。在低温胁迫下，‘伊丽莎白’中多个bZIP转录因子基因表达发生显著变化，这些bZIP转录因子可能与其他转录因子或蛋白相互作用，形成复杂的转录调控复合物，共同调节下游靶基因的表达。一些bZIP转录因子可能参与了氧化应激响应和碳水化合物代谢等过程，通过调控相关基因的表达，增强甜瓜在低温胁迫下的抗氧化能力和能量代谢水平。例如，bZIP转录因子可能通过激活抗氧化酶基因的表达，提高植物的抗氧化防御能力，减轻低温胁迫引起的氧化损伤；同时，通过调节碳水化合物代谢相关基因的表达，为植物提供足够的能量，以维持正常的生理功能。MYB转录因子家族在植物生长发育、次生代谢以及逆境胁迫响应等方面具有广泛的调控作用。在低温胁迫下，‘伊丽莎白’中部分MYB转录因子基因表达上调，这些MYB转录因子可能通过与靶基因启动子区域的MYB结合位点结合，激活或抑制相关基因的表达。一些MYB转录因子参与了苯丙氨酸代谢途径的调控，通过调节苯丙氨酸解氨酶等关键酶基因的表达，促进木质素和黄酮类等次生代谢产物的合成。木质素的合成有助于增强细胞壁的强度和稳定性，提高植物对低温胁迫的抵抗能力；黄酮类物质则具有抗氧化活性，能够清除体内过量的活性氧，减轻低温胁迫引起的氧化损伤。NAC转录因子家族在植物逆境胁迫响应中同样发挥着不可或缺的作用。在低温胁迫下，‘伊丽莎白’中多个NAC转录因子基因的表达发生显著变化。NAC转录因子含有保守的NAC结构域，能够与DNA结合，调控基因的转录。一些NAC转录因子可能参与了植物激素信号转导、氧化应激响应和细胞凋亡等过程，通过调节相关基因的表达，提高甜瓜在低温胁迫下的抗逆性。某些NAC转录因子可能通过激活抗氧化酶基因的表达，增强甜瓜的抗氧化防御能力，减轻低温胁迫引起的氧化损伤；还可能通过调控细胞凋亡相关基因的表达，维持细胞的正常生理功能，避免细胞因低温胁迫而过度凋亡。这些转录因子家族之间并非孤立作用，而是相互协作、相互影响，形成了一个复杂的转录调控网络。不同转录因子家族可能通过共同调控某些关键基因的表达，协同参与甜瓜对低温胁迫的响应过程。AP2/ERF转录因子和bZIP转录因子可能共同调节ABA信号通路相关基因的表达，增强植物对ABA的响应，从而提高耐低温能力。MYB转录因子和NAC转录因子可能在调控苯丙氨酸代谢途径和氧化应激响应相关基因表达时相互作用，共同增强植物对低温胁迫的抵抗能力。这种转录因子之间的协同作用，使得植物能够更有效地整合外界环境信号，精准调控基因表达，以适应低温胁迫环境。转录因子在甜瓜响应低温胁迫过程中具有重要的调控作用。通过对不同转录因子家族的研究，揭示了它们在低温胁迫响应中的调控机制和潜在的协同作用，为深入理解甜瓜耐低温的分子机制提供了重要线索。这些转录因子可作为潜在的分子标记和基因工程靶点，在甜瓜耐低温育种中具有重要的应用价值。未来的研究可以进一步深入探究转录因子之间的相互作用机制，以及它们与其他信号通路之间的关系，为培育耐低温甜瓜新品种提供更坚实