Page 134 - 《广西植物》2026年第5期
P. 134
8 6 6 广 西 植 物 46 卷
通过 blastx 将 116 389 个 Unigenes 序列与 NR、 表 3 NR 注释物种分布
SwissProt、KEGG 和 COG/ KOG( E value <0.000 01) Table 3 Species distribution of NR annotation
等数 据 库 进 行 同 源 序 列 比 对ꎮ 共 得 到 注 释 的 物种 数量 百分比
Percentage
Unigenes 有 80 130 个ꎬ约占总 Unigenes 的68.85%ꎬ Species Number (%)
36 259 个没有得到注释ꎬ约占总基因31.15%ꎮ 结果 阿根廷豆树 Prosopis alba 19 752 24.92
显示马大相思在 NR 数据库中注释到 79 249 个 欧洲栓皮栎 Quercus suber 10 767 13.59
Unigenesꎬ在 SwissProt 数据库中注释到 59 675 个 奥氏卡普斯孢菌(ATCC 30864 菌株)
4 050 5.11
Capsaspora owczarzaki (ATCC 30864 strain)
Unigenesꎬ 在 KEGG 数 据 库 中 注 释 到 67 790 个
银叶罗丹木 Rhodamnia argentea 1 814 2.29
Unigenesꎬ 在 KOG 数 据 库 中 注 释 到 51 083 个
蓝星睡莲 Nymphaea colorata 1 469 1.85
Unigenesꎬ有47 246个 Unigenes 在 4 大数据库都能注
费氏洪氏藻 Hondaea fermentalgiana 1 299 1.64
释到(图 1)ꎮ
微小亚得里亚共生藻
Symbiodinium microadriaticum 1 288 1.63
密花豆 Spatholobus suberectus 1 083 1.37
红车轴草 Trifolium pratense 928 1.17
其他 Other 36 799 46.43
包含的 Unigenes 基 本 持 平ꎬ 分 别 为 23 984 个 和
22 911个ꎬ占该类比的 20.61%和 19.68%ꎮ
2.2.2 差异表达基因 KEGG 分析 将差异表达基
因(DEGs)经 KEGG 数据库比对后ꎬ共有 67 790 条
Unigenes 能够在 KEGG 数 据 库 得 到 注 释ꎮ KEGG
注释被划分为代谢、遗传信息处理、环境信息处
图 1 Unigenes 的 4 大数据库注释韦恩图 理、细胞过程和有机体系统ꎮ 其中ꎬ全局和概览图
Fig. 1 Venn diagram of annotations across
谱(9 779 个ꎬ占14.43%)的 Unigenes 占总 Unigenes
four major databases of Unigenes
数最多ꎬ其次是翻译(8 962 个ꎬ占 13.22%) 和碳水
化合物代谢(3 688个ꎬ占 5.44%)ꎬ而膜转运(178
利用 NR 数据库对 79 249 个 Unigenes 进行同 个)所获得的注释最少(图 3)ꎮ
源序列对比ꎬ共匹配到 524 个物种ꎮ 其中ꎬ同源性 2.3 差异表达基因筛选与分析
最高的 是 阿 根 廷 豆 树 ( Prosopis alba)ꎬ 有 19 752 以 A 组根系和叶片为对照组ꎬ与 B、C、D、E 组
条ꎬ占 24.92%ꎻ其次为欧洲栓皮栎(Quercus suber)ꎬ有 的根系和叶片表达量进行比较ꎬ并以 | log FC | 且
2
10 767条ꎬ占13.59%ꎻ密花豆(Spatholobus suberectus)、 FDR<0.05 作为 DEGs 的筛选条件ꎬ结果在根部筛选
红车轴草(Trifolium pratense) 的同源序列较少ꎬ仅分 出有 45 584 个 DEGsꎻ在叶部筛选出 10 391 个 DEGsꎮ
别占总序列的 1.37%、1.17%(表 3)ꎮ 筛选出的 DEGs 均随着硝态氮含量的增多而增加ꎬ并
2.2.1 差异表达基因 GO 分析 由图 2 可知ꎬ功能 且根与叶的 DEGs 变化趋势基本一致(图 4)ꎮ
注释的 Unigeneꎬ分别包含在 53 个亚类当中ꎮ 生物 2.4 差异表达基因 KEGG 显著富集分析
学过程包含 24 个亚类ꎬ其中细胞过程、代谢过程 2.4.1 根差异表达基因 KEGG 显著富集分析 由
富集数量较多ꎬ分别为 32 796 个和 32 218 个ꎬ占 表 4 可知ꎬ核糖体相关通路在 4 个组比较( AR ̄vs ̄
该类 别 的 28. 18%、 27. 68%ꎬ 单 生 物 体 过 程 有 BR、 AR ̄vs ̄CR、AR ̄vs ̄DR、AR ̄vs ̄ER) 中均显著富
26 531个 Unigenesꎬ占 22.80%ꎮ 细胞组分包含 17 集ꎬ 其 中 AR ̄vs ̄ER DEGs 最 多 ( 2 016 个ꎬ P<
个功能亚类ꎬ其中细胞、细胞部分以及膜结构占比 0.001)ꎬAR ̄vs ̄BR DEGs 最少(870 个ꎬP<0.001)ꎮ
较多ꎬ分别为 35 281 个、33 815 个和 28 142 个ꎬ占 2.4.2 叶差异表达基因 KEGG 显著富集分析 由
该类别的 30.31%、29.05% 和 24.18%ꎮ 分子功能 表 5 可知ꎬ 代谢途径与次生代谢物生物合成相关
包含 12 个功能亚类ꎬ 其中催化活性与结合活性所 通路在 4 个组比较(AL ̄vs ̄BL、AL ̄vs ̄CL、AL ̄vs ̄DL、

