‘金湘玉’黄桃果实糖与类黄酮代谢及转录组测序分析

7
--
生物技术•遗传育种　引用格式：陈为峰，王春发，黄　佳，等.‘金湘玉’黄桃果实糖与类黄酮代谢及转录组测序分析[J]. 湖南农业科学，2023（12）：7-12.　DOI:
DOI:10.16498/ki.hnnykx.2023.012.002果实的甜酸风味主要由糖和有机酸的含量及比例决定，是影响水果感官特性的重要因素[1-2]。

黄酮醇类化合物是重要的次生代谢物质，不仅能影响植物的生长发育，在植株抗菌防病方面也发挥了重要作用[3]。

近年来，随着分子生物技术和生物信息学的发展，转录组学的应用越来越广泛，对果实
品质的研究也已进入转录水平阶段[4-5]。

例如：李文彬[6]在‘红阳’猕猴桃果实中发现与蔗糖相关的基因SPS 、SPP 以及与代谢相关的基因AI 、NI 、SuSy 、FK 和HK ，这些基因的表达水平都随着果实糖含量的增加而升高。

Zhang 等[7]从桃果实中克隆出蔗糖代谢酶及碳水化合物变化相关基因CINV1、CWINV1、CINV2、SUS1、SUT1及SPS1。

周兰[8]的研究表明，CHI 、F3H 、P AL 、FLS 等基因在类黄酮化合物合成中起作用，但果实不同部位这些基因的表达量不同，在栽培型苹果的果肉中，类黄酮物质的含量远远低于果皮。

这些与生物代谢相关的研究
‘金湘玉’黄桃果实糖与类黄酮代谢及转录组测序分析
陈为峰1,3，王春发2,3，黄　佳1，李　杜1，张良波1
（1.湖南省园艺研究所，湖南 长沙 410125；2.衡阳市农业科学院，湖南 衡阳 421100；3.衡
南鸿祥种养专业合作社，湖南 衡阳 421125）
摘　要：采用高效液相色谱法、RNA-seq 转录组测序等方法研究了‘金湘玉’黄桃（Prunus persica 'Jinxiangyu'）果实不同发育时期（花后20~90 d ）可溶性糖组分及类黄酮物质的含量变化，筛选出调控果实可溶性糖与类黄酮化合物形成的关键基因。

结果表明：花后60~85 d 为‘金湘玉’果实品质形成的关键阶段，果实中可溶性糖、可溶性固形物、果糖、蔗糖含量在此期间明显上升；成熟果实糖含量以果糖和蔗糖为主；山奈素苷在果实中的含量较低,槲皮素苷含量在幼果期含量较高，随果实成熟含量降低；在果实中不含有花色苷化合物；参与类黄酮代谢的相关基因（ANS 、DFR 、F3H 、FLS 、4CL1等）在果肉中表达量较低；在果实发育过程中共鉴定到181个差异表达基因参与5条糖代谢途径，其中SDH 基因表达水平较高且随着果实成熟而持续上调，主要促进果实发育后期果糖含量的积累；SPS1、SS 与SS1基因表达水平随着果实成熟而持续上调，对蔗糖调控起着关键作用；SUS3、INVA 基因在果实发育前期表达水平较高，主要促进蔗糖的降解。

关键词：黄桃；糖代谢；类黄酮化合物；差异基因表达
中图分类号：S662.1；Q811.4 文献标识码：A 文章编号：1006-060X （2023）12-0007-006
 T ranscriptome Sequencing for Sugar and Flavonoid Metabolism in Prunus persica 'Jinxiangyu '
