《Genes》：Genome-Wide Characterization of the GmCOPT Family and Its Salt-Responsive Expression Patterns in Soybean
编辑推荐：
背景：铜转运蛋白（copper transporters, COPTs）对维持植物铜（Cu）稳态至关重要，该过程对植物生长、发育和胁迫适应具有关键作用。然而，大豆（Glycine max）GmCOPT家族仍缺乏系统表征。方法：研究人员对大豆COPT基因进行了全
背景：铜转运蛋白（copper transporters, COPTs）对维持植物铜（Cu）稳态至关重要，该过程对植物生长、发育和胁迫适应具有关键作用。然而，大豆（Glycine max）GmCOPT家族仍缺乏系统表征。方法：研究人员对大豆COPT基因进行了全基因组鉴定，开展了系统发育（phylogenetic）、共线性（collinearity）和进化（evolutionary）分析，预测了蛋白基序（protein motifs）、结构（structures）及启动子顺式调控元件（promoter cis-regulatory elements），并分析了GmCOPT家族的表达谱。结果：研究人员鉴定到10个GmCOPT基因。对五个植物物种COPT蛋白的系统发育分析显示，大豆COPT与蒺藜苜蓿（Medicago truncatula）COPT亲缘关系更近。共线性分析提示，片段复制（segmental duplication）是GmCOPT家族扩张的主要驱动力，所有复制基因对的Ka/Ks比值均低于1，表明其受到纯化选择（purifying selection）。表达谱分析揭示了GmCOPT基因间的组织偏好性表达模式。GmCOPT启动子区存在ABRE、MBS和LTR顺式作用调控元件（cis-acting regulatory elements），提示其可能参与非生物胁迫响应。RT-qPCR验证表明，GmCOPT基因在大豆根和叶中对盐胁迫的响应存在差异，一些基因显著受抑制，另一些基因明显受诱导。值得注意的是，GmCOPT7和GmCOPT8在两种组织中均受盐胁迫显著诱导。结论：这些发现提示GmCOPT基因经历了复制驱动的扩张和功能多样化，具有不同的表达模式和盐诱导响应。该工作为未来研究GmCOPT基因在非生物胁迫耐受（abiotic stress tolerance）中的功能提供了候选基因。
土壤盐渍化是限制全球可持续农业的重要因素。盐胁迫下，过量Na
+破坏细胞离子稳态，触发活性氧（reactive oxygen species, ROS）过量产生并损害光合作用，导致生长受抑、叶片早衰和产量损失；盐胁迫还抑制根系发育和养分吸收，并扰动包括铜（Cu）在内的必需微量元素稳态。Cu是植物必需微量元素，作为超氧化物歧化酶、质体蓝素和细胞色素c氧化酶等关键酶辅因子，参与抗氧化防御、光合电子传递和线粒体呼吸；但过量游离Cu可触发Fenton反应产生毒性羟自由基，损伤DNA、蛋白质和脂质。植物通过铜运输和区隔系统维持Cu稳态。铜转运蛋白（copper transporters, COPTs）是负责细胞Cu
+吸收和分配的整合膜载体。COPT在拟南芥、番茄、蒺藜苜蓿、水稻和玉米中已有功能研究，且在铜充足条件下，COPT可能通过介导Cu稳态并调节脱落酸（ABA）生物合成与信号，正调控盐胁迫响应。然而，大豆（Glycine max）GmCOPT家族研究不足，其在盐胁迫中的潜在作用尚未系统研究。因此，研究人员对GmCOPT家族进行全基因组鉴定，并整合蛋白性质、种间系统发育、基因复制、启动子顺式元件、组织表达谱和盐胁迫表达模式分析，以填补这一知识空白，并为耐盐大豆育种提供候选基因。
关键技术方法概括：研究人员从Phytozome v14获取G. max Wm82.a4.v1基因组、蛋白和GFF3注释，从Pfam获取PF04145 HMM，用TBtools-II在五物种中搜索（E值1 × 10
?5），经NCBI CDD验证。用MEGA 7.0最大似然法（ML）构建系统发育树（JTT+G，1000 bootstrap）。用EnsemblPlants大豆Glycine_max_v2.1和拟南芥TAIR10，经TBtools-II MCScanX做共线性及Ka/Ks。用MEME鉴定基序，SWISS-MODEL预测三维结构，PlantCARE分析2000 bp启动子顺式元件。组织表达数据来自Phytozome JGI Plant Gene Atlas（G. max cv. Williams 82；RNA-seq BioProject PRJNA372533），含根、根瘤、茎、叶、花、荚、种子。盐处理用大豆品种Tianlong 1，种子由国家大豆改良中心、南京农业大学提供，0或200 mM NaCl处理6、12、24 h后取根叶，RT-qPCR以Gm60S为内参，2
?ΔΔCt计算，双尾Student’s t检验。
3.1 大豆GmCOPT基因家族的全基因组鉴定与系统发育分析：研究人员利用保守COPT结构域PF04145进行全基因组HMM筛选，鉴定到10个大豆COPT基因，命名为GmCOPT1至GmCOPT10。其编码蛋白长度为135–155个氨基酸，分子量为14,849.64–16,927.20 Da，均为碱性蛋白。亚细胞定位预测显示多样化分布，包括细胞膜、叶绿体和/或高尔基体。基于五个物种COPT蛋白的多序列比对构建系统发育树，发现GmCOPT蛋白分布于多个进化枝而非单一豆科特异枝，且与蒺藜苜蓿COPT最近，提示GmCOPT家族可能发生功能分化。
3.2 GmCOPT基因家族的共线性与进化分析：全基因组共线性分析鉴定到9对片段复制基因，未检测到串联复制，表明片段复制是GmCOPT家族扩张的主要驱动力。所有旁系同源基因对Ka/Ks<1，提示纯化选择，其中GmCOPT2/GmCOPT3进化约束最强。大豆与拟南芥跨物种共线性显示，GmCOPT1、GmCOPT4、GmCOPT6、GmCOPT9和GmCOPT10与AT2G37925、AT3G46900和AT5G59030直系同源，提示祖先COPT基因在双子叶植物分化前已出现，后续古多倍化事件促进基因拷贝扩张；可计算Ka/Ks的直系同源对也支持纯化选择。
3.3 GmCOPT蛋白保守基序与预测三维结构：MEME分析鉴定到4个保守基序，Motif 1–3存在于所有GmCOPT成员中，Motif 4仅见于系统发育关系较近的GmCOPT2、GmCOPT3、GmCOPT7和GmCOPT8。SWISS-MODEL预测三维结构显示，近缘旁系同源蛋白结构更相似，而系统发育分化群体间结构不同。这些差异提示GmCOPT蛋白可能发生功能分化，但仍需实验验证。
3.4 GmCOPT基因启动子区顺式作用元件分析：研究人员分析10个GmCOPT基因起始密码子上游2000 bp启动子区，PlantCARE预测到激素响应、胁迫响应和发育/光响应三类元件。激素响应元件包括ABRE、TGACG-motif、TCA-element、GARE-motif等；胁迫相关元件包括干旱响应MBS、厌氧响应ARE和冷诱导LTR；还鉴定到光响应元件如G-box和发育相关O2-site。这些预测元件在各启动子中随机分布，无固定保守位置，且每个GmCOPT基因具有独特顺式元件特征，提示其可能受特定环境或发育信号差异调控，但需功能验证。
3.5 GmCOPT基因组织表达模式：基于七种大豆组织表达谱，GmCOPT基因呈现组织偏好性。GmCOPT1在花中高表达，其次为根瘤和叶；GmCOPT2主要在花中表达，其次为茎、种子和根；GmCOPT6在叶中表达最强；GmCOPT4在多数组织中几乎检测不到；其余GmCOPT3、GmCOPT5、GmCOPT7–GmCOPT10呈低至中等表达。这些结果提示不同GmCOPT成员可能在不同大豆器官中具有功能特化。
3.6 GmCOPT基因对盐胁迫的表达响应：RT-qPCR检测盐处理6、12和24 h后根和叶中GmCOPT转录丰度。根中，GmCOPT2和GmCOPT10在6、12和24 h均显著下调，GmCOPT7在三个时间点均被诱导。叶中，GmCOPT7在6 h和24 h显著诱导，GmCOPT8在三个时间点均被诱导。GmCOPT4和GmCOPT9在根中低表达或检测不到；GmCOPT4、GmCOPT5、GmCOPT9和GmCOPT10在叶中低表达或检测不到。结果表明GmCOPT家族成员对盐处理响应存在差异，GmCOPT7和GmCOPT8为关键盐胁迫响应候选基因，其功能仍需研究。
讨论部分总结：土壤盐渍化是大豆产量和品质提升的重要非生物胁迫限制因子。Cu参与光合作用和抗氧化代谢，Cu稳态对植物生长和胁迫适应至关重要；盐胁迫诱导的微量营养失衡会加剧氧化损伤并抑制生长。尽管COPT已在拟南芥、水稻和番茄等模式植物和主要作物中表征，但大豆GmCOPT家族仍了解不足。研究人员系统表征了GmCOPT家族的系统发育关系、基因扩张模式、进化选择压力、组织表达谱及盐胁迫时空转录响应，填补了大豆研究空白。片段复制和进化选择是塑造植物基因家族扩张、结构保守和功能分化的核心驱动力；与53个绿藻和陆生植物COPT家族分析一致，大豆GmCOPT扩张主要由片段复制驱动，旁系和直系同源基因对均受纯化选择。与拟南芥和禾本科作物COPT蛋白保守的单一亚细胞定位不同，大豆GmCOPT成员预测定位于细胞膜、叶绿体和高尔基体，提示其可能具有分化功能。AtCOPT主要在缺铜条件下介导高亲和Cu吸收和维管Cu转运，而水稻COPT基因可受重金属胁迫特异诱导。本研究发现大豆GmCOPT启动子区富含激素和胁迫响应顺式元件，提示其可能与激素信号和盐胁迫响应相关。玉米COPT基因在多种组织中组成型表达，而大豆GmCOPT基因呈明显组织偏好性：GmCOPT1主要在花中表达，其次为根瘤和叶；GmCOPT2主要在花中表达，其次为茎、种子和根；GmCOPT6在叶中最高。GmCOPT在根和叶中对盐胁迫响应不同：根中GmCOPT7和GmCOPT8显著诱导，尤其在12和24 h，而多数其他可检测GmCOPT受抑制；叶中多个GmCOPT在盐胁迫6 h诱导，GmCOPT7在6和24 h显著上调。这些结果提示GmCOPT基因以时空方式响应盐胁迫，与COPT在不同器官介导铜再分配的功能一致，如AtCOPT1（AT5G59030）参与根Cu获取，AtCOPT6（AT2G26975）参与向地上部Cu分配，AtCOPT5（AT5G20650）动员细胞器储存Cu。由于铜是Cu/Zn-SOD必需辅因子，维持Cu稳态对抗氧化防御至关重要，可能参与盐响应调节，但相关假设仍需实验验证。GmCOPT不在SoyBase已知耐盐相关QTL内，提示其可能不是这些QTL的因果基因。研究人员在全基因家族水平验证了GmCOPT对盐胁迫的转录响应，但各成员具体功能和分子调控机制仍待阐明。未来可验证关键候选基因的Cu转运活性、Cu含量、过表达或敲除株系盐胁迫表型，并探索上游调控因子和下游互作蛋白。
研究结论翻译：总之，研究人员系统表征了大豆GmCOPT基因家族，揭示其片段复制驱动扩张和纯化选择。启动子和表达分析提示复杂转录调控、组织偏好性表达模式和盐诱导响应，GmCOPT7和GmCOPT8被鉴定为关键盐胁迫响应候选基因，暗示GmCOPT旁系同源基因间潜在功能分化。该研究为未来GmCOPT基因在Cu稳态和胁迫耐受中的功能研究奠定基础，并可能对耐盐大豆品种选育具有潜在意义。
