编辑推荐：
作为植物E3泛素连接酶关键家族，U-box基因在盐胁迫应答中发挥核心作用。本文首次对毛果杨103个U-box基因（PtrPUBs）开展全基因组鉴定，通过系统发育、启动子顺式元件（cis-acting elements）及表达谱分析，揭示其通过蛋白泛素化调控盐胁迫响应的分子机制，为林木抗逆育种提供新靶点。
引言：U-box基因家族在植物盐胁迫应答中的核心地位
植物长期面临盐胁迫等非生物胁迫，高盐环境会导致细胞脱水、膜损伤及离子稳态失衡。盐胁迫通过干扰Na+/K+平衡、诱导氧化应激和代谢紊乱抑制植物生长。为应对胁迫，植物进化出复杂调控网络，其中蛋白泛素化作为关键翻译后修饰机制，由E3泛素连接酶介导底物特异性降解。U-box基因家族编码的E3连接酶（如PUB25/26调控拟南芥花器官发育、OsPUB33促进水稻光合作用）在植物发育和胁迫应答中发挥重要作用。已有研究表明，毛果杨PalPUB79通过泛素化降解PalWRKY77增强ABA信号通路介导的耐旱性，而PUB30通过调控HB24降解抑制SWEET11表达，限制根系生长以应对盐胁迫。然而，毛果杨U-box基因家族的系统鉴定及其在盐胁迫中的调控网络仍待阐明。
材料与方法：多维度解析U-box基因家族特征
研究采用毛果杨v3.1基因组数据，通过HMMER（E值<1×10-5）筛选含U-box结构域（PF04564）的蛋白，经SMART和Pfam验证后获得103个PtrPUBs。理化性质分析显示蛋白分子量9.5-129.3 kDa，等电点4.31-9.08，87个蛋白定位于细胞核。系统发育分析（MEGA11邻接法）将基因分为6个亚群（Group I-IV），其中Group IV（37个成员）和Group VI（26个成员）规模最大。基因结构分析揭示36个基因无内含子，Group V内含子数量最多，而Group VI结构最精简。启动子分析（PlantCARE）发现2000 bp区域富含ABA、MeJA激素响应元件及MYB结合位点等胁迫相关顺式元件。染色体定位显示基因分布不均，第2染色体含13个基因最多，第18染色体无分布。复制事件分析识别4个串联重复和30个节段重复，表明节段重复是家族扩张主要机制。物种间共线性分析显示毛果杨与桉树、拟南芥等同源基因数量显著多于单子叶植物。
结果：U-box基因家族的结构与表达特征
- 1.
系统发育与结构特征：Group VI成员内含子数量少（利于快速表达），Group V结构复杂（支持功能多样性）。Motif 1为保守U-box结构域，其余9个基序呈亚群特异性分布。
- 2.
启动子调控元件：鉴定到干旱响应MYB位点、低温响应元件及ABA/MeJA激素响应元件，表明基因参与胁迫应答与发育调控。
- 3.
染色体分布与进化：基因密度热图显示第2、5、10染色体为富集区域。共线性分析揭示Group IV和VI为扩张核心，与桉树、拟南芥同源基因数量最多。
- 4.
组织特异性表达：热图聚类显示5类表达模式，如Cluster A（叶特异）、Cluster D（根特异），表明基因功能分化。
- 5.
盐胁迫响应：根中40个差异表达基因（23个上调，17个下调），叶中21个（12个上调，9个下调），14个基因在根叶共表达。qRT-PCR验证9个关键基因（如Potri.006G202600.1、Potri.009G165900.1）在150 mM NaCl处理12/24小时显著表达。
讨论：U-box基因家族在盐胁迫中的调控网络
- 1.
结构功能关联：Group VI精简结构利于快速响应胁迫，Group V复杂结构支持调控多样性。胞外定位的U-box蛋白（如Potri.005G148900.1）暗示其参与细胞外信号传导。
- 2.
进化驱动机制：节段重复（30对）是家族扩张主要方式，与桉树等高共线性反映双子叶植物间保守性。
- 3.
顺式元件调控：MYB结合位点、ABA响应元件富集表明基因通过激素信号通路调控胁迫应答。
- 4.
组织特异性功能：根叶差异表达模式（如Cluster D根特异表达基因）反映组织适应性调控。
- 5.
胁迫响应机制：同源基因功能注释显示，Potri.010G103100.1（同源AT3G52450）参与水胁迫响应，Potri.006G202600.1（同源AT2G35930）负调控PTI免疫，PUB18同源基因通过抑制盐胁迫响应调控耐盐性。
结论：U-box基因家族是毛果杨耐盐育种的关键资源
本研究首次系统鉴定了毛果杨103个U-box基因，揭示其通过结构多样性、顺式元件调控及组织特异性表达参与盐胁迫应答。40个根差异表达基因和21个叶差异表达基因的发现，为解析林木耐盐分子机制提供新视角。未来可通过基因编辑技术验证关键PtrPUBs功能，推动耐盐杨树新品种选育。