C HEN Wei-feng 1,3, W ANG Chun-fa 2,3, HUANG Jia 1, LI Du 1, ZHANG Liang-bo 1
（1. Hunan Horticultural Research Institute, Changsha 410125, PRC; 2. Hengyang Academy of Agricultural Sciences, Hengyang 421100,
PRC; 3. Hengnan Hongxiang Breeding Professional Cooperative, Hengyang 421125, PRC ）Abstract ：High performance liquid chromatography (HPLC) and RNA-seq were employed to study the changes of soluble sugars and flavonoids in the fruits of Prunus persica 'Jinxiangyu' at different developmental stages (20–90 d post flowering) and screen out the key genes regulating the formation of soluble sugars and flavonoids in fruits. The results showed that 60–85 d after flowering was the key stage for the fruit quality formation, and the content of soluble sugars, soluble solids, fructose, and sucrose in the fruits increased significantly during this stage. The sugars in mature fruits were mainly fructose and sucrose. The content of kaempferol glycosides was low in fruits. The content of quercetin glycosides was higher in young fruits and decreased as fruits became mature. There was no anthocyanins in the fruits. The genes (e.g., ANS, DFR, F3H, FLS , and 4CL1) involved in flavonoid metabolism showed low expression levels in fruit flesh. A total of 181 differentially expressed genes during fruit development were identified, and these genes participated in 5 glucose metabolic pathways. Among them, SDH showed a high expression level and up-regulated expression as the fruits became mature, which promoted the accumulation of fructose in late fruit development. In addition, the fruit maturity was accompanied by the up-regulated expression of SPS1, SS , and SS1, which played a key role in sucrose regulation. The high expression levels of SUS3 and INVA in the early stage of fruit development promoted the degradation of sucrose.
Key words ：yellow peach; glucose metabolism; flavonoids; differential gene expression
收稿日期：2023–10–19
基金项目：湖南省农业科技创新资金项目（2022CX113；2022CX84-09）；湖南创新型省份建设专项（2022NK4226）
作者简介：陈为峰（1988—），男，湖南安乡县人，助理研究员，主要从事果树研究。

通信作者：李　杜，张良波
　湖南农业科学（HUNAN AGRICULTURAL SCIENCES ）2023年12月8--生物技术•遗传育种
结果为探明果实生长发育过程有机物积累的机理奠定了坚实基础。

‘金湘玉’黄桃（Prunus persica 'Jinxiangyu'）是湖南省林业科学院最新选育出的湖南省首个自主知识产权的早熟黄桃新品种，为了探明该品种黄桃果实的糖成分及类黄酮物质的积累代谢规律，笔者分析了相关代谢关键基因的表达差异，以期挖掘出调控果实优良品质形成的关键基因，为实现黄桃分子育种，加快黄桃新品种更新提供参考。

1　材料与方法1.1　试验材料
供试‘金湘玉’黄桃取自湖南省沅陵县太常安置区果园（110°21'8.09"E ，28°28'7.15"N ，海拔120 m ）。

采样是从桃树谢花后20 d 开始，每隔10 d 采摘1次，果实开始成熟后5 d 采摘1次，直到果实完全成熟。

从树势健壮、栽培管理基本一致的3棵‘金湘玉’果树树冠中部位置上分别采摘5个大小均匀、无损伤、无虫害的果实，去皮切碎后用液氮速冻，置于－70℃超低温冰箱保存备用。

1.2　试验方法
1.2.1　果实可溶性糖含量的测定　可溶性固形物含量参照NY/T2637—2014采用数显手持折光仪进行测定；可溶性总糖含量采用蒽酮硫酸法测定[9]；果实中蔗糖、葡萄糖、果糖、山梨醇含量采用索莱宝科技有限公司的相关试剂盒进行测定。

1.2.2　果实黄酮类化合物含量的测定　称取0.25 g 果实样品，用1.5 mL 酚类提取液（水∶甲醇∶甲酸= 25∶24∶1，体积比，下同）研磨混匀，10 400 r/min 离心20 min ；经0.22 μm 滤膜过滤上清液，然后用液相色谱检测类黄酮含量。

流动相A （乙腈∶甲酸= 9∶1）与流动相B （水∶甲酸=9∶1）的梯度洗脱程序设置为：95%A 0 min 、85%A 25 min 、78%A 42 min 、64%A 60 min 、95%A 65 min ，最后平衡5 min ，共70 min ；流速为1 mL/min ；分析柱和保护柱为ODS-3色谱柱（5.0 μm ，4.6 mm×250 mm ），柱
温35℃；在365 nm 波长下检测黄酮醇类化合物的含量。

1.2.3　转录组测序　根据‘金湘玉’果实发育过程糖分的积累规律，选取花后30 d （A2）、花后60 d （A5）、花后90 d （A9）3个不同发育期的果实，作为RNA 提取和转录组测序的材料。

采用Nano Drop 2000分光光度计检测RNA 浓度与纯度，RNA 条带完整性质量结果用1%琼脂糖凝胶电泳进行检验；RNA 质检合格后，用 Illumina Hiseq 测序平台进行高通量测序，该步骤委托广州基迪奥生物科技有限公司完成。

