整合转录组与脂质组学揭示工程化甘蓝型油菜种子中神经酸生物合成的分子机制

《Oil Crop Science》：Transcriptomic and lipidomic analysis reveal molecular mechanisms in the biosynthesis of Nervonic acid in engineered Brassica napus seeds

　　

神经酸（Nervonic acid, NA）是一种具有显著药用价值的超长链单不饱和脂肪酸，但其进一步开发与应用受限于资源稀缺。在研究人员前期研究中，利用多基因载体（Napin-3和Napin-5）对甘蓝型油菜（Brassica napus, B. napus）进行代谢工程改造，成功实现了NA的高水平生产，转基因株系中NA积累量可达总脂肪酸的45%。本研究在开花后20天和40天（days after flowering, DAF）两个时间点，利用整合转录组学与脂质组学分析，阐明表达CgKCS基因的转基因B. napus种子中NA积累的分子基础。整合转录组与脂质组分析揭示，NA生物合成主要发生在种子发育后期，并被分隔至三酰甘油（triacylglycerol, TG）中，从而促成TG（24:1_18:2_18:3）含量的增加。转录组分析鉴定出大量与脂质代谢和种子发育相关的差异表达基因（differentially expressed genes, DEGs）。值得注意的是，参与脂肪酸延长和甘油酯生物合成的多个内源基因在转基因种子中的表达发生显著改变。其中，与超长链脂肪酸延长相关的内源KCS6和KCS19表达下调，而甘油酯组装关键酶GPAT6和LPAAT1则表达上调。研究结果表明，B. napus中NA的高水平生产不仅依赖于异源基因的活性，还依赖于内源酰基转移酶的协同响应以及对竞争性延长途径的抑制。这些结果为油料作物的代谢可塑性提供了新的机制见解，并确定了进一步优化稀有脂肪酸生产的关键遗传靶点。

论文解读

一、研究背景与意义

神经酸（Nervonic acid, NA；C24:1 ω-9）是一种超长链单不饱和脂肪酸，在缓解和预防神经退行性疾病方面显示出显著潜力。然而，NA的规模化生产和应用长期受制于资源稀缺。目前，NA的主要天然来源为元宝枫（Acer truncatum）种子油，但其作为木本植物存在成熟周期长（8-10年）和产量低的问题，严重制约了NA的工业化应用和科学研究。近年来，植物代谢工程的发展为NA生产提供了新思路。NA在植物中的生物合成起始于质体中的de novo脂肪酸合成，随后在内质网中经脂肪酸延长酶（fatty acid elongase, FAE）复合体催化完成碳链延长。其中，3-酮脂酰-CoA合酶（3-ketoacyl-CoA synthase, KCS）是决定底物和产物特异性的限速酶。CgKCS（来源于Cardamine graeca的KCS基因）是工程化植物和酿酒酵母生产NA的主要改造基因，偏好将C22:1和C20:1转化为C24:1。NA主要通过Kennedy途径组装并储存于三酰甘油（triacylglycerol, TG）中，通常由甘油-3-磷酸酰基转移酶（glycerol-3-phosphate acyltransferase, GPAT）和二酰甘油酰基转移酶（diacylglycerol acyltransferase, DGAT）分别将其掺入TG的sn-1和sn-3位。由于大多数植物仅表达偏好将C18脂肪酸掺入sn-2位的Class-A溶血磷脂酸酰基转移酶（lysophosphatidic acid acyltransferase, LPAAT），NA在sn-2位的积累潜力尚未被开发。在前期研究中，研究人员利用Napin-3（CgKCS::SLC1-1::AtDGAT1）和Napin-5（CgKCS::SLC1-1::BnFAE1::LdLPAAT::AtDGAT1）两个多基因载体改造甘蓝型油菜中双6号（ZS6），获得转基因株系N3Z和N5Z，NA含量从不可检测提升至43.8%（N3Z）和39.6%（N5Z），最佳株系达到44.9% NA，同时保持亲本油含量（39%）和产量（2,609 kg/ha），展现出商业应用潜力。为深入解析其分子机制，研究人员开展了本项整合转录组与脂质组学研究。

二、主要技术方法

