甲羟戊酸(mevalonate,MVA) 为甲羟戊酸途径(mevalonate pathway,MVP)的中间产物,是酵母细胞代谢过程中一种重要的有机酸和合成具有高附加值的次生代谢产物的前体化合物[1-2]。甲羟戊酸激酶(mevalonate kinase,MVK)是MVA代谢中关键催化酶和控制整个代谢途径的限速酶之一,也是萜类和醌类物质生物合成过程中的重要调控位点[3-5]。
在前期研究中利用低能离子注入介导药用植物乌拉尔甘草(Glycyrrhiza uralensis)基因组DNA 随机转化酿酒酵母菌,获得了遗传稳定的产萜类和醌类化合物的重组酵母菌N6076[6-7]。本文以重组菌株N6076及其原始菌株Kh08 为研究对象,通过转录组测序和MVA的液相色谱-串联质谱法(liquid chromatograph-mass spectrometer/mass spectrometer,LC-MS/MS)定量测定,以期获得2株酵母菌的差异表达基因的功能,并对MVK的表达与MVA的代谢进行关联分析,为进一步理解酵母菌MVA的代谢及调控提供依据,为后续的外源基因在酵母菌基因组中的整合研究奠定基础。
酿酒酵母(Saccharomyces cerevisiae)N6076来源于低能氮离子注入介导外源DNA随机转化S.cerevisiae Kh08获得遗传稳定的重组酵母菌,由本中心保存。
酵母菌株的培养基与培养方法参考文献[8]。
用于转录组测序(RNA-seq)的样本制备及其RNA-seq方法、差异表达基因的表达量及其GO功能和KEGG代谢通路分析方法均参考文献[9-10]的方法。
酵母菌胞内MVA测定样本的制备方法与MVA的LC-MS/MS测定方法参考文献[8]。
2株酵母菌各发酵时间点及相同发酵时间点的差异表达基因的上调和下调分析结果如图1所示。与 0 h相比,重组菌株 N6076发酵至48和 96 h时的上调基因数量分别为 1186 和 1069个,而下调基因的数量分别为 1681和 1377个;与48 h相比,发酵至96 h时上调与下调基因数量分别为 499和291个,而原始菌株 Kh08分别为523 和588个,说明重组菌株N6076在96 h/48 h时的基因表达的总体活跃程度低于原始菌株Kh08。
t1-处理组;t2-对照组
图1 重组酵母菌N6076 与其原始菌株Kh08不同发酵时间点的差异表达基因
Fig.1 Differentially expressed genes between different time points of fermentation of recombinant yeast strain N6076 and original strain Kh08
使用GOSeq方法对2株酵母菌的差异表达基因进行GO功能分析结果表明,与原始菌株Kh08相比,重组菌株N6076差异表达基因的GO 功能中的生物学过程(biological process,BP)新增14个,缺失1个;细胞组分(cellular component,CC)新增2个,缺失1个;分子功能(molecular function,MF)新增2个,缺失1个(图2)。重组菌株N6076新增加的14个BP功能涉及213个差异表达基因(表1),缺失的1个BP功能为NADH再生(GO:0006735),涉及9 个差异表达基因;在CC功能方面,重组菌株N6076新增了染色质沉默复合物(GO:0005677)和U6 snpRNP(GO:0005688)功能,分别涉及357和372 个差异表达基因;缺失的CC功能为Cul7环泛素连接酶复合物(GO:0031467),涉及6个差异表达基因;在MF功能方面,重组菌株N6076新增了肌动蛋白丝结合(GO:0051015)和蛋白酶活性(GO:0004291)功能,分别涉及22和176个差异表达基因,缺失了1个ATP酶活性功能(GO:0016887),涉及10个差异表达基因。
表1 重组菌株N6076新增的BP功能
Table 1 New added BP functions of recombinant yeast strain N6076
功能ID号描述涉及的差异表达基因数量/个GO∶0099518细胞骨架运输7GO∶0008085光适应15GO∶0000056从细胞核输出核糖小体14GO∶1901528过氧化氢介导的气孔运动信号通路14GO∶00061306-磷酸果糖激酶(还原型)16GO∶0001764迁移12GO∶0046903分泌17GO∶0030635衍生物代谢过程11GO∶0051131伴侣蛋白复合物装配24GO∶0000468三顺反子rRNA转录体产生成熟的3′-端LSUrRNA12GO∶0006495末端O-糖基化14GO∶0042421生物激素合成过程24GO∶0000482四顺反子rRNA转录物中5srRNA的成熟16GO∶0000474四顺反子rRNA转录物中ssurRNA的成熟17
图2 重组菌株N6076与其原始菌株 Kh08差异表达基因的GO功能
Fig.2 The most enriched GO Terms of differential expressed genes of recombinant yeast strain N6076 and original strain Kh08
对相同发酵时间点2株菌株的差异表达基因功能进行分析,结果表明,发酵至48 h,2株酵母菌均包含8种上调表达基因功能,但在下调表达基因功能中,重组菌株N6076与其原始菌株Kh08分别有36和23种;发酵至96 h,2株酵母菌均包含1种上调表达基因功能(DNA ADP核糖基化),但重组菌株N6076 未富集到下调表达差异基因,而原始菌株Kh08 的下调表达基因为 NADH再生功能。
上述差异表达基因GO功能分析结果表明,重组菌株N6076在BP、CC和MF功能方面均比原始菌株Kh08丰富,这有利于重组菌株次生代谢产物的合成。
KEGG代谢通路分析结果表明,2株酵母菌的KEGG注释信息主要包括代谢、遗传信息过程、细胞过程和有机体系等功能,其中,重组菌株 N6076差异表达基因涉及111条KEGG代谢途径,比原始菌株Kh08新增了13 条代谢途径(表2),涉及77个差异表达基因。
表2 重组菌株N6076新增KEGG代谢途径
Table 2 New added KEGG metabolic pathway of recombinant yeast strain N6076
编号通路差异表达基因数/个通路ID号1脂肪酸降解9ko000712丙氨酸、天冬氨酸和谷氨酸代谢7ko002503半胱氨酸与蛋氨酸代谢7ko002704赖氨酸生物合成4ko003005氨基乙磺酸和氨乙基亚磺酸代谢3ko004306甘露糖型o-聚糖的生物合成3ko005157甲醇代谢6ko006808固碳途径7ko007209萜类骨架生物合成5ko0090010MAPK信号通路4ko0401611自噬-其他11ko0413612有丝分裂8ko0413713转录失调3ko05202
萜类骨架生物合成(KEGG ID:ko00900)是重组菌株 N6076新增的一条重要的代谢途径,为其萜类和醌类等次生代谢产物的生物合成提供了基本的C骨架[11]。在自然界生物体内,萜类物质是生物次生代谢产物的一大类别,并且对于大多数真核细胞、少数原核细胞和高等植物而言,其MVP是合成萜类物质的主要途径。1958年CONRAD和 FEODOR在动物和酵母中发现MVP,且此途径主要存在于细胞质中,又称为细胞质途径[12]。2分子乙酰辅酶 A 在乙酰乙酰辅酶 A 硫解酶的作用下形成乙酰乙酰辅酶 A,然后在 3-羟基-3-甲基戊二酰辅酶 A 合成酶(HMGS)的作用下形成 HMG-CoA[13-14],进而在 3-羟基-3-甲基戊二酰辅酶 A 还原酶 (HMGR)的催化下形成不可逆的MVA[14-16],NADPH 为该反应的供氢体;MVA在MVK、磷酸甲羟戊酸激酶(MPK)作用下磷酸化,形成甲羟戊酸5-焦磷酸(MVAP)和甲羟戊酸焦磷酸脱羧酶(MVAPP);最后,MVAPP在 5-焦磷酸甲羟戊酸脱羧酶(MDC)作用下脱羧形成异戊二烯焦磷酸(IPP),从而进入次生代谢途径,合成相关的萜类、醌类及甾体类物质,与表2中的萜类骨架生物合成代谢途径相符合。
乙酰辅酶A是MVP的前体物质,对于萜类和醌类等次生代谢物质的生物合成起到总量控制的作用,改变乙酰辅酶A的浓度会直接影响次生代谢物质的产量。近年来有研究证实,通过调控乙酰辅酶A的基因表达可以增加MVP中生成乙酰辅酶A的含量,进而提高次生代谢物质在目标载体中的产量[17]。
原始菌株Kh08中调控乙酰辅酶A基因的表达量在0~48 h逐渐增加,48~96 h 降低。而重组菌株N6076中调控乙酰辅酶A基因的表达量在整个发酵过程中呈上升态势,其表达量在 0 h 最低,之后逐渐增加,这为其新增的萜类骨架生物合成提供了前体物质,有利于次生代谢产物的合成。
由于MVA的生成是一个不可逆过程,因此,HMGR (EC∶1.1.1.34)被认为是 MVP中的第1个限速酶[18]。HMGR这个酶的基因的差异表达分析结果表明,重组菌株N6076的HMGR这个酶的基因在48~96 h的KEGG表达中呈下调趋势,涉及差异表达基因5个;而在原始菌株Kh08的KEGG代谢途径中未富集到相关酶基因。这进一步表明,重组菌株N6076中MVA的代谢能力有所增强。
MVK是MVP中的限速酶之一[19-20],通过对2株酵母菌MVK基因的表达分析,获得了MVK基因随发酵时间的差异表达趋势(图3-a)。同时利用LC-MS/MS对2株酵母菌不同发酵时间的胞内MVA进行了定量测定,结果如图3-b所示。
图3 重组菌株N6076与其原始菌株Kh08各发酵时间点MK基因表达矩阵与MVA产量
Fig.3 MK gene expression matrix and MVA yield of recombinant yeast strain N6076 and original strain Kh08 at different fermentation time points
MVK将ATP-γ位上的一个磷酸基团转移到MVA第5位的羟基上形成MVAP,并释放ADP[21-22],其酶活性对次生代谢产物的合成速率具有重要的影响作用,并对次生代谢产物的产量有重要影响[23]。重组菌株 N6076 的MVK基因表达量在整个发酵过程中呈起伏上升趋势,0~48 h降低,48~96 h逐渐增加,MVA产量则逐渐积累,至96 h达到峰值。原始菌株Kh08的MVK基因表达量在48~96 h达到峰值,MVA产量在96 h达到峰值。2株酵母菌的MVK基因表达量与MVA产量的变化趋势一致。
重组菌株N6076新增1条萜类骨架生物合成途径(表2),势必要消耗更多的MVA,因此,发酵48和96 h时,MVA的积累量均低于原始菌株Kh08(图3-b),说明MVA的代谢速度略快于原始菌株Kh08。
差异表达基因的分析结果表明,重组酵母菌N6076发酵至48和96 h时的上调基因数量均低于下调基因的数量;但96 h/48 h时,上调基因数量>下调基因数量。
GO功能分析结果表明,重组酵母菌N6076新增了14项BP功能、2项CC功能和2项MF功能,涉及1 140个差异表达基因。
KEGG代谢通路分析结果表明,重组酵母菌N6076新增了13条代谢通路,涉及77个差异表达基因,其中新增的萜类骨架生物合成途径涉及5个差异表达基因,为其萜类和醌类等次生代谢产物的生物合成提供了基本的C骨架。
MVK基因表达与MVA测定结果表明,在整个发酵过程中MVK基因表达量与MVA产量的变化趋势一致。本研究为进一步认识重组酵母菌N6076的基因表达与MVA代谢调控提供了依据。
[1] 张园园. 大肠杆菌甲羟戊酸途径的构建与调控[D].济南:山东大学,2009.
[2] ENDO A.The discovery and development of HMG-CoA reductase inhibitors[J].Journal of Lipid Research,1992,5(3):1 569-1 582.
[3] 赵乐, 马利刚,杨方方,等.独行菜磷酸甲羟戊酸激酶LaPMK基因克隆、生物信息学分析及原核表达[J].中草药,2016,47(17):3 087-3 093.
[4] GARCIA D E,KEASLING J D,MANDE S C.Kinetics of phosphomevalonate kinase from Saccharomyces cerevisiae[J].Plos One,2014,9(1):e87 112.
[5] CHAPPELL J.Biochemistry and molecular biology of the isoprenoid biosynthetic pathway in plants[J].Annual Review of Plant Physiology and plant Molecular Biology,1995,46:521-547.
[6] FAN Y H,JIN X,MAO P H.Preparation of red component in the secondary metabolites of recombinant N6076 and its GC-MS detection[J].Agricultural Biotechnology,2012,1(4):58-60
[7] MAO P H,LYU J,CAI C L.Proteomic changes of recombinant yeast:Pharmaco-industrial potential[J].Journal of Biological Regulators and Homeostatic Agents,2015,29(3):569-578.
[8] 欧科, 陈福欣,王婷,等.基于代谢组学方法的2 株酿酒酵母菌胞内差异代谢产物分析[J].食品与发酵工业,2020,46(4),39-44.
[9] 唐朝,张寒玉,王婷,等.基于转录组测序的异常汉逊酵母菌不同发酵时期差异表达基因功能分析[J].食品与发酵工业,2019,45(4):1-6.
[10] 张寒玉, 唐朝,冯光文,等.离子束重组异常汉逊酵母菌 Ar_Han045的转录组学分析及谷胱甘肽代谢途径富集[J].食品与发酵工业,2019,45(10):10-15.
[11] YAMAZAKI M, MOCHIDA K, ASANO T, et al.Coupling deep transcriptome analysis with untargeted metabolic profiling in ophiorrhiza pumila to further the understanding of the biosynthesis of the anti-cancer alkaloid camptothecin and anthraquinones[J].Plant and Cell Physiology,2013,54(5):686-696.
[12] LYNENF F,EGGERER H,HENNING U,et al.Farnesyl-pyrophat and 3-Methyl-Δ3-butenyl-1-pyrophosphat,die biologischen Vorstufen des Squalens.Zur Biosynthese der Terpene,Ⅲ[J].Angewandte Chemie,2015,70(24):378-742.
[13] RUDNEY H.The biosynthesis of beta-hydroxy-beta-methylglutaric acid[J].Journal of Biological Chemistry,1957,227(1):363-377.
[14] 王志标, 王伟,程克棣,等.甲羟戊酸途径代谢酶 HMGL 研究概述[J].中国医药生物技术,2014(1):48-52.
[15] DURR IF,RUDNEY H.The reduction of beta-hydroxy-beta-methyl-glutaryl coenzyme A to mevalonic acid[J].Journal of Biological Chemistry,1960,235:2 572-2 578.
[16] 李欢, 张娜,李依民,等.利用转录组测序挖掘掌叶大黄蒽醌类生物合成相关基因[J].药学学报,2018,53(11):152-161.
[17] 汪婉宜,申业,黄璐琦,等.萜类合成甲羟戊酸途径的关键基因的克隆及表达分析[J].中国食物与营养.2017,17(9):24-28.
[18] 罗永明, 刘爱华,李琴,等.植物萜类化合物的生物合成途径及其关键酶的研究进展(续完)[J].江西中医学院学报,2003,15(2):46-49.
[19] 张雯雯, 曾日中,杨礼富,等.异戊二烯生物合成研究进展[J].微生物学杂志,2016,36(6):98-103.
[20] LIAO Z H,CHEN M,GONG Y F,et al.Isopenoid biosynthesis inplants:pathway,genes,regulation and metabolic engineering[J].Journal of Biological Sciences,2006,6(1):209-219.
[21] 王宝莲, 樊庆琦,李永波,等.甲羟戊酸激酶基因研究进展[J].中国农业科技导报,2011,13(3):17-25.
[22] 张晓东, 李彩霞,赵静,等.滇龙胆甲羟戊酸激酶基因的克隆与表达分析[J].2015,19(3):196-202.
[23] 李亮, 尚晓东,谭琦.担子菌中甲羟戊酸途径影响萜类物质合成的研究[J].中国农业科技导报,2014(3):53-61.