参考基因组版本为Prunus_persica_genome_v2.0.a1（https:///species/prunus_persica/genome_v2.0.a1）。

1.2.4　qRT-PCR 验证　取1μL 质检合格的RNA 样品按照M5 Super plus qPCR RT kit with gDNA remover 试剂盒说明书反转录合成cDNA 。

利用Primer 3 Plus 软件对7个基因的CDS （Coding sequence ）序列设计特异性引物，选择Pp RPL13作为内参基因。

引物序列见表1，所有引物均由北京擎科生物科技有限公司合成。

qRT-PCR 反应体系（20 μL ）：Template DNA 2 μL ，2×M5 HiPer SYBR Premix EsTaq （withTli RNaseH ） 10 μL ，上下游引物各0.4 μL ，ROX Refer-ence Dye （50×）or ROX Reference Dye II 0.4 μL ，ddH 2O 6.8 μL 。

qRT-PCR 反应程序：95℃ 30 s 预变性；95℃ 5 s ，60℃ 30 s ，40个循环。

按照2-ΔΔct 法计算基因相对表达水平，每个样品进行3次生物学重复。

1.3　数据处理
采用Microsoft Excel 2010软件处理原始数据，SPSS19.0软件进行相关性分析，Primer3Plus 在线软件设计引物，Origin 2021软件作图或进行图片处理。

2　结果与分析
2.1　‘金湘玉’黄桃果实糖积累动态变化规律
由图1可知，在花后20~90 d 期间，随着黄桃果实的不断发育，可溶性总糖含量总体呈上升趋
表1　基因引物序列
基因ID
基因注释上游引物（5'—3'）下游引物（5'—3'）Prupe.8G264300_v2.0.a1SUS3ACCGAGCACGTAATGGTGAG TCATGGCCTCCACAACTGTC Prupe.7G192300_v2.0.a1SS CGCACTTGTTGTTGAGGAGC GCTTTGCAGAAAGGTGACGG Prupe.2G288800_v2.0.a1SDH ACCCTGCAGATCTGTGCTTC CGCCGAGAGACTTTGCAATG Prupe.3G056600_v2.0.a1PFK3TCAACGAGCCATAAATGCAG CCACCAGGTCCTTCCAGATA Prupe.2G309400_v2.0.a1SS1ACCTTGCACATCAGGGAGTG TGTGTCAAGAGCATGGGTCC Prupe.1G268300_v2.0.a1AMY3CTCTACCGACGCAGCTACAG ACCACTTCAGTCCGCTTGAG Prupe.1G483200_v2.0.a1
SPS1
TCTGTAACAGTGGCGGTGAG AGAAACCCAACGAACCAGTG PpRPL13
60S.ribosomal protein L13
GCAGCGACTGAAGACATACAA
GGTGGCATTAGCAAGTTCCTC
陈为峰等：‘金湘玉’黄桃果实糖与类黄酮代谢及转录组测序分析
9--
生物技术•遗传育种势，可溶性固形物含量总体呈先降后升的趋势。

花后20~50 d ，黄桃果实的葡萄糖、蔗糖、果糖含量均呈先升后降的趋势；花后50~90 d ，蔗糖和果糖含量
均迅速上升，果糖含量在花后85 d 达到峰值，含量为54.38 mg/g ，蔗糖含量在花后90 d 达到峰值；山梨醇含量随果实成熟整体呈现下降趋势，在花后85 d 含量最低。

图1　‘金湘玉’黄桃果实发育过程糖含量的变化
2.2　‘金湘玉’黄桃果实类黄酮物质的积累
动态变化
由图2可知，‘金湘玉’黄桃果实中，槲皮素苷含量在花后20~50 d 含量较高，随着果实逐渐成熟其含量有所下降；山奈素苷含量总体上呈先升后降的变化趋势，但波动较为平稳，在花后80 d 达到最
大值。

图2　‘金湘玉’黄桃果实发育过程类黄酮物质含量的变化
2.3　‘金湘玉’黄桃果实转录组测序数据的
分析
对处于发育关键期（花后30、60、90 d ）的‘金湘玉’果实样品进行转录组测序，结果如表2所示，共获得3 612万~ 4 489万条Raw reads ，单个样品平均获得40 173 908条Raw reads ，经质量过滤后平
均得到39 969 440个Clean reads ，且碱基Q20百分比在96.59%~97.97%之间，均超过了95%，Q30的百分比在91.09%~93.97%之间，均超过了90%，核苷酸GC 含量在45.81%~47.56%之间，且数据整体测序错误率小于0.03%，说明转录组测序质量较好，测序结果可进行后面的数据组装。

