留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

通过整合转录组与代谢组解析不同类型椰子的脂肪酸调控网络

邓渊 赖军 毛梦迪 张越冉 李淳 杨君 王守创 罗杰

邓渊,赖军,毛梦迪,等. 通过整合转录组与代谢组解析不同类型椰子的脂肪酸调控网络[J]. 热带生物学报,2022, 13(5):478−487. DOI:10.15886/j.cnki.rdswxb.2022.05.008 doi:  10.15886/j.cnki.rdswxb.2022.05.008
引用本文: 邓渊,赖军,毛梦迪,等. 通过整合转录组与代谢组解析不同类型椰子的脂肪酸调控网络[J]. 热带生物学报,2022, 13(5):478−487. DOI:10.15886/j.cnki.rdswxb.2022.05.008 doi:  10.15886/j.cnki.rdswxb.2022.05.008
DENG Yuan, LAI Jun, MAO Mengdi, ZHANG Yueran, LI Chun, YANG Jun, WANG Shouchuang, LUO Jie. Elucidation of the fatty acid regulatory network for different coconut varieties by integrating transcriptome and metabolome[J]. Journal of Tropical Biology, 2022, 13(5): 478-487. doi: 10.15886/j.cnki.rdswxb.2022.05.008
Citation: DENG Yuan, LAI Jun, MAO Mengdi, ZHANG Yueran, LI Chun, YANG Jun, WANG Shouchuang, LUO Jie. Elucidation of the fatty acid regulatory network for different coconut varieties by integrating transcriptome and metabolome[J]. Journal of Tropical Biology, 2022, 13(5): 478-487. doi: 10.15886/j.cnki.rdswxb.2022.05.008

通过整合转录组与代谢组解析不同类型椰子的脂肪酸调控网络

doi: 10.15886/j.cnki.rdswxb.2022.05.008
基金项目: 国家自然科学基金(31960213)
详细信息
    第一作者:

    邓渊(1996−),男,海南大学热带作物学院2019级硕士研究生. E-mail:yuandeng@hainanu.edu.cn

    通信作者:

    王守创(1989−),男,教授. 研究方向:植物次生代谢调控及代谢组学. E-mail:shouchuang.wang@hainanu.edu.cn

    罗杰(1971−),男,教授. 研究方向:植物次生代谢调控及代谢组学. E-mail:jie.luo@hainanu.edu.cn

  • 中图分类号: S 667.4;Q 943

