《Plant Stress》：Physiological and transcriptomic insights into cauliflower seedlings’ responses to multiple abiotic stresses
编辑推荐：
语
本研究针对花椰菜幼苗对低温、高温、干旱、盐渍及淹水五种常见非生物胁迫的响应机制不明问题，通过整合生理指标测定与转录组测序（RNA-seq），系统解析了不同胁迫下的生长表型、抗氧化酶活性（SOD/POD/CAT）、差异表达基因（DEGs）及可变剪接事件。研究发现热胁迫诱导最多DEGs（15,276个），盐胁迫最少（3,841个）；WGCNA分析鉴定出胁迫特异性共表达模块及枢纽基因（如HSP20、MYB32）；30个共同上调基因（如MIPS1、NAC032）经qRT-PCR、亚细胞定位及酵母功能验证证实其跨胁迫耐受性。该研究为花椰菜抗逆育种提供了关键靶点与理论支撑。
论文解读
花椰菜（Brassica oleraceavar. botrytis）作为全球重要的蔬菜作物，其苗期易受低温、高温、干旱、盐渍和淹水等非生物胁迫的严重影响，导致产量和品质下降。尽管已有研究对单一胁迫响应机制有所探索，但花椰菜如何通过分子网络协调多种胁迫的适应性仍不清楚。尤其在Brassica作物中，多数研究局限于单一或双胁迫分析，缺乏对苗期多胁迫系统性比较的转录组研究。为此，宋梦菲团队在Plant Stress发表论文，通过生理与转录组整合分析，揭示了花椰菜幼苗应对五种常见非生物胁迫的特异性与共享调控机制。
研究采用花椰菜自交系“3201”为材料，在可控环境下分别施加冷（4°C）、热（40°C）、干旱（断水）、盐（1.5% NaCl灌溉）和淹水（浸没茎基部）处理，测定表型指标（株高、叶长、叶宽）和生理参数（叶绿素含量、SOD、POD、CAT活性及MDA含量）。转录组分析通过RNA-seq技术，结合差异表达基因（DEGs）筛选、GO/KEGG富集、WGCNA共表达网络构建、转录因子家族鉴定及可变剪接事件分析，系统解析了分子响应特征。关键基因（HSP20、MIPS1、MYB32、NAC032）通过qRT-PCR、亚细胞定位和酵母胁迫耐受实验验证功能。
3.1 表型与生理响应特征
热胁迫导致幼苗严重萎蔫，叶绿素降解且MDA含量极显著上升，虽激活CAT和POD活性，但植株2天内死亡，表明其急性氧化损伤不可逆。冷胁迫反而促进株高增长，通过协同激活SOD和CAT维持氧化稳态。干旱、盐胁迫抑制叶片发育并引发轻度膜脂过氧化，而淹水胁迫通过强烈诱导SOD、POD、CAT活性缓解氧化损伤，但抑制整体生长。
3.2 转录组整体响应格局
PCA与样本相关性分析显示热胁迫转录谱最独特（与对照组相关性仅0.42），盐胁迫与淹水胁迫转录谱相似性较高（相关性0.93）。热胁迫诱导DEGs数量最多（15,276个），盐胁迫最少（3,841个）。GO富集揭示胁迫特异性通路，如热胁迫富集“蛋白质折叠”，冷胁迫富集“响应低温”；KEGG分析发现“碳水化合物生物合成”等为多胁迫共享通路。30个共同上调基因中，HSP20在热胁迫下表达最显著，MIPS1在冷和淹水胁迫中高表达。
3.3 WGCNA揭示胁迫特异性模块
WGCNA鉴定出与各胁迫强相关的共表达模块：MEblue（热胁迫）富集热激蛋白与泛素化相关基因，枢纽基因为BolK_8G50810（RING/U-box蛋白）；MEbrown（冷胁迫）关联自噬通路，枢纽基因为BolK_6g18700（HB转录因子）；MEred（干旱胁迫）与ABA合成相关，枢纽基因为BolK_4g01490（CYP450）；MElightcyan（盐胁迫）和MEyellow（淹水胁迫）分别富集细胞壁生物合成和饥饿响应通路。
3.4 共同上调基因的功能验证
qRT-PCR证实HSP20、MIPS1、MYB32、NAC032在五类胁迫中均显著上调。亚细胞定位显示HSP20全细胞分布，MIPS1定位于核膜与细胞膜，MYB32和NAC032位于细胞核。酵母转化实验表明这些基因显著提升酵母在高盐（500 mM NaCl）、干旱模拟（1 M甘露醇）、高低温下的存活率，证实其广谱抗逆功能。
3.5 转录因子家族表达模式
MYB家族在各胁迫中差异表达比例最高，NAC家族在冷胁迫中特异性响应，bHLH和AP2/ERF分别主导干旱/热和盐胁迫调控。相关性网络揭示bHLH与C2C2协同表达，而与NAC、bZIP负相关，呈现家族间调控复杂性。
3.6 可变剪接的胁迫特异性调控
热胁迫诱导最多可变剪接事件（2,629个基因），44个基因在五类胁迫中均发生剪接，富集于mRNA加工通路。例如BolK_3g43910在热胁迫下产生特异性剪接变体，提示剪接 plasticity 在转录组重塑中的关键作用。
结论与意义
本研究系统揭示了花椰菜幼苗通过胁迫特异性生理调整（如抗氧化酶协同激活）与分子重编程（如DEGs、可变剪接、TF网络）适应多类逆境。核心基因（如HSP20、MYB32）的跨胁迫功能为多抗性育种提供新靶点；WGCNA模块与枢纽基因为精准改良特定胁迫耐受性奠定基础。研究首次在花椰菜中实现五类常见胁迫的跨组学比较，填补了Brassica作物多胁迫响应机制的空白，对应对气候变化下的蔬菜安全生产具有重要理论价值。
