大豆GmPEAMT1/GmPLMT双敲除生殖表型与种子营养重塑的多组学分析

《Agronomy》：Multi-Omics Analysis of Reproductive Phenotypes and Seed Nutritional Remodeling in Soybean GmPEAMT1/GmPLMT Double Knockouts

　　

研究人员构建了大豆GmPEAMT1/GmPLMT双敲除突变体（Dko），以研究在实验条件下磷脂相关基因的破坏如何影响生殖性状、种子营养成分及相关分子响应。Dko植株在光合参数和生殖结构方面表现出协同变化：花粉粒和胚珠数量减少，而花粉和胚珠大小相关性状增加，同时成熟种子数和百粒重得以维持。发育中的Dko种子表现出较高的蛋白质含量，且所测18种氨基酸中有12种水平升高，尤其是丙氨酸和半胱氨酸。转录组学和实时荧光定量PCR（qRT-PCR）分析进一步表明，氨基酸、泛酸、辅酶A、甘油脂（glycerolipid）和甘油磷脂（glycerophospholipid）相关通路发生协同变化。脂质组学（lipidomics）提供了磷脂相关及其他脂质分子特征变化的补充证据，而种子总脂质含量与对照相比无显著差异。综上，这些结果表明，GmPEAMT1/GmPLMT的破坏与所检查发育阶段的生殖性状变异和种子养分分配变化相关，其中脂质特征变化是更广泛多组学响应的一部分。由于直接代谢通量、蛋白质丰度和酶活性均未被测定，所提出的大豆GmPEAMT1/GmPLMT双敲除种子脂质和氨基酸代谢重塑应被视为一种工作假说。

**研究背景与科学问题**

大豆（*Glycine max* (L.) Merr.）是全球重要的油料与经济作物，油分含量和营养组成决定其应用价值。传统育种和基因工程多采用“数量优先”模式提高总脂含量，但单纯追求油脂积累常导致脂肪酸组成失衡、脂质储存效率低等问题。脂质代谢并非孤立通路，而是与氨基酸代谢等初级代谢密切互作的网络。磷脂酰乙醇胺（PE）和磷脂酰胆碱（PC）是两类最丰富的膜磷脂，其生物合成与甘油二酯（DG）、三酰甘油（TG）的合成共享前体DG；磷酸乙醇胺N-甲基转移酶（PEAMT）和磷脂甲基转移酶（PLMT）是PE甲基化为PC的两条途径中的关键酶。已有研究表明，PE和PC不仅参与维持膜稳态，还可能影响乙酰辅酶A（acetyl-CoA）和TG合成。然而，GmPEAMT1和GmPLMT同时缺失对大豆磷脂相关脂质特征、生殖性状、种子成分及相关分子响应的联合效应仍不清楚。为此，研究人员利用CRISPR/Cas9技术构建双基因敲除大豆，以表征磷脂头基代谢扰动与上述性状的关联。

**研究内容与总体结论**

研究人员在大豆品种“东农50”中同时敲除GmPEAMT1（Glyma_05G246500v4）和GmPLMT（Glyma_08G218900v4），获得3个独立纯合双敲除株系Dko-A、Dko-B和Dko-C。T₂代材料于2024年5月至9月种植于中国长春（43°48′N, 125°24′E）田间；Dko-A用于脂质组学和转录组学分析。研究通过表型观察、脂质组学、转录组学和种子营养测定，发现双敲除与光合参数、花粉和胚珠性状、种子营养组成以及脂质和转录谱变化相关，但成熟种子数和百粒重无显著差异，总脂质含量亦未改变。该研究提出脂质和氨基酸代谢重塑的工作假说，但强调需通过代谢通量、蛋白质丰度和酶活性测定进一步验证。

**主要技术方法**

研究人员采用CRISPR/Cas9双靶点载体pYL-GmPEAMT1/GmPLMT，经农杆菌介导的大豆子叶节转化获得突变体；使用LI-6800光合荧光分析仪测定光合参数；采用I₂-KI染色和扫描电镜观察花粉；通过卵巢切片观察胚珠；使用近红外（NIR）分析仪测定种子水分、蛋白、脂质和氨基酸；采用超高效液相色谱-串联质谱（LC-MS/MS）进行脂质组学；利用RNA-seq和qRT-PCR进行转录组测序与验证，并结合PCA、OPLS-DA、KEGG/GO富集等生物信息学方法分析。

**研究结果**

**3.1 CRISPR/Cas9诱导的GmPEAMT1和GmPLMT定向诱变**

Sanger测序确认3个独立纯合双敲除株系在GmPEAMT1和GmPLMT靶位点均携带独立的编辑事件，包括碱基缺失和替换；突变位于两个基因的编码区，预测会导致移码或局部氨基酸缺失，从而改变编码蛋白的结构。

**3.2 GmPEAMT1/GmPLMT双敲除对大豆生殖发育的影响**