2.4　‘金湘玉’黄桃果实糖代谢关键基因的筛选
对3个发育关键时期样品的转录组数据进行筛选，结果如图3所示，共发现181个基因参与果实果糖和甘露糖代谢（ko00051）、糖酵解/糖异生代谢（ko00010）、戊糖和葡萄糖醛酸的相互转化（ko00040）、半乳糖代谢（ko00052）、淀粉和蔗糖代谢（ko00500）。

其中，有26个差异基因参与果实果糖和甘露糖代谢，有34个差异基因参与果实糖酵解/糖异生代谢，有30个差异基因参与戊糖和葡萄糖醛
表2　‘金湘玉’黄桃果实转录组测序数据统计
样本名称Raw reads （条）Clean reads （条）过滤后碱基数（bp ）
Q20（%）Q30（%）GC （%）A2-139 018 61838 828 028 5 783 794 99197.7493.6647.01A2-238 154 49237 969 664 5 666 865 78597.6593.3946.94A2-337 562 91237 395 360 5 580 009 23497.8193.7946.92A5-136 750 28836 578 308 5 463 934 32897.6493.2846.79A5-244 891 31444 677 610 6 677 184 79797.3992.6547.56A5-344 059 29243 862 328 6 553 007 98297.7293.3546.73A9-136 121 51235 868 706 5 354 698 27496.5991.0945.81A9-244 626 22044 329 026 6 621 931 71896.7791.3745.93A9-3
40 380 526
40 215 932
6 00
7 467 848
97.97
93.97
45.83
注：A2、A5、A9分别表示花后30、50、90 d 的黄桃样品。

　湖南农业科学（HUNAN AGRICULTURAL SCIENCES ）2023年12月10--生物技术•遗传育种
酸的相互转化，有22、69个差异基因分别参与半乳糖代谢、淀粉和蔗糖代谢。

由此可知，在‘金湘玉’黄桃果实中参与淀粉和蔗糖代谢的差异基因数量最
多，参与半乳糖代谢的差异基因数量最少。

图3　糖代谢途径差异表达基因个数统计
（横坐标通路编号从左到右依次代表：戊糖和葡萄糖醛酸的相互转化、果糖和甘露糖代谢、半乳糖代谢、淀粉和蔗糖代谢、糖酵解/糖异生代谢途径）
2.5　糖代谢关键基因在‘金湘玉’果实3个发育关键时期的表达分析
进一步对‘金湘玉’果实3个发育关键时期转录组数据进行两两对比，以3个生物学重复的基因表达量求平均值，筛选满足|log2(fc)|＞1、P ＜0.05、FDR ＜0.05的显著表达的差异基因，结果如图4所示，在淀粉与蔗糖代谢过程中INVA 、INV*DC4和SUS3在果实发育的过程中表现为显著下调，SUS3、INVA 基因在花后30 d 的表达量较高，随果实的成熟在花后90 d 几乎不表达；
SS 、SS1、SPS1基因在果实发育的过程中表现为显著上调，花后90 d 表达量较高，这与果实发育过程中蔗糖逐渐积累的结果一致。

在果糖和甘露糖代谢过程中，PMI1在果实发育的过程中表现为显著下调；
SDH 、PFK3基因在果实发育的过程中表现为显著上调，其中SDH （Prupe.2G288800_v2.0.a1）基因在3个时期的表达量都很高，FPKM 值分别为697.63、1 140.78、1 853.14。

2.6　类黄酮代谢关键基因在‘金湘玉’果实
3个发育关键时期的表达分析
由图5可知，在类黄酮生物合成过程中，ANS 、DFR 、F3H 及FLS 基因的表达量在花后30~60 d 呈显著上调趋势，花后60~90 d 呈显著下调趋势，基因表达量的FPKM 值在花后90 d 很低；
CHI 基因的表达量在花后30~90 d 表现为显著下调。

苯丙烷类物质生物合成中，4CL3基因的表达量在花后30~90 d 呈显著下调趋势，4CL1、P AL1基因的表达量在花后30~60 d 呈显著上调趋势，在花后60~90 d 呈显著
下调趋势。

