全基因组解析大麦属耐盐碱新靶点：短茎大麦WRKY II类基因家族鉴定与HbWRKY11功能验证

《Plants》：Genome-Wide Identification of WRKY Group II Genes and the Role of HbWRKY11 in Hordeum brevisubulatum Under Saline-Alkali Stress Sihan Chen, Yicheng Yin, Guangyao Qi, Yanda Li, Bello Hassan Jakada, Dan Sun, Xinying Liu and Xingguo Lan

本研究针对日益严重的土壤盐碱化问题，聚焦在耐盐碱的野生大麦——短茎大麦(Hordeum brevisubulatum)中，首次在全基因组水平系统鉴定了WRKY II类转录因子基因家族，并深入探究了其在盐碱胁迫下的响应机制。研究人员通过比较基因组学、转录组学（RNA-seq）、实时荧光定量PCR（qRT-PCR）及酵母异源表达等多种技术，发现多个HbWRKY基因在NaHCO3胁迫下表达上调，其中HbWRKY11表现尤为突出，能显著增强酵母对NaHCO3胁迫的耐受性。这项工作不仅揭示了短茎大麦适应盐碱环境的潜在分子基础，也为禾谷类作物的耐盐碱遗传改良提供了关键候选基因。

　　

土壤盐碱化是威胁全球农业生产的重大环境挑战之一。过高的盐分和碱性条件会破坏植物的水分吸收、离子平衡和代谢过程，最终导致生长抑制和产量下降。据估计，全球超过7%的土地受到盐胁迫影响，而且由于气候变化和不合理的土地利用，受影响的面积还在不断扩大。因此，寻找能够在盐碱地顽强生存的植物“勇士”，并挖掘它们体内赋予其“抗压”能力的“王牌基因”，对于退化土地的生态修复和农作物的遗传改良都具有至关重要的意义。

在禾本科植物家族中，大麦属（Hordeum）为我们提供了一个绝佳的研究“素材库”。这里不仅有我们熟悉的栽培大麦（Hordeum vulgare），还有一些“身怀绝技”的野生亲戚，它们通常蕴含着丰富的、可耐受非生物胁迫的遗传多样性。短茎大麦（Hordeum brevisubulatum）就是其中一位佼佼者，它是一种多年生草本植物，自然分布于中亚北部的盐碱草原。它展现出了对盐碱胁迫的非凡耐受性，能在足以杀死大多数作物的极端土壤条件下生存和繁衍。那么，这位“盐碱地勇士”体内到底藏着什么样的“秘密武器”，帮助它应对如此严酷的环境呢？

科学家们将目光投向了植物体内一类重要的“基因开关”——WRKY转录因子。这是一大类植物特有的转录因子，它们在调控植物生长发育、激素信号传导、病原防御以及对非生物胁迫（如干旱、盐碱、极端温度等）的响应中发挥着核心作用。其中，WRKY II类亚家族又是成员最多、功能最多样化的一个分支，已被报道在多种植物中参与增强耐盐性。然而，在短茎大麦这种具有突出耐盐碱能力的野生植物中，WRKY II类基因家族的全貌及其在盐碱胁迫响应中的具体角色，仍然是一个未被深入探索的“黑匣子”。

为了揭开这个“黑匣子”的秘密，陈思涵（Sihan Chen）等研究人员在《Plants》杂志上发表了一项系统性研究。他们利用最新发布的短茎大麦染色体级别基因组，首次在全基因组水平上鉴定和分析了其WRKY II类基因家族，并重点关注了它们在应对盐碱胁迫中的作用，成功锁定了一个关键候选基因。

为开展此项研究，作者主要运用了以下关键技术方法：首先，利用比较基因组学（BLASTp比对、系统进化分析）在短茎大麦、海大麦、栽培大麦和水稻中鉴定WRKY II类基因，并分析其染色体分布、基因结构和保守基序。其次，对基因上游启动子区域进行顺式作用元件分析，预测其潜在功能。再次，通过生物信息学方法（MCScanX）进行基因共线性和Ka/Ks（非同义替换率/同义替换率）分析，探讨基因家族的进化历史。此外，通过RNA-seq（转录组测序）和qRT-PCR（实时荧光定量PCR）技术，分析了短茎大麦幼苗在200 mM NaHCO₃胁迫下不同时间点（0、3、6、12小时）根和叶中WRKY II类基因的表达模式。最后，利用酵母（Saccharomyces cerevisiae）异源表达系统，对关键候选基因HbWRKY11在NaHCO₃胁迫下的功能进行了初步验证。研究所用的短茎大麦种子采集自中国黑龙江省安达市的野生居群。

研究结果

2.1. 系统进化与理化性质分析

