火龙果营养丰富,富含糖、有机酸、蛋白质、矿质元素、维生素、膳食纤维等[1-2],是我国南方重要的果树树种。近年来,随着火龙果产业迅速发展,用于栽培的火龙果良种较少,通过分子生物学手段进行品种遗传改良研究具有必要性。截至目前,火龙果的全基因组测序尚未完成,遗传信息还不清楚,与果实发育相关的基因种类和功能尚不明确,这些均阻碍了火龙果果实品质形成机理等相关基础研究的进行。淀粉和蔗糖是果实内在品质形成的主要因子。目前,对火龙果果实发育过程中淀粉和蔗糖代谢途径研究较少,参与淀粉和蔗糖代谢合成的相关酶及酶基因尚未见报道,淀粉和蔗糖代谢过程是如何进行催化反应还不清楚。探究火龙果果肉中淀粉和蔗糖合成代谢途径的分子调控机制,对改善果实品质具有重要的意义。
RNA测序[3]（RNA Sequencing,简称RNA- Seq）又称转绿组测序技术,具有灵敏度高、定量准确、可重复性高、检测范围广、分析更可靠等特点,被广泛应用。该技术具有获取全部mRNA转录本的丰度信息,从中挖掘功能基因和未知基因,解析基因调控状况,揭示不同生物学过程与基因表达之间关系等方面具有独特的优势[4-5]。近年来,火龙果在转录水平上的研究取得了一系列进展,主要集中在逆境（低温[6]、干旱[7]、水淹[8]、盐胁迫[9]）、火龙果果肉颜色[10-11]转录调控、甜菜碱[12-13]生物合成相关基因表达、病虫害等[14]方面,然而基于火龙果果肉不同发育时期转录组测序方面的研究尚未见报道。
本研究以3个不同发育时期的红肉火龙果果肉为研究对象,进行转录组测序研究,采用生物信息学技术方法,对得到的Unigene进行分类和功能注释,分析差异表达基因富集的代谢通路,对参与火龙果果实发育过程中淀粉和蔗糖代谢过程的差异表达的酶基因进行筛选及qPCR验证,旨在探索火龙果果实发育过程中关于淀粉和蔗糖代谢过程中酶的催化情况及酶基因的表达规律,为进一步实现人为调控火龙果果实糖分积累、提升果实品质提供理论依据。
1 材料与方法
1.1 材料
研究材料采自山西农业大学（山西省农业科学院）棉花研究所南花农场的15株红肉火龙果树,树龄为3 a,品种为‘蜜宝’,生长状况良好。通过前期对15株红肉火龙果树物候期观察,于2019年7月15日开始采样,每5 d采1次样,共采5次样,选取3个关键时期进行研究,即：2019年7月15日（幼果期,Y）、2019年7月26日（转色期,C）、2019年7月31日（成熟期,M）,每次分别在3株树的相同部位采集果实,采样完成后迅速返回实验室,将3个果实逐一去皮,切片（从果实中部横切约0.3 cm）,并立即将果肉放入液氮中冷冻,在‒80 ℃保存备用。
1.2 总RNA提取及检测、cDNA文库构建和Illumina测序
1.2.1 总RNA提取及检测 使用TRNzol-A* Reagent提取火龙果果肉中的总RNA。采用琼脂糖凝胶电泳分析火龙果果肉中RNA完整性及是否存在DNA污染,NanoPhotometer spectrophotometer检测RNA纯度,Agilent 2100 bioanalyzer检测RNA完整性,以备合格的火龙果果肉RNA进行转录组测序。
1.2.2 cDNA文库构建、质检及测序 火龙果果肉RNA检测合格后,构建cDNA文库。纯化后的双链cDNA,经过末端修复、加A尾并连接测序接头,用AMPure XP beads筛选250~300 bp左右的cDNA,进行PCR扩增并再次使用AMPure XP beads纯化PCR产物,最终获得文库。使用Agilent 2100 bioanalyzer对文库的insert size进行检测。库检合格后,把不同文库按照有效浓度及目标下机数据量的需求pooling后进行Illumina测序。
1.3 测序数据分析
1.3.1 数据质控 为了保证数据分析的质量及可靠性,需要对原始数据进行过滤。主要包括去除带接头的reads、含N的reads、低质量reads。同时,对clean data进行Q20、Q30和GC含量计算。后续所有分析均是基于clean data进行高质量分析。
1.3.3 参考序列比对 将Trinity拼接得到的转录组作为参考序列（Ref）,将每个样品的clean reads往Ref上做mapping,分别过滤掉比对质量值低于10的reads,非成对比对上的reads,比对到基因组多个区域的reads。比对采用RSEM[18]软件,RSEM中使用到的bowtie2参数mismatch0 （bowtie2默认参数）。
1.3.4 基因功能注释 基于Nr、Nt、Pfam、Swiss-Prot进行序列相似性检索蛋白质功能进行注释。利用GO、KEGG、KOG数据库分别进行分子功能及相关代谢通路的比对预测分析。
1.4 差异基因富集分析
1.5 荧光定量PCR分析
根据qPCR引物设计原则,采用Beacon Designer 7.9软件进行引物设计,设计各候选基因的qPCR引物（表1）。引物由魔獾生物科技（上海）有限公司合成。本研究以火龙果3个发育时期的果肉为材料,YLS8为内参基因[21],采用SYBR GREEN I检测进行qPCR反应,利用2‒ΔΔCT计算各个发育时期果肉的基因相对表达量。反应总体系为10 μL,其中包括2×SYBR® Green Supermix 5 μL、Forward prime 0.5 μL、Reverse primer 0.5 μL、cDNA 1 μL、ddH 2O 3 μL。反应程序为：95 ℃ 3 min;95 ℃ 10 s,60 ℃ 30 s,39个循环。扩增完毕后,以1 ℃/4 s的速率从60 ℃逐步递增到95 ℃,获得溶解曲线。反应过程中,每个样品设置3个生物学重复。
表1 引物信息Tab. 1 Primer information |
| 基因Gene | 引物序列（5°‒3°）Primer sequence (5°‒3°) | 退火温度Annealing temperature/℃ |
|---|---|---|
| Cluster-12747.21987 | F: TGATGACACCAGTAAGATTG R: AAGACGAGGGTATTGACT | 77.0 |
| Cluster-12747.10701 | F: TATGCCACAGCCAATGAC R: CGGAACAAGTATGCCTGAA | 77.0 |
| Cluster-12747.18766 | F: GCAAGCAATATGTCCAAAGG R: GAGCAACCAACCGTGTAA | 78.5 |
| Cluster-12747.5831 | F: ATTGGATGACACTTAATGAAC R: TGACTTAGGACACTGATTG | 78.0 |
| Cluster-12747.12274 | F: GGCTCCGCTACCATCTTG R: AGAAGTGGGTCCGAGAGG | 82.0 |
| Cluster-12747.16617 | F: AAGGAGTTATGATGATGTTGAA R: GTAAGGAGTATCTGGTGTGA | 74.0 |
| YLS8-1-F | F: GCTCGGTTATGTCCACAAGGTA R: CGGAAGAAGAGCGTCTCGTTA | 80.5 |
2 结果与分析
2.1 火龙果果肉转录组测序数据评估统计及拼接
表2 ‘蜜宝’火龙果果肉发育过程测序数据评估统计Tab. 2 Evaluation and statistics of sequencing data in process of pulp development of ‘Mibao’ dragon fruit |
| 样品名称 Sample | 原始数据中的 reads数 Raw_reads | 原始数据过滤后的reads数 Clean_reads | 原始数据过滤后的碱基数 Clean_bases/Gb | 数据整体测序 错误率 Error_rate/% | Q20/% | Q30/% | GC/% |
|---|---|---|---|---|---|---|---|
| Y_1 | 58 947 082 | 57 780 354 | 8.67 | 0.03 | 97.73 | 93.62 | 49.35 |
| Y_2 | 57 772 806 | 56 719 570 | 8.51 | 0.03 | 97.7 | 93.47 | 47.88 |
| Y_3 | 46 584 304 | 45 733 802 | 6.86 | 0.02 | 98.13 | 94.57 | 47.91 |
| C_1 | 57 990 928 | 56 993 728 | 8.55 | 0.02 | 98.08 | 94.25 | 50.65 |
| C_2 | 58 081 526 | 57 021 002 | 8.55 | 0.03 | 97.98 | 94.08 | 50.68 |
| C_3 | 63 609 790 | 62 405 802 | 9.36 | 0.03 | 98.00 | 94.09 | 50.71 |
| M_1 | 54 657 236 | 53 659 912 | 8.05 | 0.03 | 97.97 | 94.00 | 50.01 |
| M_2 | 55 216 912 | 54 398 678 | 8.16 | 0.03 | 98.01 | 94.05 | 49.92 |
| M_3 | 64 994 344 | 64 080 816 | 9.61 | 0.03 | 97.97 | 94.04 | 50.25 |
| 总计 | 517 854 928 | 508 793 664 | 76.32 | ‒ | ‒ | ‒ | ‒ |
| 平均 | ‒ | ‒ | ‒ | ‒ | 97.95 | 94.02 | 49.71 |
注：Y. 幼果期;C. 转色期;M. 成熟期。 | |
Note: Y. Young fruit stage; C. Color conversion period; M. Mature stage. |
表3 拼接长度频数分布情况Tab. 3 Frequency distribution of splicing length |
| 组别Group | 300~500 bp | 500~1000 bp | 1000~2000 bp | >2000 bp | 总数Total |
|---|---|---|---|---|---|
| Transcripts | 24 964 | 27 699 | 32 908 | 29 622 | 115 193 |
| Unigenes | 14 713 | 13 201 | 10 431 | 9510 | 47 855 |
表4 拼接长度分布情况Tab.4 Distribution of gene splicing length |
| 组别 Group | 最小长度 Min_length | 平均长度 Mean_length | 中间长度 Median_length | 最大长度 Max_length | N50 | N90 | 核苷酸总数 Total_nucleotides |
|---|---|---|---|---|---|---|---|
| Transcripts | 301 | 1483 | 1121 | 16 808 | 2164 | 682 | 170 881 251 |
| Genes | 301 | 1260 | 771 | 16 808 | 2000 | 508 | 60 309 167 |
2.2 火龙果3个不同果肉发育时期差异表达基因比较
‘蜜宝’火龙果3个不同果肉发育时期基因表达差异,两两比对,结果见表5。结果表明：转色期与幼果期相比,有7226条Unigene差异基因表达,其中3371个基因上调表达,3855个基因下调表达;成熟期与幼果期相比,共有7273个差异表达基因,其中上调基因3513个,下调基因3760个;成熟期与转色期相比,共有1635个差异表达基因,其中805个基因上调,830个基因下调。上述结果表明,火龙果果肉在成熟期的差异表达基因所占比例高于其他2个时期,说明DEGs在火龙果成熟期发育阶段起着非常重要的作用。
表5 ‘蜜宝’火龙果果肉差异表达基因数目统计Tab. 5 Statistical table on the number of differentially expressed in flesh of ‘Mibao’ dragon fruit |
| 差异基因比较组 DEG set | 所有差异基因 All DEG | 上调 Up-regulated | 下调 Down-regulated | 筛选软件及阈值 Threshold |
|---|---|---|---|---|
| C vs Y | 7226 | 3371 | 3855 | DESeq2 padj<0.05|log2FoldChange|>1 |
| M vs Y | 7273 | 3513 | 3760 | DESeq2 padj<0.05|log2FoldChange|>1 |
| M vs C | 1635 | 805 | 830 | DESeq2 padj<0.05|log2FoldChange|>1 |
注：Y. 幼果期;C. 转色期;M. 成熟期。 | |
Note: Y. Young fruit stage; C. Color conversion period; M. Mature stage. |
2.3 基因注释和功能分类
2.3.1 基因注释 获得转录本后,需对转录本进行功能注释,以获得gene的功能信息。为了解火龙果果肉转录本的注释情况,与7大数据库进行比对,共有47 855个Unigene成功注释（表6）,由于火龙果的基因组信息尚未公布,缺少参考基因组信息,部分序列未能注释。与NR数据库进行比对获得的基因注释信息最高,占全部注释序列的49.2%。与PFAM和GO数据库比对获得的基因注释信息次之,占全部注释序列均为38.82%。与KOG数据库比对获得的基因注释信息最少,占全部注释序列的12.96%。
表6 基因注释成功率统计Tab. 6 Statistics of success rate of gene annotation |
| 数据库 Database | 基因数 Number of unigenes | 百分比 Percentage/% |
|---|---|---|
| NR库中 | 23 548 | 49.20 |
| NT库 | 15 488 | 32.36 |
| KO库 | 8593 | 17.95 |
| SwissProt库 | 17 724 | 37.03 |
| PFAM库 | 18 582 | 38.82 |
| GO库 | 18 582 | 38.82 |
| KOG库 | 6203 | 12.96 |
| 所有数据库 | 3617 | 7.55 |
| 至少在一个数据库注释 | 26 866 | 56.14 |
| Unigenes总数 | 47 855 | 100.00 |
2.3.2 火龙果果肉发育过程中差异表达基因GO功能富集 ‘蜜宝’火龙果3个不同果肉发育时期差异表达基因经两两比对,18 582个Unigene在GO功能注释中被分为3大类,且每一类注释种类和数量各不相同。由 图1可知,生物过程主要集中在细胞过程、代谢过程和单一有机体过程。在生物过程中,多数差异基因主要参与到细胞过程、代谢过程、单一生物过程;细胞组件中,多数差异基因主要参与了细胞、细胞组成、细胞器、细胞器组成、大分子复合物、膜、膜组成;分子功能中,多数差异基因富集到绑定、催化活性。
2.3.3 火龙果果肉发育过程中差异表达基因KOG注释分类统计 对火龙果果肉发育过程中差异表达基因进行KOG数据库比对分析,6203个Unigene在KOG注释中被分为25个group,结果如图2所示。翻译后修饰、蛋白质周转、伴侣包含的Unigene最多,共854个;其次是一般功能预测,Unigene 803个;翻译,核糖体结构与生物发生,Unigene 538个;细胞内运输,分泌和囊泡运输,Unigene 483个;RNA加工和修饰,Unigene 450个;未知功能,Unigene 426个;信号转导机制,Unigene 420个;能量生产和转换,Unigene 352个;转录,Unigene 351个;脂质转运与代谢,Unigene 307个;碳水化合物转运与代谢,Unigene 298个;氨基酸转运与代谢,Unigene 289个;最少的是细胞运动性,注释到2个Unigene。可见,在火龙果果肉发育过程中,除了转录、翻译、修饰、运输等最基本的生命活动外,信号转导,脂质、碳水化合物及氨基酸转运与代谢占重要地位,说明信号转导,脂质、碳水化合物及氨基酸转运在火龙果果肉发育过程中起到重要作用。
图2 ‘蜜宝’火龙果果肉发育过程中差异表达基因KOG注释分类统计A：RNA加工和修饰;C：能源生产和转换;B：染色质的结构和动力学;E：氨基酸的运输和代谢;D：细胞周期调控、细胞分裂、染色体分割;G：碳水化合物的运输和代谢;F：核苷酸转运和代谢;I：脂质运输和代谢;H：辅酶运输和代谢;K：转录;J：翻译,核糖体结构与生物发生;M：细胞壁/膜/被膜生物发生;L：复制、重组和修复;O：翻译后修饰、蛋白质周转、伴侣;N：细胞运动;Q：次生代谢物生物合成、运输和分解代谢;P：无机离子运输和代谢;S：功能未知;R：仅一般功能预测;U：U细胞内运输、分泌和囊泡运输;T：信号转导机制;W：细胞外结构;V：防御机制;Y：核结构;Z：细胞骨架。 Fig. 2 Eukaryotic Ortholog Groups classification of differentially expressed genes in process of pulp development of ‘Mibao’dragon fruit A: RNA processing and modification; C: Energy production and conversion; B: Chromatin structure and dynamics; E: Amino acid transport and metabolism; D: Cell cycle control, cell division, chromosome partitioning; G: Carbohydrate transport and metabolism; F: Nucleotide transport and metabolism; I: Lipid transport and metabolism; H: Coenzyme transport and metabolism; K: Transcription; J: Translation, ribosomal structure and biogenesis; M: Cell wall/membrane/envelope biogenesis; L: Replication, recombination and repair; O: Posttranslational modification, protein turnover, chaperones; N: Cell motility; Q: Secondary metabolites biosynthesis, transport and catabolism; P: Inorganic ion transport and metabolism; S: Function unknown; R: General function prediction only; U: Intracellular trafficking, secretion, and vesicular transport; T: Signal transduction mechanisms; W: Extracellular structures; V: Defense mechanisms; Y: Nuclear structure; Z: Cytoskeleton. |
2.4 火龙果果肉发育过程中差异表达基因的KEGG分类
对基因做KO注释后,根据它们参与的KEGG代谢通路进行分类（图3）,共分为5类,分别为：细胞过程、环境信息处理、遗传信息处理、代谢及有机系统,其中代谢占的比例最大。在细胞过程中,以运输和分解代谢为主;在环境信息处理中,以信号转导所占比例最高,膜运输所占比例最小;遗传信息处理中,以翻译所占比例最高,复制和修复所占比例最小;代谢中,以碳水化合物代谢所占比列最高,多糖生物合成所占比例最小;在有机系统中,以环境适应为主。
2.5 火龙果3个不同果肉发育时期差异基因(DEGs)富集的代谢通路分析
KEGG是有关Pathway的主要公共数据库,通过Pathway显著性富集,能够确定火龙果果肉在3个不同发育时期的差异表达基因参与的主要代谢途径。
2.5.1 转色期对比幼果期 转色期与幼果期共有7226个差异基因（padj<0.05|&log2FoldChange|> 1）,其中上调基因3371个,下调基因3855个（图4A）。转色期对比幼果期的差异基因共参与到119种代谢通路,其中富集于前20的KEGG代谢通路见图5。转色期上调的差异基因主要富集在植物激素信号转导、丙酮酸代谢、光合生物的固碳作用、半乳糖代谢、过氧化物酶体、类黄酮生物合成、类胡萝卜素生物合成中,其中与类胡萝卜素生物合成相关的酶基因表达量均上调,包括胡萝卜素去饱和酶、玉米黄质环氧化酶、9-顺-环氧羰基双加氧酶、花椒毒素脱氢酶,推测转色期果肉颜色形成与类胡萝卜素生物合成有关（图5A）。转色期下调的差异基因主要富集在核糖体、淀粉和蔗糖代谢、苯丙素类生物合成、糖酵解/糖异生、氨基糖和核苷酸代谢、α-亚麻酸代谢、柠檬酸循环、脂肪酸生物合成、苯丙氨酸、酪氨酸和色氨酸生物合成等（图5B）。
图4 差异基因火山图A：转色期对比幼果期;B：成熟期对比幼果期;C：成熟期对比转色期。 Fig. 4 Differental gene volcano map A: Color conversion period vs. young fruit stage; B: Mature stage vs. young fruit stage; C: Mature stage vs. color conversion period. |
2.5.2 成熟期对比幼果期 成熟期与幼果期共有7273个差异基因（padj<0.05|&log2FoldChange|> 1）,其中上调基因3513个,下调基因3760个（图4B）。成熟期对比幼果期的差异基因共参与到119种代谢通路,其中富集于前20的KEGG代谢通路见图6。成熟期上调的差异基因主要富集在内质网中的蛋白质加工、淀粉和蔗糖代谢、剪接体、丙酮酸代谢、光合生物的固碳作用等（图6A）。成熟期下调的差异基因主要富集在核糖体、植物激素信号转导、淀粉和蔗糖代谢、苯丙素类生物合成、糖酵解/糖异生等。随着火龙果果实的生长发育,成熟期较幼果期的果实口感甜,这与成熟期上调的差异基因富集在淀粉和蔗糖代谢、光合生物的固碳作用、半乳糖代谢有关（图6B）。
2.5.3 成熟期对比转色期 成熟期与转色期共有1635个差异基因（padj<0.05|&log2FoldChange|> 1）,其中上调基因805个,下调基因830个（图4C）。成熟期对比转色期的差异基因共参与到95种代谢通路,其中富集于前20的KEGG代谢通路见图7。成熟期上调的差异基因主要富集在内质网中的蛋白质加工、苯丙素类生物合成、淀粉和蔗糖代谢、光合生物的固碳作用、半胱氨酸和蛋氨酸代谢、类黄酮生物合成、果糖和甘露糖代谢、甘氨酸、丝氨酸和苏氨酸代谢、花青素生物合成等（图7A）。成熟期下调的差异基因主要富集在内质网中的蛋白质加工、植物激素信号转导、半乳糖代谢、糖酵解/糖异生、光合生物的固碳作用等（图7B）。
2.6 火龙果果肉中淀粉和蔗糖代谢途径相关酶的基因统计及表达模式
2.6.1 火龙果果肉中淀粉和蔗糖代谢途径相关的差异基因统计 本研究根据火龙果在3个不同果肉发育时期基因表达水平的显著差异,统计了参与淀粉和蔗糖代谢途径相关的酶基因个数。结果表明：转色期与幼果期比对,参与淀粉和蔗糖代谢途径的差异表达基因共有80个,其中上调基因25个,下调基因55个,下调基因占69%,成熟期与幼果期比对,参与淀粉和蔗糖代谢途径的差异表达基因共有76个,其中上调基因28个,下调基因48个,下调基因占63%;成熟期与转色期比对,参与淀粉和蔗糖代谢途径的差异表达基因共有15个,其中上调基因7个,下调基因8个,下调基因占53%;通过对火龙果3个不同果肉发育时期差异基因比对整体来看,淀粉和蔗糖途径为下调基因富集的通路。
2.6.2 火龙果果肉中淀粉和蔗糖代谢途径相关的酶基因的表达模式 以DESeq2padj<0.05|log2FoldChange|>1为筛选条件,结合KEGG富集通路,在火龙果3个不同果肉发育时期筛选参与淀粉和蔗糖代谢途径的差异表达酶基因5个（表7）。
表7 火龙果果肉不同发育时期淀粉和蔗糖代谢途径相关酶基因的表达模式Tab. 7 Expression patterns of enzymes related to starch and sucrose metabolism pathway in different development stages of dagon fruit |
| KEGG编号KEGG number | EC编号 EC number | 基因ID Gene ID | 蛋白名 Protein name | 基因表达模式 Gene expression pattern | ||
|---|---|---|---|---|---|---|
| C vs Y M vs Y M vs C | ||||||
| K00963 | 2.7.7.9 | Cluster-12747.21987 | UTP--葡萄糖-1-磷 酸尿苷酰转移酶(UGP2) | 下调 | 下调 | 下调 |
| K00975 | 2.7.7.27 | Cluster-12747.10701 | 葡萄糖-1-磷酸腺苷 酸转移酶(glgC) | 上调 | 下调 | 下调 |
| K00703 | 2.4.1.21 | Cluster-12747.18766 | 淀粉合酶(glgA) | 下调 | 上调 | 上调 |
| K01188 | 3.2.1.21 | Cluster-12747.5831 | β-葡萄糖苷酶 | 下调 | 下调 | 上调 |
| K00695 | 2.4.1.13 | Cluster-12747.16617 | 蔗糖合酶 | 上调 | 下调 | 下调 |
2.7 qPCR验证
为了验证转录组测序数据的准确性,以YLS8为内参基因,筛选出参与淀粉和蔗糖代谢途径的关键酶基因UGP2、glgC、glgA、Cluster- 12747.5831、Cluster-12747.16617等5个,并采用RT-PCR技术对5个关键的酶基因在火龙果3个不同果肉发育时期进行验证分析。结果表明,RT-PCR验证的5个与淀粉和蔗糖代谢途径相关的酶基因在火龙果3个不同果肉发育时期上调或下调的变化倍数与转录组测序数据不一致,但基因表达量变化趋势与转录组测序结果中基因的表达量趋势基本一致（图8）。即UGP2在幼果期表达量较高,随着果肉发育表达量持续下降（图8A）;glgC在转色期表达量较高,幼果期次之,成熟期表达量较低（图8B）;glgA在成熟期表达量较高,幼果期次之,转色期表达量较低（图8C）;Cluster-12747.5831在幼果期表达量较高,转色期表达量较低,成熟期较转色期的表达量略有回升（图8D）;Cluster-12747.16617在幼果期和转色期的表达量较高且一致,成熟期表达量最低（图8E）。以上结果表明,转录组测序数据可靠性较高。同时发现,火龙果果肉中UGP2、glgC、glgA、Cluster-12747.5831、Cluster-12747.16617表达受到明显的发育调控。
3 讨论
本研究针对火龙果3个不同时期的果肉组织进行转录组测序,共得到76.32 Gb测序数据,组装后共得到115 193条转录本,对获得的Unigene进行七大数据库（NR、NT、KO、Swiss-prot、PFAM、GO、KOG）比对,共有47 855个Unigene成功注释。KEGG代谢通路注释结果表明,6554个Unigene获得功能注释,共有119条代谢途径,其中碳水化合物代谢所占比例最高。
结合KEGG分析,关于淀粉和蔗糖代谢途径相关酶发现,合成淀粉和蔗糖的代谢途径主要2条,一条代谢途径为：D-葡萄糖-1-磷酸在UTP-葡萄糖-1-磷酸尿苷酰转移酶（UGP2）催化作用下,生成UDP-葡萄糖,UDP-葡萄糖在糖原合酶（GYS）和颗粒结合型淀粉合成酶（WAXY）的共同催化作用下,生成直链淀粉,直链淀粉在1,4-α-葡聚糖分支酶（GBE1）催化作用下生成淀粉;另一条代谢途径为：D-葡萄糖-1-磷酸在葡萄糖-1-磷酸腺苷酸转移酶（glgC）催化作用下,生成ADP-葡萄糖,ADP-葡萄糖在淀粉合酶（glgA）和颗粒结合型淀粉合成酶（WAXY）共同催化作用下,生成直链淀粉,直链淀粉在1,4-α-葡聚糖分支酶（GBE1）催化作用下生成淀粉。
本试验筛选出了5个参与淀粉和蔗糖代谢途径的酶基因。其中,UTP-葡萄糖-1-磷酸尿苷酰转移酶（UGP2）在一条淀粉和蔗糖代谢途径的早期阶段发挥作用,以D-葡萄糖-1-磷酸为底物,催化形成UDP-葡萄糖。UGP2的表达影响着UGP2的活性及UDP-葡萄糖生成。本研究结果显示UGP2在火龙果果肉幼果期的表达量最高,随着果肉发育表达量持续下降。推测在火龙果果肉幼果期生成的UDP-葡萄糖最多,随着果肉的不断发育成熟,生成的UDP-葡萄糖逐渐减少,成熟期最少。葡萄糖-1-磷酸腺苷酸转移酶（glgC）在另一条淀粉和蔗糖代谢途径的早期阶段发挥作用,催化D-葡萄糖-1-磷酸生成ADP-葡萄糖。glgC的表达影响着glgC的活性及ADP-葡萄糖的生成。本研究中,glgC在火龙果果肉转色期的表达量最高,幼果期的表达量次之,成熟期的表达量较幼果期的略低。推测ADP-葡萄糖在火龙果果肉转色期生成最多,幼果期次之,成熟期ADP-葡萄糖的生成略低于幼果期ADP-葡萄糖的生成。
蔗糖合酶具有分解和合成蔗糖的作用[27],既可催化蔗糖合成又可催化蔗糖分解,是植物蔗糖代谢的关键酶。在火龙果果肉淀粉和蔗糖代谢途径中,蔗糖在蔗糖合酶的分解作用下生成UDP-葡萄糖、D-果糖;UDP-葡萄糖在蔗糖合酶的合成作用下生成蔗糖。Cluster-12747.16617的表达影响着蔗糖合酶的活性及蔗糖的分解和合成。张春华等[28]研究桃果实不同发育时期蔗糖合成酶基因表达情况,结果表明：花后45~65 d,蔗糖合成酶基因表达量略微提高,85 d时,表达量下降至最低,105 d时,表达量迅速提高。本研究表明,Cluster-12747.16617在火龙果果肉幼果期和转色期的表达量最高且一致,成熟期降低。这可能是由于不同种或火龙果为非跃变型果实的的原因。推测幼果期和转色期的果肉中以蔗糖合成为主,成熟期的果肉中以分解蔗糖为主。
综上所述,在火龙果果肉淀粉和蔗糖代谢途径中,5个基因在3个果肉发育时期的表达情况各不相同。UGP2、glgC在淀粉和蔗糖代谢途径的早期阶段发挥作用,为淀粉和蔗糖的合成提供中间产物;glgA调控淀粉的合成;Cluster-12747. 5831调控葡萄糖的合成;Cluster-12747.16617调控蔗糖的分解和合成。