2.7　qRT-PCR 验证测序结果
为了进一步验证‘金湘玉’果实转录组测序结果的准确性，从RNA-Seq 测序数据中挑取与糖代谢相关的7个差异表达基因（SUS3、SS 、SDH 、PFK3、SS1、AMY3、SPS1），通过qRT-PCR 对基因的表达量进行检验，结果如表3所示，这7个基因的qRT-PCR 相对表达量与转录组测序的FPKM 值相对表达量的动态变化趋势基本一致。

这表明转录组测序表达谱结果是可靠的。

表3　不同方法对‘金湘玉’黄桃果实各发育时期差异表达
基因表达量的检测结果
基 因
检测方法A2A5A9SDH 转录组测序697.63 1 140.78 1 853.14qRT-PCR 2.718 113.722 318.981 2SPS1转录组测序0.46 4.5236.67qRT-PCR 0.000 70.032 20.208 1SS 转录组测序167.24263.16 1 750.87qRT-PCR 1.638 3 4.338 8 2.866 1SUS3转录组测序439.96296.83160.16qRT-PCR 2.984 2 1.476 20.271 9PFK3转录组测序54.5471.91141.51qRT-PCR 0.064 60.431 70.962 9SS1转录组测序21.244.4180.82qRT-PCR 0.004 30.007 30.009 9AMY3
转录组测序48.913176.26qRT-PCR
0.004 3
0.039 9
0.072 1
注：转录组测序结果为FPKM 值，qRT-PCR 结果为2-ΔΔct 值。

3　讨　论
果实的内外在品质是水果质量评估的重要指标，是影响果农收益和商品市场竞争力的重要因素[10-11]。

前人研究发现，不同植物果实中的糖组分及含量不同，甚至同一种果实不同品种间的糖组分及含量也存在较大差异。

张颖[12]的研究发现，4种不同桃品种间果糖含量和蔗糖含量存在差异，花后90 d ，‘菁香’桃果实蔗糖含量较高，‘锦绣’黄桃果糖含量最高。

前人研究还发现，桃果实栽培种的糖类型大多以蔗糖为主。

该试验结果显示，‘金湘玉’黄桃幼果发育期的蔗糖含量极低，该时期糖分主要以果糖和山梨醇为主，山梨醇含量在整个果实生长发育过程中呈现明显的下降趋势，说明山梨醇进入果实后不是一直以山梨醇的形式存在，而是会逐渐转化为果糖和葡萄糖；成熟的‘金湘玉’黄桃果实中果糖的含量最高，使其具有较高的口感甜度。

大多数桃品种的果实生长期为110~140 d ，而‘金湘玉’的果实生长期较短，为85 d 左右，为早熟品种，由于‘金湘玉’果实发育时间短，因此成熟果实中蔗糖积累较少。

陈为峰等：‘金湘玉’黄桃果实糖与类黄酮代谢及转录组测序分析
11--
生物技术•遗传育种图4　‘金湘玉’果实3个发育关键时期糖代谢基因的差异表达分析
（A2、A5、A9分别表示花后30、50、90 d 的黄桃样品，下同）
图5　‘金湘玉’果实3个发育关键时期类黄酮代谢基因的差异表达分析
糖代谢是在相关酶基因的调节下完成的，果实生长发育不同阶段相关代谢基因表达的差异可能导致不同阶段糖积累出现差异[13-14]。

该研究对‘金湘玉’果实发育3个不同时期的转录组数据进行差异基因分析发现，具有显著表达差异的基因中SS 与SS1基因的表达量水平相差较大，SS 基因的表达量远远高于SS1，这2个基因都可以调控果实糖分的积累，其中SS 基因的调控作用更明显；在果实发育过程中，SPS1与SPS4的表达量变化趋势不同，其中SPS1基因在花后30和60 d 几乎不表达，在花后90 d 表达量升高，SPS4花后30和90 d 表达量较低，在花后60 d 表达量稍有升高；
INV*DC4和INVA 均只在花后30 d 表达量很高，在花后90 d 几乎不表达，且INVA 基因表达量远高于INV*DC4。

前人研究表明，幼果中的SUS1更接近典型的SUS 。

在植物
中，SUS1在蔗糖的分解中发挥主要作用，为细胞壁构建或糖酵解提供底物，随着果实的成熟，SUS2活性增加，并在蔗糖的快速积累中发挥作用[15]。

