Page 126 - 《广西植物》2025年第11期
P. 126

2 0 7 0                                广  西  植  物                                         45 卷
                                                 表 1  RT ̄qPCR 引物序列
                                            Table 1  Primer sequences for RT ̄qPCR
               基因名称             基因 ID                 正向引物序列 (5′-3′)                 反向引物序列 (5′-3′)
              Gene name         Gene ID            Forward primer sequence (5′-3′)  Reverse primer sequence (5′-3′)
                PR ̄1         Rhdel11G0226100          TGCGCTGCTGGCAAGG             AAGTAATGAACCACCACCCG
                 PYR         Rhdel09G0282100        TGATGTGTCGGAGCGATTGG           CCTAAAGCAATGAGGCACGC
                PP2C         Rhdel05G0017700      GAGTCAGTGCTTACAACCCACGTA       GCTAATAGACACTAAAGGACGCTTG
                 ABF         Rhdel05G0102200          GCCCCCGGTGCGAATTT            TCCCCTCATTTGCACCAGTA
                MYC2         Rhdel01G0237700      TCCAACCCACCAACAACAAAAAC          CAGAGAACCCGAGCTGAAAG
                 JAZ         Rhdel09G0151900        CAACACAATGCGGTTTCGGT           AGGTATTGTTGCTTCGCGGT


                                              表 2  转录组测序数据统计结果

                                    Table 2  Statistical results of transcriptome sequencing data
                                                      平均错误率        Q20 值      Q30 值      GC 含量       匹配率
               样品编号         原始数据          干净序列
                                                       Average    Q20 value  Q30 value  GC content  Mapped ratio
              Sample No.   Raw reads     Clean reads
                                                     error rate (%)  (%)       (%)        (%)        (%)
                 CK1       72 796 176    72 432 168     0.02       98.24      94.59       44.11      93.87
                 CK2       66 000 982    65 724 896     0.02       98.33      94.86       44.45      94.04
                 CK3       72 965 128    71 765 216     0.02       98.50      95.26       44.41      94.34
                H3 ̄1       61 128 826    60 858 472     0.02       98.18      94.45       44.68      92.41
                H3 ̄2       63 427 652    63 171 888     0.02       98.30      94.71       44.53      92.62
                H3 ̄3       70 936 282    70 393 708     0.02       98.20      94.47       44.47      92.54
                H6 ̄1       60 341 510    59 855 082     0.02       98.19      94.45       44.30      92.64
                H6 ̄2       75 965 538    75 248 812     0.02       98.19      94.42       44.38      92.96
                H6 ̄3       69 575 994    69 258 608     0.02       98.33      94.79       44.89      92.65



            1 539个ꎬH6 vs CK 特异表达的基因有 2 960 个ꎬ                  集到分子功能的有 7 个亚类ꎬ主要包括阴离子结
            H6 特异表达的基因数量高于 H3ꎬ说明高温胁迫 6                         合(749 个 基 因)、核 苷 酸 磷 酸 盐 结 合 ( 748 个 基
            d 后马缨杜鹃可能通过调控更多基因表达以响应                             因)和核苷酸结合(748 个基因)ꎻ富集到生物学过
            逆境ꎮ 同时ꎬH3 vs CK 与 H6 vs CK 共同表达的差                  程的有 11 个亚类ꎬ主要包括大分子代谢过程(691
            异基因为 5 435 个( 图 1:B)ꎮ 这些组间共同差异                     个基因)、有机物生物合成过程(575 个基因)、有
            表达的基因ꎬ可能是马缨杜鹃响应高温胁迫的重                              机环化合物代谢过程(530 个基因) 等ꎻ富集到细

            要基因ꎬ对马缨杜鹃应对高温胁迫具有重要意义ꎮ                             胞组分 的 有 2 个 亚 类ꎬ 主 要 包 括 膜 结 合 细 胞 器
            对组间差异表达基因火山图可视化ꎬ在 H3 vs CK                         (695 个基因) 和细胞内细胞器(799 个基因) ( 图
            鉴定出 6 974 个 DEGsꎬ其中上调的 DEGs 为 3 227                2)ꎮ GO 富集分析显示了富集的前 20 条通路( 图
            个ꎬ下调的 DEGs 为 3 747 个ꎻ在 H6 vs CK 鉴定出                3)ꎬH3 富集的基因主要涉及有机物生物合成过
            8 395个 DEGsꎬ其中上调的 DEGs 为 3 797 个ꎬ下调                程、响应刺激、脂质代谢过程、碳水化合物代谢过

            的 DEGs 为 4 598 个(图 1:C、D)ꎮ                         程ꎬH6 富集基因最多的是磷代谢过程、磷酸化、质
            2.3 差异表达基因 GO 注释与功能富集分析                            膜、脂质代谢过程ꎬ并且马缨杜鹃 DEGs 中富集的
                 对 H3 vs CK 与 H6 vs CK 共 同 差 异 表 达 的           前 20 个 GO 项大部分属于生物学过程ꎮ
            5 435个基因进行 GO 注释ꎬ结果表明这些差异基                         2.4 差异表达基因 KEGG 注释与富集分析
            因被注释到分子功能、生物学过程和细胞组分 3                                 为进一步探究差异基因的生物学功能ꎬ分析
            个大类别ꎮ 这 3 个大类共包括 20 个亚类ꎬ其中富                        差异基因响应高温胁迫的代谢通路ꎬ对 H3 vs CK
   121   122   123   124   125   126   127   128   129   130   131