《Industrial Crops and Products》：Multi-genome comparison reveals a conserved 14–3-3 core and type-specific cold responses in
Brassica juncea
编辑推荐：
：芥菜（Brassica juncea (L.) Czern. & Coss.）的六个基因组比较揭示了一个保守的14-3-3蛋白核心以及可变的基因谱系。拷贝数变异集中在七个GF14-like类型中。在芥菜叶片中，仅GF14ψ-like基因受冷胁迫强烈诱
：芥菜（Brassica juncea (L.) Czern. & Coss.）的六个基因组比较揭示了一个保守的14-3-3蛋白核心以及可变的基因谱系。拷贝数变异集中在七个GF14-like类型中。在芥菜叶片中，仅GF14ψ-like基因受冷胁迫强烈诱导。BjuVB06G51590的表达与拟南芥（Arabidopsis）中更严重的冷损伤相关。研究人员对六个公共芥菜基因组进行了14-3-3基因家族的鉴定与比较分析，通过系统发育、直系同源群（orthogroup, OG）推断、序列进化、转座子（transposable element, TE）与启动子顺式元件分析、候选互作蛋白投射以及冷胁迫时间进程转录组分析，解析了该基因家族在拷贝数、亚基因组归属、编码序列和调控特征上的种内变异。结果表明，该家族包含40个直系同源群，其中25个在所有材料中保持单拷贝，拷贝数变异集中于GF14ν-like和GF14ω-like等七个类型。所有可评估的同线基因对均受到纯化选择（Ka/Ks < 1）。冷处理下，仅GF14ψ-like基因表现出强烈诱导，其中BjuVB06G51590在烟草叶片中定位于内质网（endoplasmic reticulum, ER）相关区域，其异源表达导致拟南芥在4 °C下过氧化氢（H2O2）和丙二醛（malondialdehyde, MDA）含量升高，超氧化物歧化酶（superoxide dismutase, SOD）、过氧化物酶（peroxidase, POD）和过氧化氢酶（catalase, CAT）活性降低，脯氨酸、可溶性糖和可溶性蛋白积累减少，表明该基因的表达与更严重的冷损伤相关。
论文解读文章
一、研究背景与科学问题
芥菜（Brassica juncea (L.) Czern. & Coss.）是一种重要的油料和蔬菜作物，其种子油可作为生物柴油原料，茎叶提供木质纤维素生物质，富含硫代葡萄糖苷的种子饼粕可制成生物熏蒸土壤改良剂。低温胁迫会降低芥菜幼苗可溶性蛋白含量、扰乱活性氧代谢、改变抗氧化酶活性并加剧脂质过氧化。14-3-3蛋白是一类能够识别磷酸化丝氨酸/苏氨酸残基并通过保守的两性沟槽与底物结合的调控蛋白，其结合可改变客户蛋白的活性、稳定性、亚细胞分布或复合物组装，在冷适应相关的磷酸化信号转导中发挥关键作用。
芥菜是B. rapa（AA）和B. nigra（BB）杂交形成的异源四倍体，其基因组同时保留A和B两个亚基因组，且祖先经历全基因组三倍化，导致基因家族成员复杂。单参考基因组无法区分物种内保守拷贝、材料依赖性变异以及组装注释造成的缺失，因此需要多基因组比较来解析家族组成、亚基因组保留和种内变异。此前关于芥菜14-3-3家族的研究尚未将材料水平拷贝保留、A/B亚基因组关系、局部调控变异、候选互作蛋白谱和时间冷响应综合考察。
二、研究内容与主要结论
研究人员比较了六个公共芥菜基因组的14-3-3拷贝保留、亚基因组关系、编码序列和局部调控特征，并通过T84-66叶片冷处理时间进程鉴定响应基因进行异源表达分析。主要发现包括：六个基因组共含262个14-3-3基因，所有十个GF14-like类型在每个材料中均存在；家族由40个直系同源群组成，其中25个在所有材料中保持单拷贝，拷贝数变异集中于七个GF14-like类型；GF14μ-like、GF14φ-like和GF14ψ-like的OG均属于稳定的单拷贝核心；所有可评估同线基因对的Ka/Ks均低于1，表明编码序列受纯化选择约束；冷处理下仅GF14ψ-like基因被强烈诱导，其中BjuVB06G51590在烟草中定位于内质网相关区域，其异源表达导致拟南芥在4 °C下氧化损伤标志物升高、抗氧化酶活性和渗透保护物质积累降低，表现出更严重的生理损伤。
三、关键技术方法
研究人员使用了以下主要技术方法：基于HMMER的PF00244结构域搜索和BLASTP互补检索鉴定候选蛋白，结合miniprot和GeMoMa对异常基因模型进行修正；MAFFT比对后以IQ-TREE进行最大似然系统发育推断，以SH-aLRT和UFBoot评估分支支持；OrthoFinder和Broccoli联合推断直系同源群，以MCScanX鉴定重复类型和同线关系，PAL2NAL生成密码子比对后以KaKs_Calculator估计Ka/Ks；HiTE进行转座子从头注释，PlantCARE和PLACE编译顺式元件目录；从BioGRID提取拟南芥14-3-3物理互作证据并以DIAMOND双向搜索投射到芥菜蛋白质组，BindSite预测结合位点；T84-66材料在4 °C下进行0、3、6、12、24小时冷处理，每个时间点三个生物学重复，RNA-seq数据以STAR比对、featureCounts定量，DESeq2分析差异表达，WGCNA构建共表达网络。
四、研究结果
4.1 保守核心支撑材料水平拷贝数变异
六个芥菜基因组共含262个14-3-3基因，T84-66和AU213各48个，HJ 46个，XC 43个，PM 39个，四川黄38个。所有十个GF14-like类型在每个材料中均出现。258个锚定基因中126个属于A亚基因组，132个属于B亚基因组。GF14ω-like是最大类型，每材料8-11个基因；GF14ι-like拷贝数范围最宽（2-6）。家族由40个OG组成：30个核心、8个软核心、2个外壳OG，无云OG。25个核心OG在所有材料中保持单拷贝，五个存在剂量变异。15个拷贝可变OG局限于七个GF14-like类型。T84-66和AU213在所有40个OG上拷贝数谱完全一致，四川黄和PM的较小家族源于不同模式。
4.2 系统发育解析类型特异性拷贝数变异
合并系统发育树将十个GF14-like类型解析为单系支系，八个类型根部支持强，GF14κ-like和GF14ν-like支持较弱。在包含两个祖先参考的九个类型中，118/122个A亚基因组蛋白最接近B. rapa，110/126个B亚基因组蛋白最接近B. nigra。40个OG中37个在芥菜单独树中单系。15个可变OG中，五个仅剂量变异，八个仅存在/缺失变异，GF14ε-S-OG03和GF14ν-SC-OG03兼具两种变异形式。
4.3 保守基因结构伴随部分OG的C端分化
基因结构和基序组成在大多数OG内保守。四外显子和七外显子结构占主导，分别出现在148和65个基因中。GF14κ-like基因一致为四外显子，GF14ι-like为七外显子，GF14π-like为六外显子。预测蛋白长度为216-322个氨基酸，基序1-3出现在全部262个蛋白中。35个OG的基序组成不变，变异局限于五个OG。OG中位平均成对蛋白一致性为99.81%，14个OG含相同蛋白，37个至少95%一致。C端分化更明显：九个OG末端20个氨基酸一致性低于80%，其中六个全长一致性至少95%。GF14μ-C-OG02全长一致性99.20%但C端仅34.67%。
4.4 纯化选择约束重复的14-3-3基因
全基因组或片段重复贡献240个基因（91.6%），22个为分散重复。18个分散重复属于B亚基因组，17个属于GF14ν-like或GF14ω-like。321个直接同线对中204个连接A和B亚基因组基因，317个连接相同GF14-like类型成员。所有319个有效Ka/Ks估计均低于1，中位值0.0492，范围0.0038-0.8317。GF14π-like对具有最高组内值（0.7824-0.8317），GF14κ-like对中位值0.0274。
4.5 局部转座子和顺式元件谱在GF14-like类型间变异
674个TE-区域关联中，321个位于启动子，288个位于下游，65个位于编码区内含子，无一与CDS/外显子重叠。DNA转座子占58.3%的关联和93.8%的内含子集。启动子TE出现在152个基因（58.0%），包括每个GF14π-like、GF14φ-like和GF14ψ-like成员。262个启动子含4783个顺式元件匹配。LTRE、DRE/CRT和ABRE分别出现在142、121和72个启动子中。冷相关元件频率在类型间差异显著：GF14π-like缺乏DRE/CRT且仅10.0%启动子含LTRE，而DRE/CRT出现在84.8%的GF14ν-like和LTRE出现在95.7%的GF14κ-like启动子中。TE区间与683个顺式元件匹配重叠（14.3%），但OG内顺式元件变异与TE位置变异仅弱相关。
4.6 投射候选互作蛋白谱在GF14-like类型间差异
从578个实验支持的拟南芥互作蛋白投射获得2326个非冗余芥菜候选群，A和B亚基因组代表相似。严格BindSite预测支持895个群（38.5%）。候选蛋白含5022个严格预测位点，其中54.4%缺乏经典Mode I、Mode II或Mode III样基序。跨材料比对将4867个位点归入1110个同源位置，中心丝氨酸/苏氨酸残基在95.2%的位置保守。GF14ψ-like、GF14ω-like和GF14κ-like具有最大投射集，分别为972、713和667个候选群。代谢与能量（20.8%）和转运与膜运输（16.9%）是最大功能类。GF14ε-like候选偏向代谢和RNA/蛋白生物合成，GF14κ-like和GF14ψ-like结合转运相关和代谢候选，转录调控因子占GF14λ-like的28.9%和GF14φ-like的32.0%。
4.7 冷胁迫引发类型特异性时间响应
T84-66叶片中，63,446个测试基因中52,104个显示整体时间响应，29,359个在至少一个时间点满足强响应标准。强响应基因从3小时的7057个增加到24小时的20,476个。48个14-3-3基因中36个显示整体响应，12个无响应。七个14-3-3基因满足强响应标准：两个诱导、五个抑制。两个诱导基因BjuVB06G51590和BjuVA04G11700均为GF14ψ-like，属于上升TC05轨迹和棕色共表达模块。第三个GF14ψ-like基因BjuVA07G20040在任何采样时间均未满足强响应标准。整体响应在非ε类（29/33）中比ε类（7/15）更频繁。39个同线A-B对中，九个仅一个同源基因显示强响应，六个相对0小时改变表达偏倚类，三个兼具两种情况。
4.8 BjuVB06G51590转基因拟南芥在冷处理中显示更大生理损伤
在烟草叶片中，BjuVB06G51590-EGFP荧光与内质网靶向mCherry标记重叠，另有ER标记阳性区域外的额外EGFP荧光，表明ER相关分布伴有非重叠组分。转基因拟南芥在冷处理前莲座较小、叶片较窄。所有八项生理性状在4 °C下0-48小时显示基因型依赖性时间轨迹。H2O2和MDA从3小时起在转基因植物中较高，48小时分别达到野生型的7.96倍和2.66倍。48小时时SOD、POD和CAT活性分别低65.1%、71.8%和61.8%。脯氨酸、可溶性糖和可溶性蛋白在48小时分别比野生型低41.0%、40.5%和46.8%。氧化损伤标志物增加和防护响应降低共同表明转基因材料在长时间冷处理中生理损伤更大。
五、讨论与结论
芥菜保留了一个在异源四倍体基因组中保守的14-3-3核心。25个OG在所有六个材料中保持单拷贝，包括每个GF14μ-like、GF14φ-like和GF14ψ-like OG，而GF14ν-like和GF14ω-like贡献了最多可变群。PM缺失多个非ε OG，四川黄主要在保留群内丢失额外拷贝，相似家族总数掩盖了基因含量和剂量的不同变化。GF14λ-like、GF14φ-like和GF14ψ-like偏向A亚基因组，GF14ε-like、GF14ν-like、GF14π-like和GF14ω-like偏向B亚基因组，解释了家族水平A/B平衡并不意味亚基因组贡献相当。
编码保守性在保留历史中持续存在：所有可评估同线对Ka/Ks低于1，大多数OG共享外显子-内含子结构和蛋白基序。GF14μ-C-OG02全长一致性99.20%但末端20个氨基酸仅34.67%一致性，表明保守折叠可容纳功能分化。GF14ψ-C-OG02在PM和XC中为四外显子，其他材料为三外显子但编码相同蛋白，区分了可能影响蛋白行为的末端氨基酸变异和不改变编码序列的外显子数变异。
调控变异与观察到的冷响应更直接相关。TE围绕基因但不与CDS/外显子重叠，所有GF14π-like、GF14φ-like和GF14ψ-like启动子含TE，而GF14μ-like缺乏下游TE。TE位置变异与顺式元件组成变异仅弱相关。在拷贝稳定的GF14ψ-like位点中，蔬菜材料T84-66、PM和XC的GF14ψ-C-OG01启动子共享两个DRE/CRT和一个LTRE，AU213、HJ和SY有一个DRE/CRT且无LTRE，但该作物用途分组在GF14ψ-C-OG03不持续。T84-66中BjuVA04G11700在含两种元件时被强烈诱导，BjuVB06G51590在两者均无时被诱导，BjuVA07G20040携带两者却仅显示时间响应，表明冷诱导需要超越基序计数的解释。拟南芥RCI1A同样经CBF非依赖途径冷诱导而无公认冷响应元件。
投射候选互作蛋白谱结合共享与类型偏好：约四分之三候选群映射到一个类型，GF14ε-like和GF14ω-like候选更受限，GF14ν-like、GF14φ-like和GF14λ-like共享更多。BindSite分析为895个群提供序列支持，多数比对位点中心丝氨酸/苏氨酸残基保守，超过一半为非经典位点。功能差异值得关注：GF14ε-like强调代谢和生物合成，GF14κ-like和GF14ψ-like连接代谢与转运，GF14λ-like和GF14φ-like包含更多转录调控因子，GF14ν-like强调蛋白稳态，GF14ω-like包含细胞骨架、细胞壁和细胞周期功能。
BjuVB06G51590转基因拟南芥中H2O2升高伴随MDA积累表明氧化负担增加伴随脂质过氧化，SOD、POD和CAT活性降低与长时间冷处理中抗氧化能力不足一致。脯氨酸和可溶性糖积累减弱及可溶性蛋白下降进一步提示代谢调整受损。冷处理前SOD活性和糖含量已较低表明存在组成性生理组分。该模式与拟南芥14-3-3ε或14-3-3ω过表达增强冷耐受相反，与RCI1A/GF14ψ负调控冻耐受和驯化方向一致。4 °C测定涉及冷响应而非冻敏感，该比较既不确立冻敏感也不将乙烯机制归于BjuVB06G51590。与ER标记的共定位将BjuVB06G51590置于与客户调控相关的空间背景中，ER相关组分可能影响相关客户的保留、组装或稳定性。
结论：芥菜14-3-3家族由一个保守核心和类型特异性变异组成，拷贝数变异集中于特定GF14-like类型，编码序列受纯化选择约束。冷处理下仅GF14ψ-like基因被强烈诱导，其中BjuVB06G51590的异源表达与拟南芥更严重的冷损伤相关，表明该家族成员的诱导并不必然赋予保护作用，生理结果取决于成员及其调控背景。多基因组比较将家族结构、调控变异和冷响应联系起来，为理解多倍体作物中基因家族进化与胁迫响应的关系提供了框架。