研究人员在短茎大麦中鉴定出23个WRKY II类基因，并将其与海大麦、栽培大麦、水稻的成员一起构建了系统进化树，将其分为五个亚组（IIa–IIe）。其中，IIc亚组成员数量最多，表现出较高的物种间进化保守性。理化性质分析表明，这些蛋白质在长度、分子量和等电点（pI）上存在显著差异，但总体呈亲水性。

2.2. 染色体分布

短茎大麦的23个WRKY II类基因不均匀地分布在6条染色体上，其中3号染色体上分布最多（9个）。比较分析显示，大麦属内WRKY II类基因的染色体框架基本保守，但也存在物种特异性的基因数目和定位差异。

2.3. 基因结构与基序分析

同一进化亚组的成员通常具有高度保守的基序组成和基因结构。例如，IIc亚组的基因通常含有较少的基序和外显子，结构紧凑，而IIa/IIb亚组成员的基序和外显子结构则更为复杂。这种结构上的分化可能是其功能多样性的分子基础。

2.4. 顺式作用元件分析

对启动子区的分析发现，HbWRKY基因的启动子中含有丰富的顺式作用元件，可分为四大类：与胁迫和防御响应、植物激素响应、生长发育以及光响应相关的元件。其中，与胁迫、茉莉酸甲酯（MeJA）和脱落酸（ABA）响应相关的元件广泛存在。值得注意的是，仅在HbWRKY11的启动子中发现了与昼夜节律相关的调控元件，暗示了其时间依赖性转录调控的潜力。

2.5. 共线性与Ka/Ks分析

在短茎大麦基因组内检测到两对由片段重复产生的WRKY II类基因，未检测到串联重复，表明该基因亚家族的扩张主要由片段重复驱动。与栽培大麦和海大麦的共线性分析分别鉴定出23对和21对直系同源基因。所有直系同源基因对的Ka/Ks值均小于1，表明WRKY II类基因在进化过程中主要受到纯化选择，功能较为保守。

2.6. 盐碱胁迫下基因表达模式

转录组（RNA-seq）分析显示，在NaHCO₃胁迫下，多个HbWRKY基因在根和叶中表达上调。通过多时间点交集筛选，鉴定出8个在根和叶中均持续上调的候选基因，包括HbWRKY11、HbWRKY14、HbWRKY43、HbWRKY1等。其中，HbWRKY11、HbWRKY14、HbWRKY43和HbWRKY1的表达量随胁迫时间延长而逐渐增加。

2.7. HbWRKYs在非生物胁迫下的表达模式

qRT-PCR进一步验证了四个候选基因的表达模式。结果显示，这些基因的胁迫诱导表达主要发生在根部，且具有不同的时间动态。HbWRKY11在根部的诱导表达幅度最大，且表现为持续增强的模式，凸显了其在长时间胁迫下的强烈、持续性响应。

2.8. HbWRKY11赋予耐盐碱胁迫的潜力

酵母异源表达功能验证实验表明，在正常条件下，表达HbWRKY11的酵母与对照生长无差异。但在含有8 mM NaHCO₃的胁迫条件下，表达HbWRKY11的酵母菌株的生长状况（斑点实验和OD₆₀₀测定）显著优于携带空载体的对照菌株。这表明HbWRKY11的表达能够增强酵母对NaHCO₃诱导的盐碱胁迫的耐受性。

研究结论与重要意义

综上所述，本研究首次在耐盐碱野生大麦——短茎大麦中，对WRKY II类转录因子家族进行了全基因组鉴定与系统性表征。研究鉴定出23个WRKY II类基因，并揭示了它们在系统进化、染色体分布、基因结构和启动子调控元件上的特征。最重要的是，通过整合转录组、qRT-PCR和酵母功能验证，研究人员发现多个HbWRKY基因响应盐碱胁迫，并成功将HbWRKY11锁定为一个关键的正向调控因子。该基因在根部对NaHCO₃胁迫表现出持续且强烈的诱导表达，其启动子中含有独特的昼夜节律相关元件，其在水稻中的直系同源基因OsWRKY11也被报道参与多种胁迫响应。酵母实验直接证实了HbWRKY11能够提高对盐碱胁迫的耐受性。

这项研究的意义深远。首先，它从耐盐碱植物资源中挖掘关键基因的策略，为揭示植物适应逆境的分子机制提供了新见解。其次，鉴定出的HbWRKY11等候选基因为大麦乃至其他禾谷类作物的耐盐碱遗传改良提供了宝贵的基因资源。最后，这项工作展示了如何利用比较基因组学、转录组学和简单的异源表达系统进行高效的功能基因挖掘，为其他非模式植物的功能基因组学研究提供了范例。未来，在植物体内进一步解析HbWRKY11的下游调控网络及其在离子稳态和胁迫信号通路中的具体作用，将有助于我们更全面地理解其功能，并推动其在作物育种中的应用。