甘蓝型油菜（Brassica napus）中CYP90基因家族的表征及其在非生物胁迫与发育阶段下的表达模式

《Agronomy》：Characterization of the CYP90 Gene Family and Its Expression Patterns Under Abiotic Stress and Developmental Stages in Brassica napus

　　

细胞色素P450（cytochrome P450）超家族对植物生长发育与胁迫响应至关重要，然而甘蓝型油菜（Brassica napus）中的CYP90基因家族仍缺乏系统表征。本研究以甘蓝型油菜品种“ZS11”为参考基因组，系统鉴定出分布于16条染色体上的25个BnCYP90基因，并依据系统发育关系将其划分为CYP90A、CYP90B、CYP90C和CYP90D四个分支。同线性与进化分析表明该基因家族主要受纯化选择作用。启动子分析揭示了多样化的胁迫与激素响应元件；网络分析预测BnCYP90_4、BnCYP90_6和BnCYP90_11在候选蛋白互作及miRNA介导的调控网络中为高度连接节点。此外，自然群体单核苷酸多态性（SNP）分析显示，BnCYP90_25中位于C09:66,297,283的无义错义变异（missense variant）与盐胁迫下种子发芽率及地上部生物量显著关联。转录组与RT-qPCR分析揭示了非生物胁迫下明显的时空表达特征：BnCYP90_25在盐胁迫下持续激活，BnCYP90_9和BnCYP90_20在渗透胁迫中期响应，BnCYP90_21和BnCYP90_25在碱胁迫下快速响应。加权基因共表达网络分析（WGCNA）确定BnCYP90_13为T1发育阶段的枢纽基因。综上，本研究全面表征了BnCYP90家族，为未来在耐逆甘蓝型油菜育种中的功能验证提供了有价值的候选基因。

研究背景方面，全球人口增长与气候变化导致干旱、盐渍化、重金属污染等非生物胁迫日益威胁作物产量与品质。甘蓝型油菜（Brassica napus L.）作为全球重要油籽作物，其产量潜力与栽培范围受非生物胁迫限制愈发严重。油菜经历复合异源多倍化事件，基因组高度扩张且结构复杂，为解析基因家族进化动态提供理想框架。油菜素内酯（brassinosteroids，BRs）是调控细胞伸长、维管分化及非生物胁迫响应的关键甾醇类激素，而CYP90家族（含CYP90A/B/C/D）是BR生物合成通路的核心催化装配线，其中CYP90B介导的C-22羟化是整条通路的初始限速步骤。尽管模式植物中单个CYP90基因已有研究，但甘蓝型油菜全基因组水平的CYP90系统鉴定、进化机制与胁迫表达谱仍属空白，因此开展本研究具有明确理论价值与育种应用前景。该研究发表于《Agronomy》。

研究人员基于甘蓝型油菜ZS11参考基因组，整合隐马尔可夫模型（HMM）筛选、BLASTP互证与系统发育分类，从2311份核心自然群体重测序数据及BnIR数据库八 accession 五阶段发育转录组、盐/旱时序转录组中开展分析；主要关键技术方法包括：基于Pfam PF00067结构域的HMMER候选池构建与拟南芥CYP90参照蛋白互证鉴定；最大似然法（LG+G模型）系统发育树构建；MCScanX同线性检测与ParaAT/KaKs_Calculator 2.0（Yang–Nielsen模型）选择压力计算；MEME/NCBI CDD基序与结构域注释及PlantCARE顺式元件预测；STRING蛋白互作网络与psRNATarget的miRNA靶向预测；WGCNA（β=12，royalblue模块）共表达枢纽挖掘；BnVIR库2311份自然群体CC/TT基因型分组与Welch t检验关联分析；以及ZS11幼苗根组织盐（0.17 mol L^?1 NaCl）、碱（0.083 mol L^?1 NaHCO₃）、渗（20% PEG 6000）处理的RT-qPCR验证。

研究结果部分，3.1节通过ProtParam与CELLO预测，25个BnCYP90蛋白长294–542 aa，多碱性、亲水，20个定位于质膜、5个定位于细胞质，符合细胞色素P450膜结合特性。3.2节染色体定位显示25基因不均分布于16条染色体，A亚基13个、C亚基12个，局部串联重复与A/C同源染色体对应保守并存不对称丢失。3.3节系统发育将25成员归为CYP90A（8）、CYP90B（8）、CYP90C（4）、CYP90D（5）四类，各含A/C亚基姊妹对。3.4节同线性分析获甘蓝型油菜与B. rapa 53对、与B. oleracea 47对合成关系，组内23对非冗余同线性对，Ka/Ks均<1（0.1229–0.5775），证实纯化选择主导。3.5节显示所有成员含CYP90-like域，171个内含子均为经典GT–AG剪接位点；启动子富集光响应（G-box、Box 4）、激素（ABRE、CGTCA、TCA、P-box）与胁迫元件（MBS、LTR、ARE）。3.6节STRING网络指出BnCYP90_4、BnCYP90_6、BnCYP90_11为高度连接中心节点。3.7节psRNATarget预测8个miRNA（miR156b/c/g、miR159、miR172a/d、miR396a、miR6034）靶向9个BnCYP90，以裂解与翻译抑制双重模式调控。3.8节BnVIR 2311份群体靶向C09:66,297,283错义SNP，CC基因型（428份）盐胁迫下发芽势与上重比（T/CK）优于TT基因型（1815份），TT基因型种子发芽率与BnCYP90_25表达更高（p<0.01）。3.9节转录组时序显示BnCYP90_2/15叶中盐诱导强，BnCYP90_13/25叶基线高且盐旱双诱导，BnCYP90_20根中盐旱后期上调。3.10节WGCNA royalblue模块与T1阶段最强正相关，BnCYP90_13与WRKY70、PME18、GSTU16、YSL1、RGL1等共表达。3.11节RT-qPCR证实BnCYP90_25盐处理24 h持续峰值，BnCYP90_9/20于PEG 12 h中段峰值，BnCYP90_21/25于NaHCO₃ 6 h早期峰值。

讨论部分指出，甘蓝型油菜CYP90家族扩张源于异源四倍化后A/C亚基保留，同线性与Ka/Ks<1印证纯化选择维持BR通路核心结构；四类分支与拟南芥高度共线，少数成员（BnCYP90_10）基序简化反映组织多样化；启动子多元激素/胁迫元件与PPI、miRNA网络共同构成多维调控；BnCYP90_25错义SNP揭示盐胁迫下“发芽速率–苗期生物量维持”生理权衡的遗传基础，CC型偏保生物量、TT型偏快发芽且基因高表达；BnCYP90_13作T1阶段枢纽联携细胞壁修饰（PME18）、氧化还原（GSTU16、WRKY70）与赤霉素互作（RGL1）模块。结论部分翻译如下：本研究在甘蓝型油菜全基因组水平系统鉴定25个BnCYP90基因，提供该家族进化、结构、表达特征预测分析，揭示多倍化中基因保留模式、胁迫下复杂调控网络及表达多样性；未来应聚焦关键基因基因组编辑功能验证与自然变异关联作图，定位关联农艺性状（发芽率、生物量）的SNP，从而为培育高产品质与耐逆甘蓝型油菜品种提供理论基础。