Page 54 - 《广西植物》2026年第7期

１ １ ５ ０ 广 西 植 物 ４６ 卷
表 １ 不同分枝类型马尾松生长量差异
Ｔａｂｌｅ １ Ｇｒｏｗｔｈ ｉｎｃｒｅｍｅｎｔ ｄｉｆｆｅｒｅｎｃｅｓ ｉｎ Ｐｉｎｕｓ ｍａｓｓｏｎｉａｎａ ｓｅｅｄｌｉｎｇｓ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ｂｒａｎｃｈｉｎｇ ｔｙｐｅｓ
树高 主枝长 地径
类型 侧枝数
Ｔｒｅｅ ｈｅｉｇｈｔ Ｍａｉｎ ｂｒａｎｃｈ ｌｅｎｇｔｈ Ｇｒｏｕｎｄ ｄｉａｍｅｔｅｒ
Ｔｙｐｅ Ｎｕｍｂｅｒ ｏｆ ｌａｔｅｒａｌ ｂｒａｎｃｈｅｓ
(ｃｍ) (ｃｍ) (ｍｍ)
正常分枝型 ＣＫ １５９.２±７.４Ｂ ７５.１±１９.５ｂ ２８.６７±３.２ｂ ５.０±０.９ｂ
狐尾型 Ｈ ２１５.７±１８.２Ａ １７２.５±２７.５Ａ ３３.５０±４.２ａ １.５±１.５Ｃ
多干型 Ｍ ９５.０±２６.６Ｃ ３７.７±５.３ｃ ２３.６７±３.２ｃ ７.８±１.７ａ

注: 不同大写字母表示不同样品间差异极显著(Ｐ<０.０１)ꎻ 不同小写字母表示不同样品间差异显著(Ｐ<０.０５)ꎮ 下同ꎮ
Ｎｏｔｅ: Ｄｉｆｆｅｒｅｎｔ ｕｐｐｅｒｃａｓｅ ｌｅｔｔｅｒｓ ｉｎｄｉｃａｔｅ ｅｘｔｒｅｍｅｌｙ ｓｉｇｎｉｆｉｃａｎｔ ｄｉｆｆｅｒｅｎｃｅｓ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｓａｍｐｌｅｓ (Ｐ<０.０１)ꎻ ｄｉｆｆｅｒｅｎｔ ｌｏｗｅｒｃａｓｅ ｌｅｔｔｅｒｓ
ｉｎｄｉｃａｔｅ ｓｉｇｎｉｆｉｃａｎｔ ｄｉｆｆｅｒｅｎｃｅｓ ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｓａｍｐｌｅｓ (Ｐ<０.０５). Ｔｈｅ ｓａｍｅ ｂｅｌｏｗ.

