《Frontiers in Plant Science》：Transcriptome analysis of fenugreek under mixed saline-alkali stress
编辑推荐：
本文聚焦土壤盐碱化对全球农业的严峻挑战，探讨了潜在耐盐碱作物——胡芦巴（Trigonella foenum-graecum L.）的分子响应机制。研究人员通过对四个胡芦巴种源进行幼苗期混合盐碱胁迫评估，筛选出耐受性最强的青海种源，并利用RNA-Seq（转录组测序）技术解析其叶片在胁迫下的基因表达谱。研究发现，青海种源耐受性最强，并鉴定出2041个差异表达基因（DEG），这些基因主要富集于类黄酮生物合成、植物激素信号转导、淀粉和蔗糖代谢等关键通路，并识别出多个属于AP2/ERF、bHLH、MYB和NAC家族的转录因子。该研究为解析胡芦巴耐盐碱分子机制提供了新见解，为培育抗逆品种提供了宝贵的遗传资源。
土壤盐碱化，一个听上去有些专业的名词，却正在无声地蚕食着全球的农业命脉。想象一下，近10亿公顷的土地因为盐碱化而难以耕种，这相当于全球7%的陆地面积受到影响，其中还包括20%的耕地和高达50%的灌溉农田。这不仅意味着作物减产，更威胁着粮食安全和生态的可持续发展。面对这一全球性环境难题，科学家们将目光投向了培育能够“吃苦耐劳”的耐盐碱作物，而深入了解植物在盐碱胁迫下的“内部斗争”——其分子响应机制，则是实现这一目标的关键钥匙。胡芦巴（Trigonella foenum-graecum L.），这种富含类黄酮、皂苷等多种生物活性物质的传统草本植物，以其对不良土壤和气候的强适应性，特别是表现出潜在的耐寒、耐旱特性，被视为一种有潜力的耐盐碱作物。然而，在此之前，还没有研究利用先进的转录组测序技术来探究胡芦巴如何应对盐碱混合胁迫。为了填补这一空白，解开胡芦巴耐受盐碱的分子密码，一项深入研究就此展开。相关研究成果发表于《Frontiers in Plant Science》。
为了系统探究胡芦巴的耐盐碱机制，研究人员首先选取了来自中国青海、甘肃、云南和安徽的四个胡芦巴种源材料。他们模拟了青海省典型盐碱土壤的离子组成，配制了以HCO3-为主导的混合盐碱溶液，并对幼苗期植株施加了0（对照）、50、100、150和200 mmol/L五个浓度的胁迫处理。通过对生长指标（如株高、地上部干重）和一系列生理生化指标（如丙二醛MDA、脯氨酸Pro、抗氧化酶SOD/POD/CAT活性、可溶性糖和可溶性蛋白含量）的综合测定与评价，研究人员明确了四个种源的耐盐碱性强弱顺序为：青海 > 安徽 > 云南 > 甘肃。基于此，耐受性最强的青海种源被选为后续转录组研究的材料。研究人员设置了对照（0 mmol/L）、中度胁迫（100 mmol/L）和重度胁迫（200 mmol/L）三个处理组，每组三个生物学重复，在胁迫处理20天后采集叶片样本。利用Illumina NovaSeq 6000平台进行高通量转录组测序，对获得的测序数据进行从头组装、基因功能注释、差异表达基因筛选及功能富集分析，并通过qRT-PCR对部分差异表达基因进行验证，以确保RNA-Seq结果的可靠性。
3.1 方差分析、3.2 混合盐碱胁迫对胡芦巴幼苗生长的影响及3.3 混合盐碱胁迫对胡芦巴幼苗生理指标的影响
研究人员首先通过方差分析确认了种源、盐碱胁迫浓度及其互作对胡芦巴幼苗的生长和生理指标均存在显著或极显著影响，这些指标是有效的胁迫响应指标。生长结果表明，随着盐碱胁迫浓度的升高，胡芦巴幼苗的地上部干重和株高均显著下降。在不同种源间，青海胡芦巴在胁迫下保持了最高的干重和株高。生理指标分析显示，盐碱胁迫诱导了抗氧化酶系统活性和渗透调节物质含量的变化，且存在种源差异。例如，青海胡芦巴表现出最高的POD和SOD活性，而安徽胡芦巴的CAT和脯氨酸含量最高。MDA（膜脂过氧化产物）含量则随胁迫先升后降，在重度胁迫下有所降低，这可能与激活的抗氧化防御系统减轻了氧化损伤有关。
3.4 不同种源胡芦巴幼苗盐碱胁迫影响的综合评价
为了更全面、客观地评价耐盐碱性，研究人员采用了隶属函数法对多个生长和生理指标进行综合评价。计算出的综合评价值（D值）明确显示，在所有测试浓度下，青海胡芦巴和安徽胡芦巴的综合耐受性表现最佳。最终的整体耐盐碱性排序为：青海胡芦巴 > 安徽胡芦巴 > 云南胡芦巴 > 甘肃胡芦巴。
3.5 转录组测序样品总RNA质量评估、3.6 转录组测序与组装及3.7 单基因功能注释
基于生理评价结果，选择耐受性最强的青海胡芦巴进行转录组测序。九个样品（三个处理各三个重复）的总RNA质量均符合建库要求。测序共获得64.45 Gb高质量数据，经从头组装得到47,757条单基因（Unigene）。主成分分析显示，不同处理组的样本能明显区分开，表明组内重复性好且处理效应显著。功能注释结果显示，共有27,300条单基因在各大公共数据库中得到注释，注释率为57.16%。其中，在NR数据库中最相似的物种是蒺藜苜蓿（Medicago truncatula，71%）。KOG分类显示，注释基因广泛参与“仅具一般功能预测”、“翻译后修饰、蛋白周转和分子伴侣”等多种功能。GO分类则显示，注释基因在“结合”、“催化活性”、“细胞过程”、“代谢过程”等类别中富集。
3.8 基因表达分析及3.9 胡芦巴盐碱胁迫下差异基因表达分析
3.9.1 差异表达基因的筛选
通过比较不同处理组间的基因表达，共鉴定出1620个独特的差异表达基因。其中，QH CK（对照） vs QH 200（重度胁迫）的差异基因数量最多，达1518个（799个上调，719个下调），表明胁迫强度越大，引发的转录组变化越剧烈。QH CK vs QH 100（中度胁迫）有328个差异基因，QH 100 vs QH 200有195个差异基因。有23个基因在三个比较组中均存在差异表达。
3.9.2 差异表达基因的功能富集分析
3.9.2.1 DEGs的GO分类与富集分析
GO富集分析表明，这些差异表达基因显著富集在“代谢过程”、“细胞过程”、“对刺激的反应”等生物学过程，“催化活性”、“结合”、“转运活性”等分子功能，以及“细胞解剖实体”、“细胞内”等细胞组分类别。这提示盐碱胁迫广泛影响了胡芦巴细胞的代谢、应激响应和物质转运等功能。
3.9.2.2 DEGs的KEGG注释与通路富集分析
KEGG通路富集分析揭示了响应盐碱胁迫的关键代谢和信号通路。在QH CK vs QH 200的比较中，差异基因显著富集于“植物-病原互作”、“植物激素信号转导”、“淀粉和蔗糖代谢”、“苯丙烷生物合成”、“类黄酮生物合成”、“MAPK信号通路-植物”等通路。这些通路涉及胁迫感知、信号传导、能量代谢以及次生代谢产物的合成，共同构成了胡芦巴应对盐碱胁迫的复杂调控网络。
3.9.2.3 DEGs的COG注释
COG注释将部分差异基因分类到15个功能类别中，其中“次生代谢产物的生物合成、运输和分解代谢”以及“碳水化合物的转运和代谢”类别占比较高，进一步印证了代谢重编程在胁迫响应中的核心地位。
此外，研究还从差异表达基因中鉴定出24个与非生物胁迫相关的转录因子，主要属于AP2/ERF、bHLH、MYB和NAC家族，这些转录因子是调控下游胁迫响应基因表达的关键开关。最后，通过qRT-PCR对随机选择的9个差异表达基因进行验证，结果与RNA-Seq数据高度一致，证实了转录组分析结果的可靠性。
结论与讨论
本研究通过整合生理表型与转录组分析，系统阐明了胡芦巴应对混合盐碱胁迫的响应机制。主要结论如下：首先，明确了四个不同地理来源胡芦巴种质的耐盐碱性存在差异，其中青海种源综合耐受性最强。其次，首次绘制了胡芦巴在盐碱胁迫下的叶片转录组图谱，鉴定出大量胁迫响应相关的差异表达基因。这些基因主要参与并调控了类黄酮生物合成、植物激素（如生长素、细胞分裂素、脱落酸等）信号转导、淀粉和蔗糖代谢等关键通路。这些通路的激活有助于胡芦巴通过合成保护性次生代谢产物、调节激素平衡以协调生长与防御、以及调整碳代谢为胁迫响应提供能量和渗透调节物质，来应对盐碱胁迫造成的氧化损伤、离子毒害和渗透胁迫。再者，研究鉴定出多个胁迫相关的转录因子家族（如AP2/ERF、MYB、bHLH、NAC），它们可能作为核心调控节点，整合上游胁迫信号并激活下游一系列功能基因的表达，从而形成复杂的转录调控网络。
该研究的意义在于：理论层面，它首次在转录组水平揭示了胡芦巴响应盐碱胁迫的分子蓝图，深入解析了其耐盐碱的潜在分子机制，特别是明确了关键代谢通路和转录调控因子的作用，为理解豆科植物乃至更广泛植物的耐盐碱生物学提供了新的见解和基因资源。应用层面，研究所鉴定出的关键差异表达基因和转录因子，是后续进行功能验证和分子标记辅助育种的宝贵候选基因资源。利用这些遗传资源，有望通过分子育种手段培育出具有更强盐碱耐受性的胡芦巴乃至其他作物新品种，这对于开发利用盐碱化土地、保障粮食安全和农业可持续发展具有重要的实践价值。综上所述，这项研究为挖掘和利用胡芦巴的耐盐碱遗传潜力奠定了坚实的理论基础。