该试验结果显示，SUS2在整个发育过程表达量极低（3个时期FPKM 值分别为0.30、25.78、0.01），SUS3基因在幼果期表达量较高（FPKM 值达439.96），但随着果实的成熟呈下调趋势，因而推测INVA 、SUS3基因在花后30 d 是使蔗糖分解为果糖和葡萄糖的关键基因；2个SDH 基因的表达量在花后30~90 d 均表现为上升趋势，其中SDH （Prupe.2G288800_v2.0.a1）表达水平较高且持续上调，推测该基因对果实发育后期果糖含量的积累调控作用较强。

综上所述，随着果实的成熟，蔗糖合成方向基因SPS1、SS 与SS1的表达水平持续上调，对果实中蔗糖含量起着关键作用；
SUS3、INVA 基因在果实发育前期表
　湖南农业科学（HUNAN AGRICULTURAL SCIENCES ）2023年12月12--生物技术•遗传育种
达水平较高促进蔗糖的降解，这些基因的调控影响着‘金湘玉’果实糖类含量的变化。

该研究还发现，参与类黄酮物质代谢的相关基因如ANS 、DFR 、F3H 、FLS 、4CL1和P AL1等基因在‘金湘玉’果肉中的表达量较低，这与果肉中测出槲皮素苷及山奈素苷含量较低的结果一致，而这些结构基因相互作用共同影响了类黄酮物质的合成。

参考文献：
[1] AUBERT C ，BONY P ，CHALOT G ，et al. Effects of storage
temperature ，storage duration ，and subsequent ripening on the physicochemical characteristics ，volatile compounds ，and phytochemicals of Western Red nectarine （Prunus persica L. Batsch ）[J]. Journal of Agricultural and Food Chemistry ，2014，62（20）：4707-4724.
[2] ZHANG P ，SHAO X F ，WEI Y Y ，et al. At-harvest fruit maturity
affects sucrose metabolism during cold storage and is related to chilling injury in peach[J]. Journal of Food Science and Technology ，2020，57（6）：2000-2009.
[3] 陈为峰. 苹果属植物山奈素苷的鉴定及其对花粉管伸长和种子数
量的影响[D]. 杨凌：西北农林科技大学，2020.
[4] 梁　烨，陈双燕，刘公社. 新一代测序技术在植物转录组研究中
的应用[J]. 遗传，2011，33（12）：1317-1326.
[5] 岳桂东，高　强，罗龙海，等. 高通量测序技术在动植物研究领
域中的应用[J]. 中国科学：生命科学，2012，42（2）：107-124.
[6] 李文彬. ‘红阳’猕猴桃果实发育转录组及花青素累积机理研究
[D]. 武汉：中国科学院研究生院（武汉植物园），2015.
[7] ZHANG C H ，SHEN Z J ，ZHANG Y P ，et al. Cloning and
expression of genes related to the sucrose-metabolizing enzymes and carbohydrate changes in peach[J]. Acta Physiologiae Plantarum ，2013，35（2）：589-602.
[8] 周　兰. 苹果果实发育中类黄酮含量变化及相关基因的研究[D].
北京：中国农业科学院，2013.
[9] 曹建康，姜微波，赵玉梅. 果蔬采后生理生化实验指导[M]. 北京：
中国轻工业出版社，2007.
[10] 彭少兵，郭军战，陈铁山，等. 树莓和黑莓开花结实规律的初步
观察[J]. 西北林学院学报，2002，17（2）：45-46，51.
[11] HUAN C ，ZHANG J ，JIA Y ，et al. Effect of 1-methylcyclopropene
treatment on quality ，volatile production and ethanol metabolism in kiwifruit during storage at room temperature[J]. Scientia Horticulturae ，2020，265：109266.
[12] 张　颖. 四个桃品种生物学特性及果实生长发育规律研究[D].
长沙：中南林业科技大学，2021.
[13] 廖光联，钟　敏，黄春辉，等. 果实糖代谢及其相关酶基因研究
进展[J]. 江西农业大学学报，2020，42（1）：187-195.
[14] 袁　晖，韦　云，李馨玥，等. ‘南果梨’及其芽变‘南红梨’
果实中糖分积累与相关基因表达差异分析[J]. 果树学报，2017，34（5）：534-540.
[15] 卢合全，沈法富，刘凌霄，等. 植物蔗糖合成酶功能与分子生物
学研究进展[J]. 中国农学通报，2005，21（7）：34-37，57.
（责任编辑：成　平）。