《Frontiers in Plant Science》：Genome-wide analysis of the glutathione S-transferase family in Vaccinium corymbosum and functional characterization of VcGSTU53 in response to abiotic stress
编辑推荐：
本文聚焦于高价值小果作物——蓝莓，通过全基因组水平系统鉴定其谷胱甘肽S-转移酶（GST）基因家族，并对铝胁迫响应的关键成员VcGSTU53进行深入的功能表征。研究揭示了GST家族在植物应对铝、镉、干旱等多种非生物胁迫中的重要作用，为培育抗逆蓝莓新品种提供了重要的分子基础与候选基因资源。
蓝莓作为一种在全球范围内备受推崇的高价值小浆果，以其丰富的健康促进特性而闻名。谷胱甘肽S-转移酶（Glutathione S-transferases, GSTs）是一个庞大而复杂的酶家族，在类黄酮代谢、植物生长发育以及对非生物胁迫的响应中扮演着至关重要的角色。然而，人们对蓝莓中GST基因家族的了解相对较少。本研究的核心即是通过全基因组鉴定与生物信息学分析，系统解析蓝莓GST基因家族，并对筛选出的关键候选基因VcGSTU53在非生物胁迫下的响应机制进行深入的功能研究。
1 蓝莓GST基因家族的全面鉴定与生物信息学分析
研究人员利用高丛蓝莓品种‘Draper’的参考基因组序列，成功鉴定出190个VcGST基因。这些基因被划分为14个亚类，包括tau（U, 75个）、phi（F, 28个）、EF1G（17个）、lambda（L, 10个）、MAPEG（10个）、zeta（Z, 8个）、GHR（8个）、Metaxin（M, 7个）、GST2N（6个）、DHAR（5个）、hemerythrin（H, 4个）、mPGES2（4个）、TCHQD（4个）和theta（T, 4个）；值得注意的是，未发现iota（I）类GST。VcGST基因的长度从366 bp到4386 bp不等，编码121至1461个氨基酸的多肽，预测分子量在13.86 kDa到163.02 kDa之间，理论等电点范围为4.77至9.56。亚细胞定位预测显示，约46.32%的VcGST蛋白定位于细胞质，其余则分布于内膜系统、叶绿体、细胞核、线粒体、细胞外空间和质膜。
2 系统进化、染色体分布与复制事件
构建蓝莓、拟南芥、葡萄和甘蓝GST蛋白的系统进化树显示，不同物种间相同GST亚家族的蛋白高度保守地聚类在一起，其中蓝莓GST与葡萄亲缘关系最近。190个VcGST基因不均匀地分布在蓝莓的48条染色体上，其中3号染色体包含的基因数量最多（13个），而32号和48号染色体上未检测到GST基因。复制事件分析表明，全基因组复制（WGD）或片段复制是驱动VcGST基因家族扩张的主要力量，共涉及150个基因。此外，还鉴定出27个串联复制基因和12个近端复制基因。非同义替换率（Ka）与同义替换率（Ks）的比值分析显示，大多数复制基因对的Ka/Ks值小于1，表明它们主要受到纯化选择；仅有五个基因对（如VcGSTL7/VcGSTL10）的Ka/Ks值大于1，暗示它们可能经历了正选择。
3 保守基序、结构域与基因结构
大多数VcGST蛋白包含两个经典的结构域：GST_N和GST_C。然而，GHR、MAPEG和GST2N等亚家族成员仅具有其中一个完整结构域，或拥有独特的结构域组成。保守基序的组成在不同亚家族间存在差异，例如motif 14为tau类GST所特有，而motif 8为phi类所特有。外显子-内含子结构分析揭示了三种主要类型：同一组内基因外显子数目相同或相近（I型）；拥有完整保守结构域的基因外显子数相同，而具有多个重复结构域的基因外显子数更多（II型）；以及同一组内基因外显子数变异较大（III型）。例如，tau和phi类GST通常只有2-3个外显子，而zeta类基因（如VcGSTZ5）的外显子数最多。
4 启动子顺式作用元件分析
对VcGST基因启动子区（转录起始位点上游2000 bp）的分析，鉴定出大量顺式作用元件，主要可分为四大类：光响应元件（798个）、胁迫响应元件（432个，涉及氧化胁迫、干旱、低温）、激素响应元件（245个，涉及茉莉酸甲酯、脱落酸、赤霉素、生长素）以及与生长发育相关的元件（157个）。值得注意的是，胁迫相关元件在tau和phi类VcGST基因的启动子中分布更为频繁。此外，多个VcGST启动子中含有MYB和MYC转录因子的结合位点，其中一些MYB结合基序与花青素生物合成相关，暗示GST可能参与蓝莓花青素的运输与储存过程。
5 VcGST基因在铝胁迫下的表达谱
为鉴定响应铝胁迫的关键VcGST基因，研究对两个蓝莓品种（‘Duke’和‘Draper’）经不同浓度AlCl3处理的30个根样本进行了转录组测序。差异表达分析共鉴定出51个差异表达的VcGST基因，其中tau类（27个）和phi类（8个）占主导地位。在所有差异表达基因中，VcGSTU53（VaccDscaff35-augustus-gene-220.36）的表达量最高，且在‘Duke’品种中随铝胁迫处理呈现持续上调。该基因被注释为编码一种参与氧化胁迫响应的抗氧化酶，因此被选作后续功能研究的核心候选基因。
6 VcGSTU53的分子特征与亚细胞定位
VcGSTU53的编码序列长693 bp，编码一个由230个氨基酸组成的蛋白。保守结构域分析显示其包含特征性的GST-N-Tau和GST-C-Tau结构域，并预测有7个谷胱甘肽结合位点。系统进化分析表明，VcGSTU53与山茶和猕猴桃中的同源蛋白亲缘关系最近。亚细胞定位实验通过在烟草表皮细胞中瞬时表达VcGSTU53-GFP融合蛋白证实，该蛋白定位于细胞质，与生物信息学预测结果一致。
7 VcGSTU53的原核表达、酶活性与胁迫耐受性
将VcGSTU53在大肠杆菌中异源表达并纯化，获得了带有His标签的重组蛋白。体外酶活实验证实，VcGSTU53能以谷胱甘肽和1-氯-2,4-二硝基苯为底物，表现出典型的转移酶活性，但对脱氢抗坏血酸无活性。为评估VcGSTU53赋予的胁迫耐受性，研究人员比较了表达VcGSTU53的重组菌株与空载体对照菌株在不同胁迫条件下的生长曲线。结果显示，在正常条件下，两者生长无差异。但在铝胁迫下，VcGSTU53能显著增强大肠杆菌的耐受性，特别是在4 mM AlCl3的高浓度下，只有表达VcGSTU53的菌株能够存活生长。此外，VcGSTU53的表达也提高了大肠杆菌对干旱（PEG8000模拟）和镉胁迫的耐受性，但对高盐（500 mM NaCl）胁迫的耐受性没有明显改善。
结论
本研究在蓝莓中全基因组鉴定了190个GST基因，并对其进行了系统的生物信息学分析，揭示了WGD/片段复制是该家族扩张的主要驱动力。通过表达谱分析筛选出铝胁迫响应的关键基因VcGSTU53，并证实其编码一个定位于细胞质、具有谷胱甘肽转移酶活性的蛋白。异源表达实验证明VcGSTU53能显著增强宿主对铝、镉和干旱等多种非生物胁迫的耐受性。这些发现深化了我们对植物GST家族在非生物胁迫下功能的理解，对于利用分子育种手段培育抗逆蓝莓新品种具有重要的意义。
