中国农业科学 ›› 2024, Vol. 57 ›› Issue (9): 1779-1793.doi: 10.3864/j.issn.0578-1752.2024.09.012
戚仁洁1,2(), 宁宇3(
), 刘静3, 刘之洋3, 徐海3, 罗志丹1,2(
), 陈龙正3(
)
收稿日期:
2023-09-29
接受日期:
2023-12-18
出版日期:
2024-05-01
发布日期:
2024-05-09
通信作者:
联系方式:
戚仁洁,E-mail:13956705571@163.com。宁宇,E-mail:20190022@jaas.ac.cn。戚仁洁和宁宇为同等贡献作者。
基金资助:
QI RenJie1,2(), NING Yu3(
), LIU Jing3, LIU ZhiYang3, XU Hai3, LUO ZhiDan1,2(
), CHEN LongZheng3(
)
Received:
2023-09-29
Accepted:
2023-12-18
Published:
2024-05-01
Online:
2024-05-09
摘要:
【目的】皂苷是苦瓜果实苦味的主要来源,具有降血糖、抗癌等多种药用价值。鉴定和挖掘参与调控苦瓜皂苷生物合成的代谢通路及基因,为进一步深入解析苦瓜皂苷形成的分子机制提供参考。【方法】以苦瓜高代自交系GK24为试验材料,分别在子房期(T1)、幼果期(T2)、商品果期(T3)和完全成熟期(T4)取果实组织,通过香草醛-冰醋酸法测定不同时期的苦瓜皂苷含量,利用转录组测序的方法鉴定差异表达基因。【结果】从T1到T4,苦瓜内源皂苷含量随着果实发育而呈现显著下降的趋势。转录组测序共鉴定到17 504个基因,在T1-vs-T2、T2-vs-T3和T3-vs-T4三组比较中,下调表达基因数量均高于上调表达基因数量。GO富集分析表明含磷复合物代谢过程、驱动蛋白复合体和2-琥珀酰-6-羟基-环己二烯-1-羧酸合成酶活性分别是生物过程、细胞组分及分子功能模块最显著富集的条目。KEGG分析显示全局与概述图谱、碳水化合物代谢和氨基酸代谢是3组比较中共同的最显著富集的代谢通路。根据表达模式可将基因划分为20个模块,其中profile 0中基因的表达量从T1到T4逐渐降低,与苦瓜果实发育过程中皂苷含量变化规律一致。从profile 0中鉴定出与苦瓜皂苷生物合成相关的基因14个,包括萜类骨架生物合成通路上的基因AAT1、HMG1、MVK、PMK、MVD2和FPS1,倍半萜与三萜生物合成途径中的基因SS12、SQE1和CPQ及后修饰酶基因CYP97A3、CYP71AN24、UGT94E5及UGT73C6。实时荧光定量PCR(qRT-PCR)结果与RNA-seq结果一致。【结论】苦瓜果实中皂苷含量随着果实不断成熟而逐渐降低。苦瓜皂苷生物合成首先通过上游的萜类骨架生物合成(ko00900)代谢途径完成萜类皂苷骨架的积累,随后通过倍半萜与三萜生物合成(ko00909)代谢途径以及下游的氧化还原酶、糖基转移酶的修饰最终生成皂苷,本研究通过转录组测序共鉴定到14个与苦瓜皂苷生物合成相关的基因。
戚仁洁, 宁宇, 刘静, 刘之洋, 徐海, 罗志丹, 陈龙正. 基于转录组测序的苦瓜皂苷合成相关基因的鉴定和分析[J]. 中国农业科学, 2024, 57(9): 1779-1793.
QI RenJie, NING Yu, LIU Jing, LIU ZhiYang, XU Hai, LUO ZhiDan, CHEN LongZheng. Identification and Analysis of Genes Related to Bitter Gourd Saponin Synthesis Based on Transcriptome Sequencing[J]. Scientia Agricultura Sinica, 2024, 57(9): 1779-1793.
表1
qRT-PCR引物序列"
基因名 Gene name | 正向引物Forward primer (5′-3′) | 反向引物Reverse primer (5′-3′) |
---|---|---|
CPQ | GTTGAAAGAGGGGGAAGAGG | CACCAGACCCGGAAGTAAAA |
FSP1 | GCGGATCTCAGGTCAACATT | TCCAGGCACATTGTAGTCCA |
CYP97A3 | GGGAGATATGGCGTGTAAGAAGAC | CTAAAGTCAGACGAGAGAACAGGG |
CYP71AN24 | GAACACAAGGCCAAGATGAGAG | CAGTTGTGTCAGTTCCAGCTAC |
UGT73C6 | TGCCTCCAGAGCTGAAACAT | TCCCATGAATCCAGCCACTT |
UGT94E5 | TCATCACACAACCAAAGGCC | ATCGGAATGTTGAGCGATGC |
SQE1 | AGCAAGGGACGGTTACTTCA | TGAGAAGCAGCCATCACAGA |
MVK | GCAGTCGTTCATGGATCCAC | TGGCCACGAAAACTCCAATG |
PMK | GGGACCATGAGAGGACTGAG | AGGGTCAGACTTTTGCCACT |
Actin | GGCAAACCCTAAAGTTTTCTTCG | GATGAGCCCTTGTAATGAAGTGG |
表2
转录组测序数据统计"
样品 Sample | 原始数据 Raw reads | 过滤后数据 Clean reads | Q30含量 Q30 content (%) | GC含量 GC content (%) | 单一比对率 Unique mapped (%) | 总比对率 Total mapped (%) |
---|---|---|---|---|---|---|
T1-1 | 46878594 | 46673674 | 92.28 | 46.80 | 93.74 | 96.88 |
T1-2 | 55531244 | 55302940 | 92.51 | 46.59 | 93.74 | 96.89 |
T1-3 | 65068774 | 64756962 | 91.96 | 46.30 | 93.49 | 96.61 |
T2-1 | 53923956 | 53678084 | 91.97 | 46.70 | 93.42 | 96.63 |
T2-2 | 46863574 | 46589694 | 92.55 | 47.09 | 93.61 | 96.87 |
T2-3 | 51548340 | 51304488 | 92.54 | 47.04 | 93.43 | 96.74 |
T3-1 | 44171240 | 43942622 | 92.37 | 46.84 | 93.10 | 95.88 |
T3-2 | 44849744 | 44668862 | 92.68 | 47.20 | 94.17 | 97.07 |
T3-3 | 55212924 | 54973380 | 92.20 | 47.14 | 93.68 | 96.43 |
T4-1 | 69720780 | 69450220 | 92.42 | 47.30 | 93.73 | 96.46 |
T4-2 | 66405070 | 66157562 | 92.94 | 46.67 | 94.16% | 96.82 |
T4-3 | 64248114 | 63979034 | 92.29 | 46.69 | 92.75 | 96.29 |
表3
Profile 0中显著富集的代谢通路"
通路名称 Pathway name | 通路ID Pathway ID | 基因数目 Number of genes | 所有基因数 Number of all genes | P值 P value |
---|---|---|---|---|
核糖体 Ribosome | ko03010 | 83 | 305 | 3.36E-13 |
同源重组 Homologous recombination | ko03440 | 20 | 58 | 1.03E-05 |
DNA复制 DNA replication | ko03030 | 16 | 46 | 7.03E-05 |
氨基酸的生物合成 Biosynthesis of amino acids | ko01230 | 44 | 202 | 9.41E-05 |
错配修复 Mismatch repair | ko03430 | 14 | 39 | 0.000135124 |
缬氨酸、亮氨酸和异亮氨酸生物合成 Valine, leucine and isoleucine biosynthesis | ko00290 | 8 | 18 | 0.000748137 |
真核生物核糖体的生物合成 Ribosome biogenesis in eukaryotes | ko03008 | 19 | 77 | 0.002133956 |
嘧啶代谢 Pyrimidine metabolism | ko00240 | 15 | 58 | 0.00377061 |
2-氧代羧酸代谢 2-Oxocarboxylic acid metabolism | ko01210 | 14 | 53 | 0.004117066 |
糖鞘脂生物合成 Glycosphingolipid biosynthesis | ko00604 | 3 | 4 | 0.006896205 |
精氨酸生物合成 Arginine biosynthesis | ko00220 | 10 | 35 | 0.008097428 |
碱基切除修复 Base excision repair | ko03410 | 13 | 52 | 0.009192653 |
核质运输 Nucleocytoplasmic transport | ko03013 | 20 | 98 | 0.01563779 |
赖氨酸生物合成 Lysine biosynthesis | ko00300 | 5 | 13 | 0.01582533 |
萜类骨架生物合成 Terpenoid backbone biosynthesis | ko00900 | 14 | 62 | 0.0173777 |
糖胺聚糖降解 Glycosaminoglycan degradation | ko00531 | 6 | 18 | 0.01771733 |
氨酰基tRNA生物合成 Aminoacyl-tRNA biosynthesis | ko00970 | 13 | 57 | 0.0198342 |
不同类型的N-聚糖生物合成 Various types of N-glycan biosynthesis | ko00513 | 8 | 30 | 0.02599519 |
组氨酸代谢 Histidine metabolism | ko00340 | 5 | 16 | 0.03922549 |
表4
Profile 0中ko00900和ko00909通路中的基因"
基因编号 Gene number | 基因注释 Gene function | 通路 Pathway | T1 FPKM | T2 FPKM | T3 FPKM | T4 FPKM |
---|---|---|---|---|---|---|
evm.TU.chr1.608 | 法尼基半胱氨酸裂解酶异构体 Farnesylcysteine lyase isoform | ko00900 | 38.15 | 28.42 | 17.24 | 10.65 |
evm.TU.chr10.805 | 2-甲基-D-赤藓糖醇2,4-环二磷酸合酶 2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase | ko00900 | 115.94 | 71.18 | 43.36 | 26.32 |
evm.TU.chr2.34 | 香叶基香叶基二磷酸还原酶样异构体 Geranylgeranyl diphosphate reductase-like isoform | ko00900 | 30.58 | 17.55 | 16.44 | 7.16 |
evm.TU.chr3.1930 | 羟基-3-甲基丁-2-烯-1-基二磷酸合酶 4-hydroxy-3-methylbut-2-en-1-yl diphosphate synthase | ko00900 | 433.17 | 194.37 | 106.87 | 85.58 |
evm.TU.chr3.2317 | 法尼基焦磷酸合酶1 Farnesyl pyrophosphate synthase 1 | ko00900 | 358.50 | 164.00 | 51.39 | 7.70 |
evm.TU.chr4.352 | 乙酰辅酶A乙酰转移酶 Acetyl-CoA acetyltransferae | ko00900 | 439.08 | 95.24 | 57.25 | 34.89 |
evm.TU.chr4.613 | 4-羟基-3-甲基丁-2-烯基二磷酸还原酶 4-hydroxy-3-methylbut-2-enyl diphosphate reductase | ko00900 | 1091.57 | 439.96 | 370.39 | 217.45 |
evm.TU.chr6.140 | 1-脱氧-D-酮糖5-磷酸还原异构酶 1-deoxy-D-xylulose 5-phosphate reductoisomerase | ko00900 | 239.29 | 131.64 | 82.92 | 76.05 |
evm.TU.chr6.518 | 脱氧-D-酮糖-5-磷酸合成酶Deoxy-D-xylulose-5-phosphate synthase | ko00900 | 637.07 | 42.68 | 19.10 | 0.17 |
evm.TU.chr4.378 | HMG-CoA合成酶 3-hydroxy-3-methyl glutaryl coenzyme A synthesis | ko00900 | 1274.61 | 1010.52 | 91.31 | 2.87 |
evm.TU.chr8.4201 | 甲羟戊酸激酶 Mevalonate kinase | ko00900 | 105.59 | 52.77 | 31.46 | 14.07 |
evm.TU.chr9.3672 | 焦磷酸甲羟戊酸脱羧酶Pyrophosphomethylvalerate decarboxylase | ko00900 | 671.59 | 131.25 | 53.80 | 28.02 |
evm.TU.chr10.2918 | 磷酸甲羟戊酸激酶 Phosphomevalonate kinase | ko00900 | 85.10 | 58.17 | 9.86 | 3.75 |
evm.TU.chr6.363 | 磷酸甲羟戊酸激酶 Phosphomevalonate kinase | ko00900 | 41.95 | 20.49 | 13.08 | 12.38 |
evm.TU.chr2.1021 | 角鲨烯合成酶 Squalene synthase | ko00909 | 437.66 | 285.81 | 89.01 | 48.19 |
evm.TU.chr3.2797 | 橙花内酯合成酶 Nerolidol synthase | ko00909 | 15.85 | 9.29 | 0.69 | 0.07 |
evm.TU.chr7.10 | 葫芦二烯醇合酶 Cucurbitadienol synthase | ko00909 | 981.27 | 444.87 | 79.12 | 4.09 |
evm.TU.chr8.588 | 角鲨烯环氧酶 Squalene epoxidase 1 | ko00909 | 1030.32 | 506.37 | 76.82 | 1.14 |
[1] |
丁雷, 胡玉立, 李梅, 吴丽丽, 秦灵灵, 刘桐华. 苦瓜的药食价值研究概括. 中国医药导报, 2021, 18(13): 117-120.
|
|
|
[2] |
崔宏伟, 韩汶延, 于蕾, 苏秀兰. 苦瓜化学成分及药理作用研究进展. 世界科学技术-中医药现代化, 2021, 23(5): 1712-1719.
|
|
|
[3] |
doi: 10.4155/fmc-2017-0255 pmid: 29432030 |
[4] |
李莎莎. 苦瓜皂苷降血糖及促进糖尿病创面愈合作用的研究[D]. 广州: 南方医科大学, 2022.
|
|
|
[5] |
刘馨泽, 吴新民, 蔺冬梅, 冯琳, 杨洪梅, 王欢, 王冬雪, 刘淑莹, 吴巍, 万茜淋. 苦瓜化学成分的提取分离及药理活性研究进展. 中成药, 2022, 44(1): 177-182.
|
|
|
[6] |
郑方园. 苦瓜皂苷的分离纯化、结构鉴定及苦瓜颗粒剂工艺研究[D]. 合肥: 安徽农业大学, 2021.
|
|
|
[7] |
|
[8] |
doi: 10.1248/cpb.30.4334 |
[9] |
doi: 10.1248/cpb.32.2044 |
[10] |
安昌, 秦源, 方韵颖, 屈祺, 王炳锐, 郑平. 药用植物三萜类生物合成途径及关键酶基因研究进展. 中国中药杂志, 2022, 47(24): 6560-6572.
|
|
|
[11] |
徐圆圆, 陈仲, 贾黎明, 翁学煌. 植物三萜皂苷生物合成途径及调控机制研究进展. 中国科学: 生命科学, 2021, 51(5): 525-555.
|
|
|
[12] |
魏一丁, 熊超, 张天缘, 高翰, 尹青岗, 姚辉, 孙伟, 胡志刚, 陈士林. 基于转录组数据对七叶树中三萜皂苷合成途径的研究. 中国中药杂志, 2019, 44(6): 1135-1144.
|
|
|
[13] |
doi: 10.1007/s10265-017-1004-7 |
[14] |
doi: 10.1186/1471-2164-15-746 |
[15] |
doi: 10.1038/srep45295 |
[16] |
单春苗, 王晨凯, 施圆圆, 张声祥, 赵历强, 吴家文. 多花黄精甾体皂苷生物合成途径分析及关键酶基因研究. 中国中药杂志, 2020, 45(12): 2847-2857.
|
|
|
[17] |
石博, 万新建, 关峰, 张景云, 黄长林, 张会国, 黄国东. 苦瓜皂苷最佳提取工艺优化及其对α-葡萄糖苷酶的抑制活性. 食品研究与开发, 2023, 44(11): 145-150.
|
|
|
[18] |
徐斌, 董英. 分光光度法测定苦瓜总皂甙含量. 食品科学, 2005, 26(10): 165-169.
|
|
|
[19] |
|
[20] |
薛禹宸, 雷秀娟, 张浩, 王英平. 不同生长时期西洋参花蕾中皂苷成分及最佳采收期的确定. 吉林农业大学学报, 2021, 43(6): 697-704.
|
|
|
[21] |
苏璐璐, 朱文莹, 陈旭, 王辉, 曹依林, 王富. 基于转录组学和代谢组学解析番茄果实品质形成过程. 山东农业科学, 2023, 55(2): 9-19.
|
|
|
[22] |
王红娟, 唐荣莉, 蒋晓英, 官玲, 雷开荣, 黄启中, 林清. 辣椒不同发育时期果实的转录组分析. 西南农业学报, 2022, 35(3): 517-525.
|
|
|
[23] |
doi: 10.1016/j.phytochem.2011.01.015 pmid: 21333312 |
[24] |
doi: 10.1126/science.1259215 |
[25] |
doi: 10.1016/j.csbj.2023.01.012 pmid: 36698975 |
[26] |
朱源, 李检秀, 李坚斌, 黄日波. 药用植物中四环三萜皂苷合成生物学研究进展. 广西科学院学报, 2020, 36(3): 309-316.
|
|
|
[27] |
doi: 10.1093/pcp/pcw038 |
[28] |
|
[29] |
季宏与, 丁常宏, 于丹, 都晓伟. 转录组学技术在人参皂苷生物合成途径研究中的应用进展. 食品与发酵工程, 2022, 48(16): 318-325.
|
|
|
[30] |
穆事成. 糖基转移酶的筛选及罗汉果苷糖基化修饰[D]. 天津: 天津科技大学, 2021.
|
|
|
[31] |
doi: 10.1016/j.tig.2010.07.001 pmid: 20739089 |
[32] |
doi: 10.1098/rstb.2012.0426 |
[33] |
doi: 10.1146/annurev-arplant-050312-120229 pmid: 24498976 |
[34] |
doi: 10.1007/s11240-009-9535-9 |
[1] | 林蔚, 吴水金, 李跃森. 转录组和蛋白质组关联分析解析巴西蕉幼苗响应低温的分子机制[J]. 中国农业科学, 2024, 57(8): 1575-1591. |
[2] | 高晨曦, 郝陆洋, 胡悦, 李永祥, 张登峰, 李春辉, 宋燕春, 石云素, 王天宇, 黎裕, 刘旭洋. 干旱条件下玉米转座子插入关联的表观调控分析[J]. 中国农业科学, 2024, 57(6): 1034-1048. |
[3] | 韩旭东, 杨传奇, 张擎, 李亚伟, 杨夏夏, 何佳甜, 薛吉全, 张兴华, 徐淑兔, 刘建超. 玉米氮效率QTL定位和候选基因筛选[J]. 中国农业科学, 2024, 57(21): 4175-4191. |
[4] | 马荆鄂, 熊信威, 周敏, 吴斯琪, 韩甜, 饶友生, 王樟凤, 许继国. 基于垂体全长转录组测序分析鸡睾丸性状相关基因[J]. 中国农业科学, 2024, 57(20): 4130-4144. |
[5] | 尹军良, 李婧怡, 韩硕, 杨培华, 马佳伟, 刘奕清, 胡海骏, 朱永兴. 生姜NHX基因家族成员鉴定及其在硅缓解盐胁迫中的表达特征[J]. 中国农业科学, 2024, 57(19): 3848-3869. |
[6] | 陈菲儿, 张志鹏, 蒋庆雪, 马琳, 王学敏. 紫花苜蓿MsSPL17的克隆及生物学功能验证[J]. 中国农业科学, 2024, 57(17): 3335-3349. |
[7] | 刘桐, 王志荣, 李伟, 刘洋, 王祥儒, 赖弟利, 何毓琦, 张凯旋, 赵振军, 周美亮. 苦荞bHLH93转录因子响应铝胁迫的功能研究[J]. 中国农业科学, 2024, 57(16): 3127-3141. |
[8] | 陈文杰, 陈渊, 韦清源, 汤复跃, 郭小红, 陈淑芳, 覃夏燕, 韦荣昌, 梁江. 利用高世代转录组测序挖掘控制南方大豆皱叶症候选基因[J]. 中国农业科学, 2024, 57(15): 2914-2930. |
[9] | 蔚书涵, 秦小杰, 吴其超, 李玲, 臧德奎, 马燕. 芍药响应细极链格孢侵染的转录组分析[J]. 中国农业科学, 2024, 57(15): 3035-3052. |
[10] | 徐明睿, 王晓娟, 杨亚莉, 马越飞, 刘万茂, 孙颖. 基于转录组学分析辣椒对磷营养逆境的响应[J]. 中国农业科学, 2024, 57(14): 2827-2846. |
[11] | 徐梦雨, 王佳洋, 王江波, 唐雯, 陈一恒, 上官凌飞, 房经贵, 卢素文. 不同葡萄种质果皮中香气物质含量及基因表达的差异分析[J]. 中国农业科学, 2024, 57(13): 2635-2650. |
[12] | 苏国钊, 李嫒嫒, 刘中华, 陈宇华, 张秀杰, 马莹雪, 杨旭红, 邓超, 徐振江. 苦瓜品种SSR分子标记鉴定技术体系构建与应用[J]. 中国农业科学, 2024, 57(11): 2227-2242. |
[13] | 张海青, 张恒涛, 高启明, 姚家龙, 王亚荣, 刘珍珍, 孟祥鹏, 周喆, 阎振立. 转录组分析筛选调控苹果分枝能力的关键基因[J]. 中国农业科学, 2024, 57(10): 1995-2009. |
[14] | 肖涛, 李辉, 罗韦, 叶涛, 余欢, 陈友波, 石钰仕, 赵德鹏, 吴芸. 基于转录组测序筛选鸡蛋绿壳性状相关基因[J]. 中国农业科学, 2023, 56(8): 1594-1605. |
[15] | 李慧, 张雨峰, 李晓刚, 王中华, 蔺经, 常有宏. 全基因组DNA甲基化和转录组联合分析鉴定杜梨耐盐相关转录因子[J]. 中国农业科学, 2023, 56(7): 1377-1390. |
|