整合表型、生理和比较转录组分析揭示洋葱细胞质雄性不育的调控网络

《Scientia Horticulturae》：Integrated phenotypic, physiological, and comparative transcriptome analysis reveals the regulatory network of cytoplasmic male sterility in onion

　　

洋葱（Allium cepa L.）是一种全球重要的辛辣蔬菜。雄性不育系的利用极大地提高了杂交种子生产的效率和纯度。然而，洋葱细胞质雄性不育（CMS）的分子机制仍不清楚。本研究采用整合表型观察、生理学检测和转录组分析的方法，比较了三对洋葱CMS系及其相应保持系在花药发育上的差异。花粉染色结果显示，CMS系的花粉粒表现出典型的败育表型，包括皱缩形态、不均匀的橙黄色染色以及细胞内部分物质丢失。生理学检测表明，CMS系的花药中可溶性糖含量显著降低，过氧化物酶（POD）、过氧化氢酶（CAT）和超氧化物歧化酶（SOD）活性显著下降，同时积累了过量的活性氧（ROS）和丙二醛（MDA），从而对花粉发育造成氧化损伤。比较转录组分析鉴定出2295个共有的差异表达基因（DEGs），功能富集分析表明这些共有DEGs主要富集在与花粉壁形成、碳水化合物代谢、糖转运、ROS代谢和植物激素信号转导相关的通路中。部分关键DEGs的表达模式通过实时荧光定量PCR（qRT-PCR）分析进一步验证。同时，CMS相关的线粒体基因orf725在CMS系花药中被发现异常高表达。基于这些发现和前期报道，研究人员提出了一个由线粒体到细胞核逆行信号介导的潜在调控网络，以解释洋葱雄性不育的分子基础。本研究为S型CMS的分子机制提供了新见解，并为洋葱杂交育种提供了有价值的理论基础。

**论文解读文章**

洋葱（*Allium cepa* L.）作为全球广泛栽培的辛辣蔬菜，其杂交种子生产中雄性不育系的利用极大提高了效率和纯度。目前细胞质雄性不育（CMS）主要分为S型、T型和R型，其中S型CMS广泛应用于“三系”育种体系。然而，洋葱CMS的分子机制尚不清楚，关键线粒体和核基因尚未系统鉴定，转录层面的调控网络有待阐明。为此，研究人员以三对CMS-S系（1217A、233A、245A）及其相应保持系（1217B、233B、245B）为材料，开展了整合表型观察、生理学检测和比较转录组分析的研究。研究提出由线粒体到细胞核逆行信号介导的调控网络，并报道了候选线粒体基因orf725，为洋葱CMS分子机制提供了新见解，为杂交育种奠定了理论基础。该论文发表在《Scientia Horticulturae》。

**主要关键的技术方法**（不超过250字）
研究人员基于三对来自山东省农业科学院试验基地的CMS-S系及保持系（核背景相同、细胞质不同）花蕾样本。采用1% I₂-KI染色进行花粉活力表型鉴定；利用对应试剂盒测定过氧化氢酶（CAT）、过氧化物酶（POD）、超氧化物歧化酶（SOD）活性及丙二醛（MDA）、活性氧（ROS）、可溶性糖含量。通过Trizol法提取总RNA后，使用Illumina NovaSeq? 6000平台进行转录组测序（RNA-seq）；以洋葱核参考基因组（nwpu）和线粒体参考基因组（KU318712.1）为参考，利用HISAT2比对、StringTie组装和DESeq2筛选差异表达基因（DEGs）；采用Blast2GO和KOBAS 2.0进行GO和KEGG富集分析；利用DeepTMHMM预测跨膜结构域；选取12个候选DEGs通过qRT-PCR验证RNA-seq数据。

**研究结果**
3.1 花药抗氧化生理特性和糖含量的测定
I₂-KI染色显示CMS系花粉粒皱缩、染色不均、内含物部分丢失。生理检测表明，相比保持系，CMS系花药中CAT、POD和SOD活性显著降低，MDA、O₂^-和H₂O₂含量显著升高，且总可溶性糖、葡萄糖和蔗糖含量均显著下降。结论：CMS花药ROS清除能力受损，氧化损伤加剧，营养供应不足。

3.2 转录组测序、序列比对和数据分析
共获得111.67 Gb高质量干净读长，Q30>98%，平均GC含量42.22%，>92%读长比对到参考基因组，鉴定出495,637个基因。样本间相关性R>0.99，主成分分析显示生物学重复聚类良好。结论：测序数据质量高，适合后续分析。

3.3 CMS与保持系间DEGs鉴定
三个比较组（1217B vs 1217A、233B vs 233A、245B vs 245A）分别鉴定出2727、4468和3725个DEGs。取共有DEGs得到2295个（927个三组共有，1368个任两组共有），其中488个为无蛋白注释的新基因。qRT-PCR验证12个关键DEGs，表达趋势与RNA-seq一致，相关系数R²=0.89。结论：共有的DEGs可能为核心调控因子。

