大麦耐旱性的生理与转录组解析：野生二棱大麦（Hordeum vulgare ssp. spontaneum）与栽培大麦的比较研究
编辑推荐：
干旱胁迫是作物生产面临的主要环境限制因素。随着气候变化导致干旱频率增加，增强耐旱性已成为大麦育种的主要目标。野生二棱大麦（Hordeum vulgare ssp. spontaneum）作为栽培大麦的野生近缘种，是育种和开发耐旱品种的重要遗传资源。本研究旨在比
干旱胁迫是作物生产面临的主要环境限制因素。随着气候变化导致干旱频率增加，增强耐旱性已成为大麦育种的主要目标。野生二棱大麦（Hordeum vulgare ssp. spontaneum）作为栽培大麦的野生近缘种，是育种和开发耐旱品种的重要遗传资源。本研究旨在比较耐旱野生基因型（H. vulgare ssp. spontaneum）与干旱敏感型栽培大麦基因型（‘Mona’品种）在对照和干旱胁迫条件下的生理及转录组响应。研究结果表明，在干旱胁迫下，野生基因型比‘Mona’品种表现出更轻微的生理和生长损伤。转录组分析共鉴定出耐旱野生基因型中2564个差异表达基因（differentially expressed genes，DEGs），以及敏感型栽培品种中2163个DEGs。值得注意的是，在耐旱基因型中，大多数DEGs呈上调表达，而在敏感型品种中，干旱诱导的基因以下调表达更为常见。基因本体论（Gene Ontology，GO）富集分析显示，这些DEGs主要与关键胁迫相关生物学过程相关，包括细胞组分组织、信号转导途径、离子转运、调控蛋白活性、活性氧（reactive oxygen species，ROS）清除、激素生物合成和渗透稳态。比较分析表明，耐旱基因型拥有更多功能有效的干旱响应基因，且转录本丰度高于敏感型基因型。总体而言，耐旱基因型在干旱胁迫下表现出持续优越的生理和功能表现。这些发现加深了对野生大麦耐旱性进化机制的理解，同时两个亚种间共享的功能机制揭示了可在水分受限环境下用于改良大麦的生理和转录性状。
大麦耐旱性的生理与转录组解析：野生二棱大麦与栽培大麦的比较研究
一、研究背景与科学问题
全球气候变化导致干旱事件频发，干旱胁迫已成为制约作物生产的主要非生物胁迫因素。自2000年以来，干旱的发生次数和持续时间较前二十年增加了29%。水约占植物鲜生物量的80%–95%，对植物生长发育和代谢活动至关重要。植物已进化出多种适应和耐受策略应对干旱，包括增加根系吸水、通过气孔关闭减少水分损失、以及积累渗透调节物质维持水分平衡。然而，干旱同时会通过降低蒸腾速率和改变膜转运蛋白活性，损害矿物质营养从根系向地上部的转运。理解干旱胁迫如何影响植物矿质营养，有助于制定缓解干旱损伤及营养缺乏的有效策略。
大麦（Hordeum vulgare L.）已成为分子生物学和遗传学研究中尤其是非生物胁迫耐受性的模式作物。野生二棱大麦（H. vulgare ssp. spontaneum）是栽培大麦（H. vulgare ssp. vulgare）的直接祖先，两个亚种完全生殖兼容，因此野生大麦一直被认为是将优良基因渗入现代品种以增强抗逆性的宝贵遗传资源。然而，现代栽培品种面临两个主要限制：一是长期高强度选择导致遗传背景极度狭窄，二是栽培大麦中真正耐旱种质有限。一个关键的科学问题是：驯化是否削弱了有效干旱适应所必需的关键调控通路？研究表明栽培大麦可能部分丧失了祖先种在水分亏缺下维持离子稳态的能力，尤其是通过液泡转运系统和膜转运蛋白调控高效分隔和区隔有毒离子的能力。此外，已有研究大多采用单一层次方法，仅提供干旱适应的部分视图。基于电感耦合等离子体质谱（ICP-MS）的离子组学可量化胁迫引起的元素组成变化，但无法解析调控养分转运、液泡分隔或膜转运蛋白活性的上游调控网络；而转录组学虽可识别基因表达变化，却不能表明这些调控变化是否转化为功能性营养稳态。因此，将离子组学与转录组学整合，可将元素表型与潜在调控网络联系起来，实现对干旱响应的更全面、机制性理解。
本研究针对上述不足，采用整合多组学方法，对野生和栽培大麦的早期干旱响应进行剖析，旨在揭示H. spontaneum优越耐旱性背后的调控网络，特别是营养稳态和胁迫相关转录重编程对干旱适应的贡献。研究人员假设：H. spontaneum优越的耐旱性源于水分亏缺条件下抗氧化防御、营养稳态、激素信号和胁迫响应转录网络之间的协同调控。
二、研究技术与方法
本研究选用两个耐旱性对比鲜明的大麦基因型：耐旱野生基因型H. spontaneum no. 11BP395（采集自伊朗西部Urmia湖区，37°11′34.0″ N, 45°21′43.4″ E，海拔1292 m）和来源于瑞典的早熟栽培品种‘Mona’（干旱敏感型）。在温室条件下进行对照和干旱胁迫处理（两叶期后干旱3周），于苗期采集叶片样本。采用ICP-MS测定叶片和成熟籽粒中元素浓度；利用Illumina NovaSeq 6000平台进行RNA-seq测序（8个文库，平均4900万条clean reads，双端151 bp），clean reads比对到‘Morex’参考基因组，用DESeq2鉴定差异表达基因（DEGs，FDR≤0.05且|log2FC|≥2）；采用基因型×处理互作模型解析基因型特异的干旱响应；使用g:Profiler进行GO富集分析，通过KAAS进行KEGG通路注释；最后用RT-qPCR验证5个候选基因（HvDhn1、HvbZIPx、HvBip-70、HvPht1-6和HvCBF6）的表达。
三、主要研究结果
（一）生长与生理性状
干旱胁迫对株高和鲜重产生不利影响，敏感型品种受害更重。野生基因型的相对含水量（RWC）降幅（38%）显著小于栽培基因型（63%），膜稳定性指数（MSI）在栽培大麦中下降58.52%，表明其膜损伤更严重。野生大麦的脯氨酸积累量大幅增加（达220%），并表现出更强的抗氧化响应，包括过氧化物酶（POX）、总酚含量（TPC）、总黄酮含量（TFC）和DPPH活性的显著升高。即使在对照条件下，栽培基因型也表现出更高的丙二醛（MDA）和过氧化氢（H?O?）水平，干旱下进一步急剧升高。激素分析显示，干旱下野生基因型的脱落酸（ABA）增加134%，而栽培基因型仅增加95%；吲哚-3-乙酸（IAA）在栽培基因型中下降46%，赤霉素（GB）下降44%，而野生基因型IAA略有增加（19%），GB无显著变化。
（二）ICP-MS元素分析
干旱胁迫在敏感和耐旱基因型间引发了对比鲜明的矿质元素分配模式。敏感基因型籽粒中铝积累增加，而耐候基因型籽粒和叶片中铝水平均降低。铜在敏感基因型中不变，但在耐候基因型籽粒中增加三倍（200%）。钙呈相反趋势：敏感基因型籽粒钙上升而耐候基因型稳定，叶片钙在敏感基因型中下降而在耐候基因型中上升。钾和磷在敏感基因型籽粒和叶片中均下降，而耐候基因型两者均增加。铁浓度在两种基因型中均增加，叶片中积累更多。锰在耐候基因型中显著增加，敏感基因型中不变。镁和锶在敏感基因型籽粒中下降，耐候基因型保持稳定。锌在敏感基因型籽粒和叶片中下降，耐候基因型籽粒锌增加。此外，钼、钠和硫在耐候基因型籽粒中增加。总体而言，耐旱基因型维持或增强了必需营养元素（Ca、K、P、Mn、Mg、Zn、Mo、S、Fe）的水平，同时限制有毒或非必需元素（Al、Ba、Sr）的积累。
（三）转录组分析
RNA-seq分析显示，野生大麦在干旱胁迫下有1576个基因上调、988个下调（共2564个DEGs）；‘Mona’品种有870个上调、1293个下调（共2163个DEGs）。609个DEGs为两基因型共有。耐旱基因型中83个DEGs的log2FC>8，而敏感型仅35个。基因型×处理互作分析共鉴定出955个耐旱基因型特异的干旱响应基因、554个敏感型特异的响应基因、39个在野生型中更强的共享响应基因、20个在‘Mona’中更强的共享响应基因、65个反向调控基因，以及1722个未达显著阈值但方向性不同的基因。
（四）GO与KEGG富集分析
耐旱野生大麦中上调基因显著富集于腺苷核糖核苷酸结合、碳水化合物衍生物结合、谷胱甘肽转移酶活性等分子功能，以及细胞呼吸调控、前体代谢物和能量生成调控、谷胱甘肽代谢过程、蛋白磷酸化、胁迫应答和防御反应等生物学过程。而敏感型品种上调基因主要富集于蛋白丝氨酸/苏氨酸磷酸酶活性、非生物刺激响应、热响应和温度刺激响应。KEGG分析显示，耐旱基因型上调DEGs富集于代谢通路、次生代谢产物生物合成、碳代谢、氨基酸生物合成、MAPK信号通路、植物激素信号转导、甘油磷脂代谢和苯丙烷生物合成等途径。
（五）耐旱候选基因
1. ABA非依赖性信号：多种钙结合蛋白基因（HvCML12、HvCML21、HvCML22、HvCML31）在耐旱基因型中特异且强烈上调；钙转运ATP酶在耐旱基因型中log2FC=5.9。蛋白激酶分析在耐旱野生大麦中发现98个上调PK基因，其中RLK-Pelle家族最大（82个）；3个CMGC家族MAPK基因仅在野生大麦中上调。
2. 植物激素生物合成：HvPYR/PYL在敏感型中下调、耐旱型中上调；HvPP2C在敏感型中显著高表达。耐候型中茉莉酸信号通路HvJAZ基因上调。
3. 转录因子：耐旱基因型中44个TF家族共189个基因差异表达（136上调、53下调），AP2/ERF、NAC、C2H2、bHLH、bZIP、MYB、TIFY等家族高表达。‘Mona’中125个TF差异表达（24上调、101下调）。AP2/ERF相关基因在野生耐候型中的表达增幅是敏感型的2倍以上。
4. 离子转运蛋白：HvABCG32和MDR-like ABC转运蛋白仅在耐候型中上调（后者log2FC>7.2）；MFS家族6个基因在耐候型中特异上调；硫酸盐转运蛋白HvSultr3;4在野生型中表达量高2.8倍；钾转运蛋白在耐候型中log2FC=5.7；Na?/Ca2?交换体（log2FC=10.34和4.01）、砷泵ATPase（log2FC=7.1）等仅耐候型上调。
5. ROS清除蛋白：黄酮生物合成途径多个酶基因（CYP75B1、CYP75A、HCT等）在野生型中显著上调，而敏感型中查尔酮合酶、花青素合酶等大幅下调（最大log2FC=-23）。17个谷胱甘肽S-转移酶基因和铁氧还蛋白HvFdx3在野生型中显著诱导。
6. 渗透稳态与相容性溶质：野生型中6个LEA蛋白基因被激活，敏感型仅2个上调、8个下调。Dhn1、Dhn7、HvP5cr3、HvP5CSB在两种基因型中均上调，Dhn10（log2FC=7.6）仅在耐候型中上调。sHSPs、HSP70、HSP90家族成员在干旱下显著高表达。敏感型中水通道蛋白显著下调（HvTIP1;1 log2FC=-9.8），而耐候型中无差异表达。
（六）RT-qPCR验证
RT-qPCR验证了5个候选基因（HvDhn1、HvbZIPx、HvBip-70、HvPht1-6和HvCBF6）的表达趋势，与RNA-seq数据高度一致。耐旱野生型中5个基因均上调，而敏感型中3个基因诱导极弱或（HvCBF6）下调。
四、讨论与结论
讨论部分强调了耐旱性是由调控、生理、抗氧化和激素机制组成的整合网络所决定。野生基因型在水分亏缺下能更有效地维持细胞稳态，其光合色素稳定性更强、抗氧化激活更强、酚类和黄酮积累更高。激素分析显示耐候型ABA信号整合更协调，IAA和GA维持更稳定。元素分析证实耐候型能维持必需元素并限制有毒元素积累。转录组分析显示耐候型表现更强的转录激活，与盐胁迫和热胁迫下伊朗野生大麦的转录模式一致，提示野生大麦具有保守的调控框架。基因型×处理互作揭示了两基因型间干旱响应机制的差异，野生型表现出更协调的转录调整。信号元件方面，Ca2?传感器、膜联蛋白、CDPK/SnRK和RLK-Pelle激酶在耐候型中强烈协同诱导。激素途径中，耐候型通过PYL/ABF/SnRK2高效激活ABA信号，而敏感型因PP2C升高和PYL降低而致ABA感知减弱。转录因子分析显示耐候型TF激活谱更广更强，包括ERF、MYB、WRKY、bHLH等家族。离子转运蛋白研究揭示ABC、MFS、硫酸盐转运蛋白、钾通道/转运蛋白和H?-ATPase在耐候型中的强诱导，支持离子稳态维持。ROS清除方面，耐候型通过GST、POD、ferredoxin以及黄酮合成途径的酶类实现高效ROS解毒；UGT基因的强诱导促进黄酮糖苷化稳定。渗透调节方面，LEA蛋白、脱水素、HSPs、脯氨酸合成酶（P5CS）以及水通道蛋白的差异调控共同维持细胞水分和蛋白质稳态。
研究结论指出：生长、生理和转录组数据一致表明野生大麦具有更强的干旱韧性，这是11BP395与‘Mona’基因型间差异的体现。野生基因型表现出更好的生长保持、叶绿素和类胡萝卜素稳定性、更低的氧化胁迫标志物以及关键胁迫响应基因的上调表达。其增强的分子响应涉及激素信号通路和蛋白激酶级联的激活、钙依赖网络、渗透保护剂产生增加、多酚生物合成以及水通道蛋白亚型的动态调控。这些机制支撑了耐旱性和生理性能的改善，转录组学见解为理解大麦适应性和基因型间性状转移提供了资源。鉴定出的干旱响应基因为开发功能性分子标记以高效筛选耐旱种质奠定了基础，通过杂交可将野生大麦的关键抗性位点渗入栽培背景，建立携带有利等位基因的双单倍体（DH）系将加速这些性状在育种中的固定和利用。