Elucidation of the fatty acid regulatory network for different coconut varieties by integrating transcriptome and metabolome

  • 摘要: 为了探究不同类型椰子中脂质种类和含量的差异,本研究利用非靶向代谢组学方法对海南高种和矮种椰子椰肉组织的代谢物进行了分析,共检测到12 579种代谢信号,其中鉴定了564种代谢物,包括氨基酸、脂肪酸、类黄酮等。定量分析表明高种椰子的脂质总体含量高于矮种椰子,其中绿矮椰子脂质含量最低。甘油脂、鞘脂和脂肪酰是不同类型椰子中主要的差异物质。转录组数据结果表明,不同类型椰子中的差异基因主要参与脂肪酸代谢、α-亚麻酸代谢、木质素代谢和苯丙烷代谢。基于表达模式的不同将差异基因划分为20个不同的模块,其中“Melightyellow”模块与差异积累的脂肪酸高度相关。通过分析模块中脂质途径结构基因启动子区的cis-element,发现大多数脂质途径基因的启动子区存在多个MYB家族的结合位点及大量的激素信号响应元件。结果表明,MYB家族的转录因子可能在脂质的合成调控中发挥重要作用,并且激素可能作为一种信号分子参与到脂质的合成调控中。本研究通过组合策略对不同品种椰肉组织脂质差异的分子机理进行分析,并进一步构建了转录因子−结构基因−代谢物的合成调控网络,为后续解析脂质合成调控的分子机制提供了资源,为油脂椰子品种选育提供了依据。
  • 图  2  不同类型椰子品种脂质代谢组的比较

    A:高矮种椰子中的脂质含量比较。DC表示矮种椰子,HT表示海南高种;B:高种、黄矮、绿矮脂质含量的比较。YD表示黄矮椰子,HT表示海南高种,GD表示绿矮椰子;C:基于脂质数据的PCA结果;:脂质数据的热图聚类结果;E:不同类型椰子品种间差异积累的脂质数目;F:不同类型椰子品种间差异积累的脂质类别。

    图  3  不同品种椰子差异表达基因的鉴定与功能富集

    A:基于转录组数据的PCA结果;B:不同类型椰子品种两两比较间的差异基因数目;C: 差异表达基因的KEGG富集通路;D:差异基因的GO富集通路。

    图  4  加权共表达分析鉴定与脂肪酸合成相关的基因

    A:模块基因与差异脂肪酸的相关性热图;B:“Melightyellow”模块中与脂肪酸合成相关的基因调控网络;C:LPC 18:1在LPCAT基因的作用下与油酸结合生成Phosphatidylcholine (18:1/18:1);D:不同椰子品种中20-Oxoarachidonic acid (Cnm05349)和Phosphatidylcholine (18:1/18:1)的含量分布。HT表示高种椰子,YD表示黄矮椰子,GD表示绿矮椰子;E: ADHLPCAT基因在不同椰子品种中的转录水平;F: ADHLPCAT基因的qRT-PCR结果。

    图  5  高矮种椰子脂肪酸合成途径基因表达水平及代谢产物积累水平的比较

    不同颜色块代表不同椰子品种比较间代谢物的倍数变化(以对数表示)。左侧色块代表 GD vs HT,右侧色块代表 YD vs HT。蓝色箭头表示该基因在高种椰子中的表达水平高于矮种椰子,红色箭头表示该基因在矮种椰子中的表达水平高于高种椰子。

  • [1] GUNN B F, BAUDOUIN L, OLSEN K M. Independent origins of cultivated coconut (Cocos nucifera L.) in the old world tropics [J]. PLoS One, 2011, 6(6): 1 − 8.
    [2] 吴翼. 椰子SSR分子标记的开发 [D]. 海口: 海南大学, 2008.
    [3] WANG S C, XIAO Y, ZHOU Z W, et al. High-quality reference genome sequences of two coconut cultivars provide insights into evolution of monocot chromosomes and differentiation of fiber content and plant height [J]. Genome Biology, 2021, 22(1): 2 − 25. doi:  10.1186/s13059-020-02226-6
    [4] PERERA L, BAUDOUIN L, MACKAY I. SSR markers indicate a common origin of self-pollinating dwarf coconut in South-East Asia under domestication [J]. Scientia Horticulturae, 2016, 211: 255 − 262. doi:  10.1016/j.scienta.2016.08.028
    [5] TEULAT B, ALDAM C, TREHIN R, et al. An analysis of genetic diversity in coconut populations from across the geographic range using sequence-tagged microsatellites (SSRs) and AFLPs [J]. Theoretical and Applied Genetics, 2000, 100(5): 764 − 771. doi:  10.1007/s001220051350
    [6] PERERA L, RUSSELL J, PROVAN J, et al. Studying genetic relationships among coconut varieties/populations using microsatellite markers [J]. Euphytica, 2003, 132(1): 121 − 128. doi:  10.1023/A:1024696303261
    [7] KAMARAL L C, PERERA S A, PERERA K L, et al. Genetic diversity of the Sri Lanka yellow dwarf coconut form as revealed by microsatellite markers [J]. Tropical Agricultural Research, 2015, 26(1): 131 − 139. doi:  10.4038/tar.v26i1.8078
    [8] LAURELES L, RODRIGUEZ F, REANO C, et al. Variability in fatty acid and triacylglycerol composition of the oil of coconut (Cocos nucifera L. ) hybrids and their parentals [J]. Journal of Agricultural and Food Chemistry, 2002, 50(6): 1581 − 1856. doi:  10.1021/jf010832w
    [9] PERERA L, RUSSELL J, PROVAN J, et al. Use of microsatellite DNA markers to investigate the level of genetic diversity and population genetic structure of coconut (Cocos nucifera L.) [J]. Genome, 2000, 43(1): 15 − 21. doi:  10.1139/g99-079
    [10] XIAO Y, XU P, FAN H, et al. The genome draft of coconut (Cocos nucifera L.) [J]. GigaScience, 2017, 6(11): 1 − 11.
    [11] KUMAR S N. Variability in coconut (Cocos nucifera L. ) germplasm and hybrids for fatty acid profile of oil [J]. Journal of Agricultural and Food Chemistry, 2011, 59(24): 13050 − 13058. doi:  10.1021/jf203182d
    [12] REYNOLDS K B, CULLERNE D P, El T A, et al. Identification of genes involved in lipid biosynthesis through de novo transcriptome assembly from cocos nucifera developing endosperm [J]. Plant Cell Physiology, 2019, 60(5): 945 − 960. doi:  10.1093/pcp/pcy247
    [13] BABU A S, VELUSWAMY S K, ARENA R, et al. Virgin coconut oil and its potential cardioprotective effects [J]. Postgraduate Medicine, 2014, 126(7): 76 − 83. doi:  10.3810/pgm.2014.11.2835
    [14] MA Z F, LEE Y Y. Virgin coconut oil and its cardiovascular health benefits [J]. Natural Product Communications, 2016, 11(8): 1151 − 1152.
    [15] BURDOCK G A, CARABIN I G. Safety assessment of myristic acid as a food ingredient [J]. Food Chemical Toxicology, 2007, 45(4): 517 − 529. doi:  10.1016/j.fct.2006.10.009
    [16] HENRY G, MOMIN R, NAIR M, et al. Antioxidant and cyclooxygenase activities of fatty acids found in food [J]. Journal of Agricultural and Food Chemistry, 2002, 50(8): 2231 − 2234. doi:  10.1021/jf0114381
    [17] KHALIL A S M, GIRIBABU N, YELUMALAIEL S, et al. Myristic acid defends against testicular oxidative stress, inflammation, apoptosis: Restoration of spermatogenesis, steroidogenesis in diabetic rats [J]. Life Sciences, 2021, 278: 119605. doi:  10.1016/j.lfs.2021.119605
    [18] FAMUREWA A C, MADUAGWUNA E K, FOLAWIYO A M, et al. Antioxidant, anti-inflammatory, and antiapoptotic effects of virgin coconut oil against antibiotic drug gentamicin-induced nephrotoxicity via the suppression of oxidative stress and modulation of iNOS/NF-ĸB/caspase-3 signaling pathway in Wistar rats [J]. Journal of Food Biochemistry, 2020, 44(1): 1 − 10.
    [19] ILLAM S, NARAYANANKUTTY A, RAGHAVAMENON A. Polyphenols of virgin coconut oil prevent pro-oxidant mediated cell death [J]. Toxicology Mechanisms and Methods, 2017, 27(6): 442 − 450. doi:  10.1080/15376516.2017.1320458
    [20] NAIRS S, MANALIL J, RAMAVARMA S, et al. Virgin coconut oil supplementation ameliorates cyclophosphamide induced systemic toxicity in mice [J]. Human & Experimental Toxicology, 2016, 35(2): 205 − 212.
    [21] LAW K, AZMAN N, OMAR E, et al. The effects of virgin coconut oil (VCO) as supplementation on quality of life among breast cancer patients [J]. Lipids in Health and Disease, 2014, 13: 139. doi:  10.1186/1476-511X-13-139
    [22] FAMUREWA A C, AJA P M, NWANKWO O E, et al. Moringa oleifera seed oil or virgin coconut oil supplementation abrogates cerebral neurotoxicity induced by antineoplastic agent methotrexate by suppression of oxidative stress and neuro-inflammation in rats [J]. Journal of Food Biochemistry, 2019, 43(3): e12748.
    [23] VERMA P, NAIK S, NANDA P, et al. In vitro anticancer activity of virgin coconut oil and its fractions in liver and oral cancer cells [J]. Anti-Cancer Agents in Medicinal Chemistry, 2019, 19(18): 2223 − 2230.
    [24] DEEN A, VISVANATHAN R, WICKRAMARACHCHI D, et al. Chemical composition and health benefits of coconut oil: an overview [J]. Journal of The Science of Food and Agriculture, 2021, 101(6): 2182 − 2193. doi:  10.1002/jsfa.10870
    [25] KUNST L, SAMUELS L. Plant cuticles shine: advances in wax biosynthesis and export [J]. Current Opinion in Plant Biology, 2009, 12(6): 721 − 727.
    [26] MACKOVA J, VASKOVA M, MACEK P, et al. Plant response to drought stress simulated by ABA application: Changes in chemical composition of cuticular waxes [J]. Environmental and Experimental Botany, 2013, 86: 70 − 75. doi:  10.1016/j.envexpbot.2010.06.005
    [27] MARTIN L B, ROMERO P, FICH E A, et al. Cuticle biosynthesis in tomato leaves is developmentally regulated by abscisic acid [J]. Plant Physiology, 2017, 174(3): 1384 − 1398. doi:  10.1104/pp.17.00387
    [28] XIONG C, XIE Q, YANG Q, et al. WOOLLY, interacting with MYB transcription factor MYB31, regulates cuticular wax biosynthesis by modulating CER6 expression in tomato [J]. Plant Journal, 2020, 103(1): 323 − 337. doi:  10.1111/tpj.14733
    [29] JU Y L, LIU M, ZHAO H, et al. Effect of exogenous abscisic acid and methyl jasmonate on anthocyanin composition, fatty acids, and volatile compounds of cabernet sauvignon (Vitis vinifera L. ) grape berries [J]. Molecules, 2016, 21(10): 1354. doi:  10.3390/molecules21101354
    [30] KOSMA D K, MURMU J, RAZEQ F M, et al. AtMYB41 activates ectopic suberin synthesis and assembly in multiple plant species and cell types [J]. The Plant Journal, 2014, 80(2): 216 − 229. doi:  10.1111/tpj.12624
    [31] RAFFAELE S, VAILLEAU F, LEGER A, et al. A MYB transcription factor regulates very-long-chain fatty acid biosynthesis for activation of the hypersensitive cell death response in Arabidopsis [J]. The Plant Cell, 2008, 20(3): 752 − 767. doi:  10.1105/tpc.107.054858
    [32] SEO P, PARK C. Cuticular wax biosynthesis as a way of inducing drought resistance [J]. Plant signaling & behavior, 2011, 6(7): 1043 − 1045.
    [33] LEE S B, SUH M C. Cuticular wax biosynthesis is up-regulated by the MYB94 transcription factor in Arabidopsis [J]. Plant Cell Physiology, 2015, 56(1): 48 − 60. doi:  10.1093/pcp/pcu142
    [34] WANG Z, LI Z F, WANG S S, et al. 2021. NtMYB12a acts downstream of sucrose to inhibit fatty acid accumulation by targeting lipoxygenase and SFAR genes in tobacco [J]. Plant Cell Environment, 2021, 44(3): 775 − 791. doi:  10.1111/pce.13957
    [35] SUN Y, SONG K, LIU L, et al. Sulfoquinovosyl diacylglycerol synthase impairs glycolipid accumulation and photosynthesis in phosphate deprived rice [J]. Journal Experimental Botany, 2021, 72(18): 6510 − 6523. doi:  10.1093/jxb/erab300
    [36] BI F C, LIU Z, WU J X, et al. Loss of ceramide kinase in Arabidopsis impairs defenses and promotes ceramide accumulation and mitochondrial H2O2 bursts [J]. The Plant Cell, 2014, 26(8): 3449 − 3467. doi:  10.1105/tpc.114.127050
    [37] SAUCEDO G M, GARCIA A, GONZALEZ S A, et al. MPK6, sphinganine and the LCB2a gene from serine palmitoyltransferase are required in the signaling pathway that mediates cell death induced by long chain bases in Arabidopsis [J]. New Phytologist, 2011, 191(4): 943 − 957. doi:  10.1111/j.1469-8137.2011.03727.x
    [38] LIANG H, YAO N, SONG J T, et al. Ceramides modulate programmed cell death in plants [J]. Genes Development, 2003, 17(21): 2636 − 2641. doi:  10.1101/gad.1140503
  • [1] 董晨, 郑雪文, 王弋, 全振炫, 李伟才.  荔枝GRF基因家族的全基因组鉴定及表达分析 . 热带生物学报, 2024, 15(2): 190-197. doi: 10.15886/j.cnki.rdswxb.20230099
    [2] 向丽君, 赖军, 郭昊, 杨君, 王守创.  番茄糖基转移酶基因SlUGT75C1-like的功能研究及自然变异分析 . 热带生物学报, 2024, 15(2): 210-216. doi: 10.15886/j.cnki.rdswxb.20230062
    [3] 朱昆宇, 潘成才, 吉怡颖, 王宇, 周扬, 江行玉.  盐胁迫下冰菜盐囊泡细胞的比较代谢组学分析 . 热带生物学报, 2024, 15(2): 224-231. doi: 10.15886/j.cnki.rdswxb.20230020
    [4] 高雪, 劳广术, 谭峥, 方渝凯, 刘文波, 靳鹏飞, 缪卫国.  贝莱斯芽孢杆菌HN-2次生代谢产物处理下黄单胞菌(Xoo)的转录组分析 . 热带生物学报, 2023, 14(4): 389-398. doi: 10.15886/j.cnki.rdswxb.2023.04.006
    [5] 王宇航, 孟秀利, 黄山春, 林兆威, 宋薇薇, 唐庆华, 覃伟权.  双条拂粉蚧基因组DNA提取方法的比较和优化 . 热带生物学报, 2023, 14(6): 636-641. doi: 10.15886/j.cnki.rdswxb.20220085
    [6] 赖军, 张越冉, 周海鸿, 王守创, 沈双欠, 杨君, 罗杰.  基于代谢组学的重要热带植物的研究进展 . 热带生物学报, 2023, 14(3): 248-258. doi: 10.15886/j.cnki.rdswxb.2023.03.002
    [7] 王颖, 吴红芬, 张颖, 汪刘浩, 那威, 吴科榜.  文昌鸡中脂联素受体CDH13基因的表达 . 热带生物学报, 2022, 13(6): 558-561. doi: 10.15886/j.cnki.rdswxb.2022.06.003
    [8] 全飞, 兰国玉, 魏亚情, 李明美, 孙树晴, 杜昊楠.  基于宏基因组测序的橡胶林病原微生物的分析 . 热带生物学报, 2022, 13(1): 27-35. doi: 10.15886/j.cnki.rdswxb.2022.01.005
    [9] 李歆毓, 丁梦莹, 冯尔辉, 张源源, 章可兰, 万迎朗.  白边侧足海天牛全基因组Survey分析 . 热带生物学报, 2022, 13(5): 457-463. doi: 10.15886/j.cnki.rdswxb.2022.05.005
    [10] 武映宏, 黄东益, 王永, 张大鹏, 刘蕊, 孙熹微, 肖勇.  椰子ALDH基因家族的鉴定及生物信息学分析 . 热带生物学报, 2022, 13(6): 541-549. doi: 10.15886/j.cnki.rdswxb.2022.06.001
    [11] 邓雅洁, 刘云豪, 孙琪, 邱越, 宋家明, 赖杭桂.  在干旱胁迫下有性四倍体木薯叶片的蛋白质组学研究 . 热带生物学报, 2021, 12(2): 210-218. doi: 10.15886/j.cnki.rdswxb.2021.02.010
    [12] 董晨, 魏永赞, 王弋, 郑雪文, 李伟才.  转录组荔枝Dof基因家族的鉴定及其表达 . 热带生物学报, 2021, 12(1): 7-14. doi: 10.15886/j.cnki.rdswxb.2021.01.002
    [13] 徐航, 刘贤青, 袁弘伦, 罗杰.  槟榔碱合成前体物质的空间分布及槟榔碱合成通路解析 . 热带生物学报, 2021, 12(3): 271-278. doi: 10.15886/j.cnki.rdswxb.2021.03.001
    [14] 张忠辉, 邓渊, 强奇, 陈日东, 杨君, 刘贤青, 王守创, 袁潜华, 罗杰, 罗越华.  基于高效液相色谱-质谱联用技术解析山栏稻营养品质的代谢物谱 . 热带生物学报, 2021, 12(4): 419-427. doi: 10.15886/j.cnki.rdswxb.2021.04.003
    [15] 张萍萍, 王震南, 杨蕊, 孙冰, 贾爱群.  金花茶花中木脂素类化学成分的分析 . 热带生物学报, 2020, 11(3): 296-300, 323. doi: 10.15886/j.cnki.rdswxb.2020.03.006
    [16] 陈婉姨, 李新佳, 罗素兰, 长孙东亭.  信号芋螺基因组DNA提取方法的优化 . 热带生物学报, 2020, 11(2): 238-244, 250. doi: 10.15886/j.cnki.rdswxb.2020.00.015
    [17] 王子贺, 兰坚强, 廖承红, 韩谦.  天维菌素处理埃及伊蚊幼虫的转录组学分析 . 热带生物学报, 2020, 11(3): 281-287. doi: 10.15886/j.cnki.rdswxb.2020.03.004
    [18] 代佳妮, 于靖, 郑蔚, 王健, 吴友根, 胡新文, 赖杭桂.  海南油茶2个优良单株的比较转录组学分析 . 热带生物学报, 2020, 11(4): 426-436. doi: 10.15886/j.cnki.rdswxb.2020.04.005
    [19] 贺延恒, 车海彦, 曹学仁.  橡胶树丛枝病的组织结构观察及分子发生机理 . 热带生物学报, 2020, 11(2): 223-230. doi: 10.15886/j.cnki.rdswxb.2020.02.013
    [20] 张耀元, 刘玲, 周秀娟, 房传营, 王守创, 丁元昊, 罗杰.  转录组和代谢组揭示独脚金内酯缺失影响水稻根和叶中的脂质和黄酮代谢 . 热带生物学报, 2020, 11(4): 415-425. doi: 10.15886/j.cnki.rdswxb.2020.04.004
  • 加载中