本研究采用整合转录组学与脂质组学分析策略。转录组分析方面，研究人员利用Illumina RNA-seq技术对N3Z、N5Z转基因株系及非转基因对照ZS6在开花后20天（20 DAF）和40天（40 DAF）的种子构建了18个cDNA文库进行测序，共获得153.5 Gb Clean Data，Q30碱基比例超过94%。测序数据经fastp质控后，使用HISAT2比对至ZS11参考基因组，并通过DESeq2进行差异表达基因（differentially expressed genes, DEGs）鉴定，筛选标准为|log₂Fold Change| ≥ 1且FDR < 0.05。脂质组分析方面，基于AB Sciex QTRAP 6500 LC-MS/MS平台对相同样本进行脂质检测，利用正交偏最小二乘判别分析（Orthogonal Partial Least Squares Discriminant Analysis, OPLS-DA）提取变量投影重要性（Variable Importance in Projection, VIP）值，以VIP ≥ 1和|log₂FC| ≥ 1.0为标准鉴定差异脂质（differential lipids, DLs）。此外，研究人员还通过qRT-PCR对关键DEGs进行了验证。样本来源为T₃代N3Z和N5Z转基因株系及ZS6非转基因对照，每种样本设置三个独立生物学重复。

三、研究结果

1. 转录组测序与DEGs鉴定

PCA分析显示，同一发育阶段的样本聚在一起，第一主成分（PC1）解释了49.14%的方差，但未能清晰区分转基因与非转基因株系，表明基因型可能并非主要变异来源。样本相关性热图显示，同组内重复样本高度一致（Pearson相关系数r > 0.95），而不同时间点样本相关性显著降低（0.65-0.78）。在20 DAF，N3Z12_20d vs Z_20d组鉴定出121个DEGs（104个上调，17个下调）；N5Z51_20d vs Z_20d组鉴定出648个DEGs（262个上调，386个下调）；N5Z51_20d vs N3Z12_20d组鉴定出816个DEGs（317个上调，499个下调）。在40 DAF，N3Z12_40d vs Z_40d组鉴定出93个DEGs（37个上调，56个下调）；N5Z51_40d vs Z_40d组鉴定出449个DEGs（143个上调，306个下调）；N5Z51_40d vs N3Z12_40d组鉴定出438个DEGs（167个上调，271个下调）。在相同植株不同发育阶段比较中，Z_40d vs Z_20d组鉴定出29,266个DEGs，N3Z12_40d vs N3Z12_20d组鉴定出27,631个DEGs，N5Z51_40d vs N5Z51_20d组鉴定出29,307个DEGs。

2. DEGs功能富集分析

GO富集分析显示，20 DAF时，N3Z12_20d vs Z_20d组DEGs主要富集于热应激反应和未折叠蛋白结合；N5Z51_20d vs Z_20d组主要富集于线粒体ATP合成和呼吸链复合体III活性；N5Z51_20d vs N3Z12_20d组主要富集于染色质结构和核小体组装。40 DAF时，N3Z12_40d vs Z_40d组主要富集于硫饥饿响应和硫代葡萄糖苷分解代谢；N5Z51_40d vs Z_40d组主要富集于DNA双链断裂修复和硫代葡萄糖苷分解代谢；N5Z51_40d vs N3Z12_40d组主要富集于未折叠蛋白结合和过氧化氢响应。不同发育阶段比较组中，DEGs均主要富集于光合作用和叶绿素代谢过程。KEGG富集分析显示，20 DAF时，N3Z12_20d vs Z_20d组主要富集于内质网蛋白质加工和脂肪酸生物合成；N5Z51_20d vs Z_20d组主要富集于苯丙烷生物合成和氧化磷酸化。40 DAF时，各组主要富集于氮代谢和氨基酸代谢等通路。

3. 脂质组分析与DLs鉴定

