紫花苜蓿盐胁迫下转录组动态与MsTIFY11B介导盐碱耐受性的调控网络解析

《Plants》：Comparative Phytochemical Profiles and Antioxidant and Tyrosinase-Inhibitory Activities of Wild and In Vitro-Derived Curcuma lampangensis Saensouk, Maknoi & Rakarcha

　　

土壤盐渍化严重限制饲草作物生产力，但中等耐盐豆科饲草紫花苜蓿(Medicago sativa L.)中调控盐胁迫适应的网络仍待探索。研究人员在200 mM NaCl胁迫下于3个时间点检测苜蓿叶片生理与转录组动态。盐胁迫诱导Na+/K+比渐进升高，抗氧化酶双相激活并伴随丙二醛(malondialdehyde, MDA)积累。时间序列RNA-seq鉴定出3631个差异表达基因(differentially expressed genes, DEGs)与132个核心盐响应转录因子(transcription factors, TFs)。通路注释显示DEGs显著参与细胞壁生物发生、氧化还原稳态及类胡萝卜素生物合成途径，其中类胡萝卜素积累在盐胁迫下被强烈激活。加权基因共表达网络分析(weighted gene co-expression network analysis, WGCNA)构建9个共表达簇，其中与Na+积累及Na+/K+比正相关的棕色模块显著富集于植物激素信号转导途径，其72.7%富集基因为TIFY家族。核心枢纽基因MsTIFY11B亚细胞定位于细胞核；酵母异源过表达增强对盐与碱耐受，对甘露醇模拟干旱几乎无保护。结果表明苜蓿盐响应具有时间框架，MsTIFY11B可能参与盐碱耐受，是应激适应品种开发的候选基因。

研究背景与立项依据

土壤盐渍化影响全球约10.7%陆地与20%灌溉农田，气候变迁与不合理灌溉或致2050年近半农田受胁。紫花苜蓿作为高生物量豆科饲草，具固氮与根系优势，但盐渍化抑制萌发、根系与光合，中等耐盐性仍受限。虽多组学已揭示部分苜蓿盐响应基因，但转录因子(regulatory components)及精确机制不清，亟需时间维度转录组结合网络分析锁定枢纽基因。

研究人员以紫花苜蓿品种‘中牧1号’(Zhongmu No.1)为材料，水培200 mM NaCl处理0、3、24、72 h取叶，联合生理生化、时间序列RNA-seq、GO/KEGG、TF普查、WGCNA、酵母异源表达与亚细胞定位，阐明盐响应时间框架并获候选基因MsTIFY11B。论文发表于《Plants》。

主要关键技术方法

研究以紫花苜蓿‘中牧1号’水培苗为样本队列，施加200 mM NaCl于3、24、72 h采样；采用Illumina NovaSeq 6000进行时间序列RNA-seq，DESeq2判定DEGs(|log₂FC|≥1, adj p<0.05)；PlantTFDB比对识别盐响应TFs；WGCNA(软阈值幂7)构建共表达模块并关联Na⁺、Na⁺/K⁺、SOD、POD、MDA等表型；候选基因MsTIFY11B做本氏烟(Nicotiana benthamiana)瞬时表达定位于核，酿酒酵母(Saccharomyces cerevisiae) INVSc1异源过表达评测NaCl、NaHCO₃、甘露醇胁迫生长；RT-qPCR验证9基因。

研究结果

2.1 盐胁迫下苜蓿生理生化响应

叶片Na⁺升、K⁺降，Na⁺/K⁺比上升；SOD于24 h达峰、POD于3 h达峰，MDA先升后降。显示早期氧化应激与后期生理调整并存。

2.2 盐处理后苜蓿叶时间序列转录组动态

12文库获~41.8–47.1 M raw reads，34950基因表达。3631 DEGs，各时点特有为主，仅76上、100下为三时共有；RT-qPCR与RNA-seq相关R²=0.89。

2.3 DEGs的GO与KEGNA富集

GO持续富集木葡聚糖代谢(GO:0010411)、细胞壁生物发生(GO:0042546)、氧化还原过程(GO:0055114)；KEGG中类胡萝卜素生物合成(ko00906)全时点富集，碳固定(ko00710)早显、光合天线蛋白(ko00196)迟显。实测叶绿素略降、类胡萝卜素显著积累。

2.4 苜蓿盐响应转录因子鉴定

132非冗余TFs属31家族，ERF、bHLH、NAC居多；14 TF三时一致，10持续下调(5 ERF等)，3持续上调(1 ERF、2 HSF)。

2.5 表型关联WGCNA模块

9模块；black负关POD/MDA，tan正关SOD；pink/purple/brown/green-yellow正关Na⁺与Na⁺/K⁺。brown(297基因)正关Na⁺比，富集亚麻酸代谢(ko00591)、萜类骨架(ko00900)、植物激素信号(ko04075)；blue(1111基因)负关Na⁺正关MDA，富集核糖体(ko03010)。brown网络含Cytochrome P450、AP2/ERF、TIFY、Bet v I等44节点。

2.6 MsTIFY11B亚细胞定位与酵母功能

brown中TIFY占激素通路11基因中8个，枢纽MsG0880042776.01拟南芥同源AtTIFY11B(35%一致)命名MsTIFY11B。本氏烟显示仅核定位。酵母中pYES2-MsTIFY11B在4 M NaCl、1.5 M NaHCO₃下10^?3–10^?4稀释生长优于空载体，1.5 M甘露醇无差。

讨论浓缩

离子稳态上Na⁺/K⁺失衡伴SOD/POD双相与MDA波动，氧化还原相关GO印证酶活；细胞壁木葡聚糖与生物发生富集提示壁重塑抗渗；类胡萝卜素持续富集与实测积累支持其作为非酶抗氧化物淬灭单线态氧及ABA前体角色，光合从碳固定转向天线蛋白微调防光抑。

TF层面ERF最大、bHLH/NAC/MYB/HSF次之，14个跨时TF暗示长期适应开关。

WGCNA中blue以核糖体蛋白为主映翻译机器补偿；brown正关Na⁺且富集激素与亚麻酸(茉莉酸(jasmonate, JA)前体)，其中TIFY占优。酵母无内源JA通路却仍耐盐碱，提示MsTIFY11B具JA非依赖的保守生化属性（如氧化还原维持），与MsTF2在酵母/大肠杆菌中功能类似。文章指明局限：仅在酵母验证，需苜蓿过表达/编辑补证。

结论翻译

本研究整合生理、时间序列RNA-seq与WGCNA阐明苜蓿盐响应网络。苜蓿通过生理适应、类胡萝卜素积累与时间特异转录程序抵御盐致渗透与氧化损伤。共表达网络识别核定位枢纽MsTIFY11B参与植物激素信号转导；酵母异源表达确认其作为盐碱耐受的正向特异调节因子，对干旱保护有限。结果拓展苜蓿TIFY成员功能多样性认知，MsTIFY11B是盐碱与边际土壤改良候选。