图(5)
计量
  • 文章访问数:  807
  • HTML全文浏览量:  247
  • PDF下载量:  58
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-02-11
  • 录用日期:  2022-05-17
  • 修回日期:  2022-02-24
  • 网络出版日期:  2022-05-24
  • 刊出日期:  2022-09-21

通过整合转录组与代谢组解析不同类型椰子的脂肪酸调控网络

doi: 10.15886/j.cnki.rdswxb.2022.05.008
    基金项目:  国家自然科学基金(31960213)
    作者简介:

    邓渊(1996−),男,海南大学热带作物学院2019级硕士研究生. E-mail:yuandeng@hainanu.edu.cn

    通讯作者: 王守创(1989−),男,教授. 研究方向:植物次生代谢调控及代谢组学. E-mail:shouchuang.wang@hainanu.edu.cn罗杰(1971−),男,教授. 研究方向:植物次生代谢调控及代谢组学. E-mail:jie.luo@hainanu.edu.cn
  • 中图分类号: S 667.4;Q 943

摘要: 为了探究不同类型椰子中脂质种类和含量的差异,本研究利用非靶向代谢组学方法对海南高种和矮种椰子椰肉组织的代谢物进行了分析,共检测到12 579种代谢信号,其中鉴定了564种代谢物,包括氨基酸、脂肪酸、类黄酮等。定量分析表明高种椰子的脂质总体含量高于矮种椰子,其中绿矮椰子脂质含量最低。甘油脂、鞘脂和脂肪酰是不同类型椰子中主要的差异物质。转录组数据结果表明,不同类型椰子中的差异基因主要参与脂肪酸代谢、α-亚麻酸代谢、木质素代谢和苯丙烷代谢。基于表达模式的不同将差异基因划分为20个不同的模块,其中“Melightyellow”模块与差异积累的脂肪酸高度相关。通过分析模块中脂质途径结构基因启动子区的cis-element,发现大多数脂质途径基因的启动子区存在多个MYB家族的结合位点及大量的激素信号响应元件。结果表明,MYB家族的转录因子可能在脂质的合成调控中发挥重要作用,并且激素可能作为一种信号分子参与到脂质的合成调控中。本研究通过组合策略对不同品种椰肉组织脂质差异的分子机理进行分析,并进一步构建了转录因子−结构基因−代谢物的合成调控网络,为后续解析脂质合成调控的分子机制提供了资源,为油脂椰子品种选育提供了依据。