共检测到874种脂质，涵盖FA、GL、GP、SP、ST、PR、SoG1和PK八类ClassⅠ脂质及25类ClassⅡ脂质，其中TG、DG、PC、PE和MGDG为主要脂质。脂质组PCA显示，PC1和PC2分别解释了40.15%和30.87%的方差，PC1主导不同发育阶段的变异，而PC2能显著区分转基因与非转基因株系。在20 DAF，N3Z12_20d vs Z_20d组鉴定出163个DLs（15个下调，148个上调）；N5Z51_20d vs Z_20d组鉴定出138个DLs（19个下调，119个上调）。在40 DAF，N3Z12_40d vs Z_40d组鉴定出391个DLs（242个下调，149个上调）；N5Z51_40d vs Z_40d组鉴定出388个DLs（234个下调，154个上调）。发育阶段比较中，Z_40d vs Z_20d组鉴定出503个DLs，N3Z12_40d vs N3Z12_20d组鉴定出375个DLs，N5Z51_40d vs N5Z51_20d组鉴定出369个DLs。

4. DLs功能富集分析

KEGG富集分析显示，20 DAF时，转基因与非转基因比较组的DLs主要富集于甘油酯代谢、α-亚麻酸代谢、花生四烯酸代谢和亚油酸代谢。40 DAF时，DLs仍显著富集于甘油酯代谢，同时涉及自噬和糖基磷脂酰肌醇（GPI）锚定生物合成等应激响应通路。发育阶段比较中，非转基因株系DLs主要富集于甘油酯代谢，而转基因株系DLs则富集于不饱和脂肪酸生物合成、次生代谢物生物合成和应激响应通路。

5. 转录组与脂质组整合分析

KEGG共富集分析揭示，甘油酯代谢是连接转录变化与脂质重塑的核心通路，这与Kennedy途径将脂肪酸导入TG储存的功能一致。靶向通路分析显示，在NA合成与积累途径中，转基因株系在40 DAF时游离脂肪酸FFA（C18:1）含量显著降低。PA含量在20 DAF时较高，转基因株系中PA（18:1_18:1）降低而PA（18:2_24:1）升高。DAG和TAG水平在40 DAF时较高，转基因株系中DG（18:1_18:1）、DG（16:0_18:1）、DG（18:1_18:3）和DG（18:1_18:2）显著降低，而DG（18:1_22:1）、DG（16:0_22:1）、DG（18:3_22:1）、DG（18:2_22:1）、TG（24:1_18:2_18:3）、TG（24:1_18:3_18:3）和TG（16:0_18:1_24:1）显著增加。

6. NA合成与积累途径中的关键DEGs

在所有样本中，内源KCS基因KCS6和KCS19在植物生长发育过程中均下调，但转基因株系中下调更为显著。内源GPAT6在种子发育过程中上调，且在转基因株系中上调更为显著。内源LPAAT1在转基因株系中显著上调，而在非转基因株系中下调。此外，多个细胞器RNA编辑因子（MORFs）、病程相关蛋白5（PR5）和伴侣蛋白DnaJ C76（DJC76）编码基因在转基因株系中显著上调。

7. qRT-PCR验证

研究人员对3个参与NA合成与积累途径的DEGs进行qRT-PCR验证，结果与RNA-Seq数据一致，证实了转录组分析的可靠性。

四、讨论与结论

讨论部分指出，NA积累伴随着能量代谢、应激响应、光合作用和次生代谢等多方面的协调转录调整。在底物竞争与KCS基因反馈调控方面，异源CgKCS可能对C18:1-CoA具有更高亲和力，有效竞争内源KCS酶，底物耗竭结合高VLCFA积累引发的转录反馈回路，导致内源延长机制被抑制。在酰基协同转运方面，内源GPAT6的上调可能促进包括NA在内的新合成VLCFA掺入甘油酯骨架；LPAAT1的上调则确保多不饱和脂肪酸有效掺入sn-2位，以维持高强度NA积累期间膜完整性。在稳态与应激适应方面，MORF、PR5和DJC76的协同上调表明代谢工程可能触发保护机制，帮助维持细胞稳态。

研究结论表明，转基因B. napus中NA的高水平生产不仅依赖于异源基因活性，还依赖于内源酰基转移酶（如GPAT6和LPAAT1）的协同响应以及对竞争性延长途径（如KCS6和KCS19）的抑制。这些结果为油料作物的代谢可塑性提供了新的机制见解，并确定了进一步优化稀有脂肪酸生产的关键遗传靶点，为NA的规模化生产奠定了理论基础。