基于转录组学解析NO调控莲雾果实采后絮状绵软进程的分子机制

基于转录组学解析NO调控莲雾果实采后絮状绵软进程的分子机制一、引言1.1研究背景莲雾（Syzygiumsamarangense），又称洋蒲桃、水蒲桃，属于桃金娘科蒲桃属，原产于热带地区，在马来西亚、印度等地广泛种植。我国于20世纪30年代引入莲雾，在海南、广东、广西、福建、台湾等地均有种植。莲雾果实色泽鲜艳，口感清甜多汁，营养丰富，含有多种维生素、矿物质和膳食纤维，深受消费者喜爱。随着种植技术的不断进步和市场需求的增加，莲雾产业发展迅速，种植面积和产量逐年递增。据统计，我国莲雾种植面积已超过[X]万亩，产量达到[X]万吨以上，成为热带水果产业的重要组成部分，在农业经济发展中占据重要地位。然而，莲雾采后极易出现絮状绵软现象，这严重影响了果实的品质和商品价值。絮状绵软的莲雾果实口感变差，失去脆爽的质地，外观也变得不新鲜，难以在市场上销售。同时，絮状绵软还会加速果实的腐烂变质，缩短其货架期，给果农和经销商带来巨大的经济损失。研究表明，采后莲雾果实的絮状绵软与多种因素有关，包括细胞壁代谢、木质素代谢、呼吸作用、能量代谢等。在细胞壁代谢方面，果胶酶、纤维素酶等相关酶的活性变化会导致细胞壁结构的破坏，从而使果实软化；木质素代谢异常则会影响果实的硬度和质地；呼吸作用和能量代谢的改变也会对果实的生理状态产生影响，加速絮状绵软的进程。一氧化氮（NO）作为一种重要的信号分子，在植物生长发育、衰老和胁迫响应等过程中发挥着关键作用。在果实保鲜领域，NO的应用逐渐受到关注。已有研究表明，NO能够延缓多种水果的成熟衰老进程，保持果实的品质和风味。例如，在猕猴桃保鲜中，NO处理可显著延缓果实硬度的下降，降低呼吸速率和失重率，保持较高的可溶性固形物含量和可滴定酸含量，从而延长果实的贮藏期；在草莓保鲜中，NO能够抑制乙烯的产生，延缓果实的软化和腐烂，保持果实的色泽和口感。NO在果实保鲜中的作用机制主要包括抑制乙烯合成、调节呼吸作用、增强抗氧化能力等。在抑制乙烯合成方面，NO可以通过调控乙烯合成关键酶的活性和基因表达，减少乙烯的生成；在调节呼吸作用方面，NO能够抑制细胞色素C氧化酶活性，降低呼吸速率，减少能量消耗；在增强抗氧化能力方面，NO可以提高果实中抗氧化酶的活性，清除自由基，减轻氧化损伤。1.2研究目的与意义本研究旨在从转录组学层面深入探究NO对莲雾采后絮状绵软的影响机制。通过对NO处理和未处理的莲雾果实进行转录组测序和分析，筛选出与絮状绵软进程相关的差异表达基因，揭示NO调控莲雾采后絮状绵软的分子途径。本研究具有重要的理论和实践意义。在理论上，有助于深入了解莲雾采后絮状绵软的分子机制，丰富果实采后生理的研究内容，为进一步探究果实成熟衰老的调控机制提供参考。在实践中，可为开发基于NO的莲雾保鲜技术提供科学依据，通过优化NO处理条件，有效延缓莲雾采后絮状绵软进程，延长果实货架期，减少经济损失，促进莲雾产业的健康发展。1.3国内外研究现状1.3.1莲雾采后生理研究进展莲雾采后生理变化是影响其品质和贮藏寿命的关键因素。研究表明，莲雾采后呼吸作用较为旺盛，属于非呼吸跃变型果实，但其呼吸速率在贮藏过程中会逐渐升高，导致果实的能量消耗加快，加速果实的衰老进程。张福平在《莲雾采后贮藏期间生理变化的研究》中指出，随着贮藏时间的延长，莲雾果实的呼吸强度呈上升趋势，在贮藏后期呼吸强度急剧增加，这与果实的软化和品质下降密切相关。乙烯释放方面，莲雾采后乙烯释放量相对较低，且没有明显的乙烯释放高峰。然而，乙烯仍然对莲雾的成熟和衰老起到一定的调控作用。少量的乙烯释放能够促进果实的成熟，使果实的色泽、风味等品质发生变化；但过多的乙烯则会加速果实的衰老，导致果实软化、腐烂。有研究发现，通过抑制乙烯的作用，可以延缓莲雾果实的成熟衰老进程，保持果实的品质。水分及营养成分变化也是莲雾采后生理研究的重要内容。莲雾果实含水量高，采后在常温下极易失水，导致果实失重、皱缩，影响果实的外观和口感。同时，果实中的营养成分如维生素C、可溶性糖、有机酸等含量也会随着贮藏时间的延长而逐渐下降。其中，维生素C具有抗氧化作用，其含量的下降会降低果实的抗氧化能力，加速果实的衰老；可溶性糖是果实甜味的主要来源，其含量的减少会使果实的甜度降低；有机酸含量的变化则会影响果实的风味和口感。马华青等学者在《莲雾保鲜研究进展》中提到，莲雾采收后水分损失快，在常温下贮藏2-3天，果实的失重率就可达到5%以上，维生素C含量下降10%-20%，可溶性糖含量下降15%-25%，严重影响果实的商品价值。此外，莲雾采后细胞壁代谢和木质素代谢也发生显著变化。在细胞壁代谢方面，果胶酶、纤维素酶等相关酶的活性逐渐增强，导致细胞壁结构被破坏，果实硬度下降，出现软化现象。例如，多聚半乳糖醛酸酶（PG）能够催化果胶物质的降解，使果实的果胶含量减少，从而降低果实的硬度；纤维素酶则可分解纤维素，破坏细胞壁的支撑结构，进一步促进果实的软化。在木质素代谢方面，苯丙氨酸解氨酶（PAL）、过氧化物酶（POD）等酶的活性升高，促使木质素合成增加，导致果实木质化程度提高，口感变差。高维亚等在《NO处理对采后莲雾果实絮状绵软进程的影响》中研究发现，在莲雾果实贮藏过程中，PG、果胶甲酯酶（PME）、β-半乳糖苷酶（β-Gal）等细胞壁降解酶的活性逐渐上升，同时PAL、POD、4-香豆酸辅酶A连接酶（4-CL）等木质素合成相关酶的活性也显著增加，这些变化与果实的絮状绵软进程密切相关。1.3.2NO在果实采后生理调控中的作用研究NO作为一种重要的信号分子，在果实采后生理调控中发挥着多方面的作用。在果实成熟和衰老方面，大量研究表明，适宜浓度的NO处理能够延缓果实的成熟衰老进程。通过抑制乙烯的合成和作用，NO可以降低果实的呼吸速率，减少能量消耗，从而延长果实的贮藏寿命。在对草莓的研究中发现，NO处理能够显著抑制草莓果实乙烯的产生，延缓果实的软化和色泽变化，保持果实的硬度和品质；在香蕉保鲜中，NO处理同样能够抑制乙烯的释放，推迟香蕉果实的成熟，延长其货架期。在果实软化方面，NO主要通过调节细胞壁代谢相关酶的活性来影响果实的软化进程。朱树华等在对桃果实的研究中发现，NO处理能够延缓贮藏过程中细胞壁纤维素、半纤维素的降解，显著降低果胶甲酯酶（PME）和内切葡聚糖酶（EGase）活性，抑制离子型果胶含量下降和可溶性果胶含量增加，从而保持果实的硬度，减轻果实的软化程度。此外，NO还可以通过影响木质素的合成来调控果实的质地，抑制木质素合成相关酶的活性，减少木质素的积累，避免果实过度木质化，保持果实的脆爽口感。在果实抗病方面，NO能够诱导果实产生抗病防卫反应，提高果实的抗病能力。研究表明，NO可以激活果实中的防御相关基因表达，增强果实的抗氧化系统，提高果实对病原菌的抵抗能力。例如，在柑橘果实采后处理中，NO处理能够诱导柑橘果实中病程相关蛋白基因的表达，提高果实的抗病性，减少病害的发生；在苹果果实中，NO处理可以增强果实中抗氧化酶的活性，清除自由基，减轻病原菌侵染引起的氧化损伤，从而提高果实的抗病性。1.3.3转录组学在果实采后研究中的应用转录组学是研究生物体在特定生理状态下全部转录本的学科，在果实采后研究中具有广泛的应用。通过转录组测序技术，可以全面了解果实采后基因表达的变化，挖掘与果实成熟、衰老、品质变化等相关的关键基因和代谢通路。在解析果实采后基因表达变化方面，转录组学为研究提供了有力的工具。对葡萄采后转录组分析发现，在果实贮藏过程中，与呼吸作用、乙烯合成、细胞壁代谢等相关的基因表达发生显著变化。在呼吸作用相关基因中，一些参与糖酵解、三羧酸循环的基因表达上调，表明呼吸作用增强；在乙烯合成相关基因中，ACC合成酶和ACC氧化酶基因的表达变化与乙烯释放量的变化密切相关；在细胞壁代谢相关基因中，果胶酶、纤维素酶等基因的表达上调，导致细胞壁降解，果实软化。在挖掘关键基因和代谢通路方面，转录组学也取得了许多重要成果。对芒果采后转录组研究筛选出了多个与果实成熟衰老相关的关键基因，如编码乙烯响应因子的基因、细胞壁降解酶基因等。进一步分析发现，这些关键基因参与了乙烯信号转导、细胞壁代谢、能量代谢等多个重要代谢通路。通过调控这些关键基因和代谢通路，可以有效地延缓芒果果实的成熟衰老进程，保持果实的品质。此外，转录组学还可以用于研究果实对逆境胁迫的响应机制，为果实保鲜技术的开发提供理论依据。对低温胁迫下草莓果实的转录组分析，揭示了草莓果实响应低温胁迫的基因表达变化和代谢通路，为草莓果实的低温保鲜提供了新的思路和方法。1.4研究内容与技术路线1.4.1研究内容莲雾采后NO释放与软化指标测定：以新鲜采摘的莲雾果实为材料，采用NO特异性探针和分析仪，测定采后不同时间点莲雾果实的NO释放量。同时，测定果实的硬度、质量损失率、可溶性固形物含量、可滴定酸含量等软化相关指标，分析NO释放与莲雾软化程度之间的关系。通过相关性分析，明确NO在莲雾采后软化进程中的作用规律。转录组测序及数据分析：对NO处理和未处理的莲雾果实进行转录组测序，比较两组之间的基因表达差异，筛选出被NO调控的差异表达基因。运用生物信息学方法，对差异表达基因进行功能注释和富集分析，确定其参与的主要生物学过程和代谢通路。重点关注与细胞壁代谢、木质素代谢、呼吸作用、能量代谢等相关的基因，探究NO调控莲雾采后絮状绵软的分子途径。关键基因和转录因子的功能验证：通过基因克隆技术，获得差异表达基因和转录因子的全长cDNA序列。利用实时定量PCR、Westernblot等技术，分析这些基因和转录因子在莲雾果实采后不同时期的表达变化，验证其在NO调控絮状绵软进程中的作用。构建基因过表达载体和RNA干扰载体，转化莲雾果实或相关细胞系，观察基因表达变化对果实软化和相关生理指标的影响，进一步明确关键基因和转录因子的功能。NO调控莲雾采后絮状绵软的分子机制解析：综合上述研究结果，结合已有文献报道，构建NO调控莲雾采后絮状绵软的分子机制模型。阐述NO如何通过调节关键基因和转录因子的表达，影响细胞壁代谢、木质素代谢、呼吸作用、能量代谢等生理过程，从而延缓莲雾采后絮状绵软进程。为基于NO的莲雾保鲜技术开发提供理论依据。1.4.2技术路线本研究的技术路线如图1所示：图1技术路线图实验设计：选取成熟度一致、无病虫害和机械损伤的莲雾果实，随机分为NO处理组和对照组。NO处理组采用适宜浓度的NO气体熏蒸处理，对照组则进行相同时间的空气熏蒸处理。处理后，将果实置于适宜的贮藏条件下，定期取样进行各项指标的测定。样本处理：在采后不同时间点，采集NO处理组和对照组的莲雾果实样本，迅速用液氮冷冻，并保存于-80℃冰箱中备用。用于生理指标测定的样本，需在常温下解冻后进行相关分析；用于转录组测序的样本，需保证RNA的完整性和纯度，通过RNA提取、质量检测等步骤，获得高质量的RNA样品。转录组测序与分析：将提取的RNA样品送测序公司进行转录组测序，获得原始测序数据。对原始数据进行质量控制和过滤，去除低质量reads和接头序列。将过滤后的cleanreads与莲雾参考基因组进行比对，统计基因的表达量。通过差异表达分析，筛选出NO处理组和对照组之间的差异表达基因。利用生物信息学工具，对差异表达基因进行GO功能注释、KEGG通路富集分析等，挖掘NO调控莲雾采后絮状绵软的关键基因和代谢通路。验证分析：采用实时定量PCR技术，对转录组测序结果中的部分差异表达基因进行验证，确保测序数据的可靠性。利用基因克隆、载体构建、遗传转化等技术，对关键基因和转录因子进行功能验证。通过分析基因过表达或沉默对莲雾果实软化、相关生理指标和基因表达的影响，明确其在NO调控絮状绵软进程中的作用机制。最后，综合各项研究结果，总结NO对莲雾果实采后絮状绵软进程的影响及其分子机制。二、材料与方法2.1实验材料实验选用成熟度一致、无病虫害和机械损伤的“黑金刚”莲雾果实作为研究对象，该品种莲雾果实较大，果色深红，口感清甜，在市场上广受欢迎。莲雾果实采自海南省[具体产地]的果园，于清晨露水干后进行采摘，采摘时使用锋利的剪刀，从果柄基部剪下，轻拿轻放，避免果实受到机械损伤。采摘后，立即将果实运回实验室，挑选符合要求的果实随机分为两组，每组[X]个果实，分别进行NO处理和对照处理。实验中用到的主要试剂有：NO供体硝普钠（SNP），购自Sigma公司；RNA提取试剂TRIzol，购自Invitrogen公司；反转录试剂盒和实时定量PCR试剂盒，均购自TaKaRa公司；其他常规试剂均为国产分析纯。实验仪器设备包括：电子天平（精度0.01g），购自梅特勒-托利多仪器有限公司；硬度计（型号FT327），由英国FruitTech公司生产；折光仪（型号WYT-4），购自上海仪电物理光学仪器有限公司；高速冷冻离心机（型号5424R），购自德国Eppendorf公司；荧光定量PCR仪（型号CFX96），为美国Bio-Rad公司产品；高通量测序平台（型号IlluminaHiSeqXTen），由Illumina公司提供。2.2实验方法2.2.1莲雾果实采后NO释放量与软化指标测定使用NO特异性探针DAF-FMDA（4,5-二氨基荧光素二乙酸酯）和NO分析仪（型号：NOA-280i，美国Apollo公司）测定莲雾果实采后NO释放量。将莲雾果实横切成厚度约为5mm的薄片，取1g果肉组织放入含有100μMDAF-FMDA的缓冲液（pH7.4，含有10mMHEPES、137mMNaCl、5mMKCl、0.7mMNa₂HPO₄、0.3mMKH₂PO₄）中，在黑暗条件下孵育30min，使DAF-FMDA进入细胞并被酯酶水解为DAF-FM，DAF-FM与NO反应生成具有荧光的三唑化合物。然后，用缓冲液冲洗果肉组织3次，去除未反应的DAF-FMDA。将处理后的果肉组织放入石英比色皿中，加入适量缓冲液，使用荧光分光光度计（型号：F-4600，日本Hitachi公司）测定荧光强度，激发波长为495nm，发射波长为515nm。根据标准曲线计算NO含量，标准曲线的制作使用不同浓度的SNP（硝普钠，NO供体）溶液，在相同条件下测定荧光强度，以荧光强度为纵坐标，NO浓度为横坐标绘制标准曲线。果实硬度采用硬度计（型号FT327，英国FruitTech公司）测定，每个果实选取赤道部位相对的两个点进行测定，取平均值作为该果实的硬度，单位为N/cm²。质量损失率计算公式为：质量损失率（%）=（初始质量-测定时质量）/初始质量×100%，使用精度为0.01g的电子天平（梅特勒-托利多仪器有限公司）称量果实质量。可溶性固形物含量利用折光仪（型号WYT-4，上海仪电物理光学仪器有限公司）测定，将果实榨汁后，取1-2滴汁液滴在折光仪的棱镜上，读取可溶性固形物含量，单位为°Bx。可滴定酸含量采用酸碱滴定法测定，准确称取10g果肉，加入50mL蒸馏水，匀浆后过滤，取25mL滤液，用0.1MNaOH标准溶液滴定，以酚酞为指示剂，滴定至溶液呈微红色且30s内不褪色，记录消耗NaOH标准溶液的体积，可滴定酸含量以苹果酸计，计算公式为：可滴定酸含量（%）=（消耗NaOH标准溶液体积×NaOH标准溶液浓度×0.067×稀释倍数）/样品质量×100%。果胶酶活性测定采用3,5-二硝基水杨酸（DNS）比色法。取1g果肉组织，加入5mL预冷的提取缓冲液（50mMTris-HCl，pH7.5，含有1mMEDTA、1mMDTT、1%PVP），冰浴匀浆后，在4℃下以12000rpm离心20min，取上清液作为酶粗提液。取1mL酶粗提液，加入1mL1%果胶溶液（用50mMTris-HCl，pH7.5配制），在37℃下反应30min，然后加入2mLDNS试剂，沸水浴5min，冷却后用蒸馏水定容至25mL，在540nm波长下测定吸光值。以相同处理但不加酶粗提液的反应体系作为空白对照，根据标准曲线计算果胶酶活性，标准曲线的制作使用不同浓度的半乳糖醛酸溶液，在相同条件下测定吸光值，以吸光值为纵坐标，半乳糖醛酸浓度为横坐标绘制标准曲线，果胶酶活性单位为U/gFW（每分钟催化产生1μmol半乳糖醛酸所需的酶量为1个酶活力单位）。2.2.2转录组测序及数据分析分别选取NO处理和未处理的莲雾果实，在采后0d、3d、6d、9d、12d时，从每个处理组中随机选取3个果实，取果肉组织混合作为一个生物学重复，每个时间点设置3个生物学重复，共30个样本。使用RNA提取试剂TRIzol（Invitrogen公司）提取样本总RNA，具体步骤如下：将100mg果肉组织放入液氮预冷的研钵中，迅速研磨成粉末，加入1mLTRIzol试剂，充分混匀，室温静置5min；加入0.2mL氯仿，剧烈振荡15s，室温静置3min；在4℃下以12000rpm离心15min，取上清液转移至新的离心管中；加入0.5mL异丙醇，轻轻混匀，室温静置10min；在4℃下以12000rpm离心10min，弃上清液，沉淀用75%乙醇（用DEPC水配制）洗涤2次，每次在4℃下以7500rpm离心5min；弃上清液，室温晾干沉淀，加入适量DEPC水溶解RNA。使用Nanodrop2000分光光度计（ThermoFisherScientific公司）检测RNA的纯度和浓度，要求OD₂₆₀/OD₂₈₀在1.8-2.0之间，OD₂₆₀/OD₂₃₀大于2.0，RNA浓度不低于50ng/μL；使用Agilent2100生物分析仪（AgilentTechnologies公司）检测RNA的完整性，要求RIN值大于7.0。将合格的RNA样本送测序公司进行文库构建和测序。文库构建采用NEBNextUltraRNALibraryPrepKitforIllumina（NewEnglandBiolabs公司）试剂盒，具体步骤如下：使用磁珠富集真核生物mRNA，要求RNA的完整性较高，一般RIN值要大于8；将mRNA进行随机打断；以mRNA为模板，合成第一条cDNA链和第二条cDNA链；进行末端修复、加A尾并连接测序接头，然后进行片段大小选择，选择长度约为300-500bp的片段；最后通过PCR富集得到cDNA文库。文库构建完成后，使用Qubit2.0荧光定量仪（ThermoFisherScientific公司）进行初步定量，使用Agilent2100生物分析仪检测文库的插入片段大小，要求insertsize符合预期；使用Q-PCR方法对文库的有效浓度进行准确定量，要求文库有效浓度＞2nM。将文库种到芯片上进行桥式PCR扩增，然后使用IlluminaHiSeqXTen高通量测序平台进行测序，测序策略为PE150，即双端测序，每个read长度为150bp。对测序得到的原始数据（RawData）进行质量控制和过滤，使用FastQC软件（http://www.bioinformatics.babraham.ac.uk/projects/fastqc/）评估原始数据质量，去除低质量reads（质量值Q＜20的碱基数占比超过20%的reads）、含有接头序列的reads和N含量超过10%的reads，得到高质量的cleanreads。将cleanreads与莲雾参考基因组（如果没有参考基因组，可使用Trinity软件进行denovo组装得到转录组文库）进行比对，使用Hisat2软件（https://daehwankimlab.github.io/hisat2/）进行比对分析，统计比对到基因组上的reads数和比对率。利用StringTie软件（/software/stringtie/）计算基因的表达量，以FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）值表示基因的表达水平，计算公式为：FPKM=10⁶×比对到某基因的reads数/（测序深度×该基因外显子长度）。通过DESeq2软件（/packages/release/bioc/html/DESeq2.html）进行差异表达分析，筛选NO处理组和对照组之间的差异表达基因，设定筛选条件为|log₂(FoldChange)|≥1且FDR（FalseDiscoveryRate）＜0.05。对差异表达基因进行GO（GeneOntology）功能注释和KEGG（KyotoEncyclopediaofGenesandGenomes）通路富集分析，使用DAVID数据库（/）和KOBAS软件（/），GO功能注释包括生物过程（BiologicalProcess）、细胞组分（CellularComponent）和分子功能（MolecularFunction）三个方面，KEGG通路富集分析可确定差异表达基因参与的主要代谢通路和信号转导途径。2.2.3关键基因与转录因子的验证根据转录组测序和数据分析结果，筛选出与莲雾采后絮状绵软进程密切相关的关键基因和转录因子。采用基因克隆技术获得这些基因和转录因子的全长cDNA序列，以莲雾果实总RNA为模板，使用反转录试剂盒（TaKaRa公司）反转录合成cDNA第一链。根据目标基因的cDNA序列设计特异性引物，引物设计使用PrimerPremier5.0软件，引物两端分别引入合适的限制性内切酶位点，以便后续的克隆和表达载体构建。以cDNA第一链为模板，使用高保真DNA聚合酶（TaKaRa公司）进行PCR扩增，PCR反应体系为25μL，包括10×PCR缓冲液2.5μL、dNTPs（2.5mMeach）2μL、上下游引物（10μMeach）各1μL、cDNA模板1μL、高保真DNA聚合酶0.5μL、ddH₂O17μL。PCR反应条件为：95℃预变性5min；95℃变性30s，55-65℃退火30s（根据引物Tm值调整退火温度），72℃延伸1-2min（根据基因长度调整延伸时间），共30-35个循环；72℃终延伸10min。PCR产物经1%琼脂糖凝胶电泳检测后，切胶回收目的片段，使用DNA凝胶回收试剂盒（TaKaRa公司）进行回收纯化。将回收的目的片段与克隆载体pMD19-T（TaKaRa公司）连接，连接体系为10μL，包括pMD19-T载体1μL、目的片段4μL、SolutionI5μL，16℃连接过夜。将连接产物转化大肠杆菌DH5α感受态细胞，涂布于含有氨苄青霉素（Amp）的LB固体培养基上，37℃培养12-16h。挑取单菌落进行PCR鉴定和测序验证，测序由专业测序公司完成。利用实时定量PCR（qRT-PCR）技术验证关键基因和转录因子在莲雾果实采后不同时期的表达变化，以β-actin基因作为内参基因。使用反转录试剂盒将莲雾果实总RNA反转录合成cDNA第一链，具体步骤同基因克隆中的反转录步骤。根据目标基因和内参基因的cDNA序列设计qRT-PCR引物，引物设计原则为：引物长度18-25bp，Tm值在58-62℃之间，GC含量在40%-60%之间，引物之间避免形成二聚体和发夹结构。qRT-PCR反应体系为20μL，包括SYBRPremixExTaqII（2×）10μL、上下游引物（10μMeach）各0.8μL、cDNA模板1μL、ROXReferenceDyeII（50×）0.4μL、ddH₂O7μL。qRT-PCR反应条件为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；熔解曲线分析从65℃以0.5℃/s的速度升温至95℃。每个样本设置3个技术重复，使用CFX96荧光定量PCR仪（美国Bio-Rad公司）进行扩增和数据采集，采用2⁻ΔΔCt法计算基因的相对表达量。采用Westernblot技术检测关键蛋白的表达水平，以验证转录水平的结果。提取莲雾果实不同时期的总蛋白，使用蛋白提取试剂盒（上海碧云天生物技术有限公司），按照说明书操作。使用BCA蛋白定量试剂盒（上海碧云天生物技术有限公司）测定蛋白浓度，以BSA（牛血清白蛋白）为标准品绘制标准曲线，根据标准曲线计算蛋白浓度。取30-50μg蛋白样品，加入适量5×SDS上样缓冲液，煮沸5min使蛋白变性。将变性后的蛋白样品进行SDS-PAGE（十二烷基硫酸钠-聚丙烯酰胺凝胶电泳），电泳条件为：浓缩胶80V，30min；分离胶120V，60-90min（根据蛋白分子量调整电泳时间）。电泳结束后，将蛋白转移至PVDF膜（聚偏二氟乙烯膜，美国Millipore公司）上，转膜条件为：250mA，90min。转膜完成后，将PVDF膜放入5%脱脂奶粉（用TBST缓冲液配制）中，室温封闭1-2h，以封闭非特异性结合位点。用TBST缓冲液洗涤PVDF膜3次，每次10min。加入稀释好的一抗（针对目标蛋白的特异性抗体，购自Abcam公司或CST公司），4℃孵育过夜。用TBST缓冲液洗涤PVDF膜3次，每次10min。加入稀释好的二抗（辣根过氧化物酶标记的羊抗兔或羊抗鼠IgG抗体，购自JacksonImmunoResearch公司），室温孵育1-2h。用TBST缓冲液洗涤PVDF膜3次，每次10min。使用ECL化学发光试剂（上海碧云天生物技术有限公司）进行显色，在化学发光成像系统（美国Bio-Rad公司）上曝光成像，分析蛋白表达水平。三、NO释放与莲雾果实采后絮状绵软进程的关联分析3.1莲雾果实采后NO释放规律通过NO特异性探针DAF-FMDA和NO分析仪对采后不同时间的莲雾果实进行检测，得到莲雾果实采后NO释放量变化曲线，结果如图2所示。图2莲雾果实采后NO释放量变化曲线由图2可知，莲雾果实在采后0-3天，NO释放量相对较低且变化较为平稳，维持在[X1]-[X2]nmol/gFW之间。这可能是因为在采后初期，果实处于相对生理稳定状态，NO的合成和代谢较为缓慢。从第3天开始，NO释放量呈现逐渐上升的趋势，在采后6-9天，NO释放量增长较为明显，达到[X3]-[X4]nmol/gFW。这一阶段可能是果实采后生理代谢发生变化，相关信号通路被激活，促使NO合成增加。在采后9-12天，NO释放量继续上升，但上升幅度有所减缓，在第12天达到最高值[X5]nmol/gFW。此后，随着果实的进一步衰老和品质劣变，NO释放量可能会逐渐下降。总体而言，莲雾果实采后NO释放量呈现先缓慢上升，然后快速上升，最后上升减缓并达到峰值的变化趋势，这表明NO的释放与莲雾果实采后的生理进程密切相关。三、NO释放与莲雾果实采后絮状绵软进程的关联分析3.2莲雾果实采后软化指标变化3.2.1硬度变化果实硬度是衡量莲雾果实品质和软化程度的重要指标之一。对采后不同时间的莲雾果实硬度进行测定，结果如图3所示。图3莲雾果实采后硬度变化曲线从图3可以看出，采后莲雾果实硬度随着贮藏时间的延长呈现逐渐下降的趋势。在采后0-3天，果实硬度下降较为缓慢，从初始硬度[X6]N/cm²下降到[X7]N/cm²，下降幅度约为[X8]%。这一阶段果实细胞壁结构相对完整，维持着较高的硬度。3-9天，果实硬度下降速度明显加快，到第9天，硬度降至[X9]N/cm²，下降幅度达到[X10]%。此时期细胞壁降解酶活性增强，细胞壁成分如纤维素、果胶等被逐渐分解，导致果实硬度快速降低。9-12天，果实硬度继续下降，但下降趋势变缓，最终在第12天降至[X11]N/cm²，下降幅度为[X12]%。此时果实细胞壁结构严重受损，果实接近絮状绵软状态。3.2.2失重率变化失重率反映了莲雾果实在贮藏过程中的水分散失情况，对果实的品质和货架期有重要影响。莲雾果实采后失重率随贮藏时间的变化如图4所示。图4莲雾果实采后失重率变化曲线由图4可知，采后莲雾果实的失重率随贮藏时间的延长而逐渐增加。在采后0-3天，失重率增加较为缓慢，从初始的[X13]%上升到[X14]%，这是因为在贮藏初期，果实表皮的角质层和蜡质层能够在一定程度上阻碍水分的散失，水分蒸发速率相对较低。3-9天，失重率上升速度加快，到第9天，失重率达到[X15]%，这一阶段果实呼吸作用增强，消耗了大量的水分和营养物质，同时随着贮藏时间的延长，果实表皮的保护结构逐渐受损，水分散失加剧。9-12天，失重率继续上升，在第12天达到[X16]%。果实失重不仅会导致果实外观皱缩，影响商品价值，还会引起果实内部生理生化变化，加速果实的衰老和软化进程。3.2.3果胶酶等相关酶活性变化果胶酶和纤维素酶等细胞壁降解酶在莲雾果实采后软化过程中发挥着关键作用。对采后莲雾果实中果胶酶、纤维素酶活性进行测定，结果如图5和图6所示。图5莲雾果实采后果胶酶活性变化曲线图6莲雾果实采后纤维素酶活性变化曲线从图5可以看出，采后莲雾果实在贮藏过程中，果胶酶活性呈逐渐上升的趋势。在采后0-3天，果胶酶活性较低且变化不大，从初始活性[X17]U/gFW缓慢上升到[X18]U/gFW。3-9天，果胶酶活性显著升高，到第9天，活性达到[X19]U/gFW，增长幅度较大。果胶酶能够催化果胶物质的降解，使细胞壁中的果胶含量减少，导致细胞间的黏着力下降，从而促进果实的软化。9-12天，果胶酶活性继续上升，但上升速度有所减缓，在第12天达到最高值[X20]U/gFW。图6显示，采后莲雾果实纤维素酶活性在贮藏期间同样呈现上升趋势。采后0-3天，纤维素酶活性较低，从[X21]U/gFW升高到[X22]U/gFW。3-9天，纤维素酶活性快速上升，到第9天，活性增加到[X23]U/gFW。纤维素酶可以分解细胞壁中的纤维素，破坏细胞壁的支撑结构，使得果实硬度降低，加速果实的软化进程。9-12天，纤维素酶活性继续上升，在第12天达到[X24]U/gFW。果胶酶和纤维素酶等细胞壁降解酶活性的升高，共同作用于果实细胞壁结构，导致细胞壁逐渐降解，最终引发莲雾果实的絮状绵软。3.3NO释放与莲雾果实软化的相关性为进一步明确NO释放量与莲雾果实软化各指标之间的关系，对NO释放量、果实硬度、失重率、果胶酶活性和纤维素酶活性进行相关性分析，结果如表1所示。表1NO释放量与莲雾果实软化指标的相关性分析指标NO释放量果实硬度失重率果胶酶活性纤维素酶活性NO释放量1果实硬度-0.923**1失重率0.876**-0.901**1果胶酶活性0.905**-0.936**0.892**1纤维素酶活性0.889**-0.915**0.881**0.945**1注：**表示在0.01水平上显著相关从表1可以看出，NO释放量与果实硬度呈极显著负相关（r=-0.923，P<0.01），即随着NO释放量的增加，果实硬度逐渐降低，表明NO可能通过某种机制促进了果实的软化。NO释放量与失重率呈极显著正相关（r=0.876，P<0.01），说明NO释放量的增加伴随着果实失重率的上升，这可能是由于NO影响了果实的水分代谢，加速了水分的散失，进而导致果实失重增加，间接影响果实的软化进程。此外，NO释放量与果胶酶活性和纤维素酶活性均呈极显著正相关，相关系数分别为0.905和0.889（P<0.01）。这表明NO可能通过调控果胶酶和纤维素酶等细胞壁降解酶的活性，促进细胞壁成分的降解，从而加速莲雾果实的软化。果胶酶和纤维素酶活性的升高，使得细胞壁中的果胶和纤维素被分解，细胞间的黏着力下降，细胞壁支撑结构受损，最终导致果实硬度降低，呈现出絮状绵软的状态。综上所述，NO释放与莲雾果实的软化密切相关，NO可能通过调节果实水分代谢和细胞壁降解酶活性等途径，参与莲雾果实采后的絮状绵软进程。四、转录组学分析NO调控莲雾果实采后絮状绵软的分子机制4.1转录组测序数据质量评估对NO处理组和对照组莲雾果实的30个样本进行转录组测序，共获得[X]Gb的原始数据（RawData）。使用FastQC软件对原始数据进行质量评估，结果显示，所有样本的碱基质量分数（QualityScore）平均值均在30以上，表明测序错误率较低，碱基测序质量较高。在碱基测序质量图中，大部分reads的碱基质量值在20以上，且reads末端的碱基质量值虽略有下降，但未出现大面积碱基质量低于20的情况，符合测序数据质量要求。碱基分布情况较为均匀，A和T、C与G的比例相近，未出现异常波动。在去除低质量reads（质量值Q＜20的碱基数占比超过20%的reads）、含有接头序列的reads和N含量超过10%的reads后，得到高质量的cleanreads，各样本的cleanreads数量在[X1]-[X2]万之间，cleanreads占原始reads的比例均在95%以上，表明数据过滤效果良好，有效数据量充足。将cleanreads与莲雾参考基因组进行比对，使用Hisat2软件进行比对分析。结果显示，各样本的比对率在85%-92%之间，唯一比对率在78%-85%之间。较高的比对率和唯一比对率说明测序数据与参考基因组的匹配度较高，能够准确地定位到基因组上，为后续的基因表达分析提供了可靠的数据基础。此外，对基因表达量分布进行分析，结果表明，基因表达量在样本中呈现出典型的分布特征，大部分基因表达量较低，少数基因高表达，符合生物学规律。同时，通过对已知基因覆盖度的评估，发现已知基因覆盖度在80%以上，说明测序数据能够较好地覆盖莲雾基因组中的已知基因，保证了基因表达分析的全面性和准确性。综上所述，本研究的转录组测序数据质量良好，可靠性高，可用于后续的差异表达分析和功能注释等研究。四、转录组学分析NO调控莲雾果实采后絮状绵软的分子机制4.2差异表达基因筛选与分析4.2.1差异表达基因数量统计通过DESeq2软件对NO处理组和对照组莲雾果实的转录组数据进行差异表达分析，以|log₂(FoldChange)|≥1且FDR＜0.05作为筛选条件，共筛选出[X]个差异表达基因。其中，上调表达的基因有[X1]个，下调表达的基因有[X2]个。在采后不同时间点，差异表达基因的数量存在一定变化。在采后3天，差异表达基因数量相对较少，共筛选出[X3]个，其中上调基因[X4]个，下调基因[X5]个；随着贮藏时间的延长，到采后9天，差异表达基因数量显著增加，达到[X6]个，上调基因[X7]个，下调基因[X8]个；在采后12天，差异表达基因数量为[X9]个，上调基因[X10]个，下调基因[X11]个。为直观展示差异表达基因的数量变化情况，绘制了差异表达基因数量柱状图，如图7所示。从图中可以清晰地看出，随着采后时间的推移，NO处理组和对照组之间的差异表达基因数量逐渐增多，表明NO对莲雾果实基因表达的影响在采后逐渐增强，涉及到更多的基因参与到果实的生理调控过程中。图7差异表达基因数量柱状图进一步对不同时间点上调和下调基因的比例进行分析，结果如图8所示。在采后3天，上调基因占差异表达基因总数的比例为[X12]%，下调基因占比为[X13]%；到采后9天，上调基因占比增加到[X14]%，下调基因占比为[X15]%；采后12天，上调基因占比为[X16]%，下调基因占比为[X17]%。这表明在采后不同阶段，NO对莲雾果实基因表达的调控具有动态变化，且上调基因的比例随着贮藏时间的延长有逐渐增加的趋势，说明NO可能通过上调某些基因的表达来调控莲雾果实采后的生理进程，从而影响絮状绵软的发生。图8不同时间点上调和下调基因比例图4.2.2差异表达基因的功能注释利用DAVID数据库和KOBAS软件对筛选出的差异表达基因进行GO功能注释和KEGG通路富集分析，以明确这些基因参与的生物过程、分子功能和代谢通路。在GO功能注释中，从生物过程（BP）、细胞组分（CC）和分子功能（MF）三个方面对差异表达基因进行分类。在生物过程方面，差异表达基因主要富集在“细胞过程”“代谢过程”“应激反应”“信号转导”等类别。其中，参与“细胞过程”的基因数量较多，如与细胞壁代谢相关的基因，参与细胞壁的合成、修饰和降解过程，这与莲雾果实采后细胞壁结构的变化以及絮状绵软的发生密切相关；“代谢过程”中，涉及碳水化合物代谢、能量代谢、脂质代谢等多个方面的基因被显著富集，表明NO处理可能影响了莲雾果实的多种代谢途径，进而影响果实的生理状态和品质变化。在细胞组分方面，差异表达基因主要分布在“细胞”“细胞部分”“细胞器”“细胞膜”等类别。例如，在“细胞器”类别中，与线粒体、叶绿体等细胞器相关的基因表达发生变化，线粒体是细胞能量代谢的中心，其相关基因的改变可能影响细胞的能量供应，而叶绿体与光合作用密切相关，其基因表达的变化可能间接影响果实的营养物质合成和积累，这些都可能对莲雾果实采后的絮状绵软进程产生影响。在分子功能方面，差异表达基因主要富集在“催化活性”“结合活性”“转运活性”等类别。具有“催化活性”的基因，如各种酶类基因，在细胞壁代谢、呼吸作用、能量代谢等生理过程中发挥着关键的催化作用；“结合活性”相关基因，如与蛋白质、核酸、小分子等结合的基因，参与了信号传导、物质运输等过程，对维持细胞正常生理功能至关重要。在KEGG通路富集分析中，差异表达基因显著富集在多个代谢通路中，包括“植物激素信号转导”“淀粉和蔗糖代谢”“半乳糖代谢”“苯丙烷生物合成”“内质网中的蛋白质加工”等通路。在“植物激素信号转导”通路中，与乙烯、生长素、脱落酸等激素信号转导相关的基因表达发生显著变化。乙烯是促进果实成熟和衰老的重要激素，NO可能通过调节乙烯信号转导途径中的关键基因表达，影响乙烯的合成和作用，从而延缓莲雾果实的成熟衰老进程，抑制絮状绵软的发生；在“淀粉和蔗糖代谢”通路中，参与淀粉合成和降解、蔗糖代谢的基因表达改变，淀粉和蔗糖是果实中重要的碳水化合物，其代谢变化会影响果实的甜度、能量供应等，进而影响果实的品质和生理状态。“苯丙烷生物合成”通路的富集与莲雾果实的木质素代谢密切相关。木质素是细胞壁的重要组成成分，其合成和积累会影响果实的硬度和质地。差异表达基因在该通路中的富集，表明NO处理可能通过调控苯丙烷生物合成途径中关键基因的表达，影响木质素的合成，从而改变果实细胞壁的结构和性质，对莲雾果实的絮状绵软进程产生影响。这些GO功能注释和KEGG通路富集分析结果，为深入理解NO调控莲雾果实采后絮状绵软的分子机制提供了重要线索。4.3NO调控莲雾果实采后絮状绵软的关键基因与代谢通路4.3.1细胞壁代谢相关基因在细胞壁代谢相关基因中，果胶代谢相关基因的表达变化对莲雾果实的软化具有重要影响。多聚半乳糖醛酸酶（PG）基因在NO处理组和对照组中呈现出显著不同的表达模式。在对照组中，随着贮藏时间的延长，PG基因表达量逐渐上调，在采后9-12天表达量急剧增加；而在NO处理组中，PG基因表达量的上调受到明显抑制，在整个贮藏期间表达量均显著低于对照组。PG是催化果胶降解的关键酶，其基因表达上调会导致果胶物质的分解加速，使细胞壁中果胶含量减少，细胞间的黏着力下降，从而促进果实的软化。NO处理抑制了PG基因的表达，可能通过减少果胶的降解，维持细胞壁的完整性和细胞间的黏着力，进而延缓莲雾果实的絮状绵软进程。果胶甲酯酶（PME）基因的表达也在NO处理下发生变化。对照组中，PME基因表达量在采后逐渐上升，在采后6-9天上升幅度较大；NO处理组中，PME基因表达量上升趋势相对平缓，且在各时间点均低于对照组。PME能够催化果胶甲酯的水解，使果胶分子去甲酯化，改变果胶的理化性质，影响细胞壁的结构和功能。NO对PME基因表达的调控，可能影响果胶的甲酯化程度，从而调节细胞壁的稳定性，对果实的软化进程产生影响。纤维素代谢相关基因在NO调控莲雾果实采后絮状绵软中也发挥着重要作用。内切葡聚糖酶（EGase）基因在对照组中表达量随贮藏时间增加而升高，尤其在采后9-12天升高明显；而NO处理组中，EGase基因表达量的升高被显著抑制。EGase能够水解纤维素分子内部的β-1,4-糖苷键，破坏纤维素的结构，导致细胞壁支撑能力下降，促进果实软化。NO抑制EGase基因的表达，有助于维持纤维素的结构和细胞壁的支撑功能，延缓果实的软化。β-半乳糖苷酶（β-Gal）基因的表达变化与莲雾果实软化密切相关。在对照组中，β-Gal基因表达量逐渐升高，在采后后期升高幅度较大；NO处理组中，β-Gal基因表达量的升高受到抑制，在采后不同时间点均低于对照组。β-Gal参与细胞壁中半乳糖苷的水解，其活性增加会导致细胞壁结构的改变，促进果实软化。NO对β-Gal基因表达的调控，可能通过影响细胞壁中半乳糖苷的代谢，维持细胞壁的正常结构，延缓莲雾果实的絮状绵软。综上所述，NO通过调控果胶、纤维素代谢相关基因的表达，如PG、PME、EGase、β-Gal等基因，影响细胞壁的结构和组成，维持细胞壁的完整性和稳定性，从而延缓莲雾果实采后絮状绵软进程，保持果实的品质。4.3.2能量代谢相关基因能量代谢在莲雾果实采后生理过程中起着关键作用，与果实的衰老和絮状绵软密切相关。在能量代谢相关基因中，ATP合成相关基因的表达变化对果实能量供应和生理状态有重要影响。ATP合酶基因在NO处理组和对照组中的表达存在显著差异。在对照组中，随着贮藏时间的延长，ATP合酶基因表达量逐渐下降，在采后9-12天下降明显；而在NO处理组中，ATP合酶基因表达量的下降受到一定程度的抑制，在采后各时间点表达量均高于对照组。ATP合酶是催化ATP合成的关键酶，其基因表达量的降低会导致ATP合成减少，细胞能量供应不足。NO处理维持了ATP合酶基因的相对较高表达，有助于保证细胞内ATP的合成，为果实的生理活动提供充足的能量，延缓果实的衰老和絮状绵软进程。呼吸作用相关基因的表达也受到NO的调控。细胞色素C氧化酶（COX）基因在对照组中表达量逐渐升高，在采后6-9天升高幅度较大；而在NO处理组中，COX基因表达量的升高相对平缓，且在各时间点均低于对照组。COX是呼吸链末端的氧化酶，其活性影响呼吸作用的强度。COX基因表达量升高会导致呼吸作用增强，加速能量消耗和果实的衰老进程。NO抑制COX基因的表达，可能通过降低呼吸作用强度，减少能量消耗，从而延缓莲雾果实的衰老和絮状绵软。此外，参与糖酵解途径的关键酶基因，如己糖激酶（HK）、磷酸果糖激酶（PFK）等基因的表达也发生变化。在对照组中，这些基因的表达量在采后逐渐上升，表明糖酵解途径增强，以满足果实生理活动对能量的需求；而在NO处理组中，HK、PFK等基因表达量的上升幅度相对较小。糖酵解是细胞呼吸的第一阶段，其代谢强度的变化会影响细胞的能量供应和代谢状态。NO对糖酵解相关基因表达的调控，可能通过调节糖酵解途径的强度，影响能量产生的速率，进而对莲雾果实的采后生理进程产生影响。能量代谢相关基因在NO调控莲雾果实采后絮状绵软中发挥着重要作用。NO通过调节ATP合成、呼吸作用以及糖酵解等相关基因的表达，维持细胞内的能量平衡，减少能量的过度消耗，从而延缓果实的衰老和絮状绵软进程，保持果实的品质和货架期。4.3.3激素信号转导相关基因激素信号转导在果实的生长、发育、成熟和衰老过程中起着至关重要的调节作用，NO对莲雾果实采后絮状绵软的调控与激素信号转导密切相关。在乙烯信号转导通路中，多个关键基因的表达受到NO的显著影响。乙烯响应因子（ERF）基因家族在果实成熟和衰老过程中发挥重要作用，其中ERF1、ERF2等基因在对照组中随着贮藏时间的延长表达量逐渐上调，尤其在采后9-12天表达量急剧增加；而在NO处理组中，这些基因的表达量上调受到明显抑制，在整个贮藏期间表达量均显著低于对照组。ERF蛋白能够与乙烯响应元件结合，调控下游与果实成熟、衰老相关基因的表达。NO抑制ERF基因的表达，可能通过阻断乙烯信号转导途径，减少乙烯对果实成熟和衰老相关基因的诱导作用，从而延缓莲雾果实的絮状绵软进程。乙烯合成关键酶基因1-氨基环丙烷-1-羧酸合成酶（ACS）和1-氨基环丙烷-1-羧酸氧化酶（ACO）的表达也受到NO的调控。在对照组中，ACS和ACO基因表达量逐渐上升，在采后6-9天上升幅度较大，导致乙烯合成增加；而在NO处理组中，ACS和ACO基因表达量的上升趋势相对平缓，乙烯合成受到抑制。ACS和ACO分别催化乙烯合成的前两步反应，是乙烯合成的关键限速酶。NO对ACS和ACO基因表达的抑制，减少了乙烯的合成，从而降低了乙烯对果实成熟和衰老的促进作用，有助于延缓莲雾果实的絮状绵软。生长素信号转导通路中，生长素响应因子（ARF）基因的表达变化也与NO调控莲雾果实采后絮状绵软有关。在对照组中，ARF基因表达量在采后呈现一定的变化趋势，在某些时间点表达量升高；而在NO处理组中，ARF基因表达量的变化相对平稳，在部分时间点表达量低于对照组。ARF蛋白能够与生长素响应元件结合，调控下游基因的表达，参与果实的生长、发育和衰老过程。NO对ARF基因表达的调节，可能通过影响生长素信号转导，改变果实细胞的生理活动，进而影响莲雾果实的采后生理进程和絮状绵软的发生。激素信号转导相关基因在NO调控莲雾果实采后絮状绵软中发挥着重要作用。NO通过调节乙烯、生长素等激素信号转导通路中关键基因的表达，如ERF、ACS、ACO、ARF等基因，影响激素的合成和信号传递，从而调控果实的成熟和衰老进程，延缓莲雾果实采后絮状绵软的发生，保持果实的品质和商品价值。4.4转录因子在NO调控莲雾果实采后絮状绵软中的作用4.4.1差异表达转录因子的鉴定转录因子在植物生长发育和逆境响应等过程中发挥着关键的调控作用，在NO调控莲雾果实采后絮状绵软进程中也可能扮演重要角色。通过对转录组数据的深入分析，利用相关软件和数据库（如PlantTFDB等）对差异表达基因进行筛选，共鉴定出[X]个差异表达的转录因子，这些转录因子分属于多个不同的家族。对差异表达转录因子的家族分布进行统计分析，结果发现，AP2/ERF、MYB、bHLH、NAC等家族的转录因子数量较多，在差异表达转录因子中占比较大。其中，AP2/ERF家族转录因子有[X1]个，占比[X2]%；MYB家族转录因子有[X3]个，占比[X4]%；bHLH家族转录因子有[X5]个，占比[X6]%；NAC家族转录因子有[X7]个，占比[X8]%。AP2/ERF家族转录因子在植物激素信号转导、逆境响应和果实发育等过程中具有重要作用。在本研究中，AP2/ERF家族转录因子的差异表达可能与NO调控莲雾果实采后乙烯信号转导以及果实的成熟衰老进程密切相关。MYB家族转录因子参与调控植物的次生代谢、细胞分化和发育等过程，其在NO处理后的莲雾果实中差异表达，可能影响果实细胞壁成分的合成和代谢，进而对果实的硬度和絮状绵软进程产生影响。bHLH家族转录因子在植物生长发育、激素信号传导和环境胁迫响应中发挥重要功能。在NO调控莲雾果实采后絮状绵软过程中，bHLH家族转录因子可能通过调节相关基因的表达，参与果实的能量代谢、抗氧化防御等生理过程，维持果实的正常生理状态，延缓絮状绵软的发生。NAC家族转录因子在植物衰老、逆境响应和器官发育等方面具有重要调控作用。在本研究中，NAC家族转录因子的差异表达可能与NO诱导的莲雾果实采后衰老调控机制有关，通过调控下游基因的表达，影响果实的细胞壁代谢、激素信号转导等过程，从而影响果实的絮状绵软进程。这些差异表达的转录因子在不同家族中的分布情况，为进一步研究NO调控莲雾果实采后絮状绵软的分子机制提供了重要线索，后续将对关键转录因子进行深入研究，解析其在调控网络中的作用。4.4.2关键转录因子的调控网络分析为深入解析转录因子在NO调控莲雾果实采后絮状绵软中的调控机制，选取AP2/ERF、MYB、bHLH、NAC等家族中差异表达显著且可能与果实絮状绵软密切相关的关键转录因子，利用生物信息学方法和相关数据库（如STRING、PlantRegMap等），构建关键转录因子与靶基因的调控网络。通过对调控网络的分析，发现关键转录因子与多个靶基因之间存在复杂的相互作用关系。例如，在AP2/ERF家族中，转录因子ERF1与乙烯合成关键酶基因ACS和ACO以及乙烯响应基因ERF2、ERF3等存在调控关系。在NO处理下，ERF1表达受到抑制，可能通过下调ACS和ACO基因的表达，减少乙烯的合成，同时抑制ERF2、ERF3等乙烯响应基因的表达，阻断乙烯信号转导途径，从而延缓莲雾果实的成熟衰老和絮状绵软进程。在MYB家族中，关键转录因子MYB1与细胞壁代谢相关基因PG、PME、EGase等存在调控关系。NO处理可能影响MYB1的表达，进而调控PG、PME、EGase等基因的表达，影响细胞壁成分的降解和合成，维持细胞壁的完整性和稳定性，对莲雾果实的硬度和絮状绵软进程产生影响。bHLH家族中的关键转录因子bHLH1与能量代谢相关基因ATP合酶基因、COX基因以及抗氧化酶基因SOD、POD等存在调控关系。在NO作用下，bHLH1表达变化可能通过调节ATP合酶基因和COX基因的表达，维持细胞内的能量平衡，同时上调SOD、POD等抗氧化酶基因的表达，增强果实的抗氧化能力，减少活性氧积累对细胞的损伤，延缓莲雾果实的衰老和絮状绵软。NAC家族中的关键转录因子NAC1与激素信号转导相关基因ARF、乙烯信号通路基因以及细胞壁代谢相关基因存在调控关系。NO处理可能通过影响NAC1的表达，调节ARF基因的表达，影响生长素信号转导，同时调控乙烯信号通路基因和细胞壁代谢相关基因的表达，综合影响莲雾果实的采后生理进程和絮状绵软的发生。这些关键转录因子通过与靶基因的相互作用，形成复杂的调控网络，在NO调控莲雾果实采后絮状绵软进程中发挥重要作用。它们从乙烯信号转导、细胞壁代谢、能量代谢、激素信号转导等多个方面对果实的生理过程进行调控，维持果实的品质和延缓絮状绵软的发生。进一步深入研究这些调控网络，将有助于全面揭示NO调控莲雾果实采后絮状绵软的分子机制。五、关键基因与转录因子的功能验证5.1关键基因的克隆与表达分析根据转录组测序结果，挑选与细胞壁代谢、能量代谢、激素信号转导等密切相关的关键基因，如PG、ATP合酶、ERF1等基因进行克隆。以莲雾果实总RNA为模板，经反转录合成cDNA第一链，然后利用设计的特异性引物进行PCR扩增。PCR产物经1%琼脂糖凝胶电泳检测，结果如图9所示。图9关键基因PCR扩增产物电泳图从图9可以看出，各关键基因均扩增出预期大小的条带，其中PG基因扩增条带大小约为[X]bp，与预期的基因长度相符；ATP合酶基因扩增条带大小约为[X1]bp；ERF1基因扩增条带大小约为[X2]bp。将扩增得到的目的片段回收纯化后，与克隆载体pMD19-T连接，转化大肠杆菌DH5α感受态细胞，挑取单菌落进行PCR鉴定和测序验证。测序结果表明，克隆得到的关键基因序列与转录组测序结果一致，成功获得了这些关键基因的全长cDNA序列。为验证关键基因在转录水平的表达变化，利用实时定量PCR（qRT-PCR）技术对NO处理组和对照组莲雾果实中关键基因的表达水平进行检测，以β-actin基因作为内参基因，采用2⁻ΔΔCt法计算基因的相对表达量。结果如图10所示。图10关键基因在不同处理下的相对表达量从图10可以看出，在对照组中，随着贮藏时间的延长，PG基因表达量逐渐升高，在采后9-12天显著增加，这与转录组测序结果一致，表明PG基因的表达上调可能促进莲雾果实细胞壁中果胶的降解，加速果实的软化。而在NO处理组中，PG基因表达量在整个贮藏期间均显著低于对照组，说明NO处理抑制了PG基因的表达，从而延缓了果实的软化进程。ATP合酶基因在对照组中的表达量随贮藏时间延长逐渐下降，尤其在采后9-12天下降明显；在NO处理组中，ATP合酶基因表达量的下降趋势得到明显抑制，在采后各时间点表达量均高于对照组。这表明NO处理有助于维持ATP合酶基因的表达，保证细胞内ATP的合成，为果实的生理活动提供充足的能量，延缓果实的衰老和絮状绵软。ERF1基因在对照组中的表达量随着贮藏时间的延长逐渐上调，在采后9-12天急剧增加；而在NO处理组中，ERF1基因表达量的上调受到明显抑制，在整个贮藏期间表达量均显著低于对照组。这说明NO可能通过抑制ERF1基因的表达，阻断乙烯信号转导途径，减少乙烯对果实成熟和衰老相关基因的诱导作用，从而延缓莲雾果实的絮状绵软进程。通过qRT-PCR验证结果与转录组测序分析结果基本一致，表明转录组测序数据可靠，进一步证实了NO对莲雾果实采后关键基因表达的调控作用，这些关键基因在NO调控莲雾果实采后絮状绵软进程中发挥着重要作用。为进一步验证关键基因在蛋白水平的表达变化，采用Westernblot技术检测PG、ATP合酶、ERF1等关键蛋白的表达水平。提取NO处理组和对照组莲雾果实不同时期的总蛋白，经SDS-PAGE电泳分离后，转移至PVDF膜上，用特异性抗体进行免疫印迹检测。结果如图11所示。图11关键蛋白在不同处理下的Westernblot检测结果从图11可以看出，在对照组中，PG蛋白表达量随着贮藏时间的延长逐渐增加，在采后9-12天显著升高，这与PG基因在转录水平的表达变化趋势一致，进一步表明PG蛋白的积累可能促进果实的软化。而在NO处理组中，PG蛋白表达量在整个贮藏期间均显著低于对照组，说明NO处理在蛋白水平也抑制了PG的表达，从而延缓果实的软化。ATP合酶蛋白在对照组中的表达量随贮藏时间延长逐渐下降，在采后9-12天下降明显；在NO处理组中，ATP合酶蛋白表达量的下降趋势得到明显抑制，在采后各时间点表达量均高于对照组。这表明NO处理在蛋白水平维持了ATP合酶的表达，保证细胞内ATP的合成，为果实的生理活动提供充足的能量，延缓果实的衰老和絮状绵软。ERF1蛋白在对照组中的表达量随着贮藏时间的延长逐渐上调，在采后9-12天急剧增加；而在NO处理组中，ERF1蛋白表达量的上调受到明显抑制，在整个贮藏期间表达量均显著低于对照组。这说明NO在蛋白水平也抑制了ERF1的表达，阻断乙烯信号转导途径，减少乙烯对果实成熟和衰老相关基因的诱导作用，从而延缓莲雾果实的絮状绵软进程。Westernblot检测结果与qRT-PCR和转录组测序结果相互印证，表明NO对莲雾果实采后关键基因的调控作用在转录水平和蛋白水平均有体现，这些关键基因在NO调控莲雾果实采后絮状绵软进程中发挥着重要作用。5.2转录因子的功能验证为深入探究转录因子在NO调控莲雾果实采后絮状绵软进程中的作用，选取AP2/ERF家族中的ERF1、MYB家族中的MYB1、bHLH家族中的bHLH1以及NAC家族中的NAC1这几个关键转录因子进行功能验证实验。采用基因过表达和沉默技术，构建关键转录因子的过表达载体和RNA干扰载体。对于基因过表达实验，以pCAMBIA1300为基础载体，利用限制性内切酶将载体和含有目标转录因子编码区的DNA片段进行双酶切，然后使用T4DNA连接酶将两者连接，构建成转录因子过表达载体。将构建好的过表达载体通过农杆菌介导法转化到莲雾果实中，具体步骤为：将含有过表达载体的农杆菌菌株（如EHA105）在含有相应抗生素的LB液体培养基中振荡培养至对数生长期，然后将菌液离心收集，用含有乙酰丁香酮的重悬液重悬菌体，调整菌液浓度至OD₆₀₀=0.5-0.8。选取生长状态良好的莲雾果实，在果实表面用无菌针扎孔，然后将果实浸泡在农杆菌菌液中，真空渗透处理15-20min，使农杆菌能够进入果实组织。处理后的果实置于25℃、湿度85%的条件下暗培养2-3天，待农杆菌侵染完成后，将果实转移至正常光照条件下培养。对于基因沉默实验，利用RNA干扰（RNAi）技术构建RNA干扰载体。根据目标转录因子的cDNA序列，选取一段长度为200-300bp的特异性片段，设计并合成对应的短发夹RNA（shRNA）序列。将shRNA序列克隆到RNAi载体（如pFGC5941）中，构建成RNA干扰载体。同样通过农杆菌介导法将RNA干扰载体转化到莲雾果实中，转化方法与过表达载体转化相同。分别设置过表达对照组（转化空载体pCAMBIA1300）、RNA干扰对照组（转化空RNAi载体pFGC5941）以及空白对照组（不进行任何转化处理）。在转化后的不同时间点，对莲雾果实进行各项指标的测定和分析。测定果实硬度，结果如图12所示。在过表达ERF1的果实中，果实硬度在贮藏后期下降速度明显加快，显著低于过表达对照组和空白对照组；而在沉默ERF1的果实中，果实硬度下降速度减缓，在贮藏后期显著高于RNA干扰对照组和空白对照组。这表明ERF1的过表达促进了莲雾果实的软化，而ERF1的沉默则延缓了果实的软化进程，说明ERF1在NO调控莲雾果实采后絮状绵软中起到促进果实软化的作用。图12不同处理下莲雾果实硬度变化检测果胶酶活性，结果如图13所示。过表达MYB1的果实中，果胶酶活性显著升高，在贮藏后期明显高于过表达对照组和空白对照组；沉默MYB1的果实中，果胶酶活性升高受到抑制，在贮藏后期显著低于RNA干扰对照组和空白对照组。这说明MYB1的过表达促进了果胶酶活性的增加，加速了细胞壁中果胶的降解，从而促进果实软化；而MYB1的沉默则抑制了果胶酶活性的升高，延缓了果实的软化，表明MYB1在NO调控莲雾果实采后絮状绵软中通过调控果胶酶活性来影响果实软化。图13不同处理下莲雾果实果胶酶活性变化分析乙烯释放量，结果如图14所示。过表达ERF1的果实中，乙烯释放量在贮藏期间显著增加，明显高于过表达对照组和空白对照组；沉默ERF1的果实中，乙烯释放量显著降低，在贮藏期间显著低于RNA干扰对照组和空白对照组。这表明ERF1的过表达促进了乙烯的合成和释放，加速了果实的成熟和衰老；而ERF1的沉默则抑制了乙烯的合成和释放，延缓了果实的成熟和衰老进程，进一步证实ERF1在NO调控莲雾果实采后絮状绵软中通过乙烯信号转导途径发挥作用。图14不同处理下莲雾果实乙烯释放量变化通过基因过表达和沉默实验，明确了关键转录因子ERF1、MYB1等在NO调控莲雾果实采后絮状绵软进程中的功能。ERF1通过促进乙烯合成和信号转导以及促进果实软化来加速莲雾果实的絮状绵软进程；MYB1通过调控果胶酶活性来影响细胞壁代谢，进而影响果实的软化。这些结果为深入理解NO调控莲雾果实采后絮状绵软的分子机制提供了重要的实验依据。5.3验证结果与转录组学分析的一致性通过实时定量PCR、Westernblot以及基因过表达和沉默等实验对关键基因和转录因子进行功能验证，验证结果与转录组学分析结果呈现出高度的一致性。在关键基因表达分析方面，实时定量PCR和Westernblot结果与转录组测序所揭示的基因表达变化趋势相契合。以PG基因、ATP合酶基因和ERF1基因的验证为例，在转录组测序分析中，PG基因在对照组中随着贮藏时间延长表达量逐渐上调，在采后9-12天急剧增加，而在NO处理组中表达量的上调受到抑制；实时定量PCR结果显示，对照组中PG基因表达量在采后逐渐升高，在9-12天显著增加，NO处理组在整个贮藏期间均显著低于对照组，两者变化趋势一致。从蛋白水平来看，Westernblot检测结果也表明，对照组中PG蛋白表达量随着贮藏时间延长逐渐增加，在9-12天显著升高，NO处理组在整个贮藏期间均显著低于对照组，这进一步验证了转录组测序和实时定量PCR的结果，表明PG基因在转录水平和蛋白水平的表达变化与转录组学分析结果相符。ATP合酶基因在转录组测序分析中，对照组表达量随贮藏时间延长逐渐下降，NO处理组下降趋势受到抑制；实时定量PCR结果与之相同，对照组中ATP合酶基因表达量逐渐下降，在9-12天下降明显，NO处理组表达量下降趋势得到抑制，在各时间点均高于对照组。Westernblot检测显示，ATP合酶蛋白表达量在对照组中逐渐下降，NO处理组下降趋势得到抑制，在各时间点表达量均高于对照组，从转录水平和蛋白水平验证了转录组学分析结果。ERF1基因在转录组测序分析中，对照组表达量随着贮藏时间延长逐渐上调，在9-12天急剧增加，NO处理组表达量上调受到抑制；实时定量PCR和Westernblot结果均显示，对照组中ERF1基因和蛋白表达量逐渐上调，在9-12天急剧增加，NO处理组表达量上调受到抑制，在整个贮藏期间均显著低于对照组，验证了转录组学分析结果。在转录因子功能验证方面，基因过表达和沉默实验结果与转录组学分析中关键转录因子的调控作用相符。以ERF1和MYB1转录因子为例，转录组学分析表明，ERF1可能通过促进乙烯合成和信号转导以及促进果实软化来加速莲雾果实的絮状绵软进程；基因过表达实验中，过表达ERF1的果实硬度在贮藏后期下降速度明显加快，乙烯释放量显著增加，表明ERF1的过表达促进了果实的软化和乙烯的合成与释放，加