编辑推荐：
本文推荐理由：本研究揭示了在豆科植物山黧豆中，微小RNA miR172f通过靶向负调控β-ODAP合成酶基因LsBAHD3的表达，从而精准调控神经活性物质β-ODAP（β-N-oxalyl-L-α,β-diaminopropionic acid）生物合成的新机制。敲低miR172f可显著提高山黧豆毛状根中β-ODAP含量，并引发包括蛋白质泛素化、半胱氨酸和蛋氨酸代谢等通路的全局转录组重编程。该发现为通过基因工程手段定向调控β-ODAP含量、平衡其神经毒性与药用潜力提供了关键靶点。
引言
草豌豆（Lathyrus sativus L.）因其对多种生物和非生物胁迫具有广泛抗性，在许多干旱和半干旱地区被广泛种植，用作人类食物或饲料。然而，将其作为日常饮食的主要组成部分面临挑战，因为它会积累具有神经活性的毒素β-N-草酰基-L-α,β-二氨基丙酸（β-ODAP）。β-ODAP以其剂量和应用背景依赖的双重效应而闻名：过量摄入可诱发神经退行性神经性山黧豆中毒，而在特定条件下又表现出神经保护和伤口愈合等药理作用。该化合物也存在于豌豆（Pisum sativum）和多种传统中草药如三七（Panax notoginseng）中。因此，调控β-ODAP的生物合成对于安全利用含β-ODAP物种或促进其分离纯化具有重要意义。已知包括丝氨酸乙酰转移酶（SAT）、β-氰基丙氨酸合酶（β-CAS）、BAHD3酰基转移酶（BAHD3）和酰基激活酶3（AAE3）在内的多个基因参与β-ODAP的生物合成。其中，LsBAHD3被证实是β-ODAP合成酶（BOS），其能以草酰辅酶A为供体，催化L-DAP草酰化形成β-ODAP。
微小RNA（miRNA）在调控植物发育、抗逆性和生物活性化合物合成中发挥重要作用。例如，miR160、miR167、miR169和miR171参与氮胁迫响应通路；miR395家族成员是硫代谢及其调控网络的常见组分。而miR172家族成员在调节植物生长、苯丙氨酸和酚酸等次级代谢产物的生物合成以及介导干旱、盐度和重金属暴露等环境胁迫响应中起着关键作用。鉴于β-ODAP含量与氮、硫代谢及多种胁迫密切相关，研究参与β-ODAP生物合成的miRNA具有重要意义。本研究旨在探究miR172f在β-ODAP生物合成中的调控作用。
材料与方法
本研究采用了一系列分子生物学和组学技术。首先，通过生物信息学网站psRNATarget预测了可能靶向LsBAHD3基因的miRNA。利用已构建的过表达LsBAHD3的豌豆毛状根系（OE LsBAHD3–13），通过RT-PCR验证了转基因系的成功构建，并通过RT-qPCR检测了其中miR172f的表达水平。为验证miR172f与LsBAHD3的靶向关系，将LsBAHD3或前体miR172f分别克隆到载体中，通过农杆菌介导的瞬时共表达实验在本氏烟叶片中进行荧光淬灭分析。
为了在原生系统中研究miR172f的功能，通过短串联靶标模拟物技术在山黧豆毛状根中敲低了miR172f的表达，并利用RT-qPCR验证了敲低效率及LsBAHD3基因的表达变化。使用高效液相色谱法测定了敲低系中β-ODAP的含量变化。此外，对miR172f敲低系和阴性对照进行了转录组测序分析，通过主成分分析、差异表达基因筛选、基因本体论分析和京都基因与基因组百科全书通路富集分析，系统探究了敲低miR172f引起的全局基因表达变化。进一步，通过基因集富集分析和加权基因共表达网络分析，挖掘了与β-ODAP生物合成高度相关的基因模块和生物学通路。
结果
1. 过表达LsBAHD3的豌豆毛状根中miR172f表达水平降低
对转基因豌豆毛状根系OE LsBAHD3–13的RT-PCR确认显示出一条大小约为1300 bp的清晰条带，表明LsBAHD3基因成功过表达。随后，RT-qPCR分析发现，与对照相比，OE LsBAHD3–13毛状根中miR172f的表达水平显著降低。这一现象暗示了miR172f与LsBAHD3之间可能存在反馈调控关系。
2. 证实miR172f靶向LsBAHD3基因
生物信息学预测显示，miR172f与LsBAHD3基因的mRNA序列存在六个连续互补的结合位点，强烈提示miR172f可能靶向LsBAHD3。为验证此相互作用，研究者将携带pTF486-BAHD3-eGFP（表达LsBAHD3与绿色荧光蛋白融合蛋白）和pOT2-miR172f（表达miR172f前体）载体的农杆菌共同浸润本氏烟叶片。与仅表达融合蛋白的阳性对照相比，共表达miR172f显著淬灭了GFP荧光信号，这直观地证明了miR172f在植物体内能够与LsBAHD3的mRNA发生相互作用，并可能导致其转录后水平的沉默。
3. 敲低miR172f提高山黧豆毛状根中β-ODAP含量
成功获得了转染空载体和STTM-miR172f载体的山黧豆毛状根。RT-qPCR分析证实，与对照相比，STTM-miR172f转基因毛状根中miR172f的表达水平显著降低，而LsBAHD3基因的表达水平则显著上调，这进一步证实了miR172f对LsBAHD3的靶向负调控作用。HPLC检测结果显示，在miR172f敲低系STTM-1和STTM-5中，β-ODAP含量相比对照均大幅增加。这一系列结果直接证明，通过敲低miR172f来解除其对LsBAHD3的抑制，可以有效提高β-ODAP的合成。
4. miR172f敲低导致山黧豆毛状根转录组重编程
转录组分析显示，miR172f敲低系与对照之间存在5415个差异表达基因，其中3249个上调，2166个下调。主成分分析能够将两组样本清晰区分。基因本体论分析表明，差异基因显著富集在防御反应、DNA模板转录调控、蛋白质泛素化等生物过程，以及细胞核、质膜等细胞组分，其分子功能涉及ATP结合、金属离子结合、DNA结合转录因子活性、丝氨酸/苏氨酸蛋白激酶活性等。
通路富集分析（KEGG）显示，差异基因主要富集在植物激素信号转导、MAPK信号通路、苯丙素生物合成、糖酵解/糖异生、淀粉和蔗糖代谢、以及半胱氨酸和蛋氨酸代谢等通路。值得注意的是，多个先前报道与β-ODAP生物合成相关的通路也被显著富集，包括丙氨酸、天冬氨酸和谷氨酸代谢、氮代谢、半胱氨酸和蛋氨酸代谢、多种植物次生代谢物的生物合成、甘氨酸、丝氨酸和苏氨酸代谢、硫代谢以及泛酸和辅酶A生物合成。
基因集富集分析进一步揭示，β-ODAP含量的增加与多个上调的基因集正相关，包括酶调节剂活性、蛋白质异二聚化活性、焦磷酸酶活性、生长素极性运输的调控等；同时与一些下调的基因集负相关，如苯丙素生物合成过程、次级活性硫酸盐跨膜转运蛋白活性、蛋白质磷酸酶抑制剂活性、茉莉酸代谢过程等。
加权基因共表达网络分析构建了29个共表达模块。模块-性状关联分析发现，模块ME1与“山黧豆毛状根中miR172f敲低”这一性状呈极显著正相关，而模块ME2则呈极显著负相关。有趣的是，已被确认为β-ODAP生物合成关键基因的LsBAHD3就位于ME1模块中，并且与miR172f敲低性状呈高度正相关。这表明ME1模块可能是一个与β-ODAP生物合成功能协同的基因集合。
讨论
先前研究表明，LsBAHD3-LsAAE3模块在山黧豆和豌豆的β-ODAP生物合成中发挥保守功能。在本研究中，过表达LsBAHD3的豌豆毛状根中miR172f表达下调，以及敲低miR172f的山黧豆毛状根中LsBAHD3表达上调和β-ODAP含量增加，共同证实了miR172f是一个通过靶向LsBAHD3来调控β-ODAP水平的新型调节因子。这为理解β-ODAP合成的转录后调控机制提供了新视角。
硫代谢途径是β-ODAP生物合成的核心生化支架，而miR395家族成员通过靶向ATP硫酸化酶等基因参与硫代谢调控。本研究的转录组结果与前期发现高度一致：蛋白质泛素化、半胱氨酸和蛋氨酸代谢等过程在miR172f敲低后发生显著变化。基因集富集分析指出，酶调节剂活性和蛋白质异二聚化活性与β-ODAP生物合成正相关。这具有重要生物学意义，因为参与硫代谢和β-ODAP合成的关键酶——丝氨酸O-乙酰转移酶（SAT）的活性，通常通过与O-乙酰丝氨酸(硫醇)-裂解酶相互作用形成异源寡聚的半胱氨酸调节复合物来调控。此外，钙依赖蛋白激酶LsCDPK-SK5可通过磷酸化LsSAT2来减弱半胱氨酸对LsSAT2活性的反馈抑制。这些发现暗示，miR172f的调控作用可能通过影响酶复合物的形成或修饰状态，间接而广泛地影响β-ODAP的合成网络。
结论
综上所述，本研究证实miR172f是β-ODAP生物合成的一个新型调控因子，其通过靶向β-ODAP合成酶基因LsBAHD3发挥作用。在山黧豆毛状根中敲低miR172f，通过影响蛋白质泛素化、半胱氨酸和蛋氨酸代谢等多种生物学过程和代谢通路，导致β-ODAP含量增加。同时，酶调节剂活性、蛋白质异二聚化活性和蛋白质磷酸酶抑制剂活性等与β-ODAP生物合成表现出强相关性。这些结果表明，敲低miR172f引发了山黧豆毛状根中涉及5415个基因的全局转录重编程。该研究为通过遗传操纵手段精准调控山黧豆中β-ODAP含量提供了新的分子靶点，并为深入理解β-ODAP生物合成的遗传机制，特别是其合成酶基因的表达调控网络，提供了宝贵见解。
