《Cells》：Genome-Wide Identification and Expression Analysis of the Heavy-Metal ATPase (HMA) Gene Family Reveal Their Correlation with Copper and Cadmium Stress Adaptation in Stevia rebaudiana
编辑推荐：
中文摘要翻译：重金属ATP酶（HMA）广泛参与植物对铜、锌和镉等金属元素的吸收与转运。尽管HMA已在多种植物中得到表征，但甜叶菊（Stevia rebaudiana）中HMA基因的研究仍然有限。本研究鉴定出12个SrHMA基因，其编码的蛋白质大多定位于质膜，具
中文摘要翻译：重金属ATP酶（HMA）广泛参与植物对铜、锌和镉等金属元素的吸收与转运。尽管HMA已在多种植物中得到表征，但甜叶菊（Stevia rebaudiana）中HMA基因的研究仍然有限。本研究鉴定出12个SrHMA基因，其编码的蛋白质大多定位于质膜，具有保守的E1-E2_ATPase和水解酶（hydrolase）结构域，并含有1至3个N端HMA结构域。系统发育分析将其分为两个亚家族（Zn/Co/Cd/Pb和Cu/Ag P1B-ATPase）。这些基因不均匀地分布在5条染色体上，存在两次片段重复事件和一次双子叶植物特有的基因组重复事件，其启动子区域富含多种胁迫响应元件。在铜和镉胁迫下的表达谱显示出明显的组织特异性：镉胁迫在根中上调了大多数基因的表达，但在叶片中则抑制其表达；而铜胁迫在叶片中上调了所有基因的表达，在根中上调了半数基因的表达。本研究是对SrHMA基因家族的首次系统分析，揭示了其在重金属胁迫响应中的潜在作用，并为未来功能鉴定以及培育具有增强金属耐受性的甜叶菊品种提供了候选基因。
论文解读文章
一、研究背景与科学问题
随着快速工业化、集约化农业和城市扩张，土壤重金属污染问题日益严重，对经济作物的品质与安全构成持续威胁。更为严重的是，植物来源的重金属会沿食物链发生营养级传递和生物富集，最终危及生态安全和公众健康。在众多重金属中，铜（Cu）和镉（Cd）尤为值得关注：铜在过量浓度下对植物具有植物毒性，而镉即使在低浓度下也具有高度毒性，且二者在环境中广泛分布。过量的Cu2?和Cd2?作为高强度胁迫因子，通过产生活性氧（ROS）爆发和破坏离子稳态等共同机制引发系统性生理失衡。在细胞层面，根尖分生组织细胞遭受直接生长抑制，叶肉细胞则发生光合系统崩溃；在细胞器层面，叶绿体和线粒体因其高电子传递活性而优先成为氧化损伤的靶标，导致结构降解。因此，探究植物重金属胁迫耐受机制对提高其在逆境环境中的适应能力具有重要意义。
为应对重金属毒害，植物进化出多种调控策略以防止组织中重金属过量积累，包括限制吸收、细胞外排、区室化、螯合和解毒等。这些过程大多依赖多种跨膜转运蛋白来动态维持植物体内不同金属离子的稳态。其中，重金属ATP酶（HMA），又称P1B-ATPase，是P型ATP酶大家族的一员，通过水解ATP提供能量，调控包括Cu2?和Cd2?在内的多种重金属离子跨膜吸收与转运，在植物重金属稳态维持中发挥核心作用。尽管HMA基因家族已在拟南芥（Arabidopsis thaliana）、大麦、玉米、水稻、大豆、毛白杨和绣球等多种植物中完成了全基因组水平的系统鉴定，但在重要经济作物甜叶菊中，HMA成员的鉴定及其对铜和镉胁迫的响应机制仍属空白。
甜叶菊（Stevia rebaudiana Bertoni）是菊科多年生草本植物，其活性成分甜菊糖苷作为一种天然、安全、低热量、高甜度的甜味剂，已被广泛应用于全球食品和制药工业，并获欧洲食品安全局（EFSA）和国际食品法典委员会等国际监管机构认证。然而，甜叶菊主产区多位于干旱和半干旱地区，降水少、蒸发强，重金属元素不易淋溶和迁移，易在土壤中累积，导致重金属背景值较高。铜和镉作为主要污染物，已被证明显著影响甜叶菊的生理代谢和生长。因此，在甜叶菊中开展HMA基因家族的鉴定与功能研究，对于保障原料安全、实现可持续生产至关重要。
二、研究内容与主要结论
本研究首次在甜叶菊全基因组水平上系统鉴定了HMA基因家族成员，并对其系统发育关系、保守结构域、基因结构、染色体分布、共线性关系、启动子顺式作用元件、蛋白互作网络、miRNA调控网络及表达特征进行了全面分析，重点筛选了响应铜和镉胁迫的关键候选基因。研究共鉴定出12个SrHMA基因，这些基因编码的蛋白质长度为763至1149个氨基酸，分子量介于82,105.74至127,330.27 Da之间，等电点范围为5.23至7.53，多数蛋白质为酸性且预测稳定。亚细胞定位预测显示，除SrHMA3（叶绿体）和SrHMA9（液泡膜）外，大多数蛋白质定位于质膜。系统发育分析将12个SrHMA基因分为两个亚家族：3个成员（SrHMA1、SrHMA6、SrHMA12）属于Zn/Co/Cd/Pb P1B-ATPase亚群，其余9个属于Cu/Ag P1B-ATPase亚群，这一分类与拟南芥和水稻一致，表明HMA家族在双子叶和单子叶植物分化之前就已存在。所有SrHMA蛋白均含有E1-E2_ATPase结构域、水解酶结构域和1至3个N端HMA结构域。基因结构分析显示，SrHMA基因的外显子数目为6至17个不等。染色体定位分析表明，12个SrHMA基因不均匀地分布在5条染色体上，其中3号染色体包含最多成员（4个）。种内共线性分析鉴定出两次片段重复事件（SrHMA11-SrHMA7和SrHMA4-Streb.6G009190.1），Ka/Ks值均小于1，表明SrHMA基因在进化过程中受到纯化选择。种间共线性分析显示，SrHMA基因与双子叶植物（向日葵、拟南芥、大豆）存在共线性关系，但与单子叶植物（水稻、玉米）无共线性，提示HMA基因家族在双子叶植物分化后经历了独立的基因组重复事件。启动子分析发现，所有SrHMA基因的启动子区域均富含光响应元件、激素响应元件和逆境胁迫响应元件。蛋白互作预测显示，6个SrHMA蛋白（SrHMA3、SrHMA4、SrHMA5、SrHMA6、SrHMA11和SrHMA12）均预测与铜伴侣蛋白CCS互作。miRNA调控网络预测包含99个miRNA和12个SrHMA基因，共140个互作关系，其中SrHMA4受17个miRNA靶向调控，表明其表达可能受到更严格的后转录调控。
基于公共转录组数据的组织表达模式分析显示，SrHMA基因大多在幼苗的根、茎和叶中高表达，但不同基因在不同组织和发育阶段表现出差异化的表达模式。SrHMA2、SrHMA3、SrHMA4和SrHMA6在所有组织中均表现出相对较高的表达水平，提示其可能参与甜叶菊基础金属离子转运过程；SrHMA11和SrHMA4在幼苗期叶片中特异性高表达；其余成员在幼苗期根中表达最高。
在镉胁迫（100 μM CdCl?）和铜胁迫（721 μM CuSO?）下的qRT-PCR表达分析揭示了该基因家族显著的組織特异性和时间动态。在镉胁迫下，叶片中大多数基因在1小时内出现早期快速诱导后迅速下降，24小时后几乎所有基因的表达水平均低于对照，表明长期镉胁迫对叶片中SrHMA基因整体呈现抑制作用；根中SrHMA基因响应更为强烈，其中SrHMA7和SrHMA8在3小时达到峰值（约为对照的15倍和12倍），且呈现多峰表达模式，SrHMA2则呈现独特的"U型"表达模式。在铜胁迫下，叶片中SrHMA5持续上调，从1小时的2倍上升至12小时的11倍并维持在高水平；根中SrHMA7和SrHMA8被强烈诱导，分别在24小时和12小时达到25倍和14倍的峰值表达水平。综合来看，SrHMA7和SrHMA8在铜和镉胁迫下均表现出强烈且特异的诱导表达，提示其可能是甜叶菊根部重金属解毒的核心成员。
三、主要技术方法
研究人员利用拟南芥（A. thaliana）和水稻（O. sativa）的HMA蛋白序列作为查询序列，通过BLASTP同源比对和HMM（隐马尔可夫模型）搜索策略，结合Pfam和CDD数据库对保守结构域进行验证，在甜叶菊全基因组中鉴定HMA基因家族成员。随后运用MEGA11软件进行多序列比对并构建邻接法（NJ）系统发育树；利用MEME工具预测保守基序；通过MCScanX进行基因重复事件和共线性分析；使用PlantCARE数据库预测启动子顺式作用元件；借助STRING数据库预测蛋白互作网络，并通过psRNATarget预测miRNA靶向关系。在表达分析方面，研究人员利用NCBI SRA数据库（BioProject ID: PRJNA705537）中的公共转录组数据（涵盖幼苗期根、茎、叶及营养期、现蕾期、初花期和盛花期的叶片）进行组织表达模式分析，并通过qRT-PCR对100 μM CdCl?和721 μM CuSO?处理下0、1、3、6、12、24小时六个时间点的叶片和根组织进行胁迫响应表达分析，每个处理设置三个生物学重复。
四、讨论与结论
讨论部分指出，甜叶菊中鉴定出的12个HMA基因数量略高于拟南芥、水稻和绣球，与槟榔和亚麻相近，但显著低于大豆和甘蓝型油菜等多倍体物种，表明多倍化事件可能在HMA基因家族的扩增中发挥更重要的作用。系统发育分析显示，与AtHMA5（主要在根中表达、受铜特异性诱导、对铜解毒至关重要）聚在同一分支的SrHMA7、SrHMA8和SrHMA10同样主要在根中表达，推测其在甜叶菊中可能发挥类似的铜解毒功能。亚细胞定位预测显示，SrHMA3作为AtHMA8的同源蛋白被预测定位于叶绿体，与拟南芥中AtHMA1、AtHMA6和AtHMA8参与叶绿体铜稳态的定位模式一致。值得注意的是，SrHMA7和SrHMA8在镉和铜胁迫下均表现出强烈且特异的诱导表达，这种独特的强响应模式在模式植物拟南芥和水稻的同源基因中未见报道，使其成为未来功能研究的重要候选基因。
研究结论指出，本研究首次对甜叶菊HMA基因家族进行了全面的生物信息学分析，鉴定出12个SrHMA基因，系统发育上分为Zn/Co/Cd/Pb和Cu/Ag两个亚群。同一亚群内成员具有保守的基序模式，亚群间存在差异。共线性分析表明SrHMA家族的扩增主要由片段重复事件驱动。启动子分析显示所有SrHMA基因均含有丰富的与逆境胁迫、激素响应和光响应相关的顺式作用元件。组织表达模式和胁迫响应分析揭示了SrHMA基因的组织特异性表达，以及在镉和铜胁迫下响应幅度、动态趋势和组织分布上的显著差异。这些发现为后续SrHMA基因功能研究及甜叶菊镉、铜适应分子机制的解析奠定了基础。特别是SrHMA7和SrHMA8在重金属胁迫下的强烈特异性诱导，使其成为通过转基因手段（如酵母或拟南芥异源表达、甜叶菊敲除/过表达株系构建）进行功能表征的优选候选基因，有助于阐明其在金属解毒和耐受中的精确作用。该研究发表于《Cells》期刊。
