Page 83 - 《广西植物》2025年第7期

７ 期 何斌等: 刺梨全长转录组测序分析及黄酮类化合物生物合成相关基因挖掘 １ ２ ７ ５

表 ２ ｉｓｏｆｏｒｍｓ 功能注释统计情况 进行 ＧＯ 注 释 之 后ꎬ 将 注 释 成 功 的 １３５ １１５ 个
Ｔａｂｌｅ ２ Ｓｔａｔｉｓｔｉｃｓ ｏｆ ｆｕｎｃｔｉｏｎａｌｌｙ ａｎｎｏｔａｔｅｄ ｉｓｏｆｏｒｍｓ Ｕｎｉｇｅｎｅ 序 列 按 照 ＧＯ ３ 个 大 类 [ 生 物 学 过 程
基因数量 百分比 ( ｂｉｏｌｏｇｉｃａｌ ｐｒｏｃｅｓｓꎬ ＢＰ )ꎻ 细 胞 成 分 ( ｃｅｌｌｕｌａｒ
数据库
Ｎｕｍｂｅｒ ｏｆ Ｐｅｒｃｅｎｔａｇｅ
Ｄａｔａｂａｓｅ ｃｏｍｐｏｎｅｎｔꎬ ＣＣ)ꎻ 分 子 功 能 ( ｍｏｌｅｃｕｌａｒ ｆｕｎｃｔｉｏｎꎬ
ｇｅｎｅｓ (％)
ＭＦ)]的下一层级进行分类(３２ 个功能组)ꎬ分类
ＫＥＧＧ ２０ ６２８ ８２.５０
Ｎｒ ２４ ７７８ ９９.１０ 结果如图 ５ 所示ꎮ 生物学过程包括 ２２ 个功能组ꎬ
Ｓｗｉｓｓ￣Ｐｒｏｔ ２０ ８７４ ８３.４９ 最多的是 细 胞 过 程 ( １６ １４４ 个)ꎬ 代 谢 过 程 次 之
ＴｒＥＭＢＬ ２４ ７５２ ９９.００ (１３ ６１４个)ꎬ最少的则是色素沉着( 仅 ２２ 个)ꎻ细
ＫＯＧ １５ ３４４ ６１.３７
胞组分仅有 ２ 个功能组ꎬ分别是细胞有机体的一
ＧＯ ２２ ６９２ ９０.７６
部分(１９ ８８４ 个)和含蛋白质的复合物(３ ２８１个)ꎻ
Ｐｆａｍ ２２ ３８３ ８９.５２
分子功能有 ２０ 个功能组ꎬ结合活性最多(１４ ３７６
已在至少 １ 个数据库中注释 ２４ ８５９ ９９.４２
个)ꎬ催化活性次之(１１ ９２０ 个)ꎬＲＮＡ 折叠伴侣最
Ａｎｎｏｔａｔｅｄ ｉｎ ａｔ ｌｅａｓｔ ｏｎｅ ｄａｔａｂａｓｅ
合计 Ｔｏｔａｌ ２５ ００３ １００.００ 少(仅 ２ 个)ꎮ

图 ３ Ｎｒ 注释同源物种统计
Ｆｉｇ. ３ Ｓｔａｔｉｓｔｉｃｓ ｏｆ Ｎｒ ａｎｎｏｔａｔｉｏｎ ｈｏｍｏｌｏｇｏｕｓ ｓｐｅｃｉｅｓ ｄｉｓｔｒｉｂｕｔｉｏｎ

２.３ 黄酮类化合物生物合成代谢途径相关基因的挖掘 出的转录因子进行本地 ＢＬＡＳＴꎬ找到 ５５ 个参与
刺梨作为药食两用的水果ꎬ含有丰富的黄酮 调节刺梨黄酮类化合物生物合成的潜在转录因
类化合物ꎮ 分析 ＫＥＧＧ 注释结果ꎬ共统计到 １５１ 子( 表 ４) ꎬ包括 ＷＲＫＹ 家族 ４１ 个、ＭＹＢ 家族 １２
条代谢途径ꎬ其中参与黄酮类化合物生物合成过 个及 ｂＨＬＨ 家族 ２ 个ꎮ 这些转录因子可为今后黄
程中相关转录本 ９９ 条(表 ３)ꎮ 这些转录本为深入 酮类化合物生物合成的研究提供一定的参考价
探究刺梨黄酮类化合物的潜在药用价值提供了一 值ꎮ 此外ꎬＭＹＢ￣ｒｅｌａｔｅｄ 和 ＷＲＫＹ 可 能 参 与 刺 梨
定的参考ꎮ 超 氧 化 物 歧 化 酶 的 生 物 合 成ꎻ Ｃ２Ｈ２、 ＧＲＡＳ、
２.４ 转录因子预测 ＷＲＫＹ、Ｃ３Ｈ 等可能参与刺梨的生长发育和对环
对去冗余后的序列进行转录因子预测ꎬ共鉴 境胁迫的应答ꎮ
定出 ８２ 种转录因子家族ꎬ共 １ ９３０ 个基因ꎮ 数量 ２.５ ＣＤＳ 预测
最多的前十个转录因子家族如图 ６ 所示ꎬｂＨＬＨ、 经过 ＣＤＳ 预测ꎬ共得到 ２４ ３５７ 个 ＣＤＳꎬ最大长
ＭＹＢ￣ｒｅｌａｔｅｄ 及 ＷＲＫＹ 可能参与刺梨黄酮类化合 度为 ７ ７５８ ｂｐꎬ最小长 度 为 ２８５ ｂｐꎬ平 均 长 度 为
物和萜类化合物的生物合成ꎬ将已知参与调控黄 １ ７２７ ｂｐꎮ 大部分 ＣＤＳ 长度集中在 ３００ ｂｐ 至３ ０００
酮类化合物生物合成的转录因子与本研究鉴定 ｂｐ 之间ꎬ占比 ９１.８４％(图 ７)ꎮ