棉花响应高温胁迫的转录组学分析

《Plant Physiology and Biochemistry》：Transcriptomic analysis of response to high-temperature stress in cotton

　　

高温胁迫（HT）是限制全球作物生产的主要非生物胁迫因子。棉花（Gossypium hirsutum L.）对热胁迫尤为敏感，但其潜在的分子机制尚不清楚。研究人员整合生理学、转录组学及功能分析，阐明了棉花材料中棉所087（ZMS087）对HT的响应。田间试验表明，热胁迫通过损害花粉活力、降低叶绿素合成及叶面积、改变生物量分配，从而导致显著减产。控制环境实验表明，HT破坏了氧化稳态和光合能力，表现为脂膜过氧化程度增加和渗透调节物质积累、抗氧化酶活性降低以及叶绿素含量下降。转录组测序鉴定出5,211个差异表达基因，其中棉子糖代谢过程显著激活。GhRS2和GhGOLS1作为关键上调基因出现，与棉子糖含量增加相对应。通过病毒诱导基因沉默（VIGS）进行的功验证表明，沉默任一基因均会通过降低棉子糖积累、削弱抗氧化酶活性并加剧氧化损伤而损害耐热性。此外，外源施加棉子糖（20–200 mg/L）通过改善叶片水分状况、提高存活率并促进生物量积累，缓解了热诱导的损伤。综上所述，本研究在棉花中提供了物种特异性功能证据，表明GhRS2/GhGOLS1介导的棉子糖生物合成途径正向调控耐热性，并揭示了热胁迫下棉子糖代谢与抗氧化反应之间的机制性联系。

随着全球气候变暖加剧，高温胁迫（HT）已成为农业生产面临的主要气象胁迫因子，对作物产量和品质构成严重威胁。棉花（Gossypium hirsutum L.）作为重要的经济作物和纺织原料，其生长、发育及产量形成高度依赖环境条件。由于植物无法像动物一样通过迁移躲避不利环境，其代谢、生长和发育等生理生化过程对温度变化极为敏感。持续超过最适温度范围的环境高温会导致不可逆的损伤，当温度超过42℃时，大多数植物即遭受热胁迫。热胁迫引发叶片和茎秆灼伤、加速落叶和衰老、抑制营养生长、阻碍果实发育，最终导致生产力下降。因此，解析植物响应HT的分子机制对于保障棉花生产具有重要意义。

在长期进化与环境适应过程中，植物形成了复杂的多途径调控网络以调节生理和代谢过程，从而减轻胁迫损伤。可溶性糖作为重要的渗透调节物质备受关注。在高温条件下，植物抗氧化酶活性下降，大量活性氧（ROS）自由基产生，导致氧化损伤。而棉子糖家族寡糖（RFOs）及其前体半乳糖醇能够在胁迫条件下缓解氧化损伤。研究表明，鹰嘴豆在高温和氧化胁迫下半乳糖醇和棉子糖含量显著增加，但其转录调控机制，特别是在棉花响应热胁迫中的作用，尚不清楚。基于上述背景，研究人员以棉花品种中棉所087（ZMS087）为材料，开展生理学和转录组学分析，旨在从转录组水平深入理解棉花热胁迫响应机制，为改良棉花耐热性提供科学依据。该论文发表于《Plant Physiology and Biochemistry》。

研究主要采用了以下关键技术方法：田间试验结合控制环境实验进行表型与生理学分析；Illumina NovaSeq X Plus平台双端150 bp测序策略进行转录组测序（RNA-seq）及差异表达基因（DEG）分析；液相色谱-质谱联用（LC-MS）进行代谢物定量；病毒诱导基因沉默（VIGS）进行基因功能验证；以及定量实时PCR（RT-qPCR）进行基因表达验证。

研究结果部分包括以下内容。

3.1 高温胁迫对棉花田间生长发育的影响。田间条件下，HT显著损害棉花生长发育，导致花粉活力下降、叶绿素合成受抑、叶面积减少及生物量分配失衡，最终造成显著减产。花粉活力在HT处理24小时和48小时后分别降低35.88%和45.46%；蕾铃脱落率从48.37%显著上升至70.67%；土壤植物分析发育（SPAD）值显著下降8.5，表明叶绿素降解增强、光合作用受抑；单株叶面积从1755 cm2锐减至1375.82 cm2；各器官干重均显著降低，其中果枝和棉铃干重分别减少35.62%和38.84%。

3.2 棉花对高温胁迫的表型及生理响应。控制环境HT试验显示，HT处理48小时后植株出现严重萎蔫，存活率显著降低47.83%。超氧化物歧化酶（SOD）活性在HT处理12小时后显著升高，过氧化物酶（POD）活性在24和48小时后分别显著增加103.5%和49.55%。丙二醛（MDA）含量和脯氨酸（Pro）含量在HT处理后均显著上升。与此同时，叶绿素a、叶绿素b及总叶绿素含量在48小时HT处理后分别显著降低21.35%、29.63%和24.71%。