在播种后54 d（始花后约7 d），Dko的胞间CO₂浓度低于对照，净光合速率和总能量支出高于对照，而气孔导度、边界层导度和蒸腾速率与对照相近。I₂-KI染色和扫描电镜观察显示，Dko花粉粒总数和可染色花粉数减少，花粉直径和染色特征改变；卵巢切片显示胚珠数减少，胚珠大小相关性状增加；成熟种子数和百粒重与对照无显著差异。这些结果被界定为形态和染色相关表型，而非生殖适合度的直接证据。

**3.3 通过GmPEAMT1/GmPLMT双敲除重塑大豆种子脂质代谢**

脂质组学分析中，OPLS-DA模型显示Dko-A与对照分离（负离子模式R²Y=0.998，Q²=0.905；正离子模式R²Y=0.999，Q²=0.902）。共鉴定到151个上调、83个下调脂质，其中甘油磷脂类数量最多；TG、Cer、PE和DG变化突出，Cer(t18:0/18:2)显著上调。KEGG富集显示鞘脂代谢、甘油脂代谢、糖基磷脂酰肌醇（GPI）锚定生物合成、甘油磷脂代谢和自噬等通路富集，脂质碳链长度和饱和度分布也发生改变。由于生物学重复较少，该结果被解释为探索性证据。

**3.4 GmPEAMT1/GmPLMT双敲除大豆种子的转录重塑**

转录组分析鉴定到766个上调基因和807个下调基因，PCA显示Dko与对照明显分离。GO富集突出微管形成、形态发生、神经酰胺胆碱磷酸转移酶活性、肌醇磷酸神经酰胺合酶活性、鞘磷脂合酶活性及L-丙氨酸转氨等条目；KEGG富集前三为不饱和脂肪酸生物合成、谷胱甘肽代谢和丙酮酸代谢，同时半胱氨酸与蛋氨酸代谢、氧化磷酸化、异黄酮生物合成、丙氨酸/天冬氨酸/谷氨酸代谢、脂肪酸降解、α-亚麻酸代谢和缬氨酸/亮氨酸/异亮氨酸降解等通路也显著富集。qRT-PCR验证了关键氨基酸和脂质代谢基因的表达趋势。

**3.5 GmPEAMT1/GmPLMT双敲除大豆种子的营养成分与运输功能**

近红外分析显示，Dko种子水分含量降低、蛋白质含量升高，总脂质含量无显著变化；18种氨基酸中12种相对含量显著升高，其中丙氨酸升高24.7%，半胱氨酸升高43.1%。由于氨基酸数据以总蛋白质量为基准，这表示蛋白质归一化的相对氨基酸组成变化，而非单粒种子的绝对氨基酸含量。

**讨论与结论**

讨论部分围绕表型变异、脂质重塑、氨基酸代谢及综合营养重塑展开。研究指出，Dko在单一发育时间点表现出较高的净光合速率和较低的胞间CO₂浓度，同时花粉、胚珠和种子性状发生变化，但这些属于关联观察，不能证明光合变化直接驱动种子能量分配。脂质组和转录组共同提示鞘脂代谢抑制、TG相关特征改变和膜脂相关代谢扰动；C_6:0、C_8:0等短链饱和脂肪酸水平升高，但总脂含量未变，因此是脂质组成变化而非储脂总量增加。氨基酸代谢、泛酸和辅酶A生物合成相关基因表达上调，与蛋白质含量和氨基酸比例升高一致，但该研究未直接测定acetyl-CoA和acyl-CoA通量、蛋白质丰度及酶活性，故所提出的脂质和氨基酸代谢重塑属于工作假说。

研究结论可概括为：GmPEAMT1和GmPLMT缺失与大豆代谢重塑相关，可能影响种子形成、脂质相关性状和营养组成。具体表现为：(1)在温室条件和所检查的发育阶段，Dko的光合参数、花粉/胚珠性状和种子组成存在差异，但该研究未直接检验磷脂代谢变化是否影响碳分配或生殖发育；(2)C_6:0、C_8:0、TG相关脂质和氨基酸相关性状的变化提示脂肪酸重塑、氨基酸相关代谢与种子脂质组成之间可能存在关联，但acetyl-CoA/acyl-CoA可用性、代谢通量和TG合成速率未直接测定；(3)Dko种子蛋白质含量升高、蛋白质归一化氨基酸组成改变，但因数据以总蛋白为基准，说明是相对组成变化而非每粒种子绝对氨基酸量增加。双敲除并未完全消除PC代谢，残余PC中的18:2和18:3酰基编辑在种子发育中仍活跃。该研究受限于NIR间接测定、仅以单一代表株系进行组学分析，以及未检测GmPEAMT1和GmPLMT蛋白丰度与酶活性；未来需要通过湿化学方法、多个独立株系、稳定同位素示踪、蛋白定量和酶活性分析进一步验证。同时，花粉萌发、花粉管生长、受精效率、种子萌发和后代表现未测定，因此生殖表型的适应意义仍有待证实。