English Abstract

邓渊,赖军,毛梦迪,等. 通过整合转录组与代谢组解析不同类型椰子的脂肪酸调控网络[J]. 热带生物学报,2022, 13(5):478−487. DOI:10.15886/j.cnki.rdswxb.2022.05.008 doi:  10.15886/j.cnki.rdswxb.2022.05.008
引用本文: 邓渊,赖军,毛梦迪,等. 通过整合转录组与代谢组解析不同类型椰子的脂肪酸调控网络[J]. 热带生物学报,2022, 13(5):478−487. DOI:10.15886/j.cnki.rdswxb.2022.05.008 doi:  10.15886/j.cnki.rdswxb.2022.05.008
DENG Yuan, LAI Jun, MAO Mengdi, ZHANG Yueran, LI Chun, YANG Jun, WANG Shouchuang, LUO Jie. Elucidation of the fatty acid regulatory network for different coconut varieties by integrating transcriptome and metabolome[J]. Journal of Tropical Biology, 2022, 13(5): 478-487. doi: 10.15886/j.cnki.rdswxb.2022.05.008
Citation: DENG Yuan, LAI Jun, MAO Mengdi, ZHANG Yueran, LI Chun, YANG Jun, WANG Shouchuang, LUO Jie. Elucidation of the fatty acid regulatory network for different coconut varieties by integrating transcriptome and metabolome[J]. Journal of Tropical Biology, 2022, 13(5): 478-487. doi: 10.15886/j.cnki.rdswxb.2022.05.008
  • 椰子(Cocos nucifera L.)是许多热带国家最重要的经济作物[1]。椰子属于棕榈科椰子属,基因组大小约为2.4 Gb[2-3]。根据植物形态和繁殖类型,椰子可以分为2种类型:typica (tall)和nana (dwarf)[4]。椰子根据进化起源可以分为2个亚群——印度洋种群和太平洋种群[5-6]。高种椰子异花授粉,基因型表现为一定程度的杂合,而矮种椰子自花授粉,基因型纯合。矮种椰子生长缓慢、植株矮小、树冠较小、果小而多、椰肉厚度薄、开花时间较早,通常在种植4~6 a后开花,株高10~15 m[7-9]。高种椰子耐寒,果实大,能够适应更广泛的环境,株高20~30 m,初次开花时间晚(种植8~10 a后开花)。在中国,椰子主要种植在海南和云南等地区,海南地区的高种椰子“Hainan Tall”有着36 000 hm2的种植面积[10],高种椰子耐盐碱和高温,但对低温敏感。目前普遍认为矮种椰子经高种椰子驯化而来[4],且矮种椰子更易受不利环境的影响。椰油占椰肉重量的65%[8,11],与其他植物油相比,椰油具有独特的理化性质和药用价值。椰子油富含饱和脂肪酸(93%),主要成分是己酸(C6)、辛酸(C8)、癸酸(C10)、月桂酸(C12)、肉豆蔻酸(C14)、棕榈酸(C16)、硬脂酸(C18)和花生酸(C20),其中85%是中链脂肪酸(medium-chain fatty acids, MCFA)。据报道,月桂酸在棕榈科中的主要来源是椰子和油棕,两者都仅限于在热带和亚热带气候中生长[11-14]。月桂酸能够提高氧化稳定性,具有低熔点,能产生稳定的乳液,这对食品和化学工业应用具有重要意义[11]。目前月桂酸在医药和食品工业等各个领域都具有广泛的应用,包括用于制造肥皂、洗涤剂、纺织品、油漆、清漆、化妆品、药品等。肉豆蔻酸在食品工业中常用作调味剂和食品添加剂[15]。肉豆蔻酸有一定的药理活性,相较于其他饱和脂肪酸,肉豆蔻酸有更强的抗氧化活性[16],肉豆蔻酸还可以降低血糖[17]。除了饱和脂肪酸以外,椰子油还含有少量不饱和脂肪酸,如油酸(C18:1)、棕榈油酸(C18:1)和亚麻酸(C18:3)。目前有报道使用椰子油可以降低阿尔茨海默病的患病风险[13-14],抑制肿瘤生长、抑制氧化应激和神经炎症,以及消除抗肿瘤药甲氨蝶呤引起的脑神经毒性[18-24]。植物脂肪酸代谢途径在其他物种中已有广泛报道[25-34],但目前在椰子中罕有报道。不同类型椰子品种在脂肪酸积累和合成调控途径是否存在差异,仍有待阐明。椰肉作为脂肪酸积累的主要部位,是研究脂肪酸合成调控的理想场所。本研究利用非靶向代谢组学,比较了海南高种、黄矮椰子和绿矮椰子的脂肪酸积累水平,研究表明,高种椰子的脂肪酸含量最高,黄矮椰子其次,绿矮椰子最低。结合转录组学数据,重建了椰子脂肪酸代谢途径,研究发现,相较于矮种椰子而言,脂肪酸途径的基因大多在高种椰子中高表达。WGCNA结果表明,“Melightyellow”模块中的基因与差异积累的脂肪酸有着相似的变化趋势,在此基础上结合模块中结构基因启动子区的cis-element,预测并构建了完整的脂肪酸调控网络。本研究通过组学策略探究了不同类型椰子品种脂肪酸差异的分子机理,预测并构建了完整的椰子脂肪酸调控网络,为后续椰子脂肪酸合成调控的进一步解析提供了可用资源,为高油脂椰子品种的选育提供了理论依据。

    • 本试验材料由中国热带农业科学院椰子研究所提供,包括2份高种椰子(Hainan Tall, HT),2份黄矮椰子(Yellow Dwarf,YD)和2份绿矮椰子(Green Dwarf, GD)。每份材料随机选取3株长势一致的植株,每株植株采摘3个大小和发育时期一致的椰果,取椰肉混合后作为1个生物学重复,并将其一分为二,一份作为转录样,一份作为代谢样。将样品立即置于液氮冷冻,−80 ℃超低温冰箱保存备用。

    • 椰肉代谢样品经冷冻干燥后,采用BM500球磨仪研磨成粉末(频率30 Hz,持续时间 1.5 min)。称取100 mg研磨好的干粉样品,加入1.2 mL含有0.1 mg·L−1利多卡因内标的φ = 70%甲醇提取液,涡旋30 s,低温静置10 min,重复上述操作3次后,放入4℃冰箱继续萃取代谢产物12 h,10 300 r·min−1 4 ℃离心10 min,取上清液,用0.22 μm的滤膜过滤后,放入进样瓶中即可用于LC-MS进行代谢物检测分析。

    • 对于原始的下机数据使用Fastp软件的默认参数去除低质量序列和接头序列,用HISAT2软件默认参数将其比对到矮种椰子的参考基因组,Samtools排序后使用featureCounts进行定量,用TPM (transcripts per million)对原始定量结果进行标准化。R软件中的FactoMineR 和factoextra用于主成分分析。DEseq2用于样本间差异表达基因的鉴定,同时满足|log2FoldChange| > 1和FDR (False Discovery Rate) < 0.05的基因被认为是2个样本比较间的差异表达基因。

    • 为了探究不同类型椰肉的代谢组差异,选择了6份具有代表性的椰子种质资源,包括2份高种椰子(HT)、2份黄矮椰子(YD)、2份绿矮椰子(GD),进行代谢物检测。通过 Q Exactive-Orbitrap UHPLC-ESI-MS/MS 在正离子模式下进行非靶向扫描,共检测到12 579种质谱信号。为了研究椰子果肉的代谢组构成,利用多种策略对所检测的质谱信号进行代谢物结构解析。对于能获取商业标准品的代谢物,通过比对保留时间(retention time, RT),母离子(Q1)及质谱碎片信息等因素来鉴定代谢物。如RT为9.15 min的代谢信号(Cnm04191),其母离子精确的m/z为279.230 28,并且有149.0、219.1、57.1和105.0等一系列的特征碎片。与标准品比对,代谢信号Cnm04191由于相同的质谱碎片信息被鉴定为α-Linolenate (图1-A,1-B)。对于不能获得商业标准品的代谢物,将试验获取的质谱信息与已经公开发表的文献或者代谢数据库进行比对来初步解析代谢物结构。RT为8.39 min的质谱信号(Cnm09776)的MS1谱图中检测到准确的m/z为518.322 20,二级质谱碎片中观察到[Y0]+离子184.1,这是由于形成了磷酸胆碱这一特征碎片,因此,推测质谱信号Cnm09776为LysoPC (18:3/0:0) (图1-C)。质谱信号Cnm06436被注释为Glyceryl monolinoleate (图1-D)。通过Compound Discoverer 3.1软件自动批量匹配在线代谢数据库,包括mzCloud、ChemSpider、MassBank等进行代谢物结构鉴定。最终,超过 564 种代谢物被鉴定,主要为氨基酸、黄酮类化合物、有机酸、脂质 、苯丙烷 、糖类 、多酚、生物碱、核苷酸及其衍生物、植物激素。

      图  1  利用UHPLC–HRMS对次生代谢物进行检测和鉴定

      基于检测到的97种脂质信号进行了主成分分析(principal component analysis, PCA),结果表明,每个品种的3个生物学重复都能较好地聚在一起,黄矮椰子和绿矮椰子聚集较近,但两者与高种椰子相聚较远(图2-A)。结果表明,同一样本不同重复间的代谢数据有较高的一致性,且不同类型的椰子品种在脂质积累水平有较大的差异,高矮种间的差异更大。进一步分析结果表明,高种椰子的脂质水平极显著(P < 2.22E-16)高于矮种椰子(图2-C),且绿矮椰子的脂质积累水平最低(图2-D)。热图结果表明,检测到的97种脂质在3种不同类型的椰子果肉中具有不同的积累模式(图2-B)。为进一步探究三者间脂质积累存在差异的代谢产物,利用OPLS-DA对样本两两间差异代谢物进行鉴定(YD vs HT, GD vs HT, GD vs YD)。在两两比较间脂质积累存在差异的物质共有50种,其中5种脂质在3组比较间都表现为积累差异(图2-E)。进一步对积累存在差异的代谢物进行注释,脂肪酰、甘油脂、鞘脂是主要的差异物质,其中脂肪酰占比最高达74%,鞘脂和甘油酯分别各占8% (图2-F)。糖脂主要与光合作用及磷响应相关[35],鞘脂主要与植物生长、植物育性和胁迫应答等相关[36-38]。以上结果表明,这些物质在不同类型的椰子中迥异的积累模式,可能与它们响应生物和非生物胁迫以及适应不同地区的环境条件有关。

      图  2  不同类型椰子品种脂质代谢组的比较

    • 为进一步研究不同类型椰子品种在脂质通路的转录水平是否存在差异。对代表性品种的椰肉组织进行转录组测序。为了解不同类型椰子的转录动态,基于转录组数据进行PCA (图3-A)。结果表明,绝大多数的生物学重复彼此聚集,表明了转录数据的高可靠性。此外,这些转录本根据样本来源聚为3类,其中PC1捕获了43.5%的数据方差,代表不同类型的椰子样本之间的差异,这些结果均表明三者在基因的转录水平上存在较大的差异。随后对不同类型椰子品种的转录数据进行两两比较来鉴定不同样本间的差异表达基因(differently expressed genes, DEGs)。笔者在HT vs YD,HT vs GD,GD vs YD间分别鉴定到了2 139、4 729和2 599个DEGs,在两两比较间存在差异的基因共计5 841个,在3个比较中都存在差异表达的基因有384个(图3-B)。为了进一步探究差异基因的分子功能及其参与的生物学途径,对5 841个差异基因进行Kyoto Encyclopedia of Genes and Genomes (KEGG)和Gene Ontology (GO)富集分析,KEGG结果表明这些差异基因主要与植物激素信号转导、植物与病原菌互作、氨基酸代谢相关(图3-C)。GO富集分析结果也表明,不同椰子品种的表达差异主要体现在赤霉素信号通路、脂质储存和油菜素内酯响应等途径(图3-D)。

      图  3  不同品种椰子差异表达基因的鉴定与功能富集

    • 采用加权共表达网络分析(weighted gene co-expression network analysis, WGCNA),将高连接度的基因定义为基因簇,同一簇内的基因具有较高的相关性。这些表达模式相近的基因在功能上具有相关性。根据表达模式的不同可以将5 841个差异基因划分为20个模块,其中模块“Melightyellow”中基因的表达模式与差异脂肪酸的积累模式高度相关(图4-A)。随后对“Melightyellow”中的基因进行蛋白功能注释和功能富集,其中参与“Fatty acid metabolism”,“Glycerolipid metabolism”,“Biosynthesis of unsaturated fatty acids”等脂肪酸相关通路的基因作为节点构建调控网络。基于保守功能结构域、蛋白功能注释和同源比对,笔者初步预测了13个可能参与脂肪酸代谢的基因,其中包括long chain acyl-CoA synthetaseglycerol-3-phosphate dehydrogenase1-acyl-sn-glycerol-3-phosphate acyltransferaselysophospholipid acyltransferase等。随后基于基因−基因、基因−代谢物、代谢物−代谢物间的皮尔森相关系数,结合脂肪酸代谢途径基因启动子区cis-elements的潜在结合能力,在“Melightyellow”模块中鉴定到了10个转录因子,包括ERF,MYB 和bHLH等。这些转录因子与“Melightyellow”模块中鉴定的13个候选基因,差异积累的脂肪酸都表现出较高的相关性(图4-B)。其中ADH (alcohol dehydrogenase, AZ13G0249870)基因参与20-Hydroxyarachidonic acid到20-Oxoarachidonic acid (Cnm05349)的合成,而LPC 18:1在LPCAT (lysophosphatidylcholine acyltransferase, AZ11G0210280)基因的参与下与Oleic acid结合生成Phosphatidylcholine (18:1/18:1, Cnm11965) (图4-C)。研究发现20-Oxoarachidonic acid和Phosphatidylcholine (18:1/18:1)在高种椰子中的含量显著高于矮种椰子(图4-D),转录组的数据也表明,ADHLPCAT在高种椰子中有着更高的转录水平(图4-E)。为了进一步验证转录数据的可靠性,笔者通过qRT-PCR证实,ADHLPCAT在高种椰子中的表达水平显著高于黄矮椰子和绿矮椰子(图4-F)。这些基因在高种椰子中的高水平表达可能是导致高种椰子脂质含量显著高于黄矮椰子和绿矮椰子的重要原因。

      图  4  加权共表达分析鉴定与脂肪酸合成相关的基因

    • 脂质的代谢和转录数据表明高种椰子在脂质积累水平和脂质合成基因的转录水平都显著高于矮种椰子。为了进一步探究脂肪酸合成通路在高矮种椰子间是否存在差异,构建了亚油酸相关的合成代谢通路。该通路由丙酮酸起始,在一系列酶的催化下生成亚油酸,亚油酸在酶的作用下最终生成茉莉酸甲酯(图5)。基于蛋白的同源性关系,对椰子脂肪酸合成通路的基因进行注释,包括 ACCase、FAD3、LOX2S、AOS、AOC、OPR、OPCL1、JMT等基因。结合转录组数据的定量结果和代谢组的相对物质含量,笔者发现在脂肪酸上游途径,参与亚油酸、亚麻酸合成的的基因,如ACCase在高种椰子中的表达量高于矮种,而代谢物含量变化也与转录数据相一致。亚麻酸至9,10-Epoxy-10,12Z-octadecadienoate途径和亚油酸至8,9-DiHETrE途径及亚油酸至10-Oxo-11,15-phytodienoic acid和Heptadecatrienal途径,相关的代谢物在高种椰子中也始终表现为高积累。在代谢通路下游,茉莉酸甲酯合成途径中的关键基因如AOC、OPR、JMT等均在高种椰子中有着更高的表达水平,相关代谢物也在高种椰子中表现为高积累。以上结果均表明,高种椰子椰肉中的脂肪酸含量高于矮种椰子,这一现象产生的原因应该是脂肪酸合成通路的基因大多在高种椰子中特异性高表达。

      图  5  高矮种椰子脂肪酸合成途径基因表达水平及代谢产物积累水平的比较

    • 本研究不仅成功构建了不同类型椰子椰肉组织的代谢数据库,而且发现高种椰子的脂质积累水平远高于矮种椰子,不同类型的椰子品种在脂质的积累上具有特异性。此外通过整合转录组学和代谢组学的相关数据,发现“Melightyellow”模块与脂肪酸高度相关,并从中鉴定到了大量参与脂肪酸合成的结构基因。笔者通过分析脂肪酸合成途径基因启动子区的cis-element,发现这些基因的启动子区存在大量的激素信号响应元件和MYB家族转录因子的结合位点,这表明MYB家族在脂肪酸合成调控中发挥重要作用,而且激素作为一种信号分子可能通过响应外界环境条件的变化来参与到脂肪酸的合成途径中。蜡质主要由脂肪酸和相关衍生物构成,目前大量的研究结果表明,蜡质的合成受到多种环境因素的诱导,如干旱胁迫等[26],ABA等激素可以正向调控蜡质合成基因的表达[26-27]。熊程等[28]的研究结果表明,MYB31通过调控CER6的表达,协同调控番茄表皮蜡质的合成,并且CER6的表达进一步受到ABA的诱导。鞠延仑等[29]在葡萄中发现外源施加茉莉酸甲酯(MeJA)可以调控脂肪酸代谢途径的脂氧合酶,进而对葡萄的风味造成影响。此外,研究者鉴定到大量的参与脂质途径的MYB转录因子,如MYB30MYB41能够调控超长链脂肪酸的合成来影响棉花纤维的伸长[30-31]MYB96MYB94能够激活蜡质合成基因的表达[32-33],烟草中的MYB12a能够促进脂肪酸的降解等[34],这些均与本研究的结果相契合。本研究通过结合共表达分析和分析脂质途径结构基因启动子区的cis-element,对模块中可能参与脂质途径的MYB家族转录因子进行了进一步鉴定,预测并构建了椰子脂肪酸途径的完整调控网络,为后续脂肪酸途径的进一步解析提供了有力资源。本研究不仅阐明了高矮种椰子间脂肪酸差异及相关分子机制,而且可为椰子油脂的遗传改良和高油脂椰子品种的选育提供一定的依据。

参考文献 (38)

目录

    /

    返回文章
    返回