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、
   129   130   131   132   133   134   135   136   137   138   139