3.3 棉花高温胁迫下的转录组分析。转录组测序鉴定出5,211个差异表达基因，其中3,045个上调、2,166个下调。基因本体论（GO）富集分析显示，差异表达基因主要富集于"对非生物刺激的响应""对温度刺激的响应"等生物学过程。通过设定筛选标准（FPKM ≥ 10、FDR < 0.05、倍数变化≥1.5并经RT-qPCR验证），鉴定出Gh_D01G141200和Gh_Contig00249G000600两个关键基因，分别命名为半乳糖醇合成酶（GhGOLS1）和棉子糖合成酶（GhRS2）。GO分析中，上调差异表达基因独特富集于棉子糖代谢过程，且LC-MS检测证实HT处理后棉子糖含量显著增加15.29%。

3.4 GhGOLS1——棉花耐热性的正向调控因子。通过VIGS验证，GhGOLS1沉默植株中该基因转录水平降低约76.52%–88.58%。HT处理后，GhGOLS1沉默植株存活率显著降低37.93%，棉子糖含量显著减少13.49%，硝基蓝四氮唑（NBT）和3,3'-二氨基联苯胺（DAB）染色显示活性氧积累加剧，SOD和POD活性下降，MDA含量显著升高。

3.5 棉子糖积累后棉花对高温胁迫的表型及生理响应。GhRS2沉默植株中该基因转录水平降低65.95%–82.25%。HT处理后，沉默植株存活率从83.33%降至40.54%，棉子糖含量显著降低，MDA含量在24和48小时后分别显著增加，NBT和DAB染色显示活性氧积累加剧，SOD和POD活性显著下降。这些结果表明GhRS2通过维持活性氧稳态在耐热性中发挥关键作用。

3.6 外源棉子糖对棉花植株表型的影响。叶面喷施不同浓度棉子糖（20、100、200 mg/L）可显著改善HT下棉花的生长状态。HT处理24和48小时后，处理组叶片含水量从对照的58.78%和55.83%显著提升至约72%和70%左右；处理48小时后存活率从59.53%分别提升至81.3%、76%和69.33%；总干重从0.15 g显著增加至0.22–0.25 g。

3.7 外源棉子糖田间验证对棉花产量和纤维品质的影响。田间花期喷施试验表明，20–200 mg/L棉子糖处理组籽棉产量从6341.70 kg/hm2分别提升至6627.06、6672.76和6718.46 kg/hm2，增幅达4.50%–5.94%。各处理组纤维品质指标（上半部平均长度、整齐度指数、断裂比强度、马克隆值、伸长率）无显著差异，表明外源棉子糖在提高产量的同时不影响纤维品质。

讨论部分，研究人员首先指出HT通过破坏活性氧动态平衡、引发氧化损伤和膜脂过氧化，最终导致细胞程序性死亡，本研究结果与前人发现一致。其次，研究人员系统阐述了GhGOLS1/GhRS2-棉子糖途径的调控机制：GhRS2作为关键上游调控因子促进棉子糖生物合成，积累的棉子糖家族寡糖通过三种非互斥机制发挥保护作用——作为渗透保护物质维持细胞膨压和蛋白质结构稳定、作为活性氧直接清除剂淬灭羟自由基和超氧阴离子、作为膜稳定剂与膜磷脂极性头部基团相互作用以维持膜完整性。研究人员特别指出，沉默植株中抗氧化酶活性下降与棉子糖减少的关联，暗示棉子糖可能通过类似分子伴侣的功能保护蛋白质结构，从而维持抗氧化酶活性，形成调控反馈环路。本研究首次在棉花中提供证据，表明棉子糖途径正向调控HT响应。最后，在农业应用层面，研究人员结合水稻耐热育种及外源物质应用案例，论证了遗传改良与外源喷施相结合的互补策略，并提出本研究鉴定的最佳浓度范围（20–200 mg/L）为田间实际应用提供了可行方案。

研究结论部分翻译如下：本研究系统阐明了棉花响应HT的生理和分子机制。主要结论如下：HT严重损害棉花生长，降低花粉活力、叶绿素含量和叶面积，破坏生物量分配，最终导致显著减产。生理学分析揭示，HT破坏了氧化稳态和光合能力，表现为脂膜过氧化增加、渗透调节物质积累、抗氧化酶活性改变以及叶绿素含量降低。转录组测序鉴定出5,211个差异表达基因，其中棉子糖代谢过程显著激活。GhRS2和GhGOLS1作为关键上调基因出现，与棉子糖含量增加相对应。通过VIGS进行的功能验证表明，沉默任一基因均会通过降低棉子糖积累、削弱抗氧化酶活性并加剧氧化损伤而损害耐热性，确认了它们的正向调控作用。外源施加棉子糖（20–200 mg/L）显著缓解了热诱导损伤，改善了叶片水分状况，提高了存活率，并促进了生物量积累。综上所述，本研究提供了GhRS2和GhGOLS1在热胁迫下棉花中的功能表征，并提供了将棉子糖代谢与抗氧化反应联系起来的物种特异性证据。这些发现为通过遗传或农艺手段提高棉花耐热性奠定了理论基础。