3.4 GO功能注释和KEGG通路分析
GO分析显示1807个有注释的共有DEGs显著富集于花粉生长调控、细胞壁组织、蛋白质磷酸化、果胶酯酶活性、甲基氧化酶活性等术语。KEGG分析显示719个DEGs富集于123条通路，前20包括淀粉和蔗糖代谢、戊糖和葡萄糖醛酸转换、角质/木栓质/蜡生物合成，以及ROS清除相关通路（半胱氨酸和蛋氨酸代谢、苯丙烷生物合成、类胡萝卜素生物合成）。结论：CMS与多种代谢和信号通路失调密切相关。

3.5 线粒体基因组相关DEGs鉴定
两个线粒体基因orf725和nad5在CMS系中显著上调，且其编码蛋白均预测有跨膜结构域。参与氧化磷酸化和ATP合酶生物合成的多个核基因（如质膜ATP酶、ATP合酶亚基）在CMS系中显著下调。结论：orf725和nad5为洋葱CMS候选线粒体基因，可能通过破坏ATP合酶功能引发不育。

3.6 花药和花粉结构发育相关DEGs
共鉴定26个DEGs参与花药与花粉发育，其中23个在CMS系中下调。包括16个细胞壁动态重塑基因（PME、PL、PAE、果胶酶抑制子、XTH、CSLD4）、4个细胞骨架基因（ADF、微管相关蛋白）、4个花粉Ole e 1过敏原/延伸家族蛋白基因、2个ABC转运蛋白（ABCG28）。结论：这些基因下调导致花粉壁结构和细胞骨架缺陷，引起花粉败育。

3.7 碳脂代谢、ROS稳态和激素信号通路相关DEGs
碳水化合物与能量代谢：20个DEGs参与淀粉与蔗糖代谢、糖酵解/糖异生、TCA循环等，17个在CMS系中下调（如SUS、PFK、FRK、TPP）。脂质代谢与转运：14个DEGs参与角质/木栓质/蜡生物合成、类固醇、甘油酯/鞘脂代谢，其中CYP704A2、CYP96A10、GDSL脂肪酶下调。ROS清除：17个DEGs涉及苯丙烷生物合成（过氧化物酶PER1等）、半胱氨酸/蛋氨酸代谢（SAM、SAMDC等），多数过氧化物酶基因下调。植物激素信号转导：20个DEGs涉及生长素、脱落酸、赤霉素、茉莉酸信号，所有8个生长素相关基因下调，GA信号抑制因子GAI1下调，茉莉酸信号抑制因子上调。结论：糖供应不足、脂质代谢紊乱、ROS过量积累和激素信号失调共同参与CMS。

3.8 转录因子相关DEGs
鉴定24个主要转录因子（TF）基因，包括12个MYB/MYB相关、6个bHLH、4个MADS-box、2个NAC。19个在CMS系中下调（如MYB、bHLH、MADS），仅5个上调。结论：这些TF的表达异常可能破坏绒毡层程序性细胞死亡（PCD）和花粉发育。

**讨论与结论**
讨论部分深入分析了CMS的分子机制：线粒体基因orf725编码跨膜蛋白，可能干扰ATP合酶组装，导致氧化磷酸化受阻和ROS爆发，进而通过逆行信号调控核基因表达。花粉壁缺陷（PME、PL、XTH、CESA、CYP450、GDSL下调）导致外壁和内壁形成异常；糖代谢基因（SUS、PFK、FRK、TPP下调）造成能量供应不足；ROS清除基因（过氧化物酶、SAM等下调）引发氧化应激；激素信号（生长素、GA、JA通路）失调及TF（MYB、bHLH、MADS、NAC）表达改变进一步加剧绒毡层PCD异常和花粉败育。研究结论（翻译原文结语部分）：基于当前数据和前期报道，研究人员提出了一个由线粒体到细胞核逆行信号介导的调控网络来解释洋葱CMS-S系的雄性不育机制（图10）。首先，线粒体重排导致CMS相关基因orf725过表达，其跨膜蛋白破坏线粒体ATP合酶，抑制氧化磷酸化和ATP生成，造成ROS过量积累。随后，线粒体缺陷发送逆行信号调节核基因表达：下调碳水化合物代谢基因（如PFK、TPP、FRK）和花粉壁相关基因（如PME、PL-like、XTH、CESA）导致营养缺乏和内壁缺陷；ROS清除受损导致大量ROS积累和异常绒毡层PCD；激素信号紊乱和MYB/bHLH/MADS/NAC转录因子表达降低进一步抑制孢粉素合成相关基因（如CYP704A、CYP94D、CYP96D、GDSL），引起外壁缺陷。这些级联紊乱共同导致花粉败育和洋葱雄性不育。未来计划在拟南芥或本氏烟中进行orf725的异源功能验证。本研究为S型CMS分子机制提供了新见解，并为洋葱杂交育种提供了有价值的理论依据。