《Frontiers in Plant Science》：Integrated multi-omics analysis reveals distinct molecular features of key stages in cotton fiber development in Gossypium hirsutum
编辑推荐：
引言：棉花纤维发育是纤维品质的关键生物学过程，其调控涉及转录、翻译和代谢等多个分子层面。方法：为系统阐明陆地棉（*Gossypium hirsutum*）纤维发育的分子基础，研究人员选用了两个纤维品质对比鲜明的品种，即苏棉11号（SM11）和豫棉5号（YM5）
引言：棉花纤维发育是纤维品质的关键生物学过程，其调控涉及转录、翻译和代谢等多个分子层面。方法：为系统阐明陆地棉（*Gossypium hirsutum*）纤维发育的分子基础，研究人员选用了两个纤维品质对比鲜明的品种，即苏棉11号（SM11）和豫棉5号（YM5），并对开花后15、20、25和30天（DPA）采集的纤维进行了整合的转录组学、蛋白质组学和代谢组学分析。结果：随着纤维发育进程，两个品种中差异表达基因（DEGs）、差异积累蛋白（DAPs）和差异积累代谢物（DAMs）的数量均增加，表明发育后期存在广泛的分子重编程。跨组学比较识别出植物激素信号转导、淀粉和蔗糖代谢、苯丙烷生物合成、类黄酮生物合成、角质、木栓质和蜡质生物合成以及ABC转运蛋白为纤维发育中共同涉及的核心通路。SM11在苯丙烷代谢、细胞色素P450和MAPK信号通路中表现出更强的富集，而YM5在核糖体相关过程、脂肪酸延长和氮代谢中表现出更显著的富集。蛋白质组学和代谢组学分析进一步证实了两个品种在苯丙烷和类黄酮代谢、糖代谢以及脂质代谢方面的显著差异。整合多组学分析进一步表明，苯丙烷和类黄酮生物合成构成了跨三个组学层面的关键协调模块。代谢物如花旗松素、二氢杨梅素和芥子醛与候选基因和蛋白质密切相关，共同形成一个相互调控的网络。讨论：这些发现为棉花纤维发育的分子调控提供了机制性见解，并识别出用于陆地棉（*G. hirsutum*）纤维品质改良的候选分子靶点。
棉花纤维发育是决定纤维长度、强度和细度等品质性状的关键生物学过程，其调控涉及转录、翻译和代谢等多个分子层面。尽管前人研究已阐明纤维起始、伸长、过渡和次生壁形成等阶段，但15–30 DPA这一关键窗口期（涵盖纤维伸长终止、细胞壁重塑和快速次生壁沉积）仍缺乏系统表征，尤其是整合转录组、蛋白质组和代谢组数据对陆地棉（*Gossypium hirsutum*）不同基因型进行对比分析的研究。为此，研究人员选取两个纤维品质对比鲜明的品种（苏棉11号SM11，低品质；豫棉5号YM5，高品质），在开花后15、20、25和30 DPA采集纤维，开展整合转录组学、蛋白质组学和代谢组学分析，旨在揭示该关键窗口期分子差异，鉴定与纤维品质变异相关的候选基因、蛋白质和代谢物，为未来功能验证提供基础。该研究发表在《Frontiers in Plant Science》。
研究人员主要采用以下关键技术方法：利用Illumina HiSeq 2500平台进行RNA-seq转录组测序（样品来源：新疆农科院棉花研究所，SM11和YM5品种，纤维样品采集于15、20、25、30 DPA，各3个生物学重复）；采用液相色谱-串联质谱（LC-MS/MS）进行非靶向代谢组学分析（各6个生物学重复）；采用数据非依赖性采集（DIA）模式结合Q Exactive HF-X质谱仪进行蛋白质组学分析（各3个生物学重复）。通过DESeq2、metaX和DIA-NN等软件分别鉴定差异表达基因（DEGs）、差异积累代谢物（DAMs）和差异积累蛋白（DAPs），并进行KEGG通路富集分析及Pearson相关性网络构建。
研究结果如下：
**Quality and reproducibility assessment of transcriptomic sequencing data**：对24个纤维样品进行RNA-seq，获得155.51 Gb清洁数据，各样品Q30值96.98%–97.95%，与参考基因组比对率87.75%–90.29%，生物学重复间相关性系数>0.94，PCA显示重复聚类良好，qRT-PCR验证10个基因与RNA-seq结果显著正相关（R=0.90，P<0.01），表明转录组数据可靠。
**Expression patterns and functional enrichment analysis of differentially expressed genes across developmental stages**：SM11中鉴定出25,073个DEGs，YM5中鉴定出28,599个DEGs，两个品种的DEGs数量均随发育进程递增，表明后期转录重编程广泛。KEGG富集分析显示，两个品种共同富集于植物激素信号转导、淀粉和蔗糖代谢、糖基转移酶、类黄酮生物合成、角质/木栓质/蜡质生物合成、ABC转运蛋白等通路，但SM11在苯丙烷生物合成、细胞色素P450和MAPK信号中富集更强，YM5在核糖体、脂肪酸延长和氮代谢中更显著。品种间比较鉴定出5,644个非冗余DEGs，主要富集于淀粉和蔗糖代谢、苯丙烷生物合成、类黄酮生物合成、植物激素信号转导等，提示糖代谢、次生代谢和激素信号参与品种间纤维发育差异。共享DEGs的k-means聚类分析揭示四个表达模式，涉及细胞壁相关代谢、次生代谢和激素信号。
**Proteomic expression characteristics and analysis of differentially accumulated proteins**：SM11中鉴定出878、1,802和2,393个DAPs（vs 15 DPA），YM5中鉴定出801、1,545和2,058个DAPs，均随发育递增。DAPs主要富集于苯丙烷生物合成、类黄酮生物合成、角质/木栓质/蜡质生物合成、淀粉和蔗糖代谢等通路。品种间比较鉴定出1,299个非冗余DAPs，20 DPA差异最大，富集于类黄酮生物合成、苯丙烷生物合成、黄酮和黄酮醇生物合成等，表明苯丙烷/类黄酮代谢、碳水化合物代谢和表皮/细胞壁过程与品种间发育差异密切相关。
**Analysis of differentially accumulated metabolites and metabolic pathway enrichment**：共鉴定1,995种代谢物，脂质和脂质样分子、有机酸及衍生物、有机杂环化合物最丰富。SM11中鉴定出649、1,056和1,173个DAMs，YM5中鉴定出834、935和1,064个DAMs，均随发育递增，且下调积累代谢物占主导。DAMs主要富集于类黄酮生物合成、苯丙烷生物合成、氧化磷酸化、三羧酸循环等通路。品种间比较显示20 DPA代谢差异最大，DAMs富集于类黄酮生物合成、苯丙氨酸代谢、苯丙烷生物合成等，表明苯丙烷/类黄酮代谢、氨基酸代谢和脂质相关通路贡献于品种间代谢分歧。
**Integrated transcriptomic and metabolomic analysis**：整合分析显示转录组和代谢组共同显著富集于苯丙烷生物合成、类黄酮生物合成和ABC转运蛋白。Pearson相关性分析构建了基因-代谢物关联网络，关键代谢物如阿魏酸、芥子醇、芥子醛、5-羟基阿魏酸、芥子酸、花旗松素、二氢杨梅素与多个候选基因显著相关，鉴定出GH_A05G1187、GH_A05G1439等潜在枢纽基因，表明苯丙烷和类黄酮代谢是后期纤维发育中跨组学调控的主要模块。
**Integrated transcriptomic and proteomic analysis**：共有2,995个差异分子在转录组和蛋白质组间共享，呈显著正相关（R=0.31，P<0.01），但部分基因和蛋白表达趋势不一致，提示存在转录后和翻译后调控。共享分子主要富集于类黄酮生物合成、角质/木栓质/蜡质生物合成、甘油脂代谢、脂肪酸代谢等通路，表明次生代谢、脂质代谢和表皮生物合成是转录组-蛋白质组协调变化的核心过程。
**Integrated transcriptomic, proteomic, and metabolomic analysis**：三组学整合分析显示苯丙烷生物合成和类黄酮生物合成在三个组学层面均显著富集。构建mRNA-蛋白质-代谢物关联网络，花旗松素、二氢杨梅素和芥子醛与相应候选基因和蛋白质表现出高度一致的积累模式，进一步确认苯丙烷和类黄酮代谢是棉花纤维发育中关键的跨组学调控模块。
讨论部分总结：研究结果表明，两个品种的发育差异并非源于单一通路的紊乱，而是多个调控层在关键窗口期协调重塑的结果。后期纤维发育以强烈的阶段依赖性重编程为特征，尤其是植物激素信号转导、淀粉和蔗糖代谢、苯丙烷和类黄酮代谢、脂质和表皮相关代谢及ABC转运蛋白过程。20 DPA是品种间蛋白质和代谢物差异最显著的阶段，可能对应纤维从伸长向次生壁沉积的过渡。两个品种的代谢偏好差异（SM11偏向信号响应和次生代谢，YM5偏向蛋白质合成和脂质延长）揭示了不同遗传背景下纤维发育的替代分子途径。苯丙烷和类黄酮代谢在三个组学层面一致性最强，花旗松素、二氢杨梅素和芥子醛作为关键候选代谢物，与细胞壁生物合成的前体供应和代谢通量重新分配相关，但需进一步功能验证。转录组-蛋白质组相关性中等（R=0.31），表明翻译后调控在纤维发育中起重要作用。研究结论（Conclusions）翻译如下：总之，对SM11和YM5发育中纤维的整合转录组、蛋白质组和代谢组分析揭示了在15–30 DPA关键窗口期内广泛的分子重编程，尤其是在植物激素信号转导、淀粉和蔗糖代谢、角质/木栓质/蜡质生物合成、ABC转运蛋白相关过程以及苯丙烷和类黄酮代谢中。其中，苯丙烷生物合成和类黄酮生物合成表现出最强的跨组学一致性，支持其与品种特异性纤维发育分歧的密切关联。整合分析进一步识别出花旗松素、二氢杨梅素和芥子醛为与相应基因和蛋白质相关的关键候选代谢物，为未来功能验证和纤维品质改良提供了潜在靶点。
