《Plant Stress》：Identification of genetic loci and candidate genes for salt tolerance in barley (
Hordeum vulgare L.) at the germination stage using a genome-wide association study and multi-omics analysis
编辑推荐：
为解决土壤盐渍化严重制约作物生产的问题，本研究通过整合全基因组关联研究(GWAS)、转录组和代谢组学，深入解析了大麦(Hordeum vulgare L.)萌芽期耐盐性的复杂调控网络。研究在250份大麦种质中鉴定出306个显著SNPs，揭示了盐胁迫响应中氨基酸与苯丙烷类生物合成等关键通路，并筛选出24个与抗氧化酶、膜蛋白转运、转录因子等相关的候选基因，为培育耐盐大麦品种提供了宝贵的遗传资源和理论依据。
在全球范围内，有近十亿公顷的土地受到盐渍化的威胁，不合理的灌溉和剧烈的气候变化更使盐渍地面积持续扩张。土壤盐分不仅直接抑制植物生长，更对全球粮食安全构成了严峻挑战。大麦，作为全球第四大谷类作物，因其相对较强的耐盐性而被视为开发利用盐碱地的“先锋作物”。然而，植物在种子萌发阶段对盐胁迫最为敏感，而作物的耐盐性又是一个由众多基因相互作用的复杂数量性状，其背后的遗传基础和分子机制仍有待深入解析。为了挖掘大麦萌芽期耐盐性的关键基因资源，并为分子育种提供靶点，这项研究应运而生。
为了系统解答上述问题，研究人员开展了一项综合性研究，成果发表在期刊《Plant Stress》上。该研究通过对一个包含250份不同来源大麦种质的自然群体，在连续两年内施加200 mM NaCl盐胁迫，评估了相对发芽率、相对根/芽长、相对根/芽鲜重等五个关键形态性状。随后，整合了全基因组关联研究(GWAS)、非靶向代谢组学(基于LC-MS/MS)和转录组学(RNA-seq)分析，并结合实时定量PCR(RT-qPCR)验证，系统鉴定了与耐盐性相关的遗传位点、差异代谢物、差异表达基因及核心候选基因，最终提出了大麦萌芽期响应盐胁迫的潜在分子机制模型。
本研究主要运用了以下几项关键技术方法：首先，基于250份大麦自然种质（部分材料来自中国江苏沿海地区农业科学研究所）的表型数据，利用混合线性模型(MLM)进行了全基因组关联研究(GWAS)。其次，以耐盐品种“Liu Leng Zi Da Mai”(LLZDM)的根组织为材料，通过液相色谱-串联质谱(LC-MS/MS)进行非靶向代谢组学分析。同时，对盐处理和对照条件下的LLZDM根组织进行了转录组测序(RNA-seq)。最后，通过整合GWAS和转录组数据筛选候选基因，并使用RT-qPCR在多个品种中对部分基因的表达进行了验证。
研究结果如下：
3.1. 表型变化分析与性状间相关性
盐处理显著抑制了所有大麦基因型的芽和根生长，各性状的相对值均显著降低。250份大麦种质在耐盐性上表现出显著的遗传变异。其中，地方品种如“Liu Leng Zi Da Mai”(LLZDM)比栽培品种如“Nasi Nijo”(NN)表现出更强的耐盐性。在正常条件下两者无显著差异，但在盐胁迫下，NN的性状下降幅度(50.00%至78.03%)远大于LLZDM(21.16%至39.93%)，清晰证明了LLZDM的耐盐优势。
3.2. 萌芽期耐盐性的全基因组关联分析
利用106,131个高质量SNP进行GWAS分析，以-log10P≥ 4.0为阈值，在两年间共检测到306个与五个性状显著相关的SNP。其中，有23个SNP被不同性状同时检测到，8个SNP在两年中均被稳定检测到，这些稳定SNP主要定位于2号和6号染色体。单个SNP所能解释的表型变异率(PVE)范围为0.09%至3.91%。
3.3. 大麦根响应盐胁迫的代谢组学分析
对耐盐品种LLZDM根组织的非靶向代谢组学分析表明，盐胁迫诱导了324个差异积累代谢物(DAMs)，其中167个上调，157个下调。这些代谢物主要属于脂质和类脂分子、有机酸及其衍生物等类别。京都基因与基因组百科全书(KEGG)通路富集分析显示，DAMs显著富集于生物素代谢、氨基糖和核苷酸糖代谢等通路。
3.4. 盐胁迫条件下大麦根的转录组学分析概览
RNA-seq分析在LLZDM根组织中鉴定出2351个差异表达基因(DEGs)，其中1544个上调，807个下调。基因本体(GO)富集分析显示，DEGs主要富集于离子转运、对激素响应、胁迫响应等功能类别。KEGG通路分析进一步表明，DEGs显著富集于苯丙烷生物合成、丙氨酸、天冬氨酸和谷氨酸代谢、谷胱甘肽代谢等通路。
3.5. 代谢组与转录组的整合分析
整合分析发现，DEGs和DAMs共同显著富集于“氨基酸的生物合成”和“苯丙烷生物合成”等关键通路。在氨基酸生物合成通路中，盐胁迫导致谷氨酸、脯氨酸等含量显著上升，而赖氨酸、组氨酸等含量下降，同时有12个相关基因上调。在苯丙烷生物合成通路中，盐胁迫诱导了咖啡酰辅酶A、圣草酚等代谢物的积累，并上调了包括5个肉桂酰辅酶A还原酶(COMT)基因在内的多个关键基因的表达。
3.6. 候选基因的鉴定与验证
通过整合GWAS和RNA-seq数据，最终鉴定出24个与耐盐性相关的候选基因。这些基因在染色体上分布不均，其中17个上调，7个下调。根据功能可分为几类：与抗氧化酶相关，如过氧化氢酶(CAT)、谷胱甘肽S-转移酶(GST)；与膜蛋白转运相关，如跨膜蛋白184C(TP184C)、钾转运蛋白(KT)；转录因子，如一个WRKY转录因子显著上调；以及受体蛋白，如受体样蛋白激酶(RLKs)。RT-qPCR对其中12个基因的表达验证结果与转录组数据基本一致，证实了筛选结果的可靠性。
归纳研究结论与讨论：
本研究通过多组学整合策略，系统揭示了大麦萌芽期耐盐性的复杂调控网络。结论表明，盐胁迫下大麦通过协调氨基酸代谢和苯丙烷类代谢来应对氧化胁迫和离子失衡。其中，脯氨酸、谷氨酸等渗透调节物质的积累，以及苯丙烷类途径中COMT等基因的上调和类黄酮等抗氧化物质的积累，是重要的耐受机制。
更重要的是，研究筛选出的24个核心候选基因，包括GST、TP184C、KT、WRKY和RLKs等，分别从抗氧化防御、离子稳态、胁迫信号感知与转导等不同层面发挥作用。基于这些发现，论文提出了一个耐盐机制模型：盐胁迫信号首先被细胞膜上的RLKs等受体感知，进而激活WRKY等转录因子，后者调控下游靶基因（如抗氧化酶基因、离子转运蛋白基因等）的表达，最终通过增强ROS清除能力、维持离子平衡和促进保护性代谢物合成，来改善大麦在盐胁迫下的发芽和幼苗生长。
这项研究的意义重大。它不仅为理解大麦萌芽期耐盐性的分子机制提供了新的见解，更重要的是，所鉴定的306个显著SNP和24个候选基因，特别是那些在两年中稳定出现的位点，为后续通过基因组辅助选择培育耐盐大麦品种提供了宝贵的遗传资源和分子靶点。尽管基于有限样本的转录组分析可能遗漏其他种质中的有利等位基因，且单个标记的解释率有限，但将多个位点的优良等位基因聚合，有望显著提升品种的耐盐性。未来，通过增加测序深度、扩大群体规模以及利用基因编辑等技术对候选基因进行功能验证，将进一步推动大麦耐盐分子育种的发展，对于利用广大盐碱地、保障粮食安全具有重要的实践价值。
