时间序列转录组分析鉴定MsTIFY11B作为紫花苜蓿（Medicago sativa L.）耐盐碱调控的候选基因

《Plants》：Time-Course Transcriptomic Analysis Identifies MsTIFY11B as a Promising Candidate Regulator of Salt–Alkali Tolerance in Alfalfa (Medicago sativa L.)

　　

土壤盐渍化严重限制饲草作物生产力，而中等耐盐豆科饲草紫花苜蓿中调控盐胁迫适应的网络仍大多未被探索。研究人员在200 mM NaCl胁迫下于三个时间点检测苜蓿叶片的生理与转录组动态。盐胁迫诱导Na+/K+比渐进升高，抗氧化酶双相激活并伴随丙二醛积累。时间序列RNA-seq分析鉴定出3631个差异表达基因（DEGs）和132个核心盐响应转录因子（TFs）。通路与功能注释分析表明这些DEGs显著参与细胞壁生物发生、氧化还原稳态及类胡萝卜素生物合成途径，其中类胡萝卜素积累在盐胁迫下被强烈激活。利用加权基因共表达网络分析（WGCNA）构建了9个不同共表达簇。值得注意的是，与Na+积累和Na+/K+比正相关的棕色模块在植物激素信号转导途径中显著富集，该途径富集基因中72.7%属于TIFY家族。其中核心枢纽基因MsTIFY11B被分离用于功能表征。亚细胞定位证明MsTIFY11B仅定位于细胞核。酵母异源表达显示MsTIFY11B过表达增强对盐与碱的耐受性，而对甘露醇诱导的干旱胁迫几乎无保护作用。综上，研究结果为苜蓿盐响应的转录组时间框架提供了全面视图，并提示MsTIFY11B可能参与盐碱耐受，是进一步功能表征及开发耐逆苜蓿品种的有前景候选基因。

研究背景：土壤盐渍化影响全球约10.7%陆地与20%灌溉农田，预计2050年近半农业区受盐渍化威胁。紫花苜蓿（Medicago sativa L.）作为高生物量豆科饲草具中等耐盐性，但盐害仍制约其产量与品质。尽管多组学已揭示部分苜蓿盐响应基因，额外调控组件及其精确分子机制大多未知，因此亟需借助转录组学鉴定新型盐响应候选基因与通路。该研究发表于《Plants》，结合生理生化与转录组学阐明苜蓿盐适应调控过程，对培育耐盐种质具潜在应用价值。

主要关键技术方法：研究以中国北方广栽耐盐苜蓿品种中牧一号（Zhongmu No.1）为材料，设0、3、24、72 h 200 mM NaCl水培处理，每时间点3次生物学重复；采用Illumina NovaSeq 6000进行150 bp双端RNA-seq，以HISAT2比对至苜蓿参考基因组，StringTie定量FPKM，DESeq2筛选|log₂FC|≥1且padj<0.05的差异表达基因（DEGs）；对DEGs行GO与KEGG富集；基于PlantTFDB鉴定转录因子（TFs）；使用WGCNA（软阈值幂7）构建2625个非冗余DEG共表达模块并关联Na⁺、K⁺、Na⁺/K⁺比、SOD、POD、MDA等性状；以RT-qPCR验证9个DEG（R²=0.89）；选棕色模块枢纽MsTIFY11B（MsG0880042776.01）做本氏烟（Nicotiana benthamiana）亚细胞定位（35S:GFP融合）与酵母INVSc1异源表达（pYES2载体，Gal诱导），在含NaCl、NaHCO₃、甘露醇平板评估生长。

研究结果：

2.1 苜蓿盐胁迫生理生化响应：叶片Na⁺上升、K⁺下降致Na⁺/K⁺比升高；SOD在24 h达峰、POD在3 h达峰，MDA先升后降，表明早期氧化应激与后期生理调整。

2.2 盐处理后苜蓿叶片时间序列转录组动态：12文库获约41.84–47.10 M raw reads，34950表达基因，3631 DEGs；各时间点特有DEG占多数，3 h特异上调526/下调378，24 h 265/247，72 h 260/241，三时点共有上调76/下调100；RT-qPCR验证一致。

2.3 DEG的GO与KEGNA富集：GO持续富集木葡聚糖代谢（GO:0010411）、细胞壁生物发生（GO:0042546）、氧化还原过程（GO:0055114）；KEGG中类胡萝卜素生物合成（ko00906）三时点均富集，碳固定（ko00710）在3 h特异，光合天线蛋白（ko00196）在24/72 h富集；实测叶绿素略降、类胡萝卜素显著累积。

2.4 盐响应转录因子鉴定：132个非冗余TFs跨31家族，ERF、bHLH、NAC居多；14个TF三时点持续差异，含5 ERF持续下调、1 ERF与2 HSF持续上调。

2.5 WGCNA模块与性状关联：9模块中棕色模块（297基因）与Na⁺及Na⁺/K⁺比正相关，富集植物激素信号转导（ko04075）、亚油酸代谢（ko00591）、萜类骨架生物合成（ko00900）；蓝色模块（1111基因）负关Na⁺正关MDA，富集核糖体途径。棕色模块共表达网络44节点含Cytochrome P450、AP2/ERF、TIFY、Bet v I/MLP，其中植物激素通路11基因里8个为TIFY家族。

2.6 MsTIFY11B亚细胞定位与酵母功能：MsTIFY11B（与拟南芥TIFY11B同源35%）仅核定位；酵母中过表达提升4 M NaCl与1.5 M NaHCO₃耐受，对1.5 M甘露醇无显著保护。

讨论总结：离子稳态失衡伴SOD/POD双相与MDA波动反映氧化还原重编程；细胞壁重构（木葡聚糖、生物发生）与类胡萝卜素累积（光保护与非酶抗氧化、ABA前体）构成盐适应代谢特征；ERF/bHLH/NAC等TF动态微调生长-防御权衡；WGCNA棕色模块将TIFY家族锚定于激素信号与亚油酸（茉莉酸JA前体）代谢，但酵母无内源JA通路仍见保护，提示MsTIFY11B具JA非依赖的保守生化活性（如 redox 稳态）；MsTIFY11B在植物为转录调控子，在异源系统呈基础抗逆属性，限于酵母验证尚不能直接推定苜蓿内功能，需后续过表达/编辑回补。

结论翻译：本研究整合生理谱、时间序列RNA-seq与WGCNA阐明苜蓿盐响应调控网络。结果显示苜蓿通过生理适应、类胡萝卜素累积与时间区分的转录程序抵御盐诱导渗透与氧化损伤。共表达网络鉴定并表征核定位枢纽基因MsTIFY11B，其参与植物激素信号转导；酵母异源表达确认MsTIFY11B为盐与碱耐受的特异性正调控子，对干旱保护有限。结果拓宽了对苜蓿TIFY家族功能多样性的理解，并将MsTIFY11B列为盐碱与边际土壤条件下进一步探究与改良苜蓿性能的候选基因。