藏红花谷胱甘肽S-转移酶（GSTs）的全基因组鉴定及激素介导的应激反应机制研究
引言
谷胱甘肽S-转移酶（Glutathione S-transferases, GSTs; EC 2.5.1.18）是广泛存在于植物中的多功能酶，在胁迫管理、有害化学物质解毒以及次级代谢物转运中发挥关键作用。本研究旨在对具有重要经济价值的园艺作物藏红花（Crocus sativus）进行GST基因家族的全基因组分析，以填补该领域的研究空白，并为理解其在胁迫耐受性和有价值化合物（如藏花素和藏花醛）合成中的作用奠定基础。
材料与方法
通过BLASTP搜索藏红花基因组数据，结合Pfam和SMART数据库的保守结构域验证，共鉴定出29个非冗余的CsGST基因。利用ClustalW进行多序列比对，并采用RAxML构建系统发育树，将CsGSTs划分为Phi、Tau、Theta和Lambda亚家族。MEME套件用于基序分析，Pfam数据库用于结构域验证。通过Ka/Ks分析评估选择压力，利用TBtools进行染色体分布和共线性分析。WolfPsort预测亚细胞定位，ExPASy ProtParam分析理化性质。使用PlantCARE数据库分析启动子区顺式作用元件。基于公共RNA-seq数据集（PRJNA400472, PRJNA976833）分析基因表达模式，并利用eggNOG-mapper进行GO富集分析。通过SWISS-MODEL进行同源建模以预测蛋白质三维结构。
结果
藏红花GST基因的发现
全基因组分析在藏红花中鉴定出29个GST基因（CsGST1-CsGST29），这些基因分布在除chr3之外的其余8条染色体上，并在chr1、chr5和chr6上呈现簇状分布。所有基因均含有GST特异的保守结构域，如N端硫氧还蛋白样结构域（GST_N）和C端α-螺旋结构域（GST_C）。
GST基因的系统发育分析
系统发育分析将CsGSTs主要分为四个亚家族：Phi类（13个成员，如CsGST1-3, 6-9等）、Tau类（10个成员，如CsGST4-5, 10-11等）、Theta类（3个成员，CsGST19-20, 27）和Lambda类（3个成员，CsGST8, 16-17）。Phi和Tau类是植物中最为丰富的亚家族，与它们在胁迫响应和解毒中的重要作用相符。
基序和结构域分析保守性
基序分析揭示了CsGST蛋白中存在10个保守基序，其分布模式与系统发育分类高度一致。结构域分析确认了所有CsGST蛋白均包含典型的GST_N和GST_C结构域，但也发现了一些变异形式（如GST_N_2, GST_C_2等），反映了该基因家族的结构多样性。
GST基因序列的Ka/Ks分析
对CsGST基因对的Ka/Ks分析显示，所有比值均小于1（范围0.0596-0.4868），表明该基因家族主要受到纯化选择的作用。例如，CsGST2-CsGST3对的Ka/Ks值最低（0.0596），表明其功能受到严格约束；而CsGST12-CsGST18对的比值相对较高（0.4868），暗示可能存在一定程度的功能分化。
GST基因的染色体分布与共线性分析
29个CsGST基因不均匀地分布在8条染色体上，chr1、chr5和chr6是主要富集区域。共线性分析检测到12个片段复制事件，而未发现串联复制事件，表明片段复制是CsGST基因家族扩张的主要驱动力。与番茄（Solanum lycopersicum）的种间共线性分析发现了14个保守的直系同源基因对，表明GST基因家族在单子叶和双子叶植物中具有进化保守性。
GST蛋白的亚细胞定位预测
WolfPsort预测显示，大多数CsGSTs定位于细胞质（如CsGST1, 11, 24, 25），部分定位于叶绿体（如CsGST7, 16）、细胞核（如CsGST17, 18）或液泡，这种多样化的亚细胞定位暗示了其在不同细胞区室中可能具有多种功能。
GST蛋白的理化特性分析
理化性质分析表明，CsGST蛋白的分子量范围为20.49至68.49 kDa，等电点（pI）在5.33至7.74之间，表明多为中性至弱酸性蛋白。GRAVY值在-0.281至0.143之间，提示大多数CsGSTs为亲水性蛋白，有利于其在细胞水相环境中发挥功能。
启动子区激素响应元件的分布
对29个CsGST基因启动子区（转录起始位点上游2000 bp）的分析共鉴定出844个顺式作用元件。其中，激素响应元件尤为丰富，包括236个茉莉酸甲酯（MeJA）响应元件、120个脱落酸（ABA）响应元件（ABRE）和34个赤霉素（GA）响应元件，此外还有水杨酸、生长素和乙烯响应元件。这些发现强烈暗示CsGSTs的表达受到多种植物激素的精密调控。
CsGST蛋白的结构与相互作用分析
对代表性CsGST蛋白（如CsGST1, 2, 3, 7, 9, 14）的同源建模获得了高质量的三维结构模型（QMEAN得分介于-0.93至-0.86）。模型显示其具有典型的GST结构：N端结构域（残基10-91）形成谷胱甘肽结合位点（G-site），C端结构域（残基96-221）形成疏水的底物结合位点（H-site）。蛋白质-蛋白质相互作用（PPI）网络分析显示CsGST3是一个关键枢纽，与多个其他CsGSTs（如CsGST11, 17, 20, 29）存在相互作用，提示其在协调胁迫响应通路中可能起核心作用。
藏红花中GST基因表达的激素介导调控
RNA-seq表达分析揭示了CsGSTs在柱头发育、组织特异性胁迫响应及芽休眠过程中存在差异表达，且其表达受到MeJA、ABA和GA的调控。例如，在柱头成熟过程中（Stigma_RED vs. Stigma_0_DAY），CsGST1、CsGST13和CsGST23显著上调表达。在组织特异性比较中（Apical_Bud_Nov vs. Leaf_Nov），CsGST9、CsGST21、CsGST24-27在顶芽中显著上调，而CsGST13和CsGST14则显著下调。GO富集分析进一步证实了差异表达的CsGSTs基因在“细胞对胁迫的反应”、“对脱落酸的反应”和“赤霉素介导的信号通路”等条目中显著富集。
讨论
本研究首次对藏红花GST基因家族进行了全面的基因组学分析。鉴定出的29个CsGSTs在进化上主要受到纯化选择，其家族扩张主要通过片段复制实现。系统发育分析将其分为四个主要亚家族，其中Phi和Tau类占主导地位，这与它们在植物胁迫响应中的核心作用一致。表达分析和启动子元件预测共同表明，MeJA、ABA和GA等激素在调控CsGSTs表达以应对生物和非生物胁迫中扮演重要角色。特别是在柱头发育和芽休眠等关键生理过程中，特定的CsGSTs（如CsGST9, 16, 23, 27）表现出显著的表达变化，暗示它们可能在保护有价值的次级代谢物（如藏花素）合成、以及帮助植株适应季节性环境变化方面具有特殊功能。蛋白质结构建模和互作网络分析为理解其分子功能提供了结构基础。
结论
本研究成功地在藏红花基因组中鉴定并注释了29个GST基因，系统分析了它们的系统发育关系、基因结构、染色体分布、进化选择压力、顺式调控元件、表达模式以及蛋白质结构特征。研究结果揭示了CsGST基因家族在藏红花激素介导的胁迫响应中的潜在重要作用，为了解这一珍贵作物应对环境挑战的分子机制提供了宝贵的见解，并为未来通过遗传手段提高藏红花的抗逆性和品质奠定了重要的理论基础。后续工作应侧重于通过功能基因组学实验（如基因敲除、过表达）验证关键CsGST基因的具体功能。
打赏
