





中国农业科学 ›› 2025, Vol. 58 ›› Issue (6): 1173-1194.doi: 10.3864/j.issn.0578-1752.2025.06.010
孙萍1(
), 朱文灿2,3(
), 林贤锐1, 吴嘉颀1, 曹译文1, 陈辰斐1, 王轶1, 朱建锡1, 贾惠娟4, 钱敏杰2,3(
), 沈建生1(
)
收稿日期:2024-09-07
接受日期:2024-10-20
出版日期:2025-03-16
发布日期:2025-03-25
通信作者:
联系方式:
孙萍,E-mail:sunpingzju@163.com。朱文灿,E-mail:wencanzhu@hainanu.edu.cn。孙萍和朱文灿为同等贡献作者。
基金资助:
SUN Ping1(
), ZHU WenCan2,3(
), LIN XianRui1, WU JiaQi1, CAO YiWen1, CHEN ChenFei1, WANG Yi1, ZHU JianXi1, JIA HuiJuan4, QIAN MinJie2,3(
), SHEN JianSheng1(
)
Received:2024-09-07
Accepted:2024-10-20
Published:2025-03-16
Online:2025-03-25
摘要:
【目的】 桃的果皮色泽与果实的外观品质和经济价值紧密相关,花青苷是桃果皮着色的决定性色素物质。本研究基于代谢组和转录组联合分析,探究梅雨季节寡照多雨的环境条件对桃果皮类黄酮(花青苷、黄酮醇、原花青素)含量和类黄酮合成相关基因表达的影响,鉴定和挖掘桃类黄酮生物合成通路关键基因以及转录调控因子,为改善梅雨季节桃栽培措施、促进桃果皮着色以及进一步丰富桃类黄酮生物合成分子机制提供应用依据和理论基础。【方法】 以桃品种‘中金蟠7-12号’为试材。模拟多雨环境栽培为T1试验组,模拟寡照多雨环境栽培为T2试验组以及避雨栽培(模仿正常栽培环境)为CK对照组。在不同时期(0D和24D),对不同处理的桃果实果皮样品进行代谢组测定以及转录组测序(RNA-Seq)分析。进一步通过KEGG富集通路以及加权基因共表达网络分析(WGCNA)鉴定调节类黄酮化合物生物合成的关键候选基因。【结果】 代谢组测定结果显示,矢车菊素-3-O-葡萄糖苷、原花青素B1和槲皮素-3-O-葡萄糖苷分别为桃果皮中花青苷、原花青素和黄酮醇的主要成分。其中,矢车菊素-3-O-葡萄糖苷是桃果皮红色着色的决定性物质。T1和T2处理均抑制了桃果皮的花青苷含量,并且以T2处理的抑制效果更显著。通过RNA-Seq分析,共鉴定得到8 296个差异表达基因(DEG),其中24D-T1 vs 0D组间筛选得到的DEG数目最多,为6 879个。通过WGCNA,来自于绿松石、红色、绿黄色、棕色、蓝色和紫红色模块的基因被鉴定为参与调节桃果实果皮中类黄酮化合物生物合成的候选基因。KEGG富集分析显示,代谢途径是除了绿松石模块外,其他所有模块候选基因中被富集到最多的通路。基于WGCNA,鉴定确认了15个类黄酮生物合成通路相关结构基因,此外还鉴定得到了MYB、bHLH、ERF、bZIP以及C2H2等转录因子。【结论】 梅雨季节的多雨寡照显著抑制了桃果皮中花青苷的积累以及桃果实果皮的转红,避雨栽培可以在雨季提高桃果实的外观品质和经济价值。此外,鉴定得到一些与类黄酮生物合成密切相关的关键结构基因和调控基因,可为梅雨季节改善桃果实着色提供理论指导。
孙萍, 朱文灿, 林贤锐, 吴嘉颀, 曹译文, 陈辰斐, 王轶, 朱建锡, 贾惠娟, 钱敏杰, 沈建生. 基于代谢组和转录组解析多雨寡照对桃果皮着色和类黄酮积累的影响[J]. 中国农业科学, 2025, 58(6): 1173-1194.
SUN Ping, ZHU WenCan, LIN XianRui, WU JiaQi, CAO YiWen, CHEN ChenFei, WANG Yi, ZHU JianXi, JIA HuiJuan, QIAN MinJie, SHEN JianSheng. Effects of Rainy and Low Light Conditions on Coloration and Flavonoid Accumulation in Peach Peel Based on Metabolomic and Transcriptomic Analyses[J]. Scientia Agricultura Sinica, 2025, 58(6): 1173-1194.
表1
引物序列"
| 基因编号Gene ID | 正向引物Forward primer (5′-3′) | 反向引物Reverse primer (5′-3′) |
|---|---|---|
| Prupe.7G168300_v2.0.a1 | GGACTGGACACAGAGGCATT | GGGGCATTTTGGGTAGAAAT |
| Prupe.1G376400_v2.0.a1 | GGGGTGCTAGACATCCTCAA | CTGGTGCTCTTCGACATTCA |
| Prupe.4G200500_v2.0.a1 | TGAAGGGCTGTGATGGTGTA | TCCAAGGTTCCTTTGATTGC |
| Prupe.2G324700_v2.0.a1 | TGCCTCTCCCAACACTCTCT | CCATCAGCCACATCAAACAC |
| Prupe.2G212900_v2.0.a1 | TCTGGGGAGCAGAAGAAGAA | CAGGCATTGCATAGGGTTTT |
| Prupe.8G125100_v2.0.a1 | ACGCCCTATTGACACAGTCC | CTTGGGCCCCTTTTTGTTAT |
| Prupe.7G225600_v2.0.a1 | TGTCTCAGCCGTGAAAGATG | TGCAGAAGCAGAAGCAGAAA |
| Prupe.3G240000_v2.0.a1 | AAGGGTGGACTCGTTTTGTG | TCTTCTCCTCCCCTCTCGAT |
| β-actin | GATTCCGGTGCCCAGAAGT | CCAGCAGCTTCCATTCCAA |
图1
不同处理下桃果实的果皮表型以及代谢组学分析 A:3种处理下桃果实的代表性图像Representative images of peach fruits under three treatments。B:3种处理下桃果实果皮中不同类黄酮化合物组分含量The contents of different flavonoid compounds in peach fruit peel under three treatments。每个值代表3个生物学重复的平均值±标准差,无相同字母的数值表示差异显著Each value represents the mean±standard deviation of three biological replicates. There is no significant difference between values containing the same letter (P<0.05)"
表2
RNA-Seq文库的质量和产量统计"
| 样品 Sample | 原始数据 Raw reads | 过滤后数据 <BOLD>C</BOLD>lean reads | 过滤后碱基总数 Clean base (G) | 测序错误率 Error rate (%) | 碱基质量值 Q20 (%) | 碱基质量值 Q30 (%) | GC含量 GC content (%) |
|---|---|---|---|---|---|---|---|
| 0D-1 | 46947136 | 45207546 | 6.78 | 0.03 | 97.83 | 93.78 | 46.90 |
| 0D-2 | 45181526 | 43587342 | 6.54 | 0.02 | 98.16 | 94.61 | 46.71 |
| 0D-3 | 43048874 | 41413698 | 6.21 | 0.02 | 98.12 | 94.44 | 46.82 |
| 24D-T1-1 | 45453162 | 43941034 | 6.59 | 0.02 | 98.12 | 94.49 | 46.88 |
| 24D-T1-2 | 43462974 | 42019444 | 6.30 | 0.02 | 98.16 | 94.56 | 46.90 |
| 24D-T1-3 | 51322012 | 49684760 | 7.45 | 0.02 | 98.07 | 94.34 | 47.14 |
| 24D-T2-1 | 46961188 | 45358530 | 6.80 | 0.02 | 98.14 | 94.50 | 46.56 |
| 24D-T2-2 | 51255062 | 49645836 | 7.45 | 0.02 | 98.04 | 94.26 | 46.88 |
| 24D-T2-3 | 50412392 | 48881676 | 7.33 | 0.02 | 98.16 | 94.57 | 46.82 |
| 24D-CK-1 | 47734934 | 46239584 | 6.94 | 0.03 | 97.99 | 94.12 | 46.58 |
| 24D-CK-2 | 52248104 | 50503542 | 7.58 | 0.02 | 98.15 | 94.52 | 46.88 |
| 24D-CK-3 | 56671330 | 54632786 | 8.19 | 0.03 | 98.02 | 94.22 | 47.14 |
图2
RNA-Seq结果分析 A:样品相关性分析图,数值表示皮尔逊相关系数Sample correlation analysis chart, the value represents Pearson correlation coefficient;B:组间对比筛选到的DEG数目统计Statistics of the number of DEGs selected through comparison between different groups;C:DEG层次聚类分析热图,蓝色到红色代表基因表达水平从低到高The hierarchical clustering analysis heatmap of DEGs, with blue to red representing gene expression levels from low to high"
图5
GO富集分析 1:细胞过程Cellular process;2:代谢过程Metabolic process;3:刺激反应Response to stimulus;4:生物调节Biological regulation;5:生物过程调控Regulation of biological process;6:发展过程Developmental process;7:多细胞生物过程Multicellular organismal process;8:本地化Localization;9:生殖Reproduction;10:信号Signaling;11:生殖过程Reproductive process;12:物种间的生物过程Biological process involved in interspecies;13:生物过程的负向调节Negative regulation of biological process;14:生物过程的正向调节Positive regulation of biological process;15:发育Growth;16:细胞解剖实体Cellular anatomical entity;17:含蛋白质复合物Protein-containing complex;18:结合Binding;19:催化活动Catalytic activity;20:转运活动Transporter activity;21:ATP依赖性活动ATP-dependent activity;22:转录调节活动Transcription regulator activity;23:结构分子活动Structural molecule activity;24:分子功能调节活动Molecular function regulator activity;25:分子换能器活动Molecular transducer activity;26:抗氧化活动Antioxidant activity;27:翻译调节活动Translation regulator activity;28:分子衔接子活动Molecular adaptor activity;29:蛋白质折叠伴侣Protein folding chaperone;30:分子载体活动Molecular carrier activity;31:细胞骨架运动活动Cytoskeletal motor activity;32:小分子传感器活动Small molecule sensor activity;33:蛋白标签Protein tag"
| [1] |
刘冠儒. 湖北省孝感市桃产业现状与发展研究[D]. 武汉: 武汉轻工大学, 2023.
|
|
|
|
| [2] |
李晓颖, 谭洪花, 房经贵, 韩键, 宋长年. 果树果实的风味物质及其研究. 植物生理学报, 2011, 47(10): 943-950.
|
|
|
|
| [3] |
doi: 10.1093/jxb/erq442 pmid: 21278228 |
| [4] |
doi: 10.1039/b311404j pmid: 15282635 |
| [5] |
|
| [6] |
|
| [7] |
pmid: 16669768 |
| [8] |
|
| [9] |
|
| [10] |
|
| [11] |
|
| [12] |
|
| [13] |
doi: 10.1093/oxfordjournals.aje.a009257 pmid: 9247006 |
| [14] |
pmid: 11121513 |
| [15] |
doi: 10.1016/s0031-9422(00)00235-1 pmid: 11130659 |
| [16] |
doi: 10.1016/j.tplants.2013.06.003 pmid: 23870661 |
| [17] |
|
| [18] |
|
| [19] |
|
| [20] |
|
| [21] |
|
| [22] |
|
| [23] |
|
| [24] |
|
| [25] |
孟宪宸, 崔彦志, 李星宇, 张秀玲, 马振伟, 柴菊华. 成熟期降水对碣石山产区‘马瑟兰’果实品质的影响. 中外葡萄与葡萄酒, 2023(3): 24-31.
|
|
|
|
| [26] |
姜琳琳, 王静, 张晓煜, 陈仁伟, 胡宏远. 成熟期降水对贺兰山东麓酿酒葡萄果实品质的影响. 中国农业气象, 2020, 41(3): 156-161.
|
|
|
|
| [27] |
沈玉英, 熊彩珍, 贾惠娟. 盖棚对湖景蜜露桃着果及果实品质的影响. 中国南方果树, 2006(5): 67-68.
|
|
|
|
| [28] |
李玲, 王慧, 张梅, 陈修德, 高东升. 设施栽培桃果实芳香成分的GC-MS分析. 山东农业大学学报(自然科学版), 2011, 42(1): 41-48.
|
|
|
|
| [29] |
|
| [30] |
成梦园. 不同避雨膜对鲜食葡萄生长发育和果实品质的影响研究[D]. 太谷: 山西农业大学, 2022.
|
|
|
|
| [31] |
|
| [32] |
|
| [33] |
|
| [34] |
熊彩珍, 凌柏芳, 李洁, 景筱荣. 避雨设施栽培对桃果实生长发育及糖分的影响. 浙江林业科技, 2012, 32(3): 46-49.
|
|
|
|
| [35] |
|
| [36] |
徐长宝, 陈国强. 江淮梅雨对梨树生长结果及果实品质的影响. 中国南方果树, 2007(6): 74-75.
|
|
|
|
| [37] |
马艳儿, 杨丽, 何玉云, 李华, 王华. 降水对葡萄多酚类物质的影响. 江苏农业科学, 2016, 44(5): 221-224.
|
|
|
|
| [38] |
|
| [39] |
|
| [40] |
李舒婷, 李小溪, 高媛, 潘秋红. 简易避雨栽培对‘赤霞珠’葡萄果实类黄酮物质表征的影响. 果树学报, 2018, 35(7): 802-816.
|
|
|
|
| [41] |
曹锰, 郭景南, 魏志峰, 高登涛, 程大伟, 孙晓文. 避雨栽培对‘金手指’葡萄果实生长及香气物质组分的影响. 果树学报, 2015, 32(5): 894-902.
|
|
|
|
| [42] |
|
| [43] |
杨乐, 李文杨, 赵师成, 孙耀清, 申洁梅. 避雨栽培措施对信阳五月鲜桃果实生长的影响. 经济林研究, 2019, 37(1): 205-210.
|
|
|
|
| [44] |
薄流齐洋. 避雨栽培对刺葡萄生长发育及果实品质的影响[D]. 长沙: 湖南农业大学, 2021.
|
|
|
|
| [45] |
刘亚妮. 避雨栽培对葡萄生长发育及果实品质的影响[D]. 秦皇岛: 河北科技师范学院, 2019.
|
|
|
|
| [46] |
赵云. 两品种桃果皮花青苷积累差异及光照调控机制[D]. 杭州: 浙江大学, 2019.
|
|
|
|
| [47] |
|
| [48] |
|
| [49] |
doi: 10.1104/pp.105.064238 pmid: 16169968 |
| [50] |
doi: 10.3389/fpls.2015.00243 pmid: 25914714 |
| [51] |
|
| [52] |
|
| [53] |
|
| [54] |
doi: 10.1111/nph.17669 pmid: 34351627 |
| [55] |
doi: 10.1111/pbi.14024 pmid: 36752420 |
| [56] |
|
| [57] |
doi: 10.3389/fpls.2015.00908 pmid: 26579158 |
| [58] |
|
| [59] |
|
| [60] |
|
| [61] |
|
| [62] |
|
| [63] |
|
| [64] |
|
| [65] |
|
| [66] |
doi: 10.1111/tpj.12792 pmid: 25688923 |
| [1] | 岳丽昕, 王清华, 王振宝, 尼玛琼吉, 刘泽洲, 孔素萍, 张立峰, 高莉敏. 基于广靶代谢组学分析藏葱和细香葱营养品质及黄酮类代谢产物差异[J]. 中国农业科学, 2026, 59(5): 1070-1086. |
| [2] | 王忠妮, 雷月, 李佳丽, 宫彦龙, 朱速松. ABC转运蛋白OsARG1调控水稻抽穗期的功能[J]. 中国农业科学, 2026, 59(1): 1-16. |
| [3] | 王思琪, 邹利人, 白瑞雯, 闫可, 王思洋, 齐晓光, 申海林, 温景辉. 赤霉素调控‘蜜汁’葡萄穗轴硬化关键基因的挖掘[J]. 中国农业科学, 2026, 59(1): 179-189. |
| [4] | 张义茹, 韩雪, 姚鑫杰, 冯军, 魏爱丽, 李文超, 张彬, 韩渊怀, 李红英. 基于多组学解析谷子后熟米色变化的分子机制[J]. 中国农业科学, 2025, 58(9): 1702-1718. |
| [5] | 谭西北, 兰徐颖, 刘崇怀, 樊秀彩, 姜建福, 孙磊, 李鹏, 余书鑫, 张颖. 不同抗性葡萄响应白腐病侵染的次生代谢物变化[J]. 中国农业科学, 2025, 58(9): 1767-1778. |
| [6] | 吴郁, 曲翔汝, 杨丹, 伍芩, 陈国跃, 江千涛, 魏育明, 许强. 广泛非靶向代谢组学解析小麦抗条锈病反应中叶绿体代谢产物[J]. 中国农业科学, 2025, 58(7): 1333-1343. |
| [7] | 杨彩丽, 李永洲, 贺亮亮, 宋银花, 章鹏, 刘肇先, 李鹏慧, 刘三军. 葡萄TPS基因家族全基因组鉴定及VvTPS4在单萜形成中的功能验证[J]. 中国农业科学, 2025, 58(7): 1397-1417. |
| [8] | 邹晓威, 夏蕾, 朱晓敏, 孙辉, 周琦, 齐霁, 张亚封, 郑岩, 姜兆远. 基于转录组测序的玉米瘤黑粉菌UM01240过表达菌株诱导玉米抗病性分析[J]. 中国农业科学, 2025, 58(6): 1116-1130. |
| [9] | 谢露露, 李福, 张思远, 高建昌. 基于跨物种转录组解析影响不定根发生的保守基因[J]. 中国农业科学, 2025, 58(6): 1195-1209. |
| [10] | 高岩浩, 王婷婷, 白卫卫, 杜兴杰, 刘贤, 秦本源, 付彤, 孙宇, 高腾云, 张天留. 脂质组与转录组联合揭示南阳牛不同肌肉组织脂质特征的差异表达模式[J]. 中国农业科学, 2025, 58(6): 1239-1258. |
| [11] | 郭奥琳, 林俊璇, 赖恭梯, 贺丽媛, 车建美, 潘若, 杨方学, 黄玉吉, 陈桂信, 赖呈纯. 过表达VdF3′5′H2对刺葡萄细胞花青素组分积累的影响[J]. 中国农业科学, 2025, 58(4): 802-818. |
| [12] | 罗朝丹, 冯春梅, 李建强, 黎新荣, 韦勇, 杨李益, 刘筱瑾, 谭和, 任二芳, 罗小杰. 基于GC-MS非靶向代谢组学分析不同成熟度‘台农一号’芒果香气物质[J]. 中国农业科学, 2025, 58(3): 564-581. |
| [13] | 王凡, 刘陈玮, 陆红臣, 徐仁超, 卞晓春. 蚕豆响应交链格孢侵染的转录组分析及VfPR4的抗病功能验证[J]. 中国农业科学, 2025, 58(22): 4656-4672. |
| [14] | 穆赢通, 路景诗, 张雨桐, 石凤翎. 基于转录组和WGCNA的直立型花苜蓿抗旱关键基因识别[J]. 中国农业科学, 2025, 58(21): 4528-4543. |
| [15] | 潘媛, 王德, 刘楠, 孟祥龙, 戴蓬博, 李波, 胡同乐, 王树桐, 曹克强, 王亚南. 两种高通量测序技术鉴定苹果病毒效果评价及两种新病毒的鉴定[J]. 中国农业科学, 2025, 58(2): 266-280. |
|
||