２.３.２ 差异表达基因筛选 将各处理两两组合ꎬ筛
选 不 同 分 枝 类 型 马 尾 松 的 差 异 表 达 基 因
(ｄｉｆｆｅｒｅｎｔｉａｌｌｙ ｅｘｐｒｅｓｓｅｄ ｇｅｎｅｓꎬ ＤＥＧｓ)ꎮ 如 图 ５ 所
示ꎬＨ 与 ＣＫ 之间有 ２ ８９４ 个 ＤＥＧｓꎬ其中 １ ２０３ 个
上 调ꎬ １ ６９１ 个 下 调ꎻ Ｍ 与 ＣＫ 之 间 有 ３ ９５６ 个
ＤＥＧｓꎬ其中１ ８１８个上调ꎬ２ １３８ 个下调ꎻＨ 与 Ｍ 之
间有 １ １２９ 个 ＤＥＧｓꎮ 其中ꎬ在 Ｈ＿ｖｓ＿ＣＫ 和 Ｍ＿ｖｓ＿
ＣＫ 中均呈现差异表达的基因有１ ３６７个(图 ６)ꎮ
２.３.３ 差异表达的基因功能注释 将 ＤＥＧｓ 进一步
进行 ＫＥＧＧ 通路分析ꎬ结果( 图 ７) 发现ꎬＨ 和 ＣＫ
被富集到戊糖和葡萄糖醛酸酯的相互转化、赖氨
Ｈ. 狐尾型ꎻ Ｍ. 多干型ꎻ ＣＫ. 正常分枝型ꎮ 不同小写字母
酸生物合成、肽聚糖生物合成等通路ꎻＭ 和 ＣＫ 富
表示不同样品间差异显著(Ｐ<０.０５)ꎮ 下同ꎮ
集到 β￣丙氨酸代谢、赖氨酸生物合成等ꎻＨ 和 Ｍ 的
Ｈ. Ｆｏｘ￣ｔａｉｌ ｔｙｐｅꎻ Ｍ. Ｍｕｌｔｉ￣ｓｔｅｍ ｔｙｐｅꎻ ＣＫ. Ｎｏｒｍａｌ ｂｒａｎｃｈｉｎｇ
ｔｙｐｅ. Ｄｉｆｆｅｒｅｎｔ ｌｏｗｅｒｃａｓｅ ｌｅｔｔｅｒｓ ｉｎｄｉｃａｔｅ ｓｉｇｎｉｆｉｃａｎｔ ｄｉｆｆｅｒｅｎｃｅｓ
ＤＥＧｓ 富集到丁酸代谢、玉米素生物合成、精氨酸
ｂｅｔｗｅｅｎ ｄｉｆｆｅｒｅｎｔ ｓａｍｐｌｅｓ (Ｐ<０.０５). Ｔｈｅ ｓａｍｅ ｂｅｌｏｗ.
生物合成等ꎮ 此外ꎬＤＥＧｓ 富集到部分与植物激素
图 ２ 不同比较组植物激素含量比较
有关的通路ꎬ如玉米素生物合成、植物激素信号转
Ｆｉｇ. ２ Ｃｏｍｐａｒｉｓｏｎ ｏｆ ｐｌａｎｔ ｈｏｒｍｏｎｅ ｃｏｎｔｅｎｔｓ
导等ꎮ 结合分析ꎬ推测马尾松差异表达基因可能
ｉｎ ｄｉｆｆｅｒｅｎｔ ｃｏｍｐａｒｉｓｏｎ ｇｒｏｕｐｓ
通过玉米素生物合成、氨基酸代谢、糖代谢等显著
则是一种天然的 ＣＴＫꎮ 此外ꎬ在 Ｈ＿ｖｓ＿ＣＫ 比较组中 富集通路参与调控顶芽生长发育ꎮ
显著富集的色氨酸代谢途径ꎬ负责生成 ＩＡＡ 的合成 ２.３.４ 植物激素相关的差异表达基因 生长素有
前体ꎮ 综上所述ꎬ马尾松顶芽对分枝的调控涉及多 关途径的生长素输入载体 ＡＵＸ１( ｃｌｕｓｔｅｒ￣６６３７５.０、
种植物激素的协同作用ꎻＧＡ 与 ＡＢＡ 的合成与多干 ｃｌｕｓｔｅｒ￣６６３７５. ９、 ｃｌｕｓｔｅｒ￣５５６４４. ５ 等) 和 ＡＵＸ / ＩＡＡ
型、狐尾型马尾松的形成密切相关ꎬ而 ＩＡＡ 参与调 ( ｃｌｕｓｔｅｒ￣５５６４４.０、ｃｌｕｓｔｅｒ￣４７１６４.１、ｃｌｕｓｔｅｒ￣４９６２５. ０、
ｃｌｕｓｔｅｒ￣５０９６３.０) 等 ＤＥＧｓ 在 Ｈ 中呈上调趋势( 图
控狐尾型马尾松的顶芽生长和分枝过程ꎮ
２.３ 转录组测序结果与分析 ８)ꎻ细胞分裂素信号转导途径的 ２ 个 ＤＥＧｓ 分别被
２.３.１ 转录组测序和数据组装 不同样品间原始 注释为细胞分裂素受体 ＣＲＥ１( ｃｌｕｓｔｅｒ￣２１５８８.２) 和
数据分布范围为４２ ５６２ ０１２ ~ ６５ ４６４ ４８２ꎬ过滤后 磷酸转运蛋白 ＡＨＰ( ｃｌｕｓｔｅｒ￣３４０３４.３) 并在 Ｍ 中上
获得 不 同 样 品 的 过 滤 数 据 读 段 数 量ꎬ 范 围 为 调ꎻ赤霉素信号 转 导 途 径 上 的 赤 霉 素 受 体 ＧＩＤ１

４１ ７６５ ６３６ ~ ６４ ５４５ ８６２ꎬＱ２０ 不低于 ９８.１８％ꎬＱ３０ (ｃｌｕｓｔｅｒ￣３４５３４. ２ ) 和 核 心 负 调 控 因 子 ＤＥＬＬＡ
高于 ９４.８４％ꎬＧＣ 含量大于 ４４.１６ ％( 表 ３)ꎮ 这表 (ｃｌｕｓｔｅｒ￣５７８２２.１)在 Ｍ 中上调ꎬ在 Ｈ 中下调ꎻ 脱落
明样品转录组测序质量较好ꎬ可继续进行后续分 酸信号通路上的负调控因子蛋白磷酸酶 ＰＰ２Ｃ 注释
析验证ꎮ 到ｃｌｕｓｔｅｒ￣５３５１９.７、 ｃｌｕｓｔｅｒ￣６５２８７.０、ｃｌｕｓｔｅｒ￣５３５１９.６等