Page 180 - 《广西植物》2025年第3期
P. 180
５ ６ ０ 广 西 植 物 ４５ 卷
最终 获 得 长 ４４ ０２５ ｂｐ 的 序 列ꎮ １２ 个 物 种 共 ２.２.２ 遗传多样性 １０ 个微卫星位点在 １２ 个物种
４４ ０２５ ｂｐ 长度的序列上只有 ４８ 个变异位点( 占 中共检测到 １７９ 个等位基因ꎬ平均每位点等位基
总序列长度的 ０.１１％)ꎬ其中单变异位点为 ２６ 个 因数为 ７.８ꎬ建润楠种群最少(平均为 ５.０)ꎬ红楠种
(占总序列长度的 ０.０６％)、简约信息位点为 ２２ 个 群最多(１０.４)ꎮ 物种的独有等位基因数都很少ꎬ
(占 ０.０５％)ꎮ 在 ７２ 个叶绿体 ＤＮＡ 序列片段中ꎬ 其中红楠种群的独有等位基因最多( 平均为每位
只有 １３ 个片段具多态性ꎮ 因此ꎬ后续使用这 １３ 个 点 ０.７ 个)ꎬ基脉润楠和建润楠中没有独有等位
叶绿体 ＤＮＡ 序列片段组成长度为 １０ ７５６ ｂｐ 的序 基因ꎮ
列矩阵进行系统进化分析ꎮ 对 １２ 个润楠属物种 １２ 个物种中观察杂合度和期望杂合度中等ꎬ
１３ 个序列的异质性检验表明ꎬ序列间替换型式异 期望杂合度为 ０.５０１( 粗壮润楠) ~ ０.８３８ ( 琼桂润
质性在多数物种对间不显著ꎬ不会影响系统进化 楠)ꎬ观察杂合度为 ０.４４２ (粗壮润楠) ~ ０.７１８ ( 红
分析ꎮ 对已测序的红楠、刨花润楠、锈毛润楠和滇 楠) (表 １)ꎮ 近交系数都为正值ꎬ最大的为 ０.２７７
润楠的叶绿体基因组比较发现ꎬ叶绿体基因组之 (刨花润楠)ꎬ最小的为 ０.１０３( 薄叶润楠)ꎬ表明这
间序列的相似度为 ９９.９％ꎬ远高于其他属内物种 些物种存在较高的近交(表 １)ꎮ
间叶绿体基因组的相似度ꎮ 所有物种总的遗传分化系数 Ｆ 为 ０.１５６ꎮ 两
ＳＴ
２.１.２ 润楠属系统进化关系 由于物种间很小的 两间的遗传分化系数范围为 ０.００６ ~ ０.３０４ꎬ其中华
序列差异ꎬ因此基于多态叶绿体 ＤＮＡ 片段采用不 润楠和琼桂润楠间最小ꎬ粗壮润楠和短序润楠间
同的构建进化树方法未得到具一致拓扑结构的系 最大ꎮ
统进化树(图 ２)ꎮ ＵＰＧＭＡ 树显示红楠、建润楠和 ２.３ 不同分子标记的物种进化关系及其比较
狭叶润楠聚为一支( 支持率为 ０.９９)ꎬ其他的润楠 通过叶绿体基因组序列分析发现ꎬ在这 １４ 个
属物种聚为一支( 支持率为 ０.９７)ꎬ除这两大支以 润楠属物种组成的突变网络图中ꎬ共有 ５ 个推断
外ꎬ其他支系的自展支持率都很低ꎮ ＭＬ 树表明建 的共同祖先节点ꎬ在图中分别为 Ｖ１、Ｖ２、Ｖ３、Ｖ４ 和
润楠、红楠和狭叶润楠聚为一支ꎬ剩下的物种聚为 Ｖ５ (图 ３:Ａ)ꎮ 锈毛润楠和短序润楠具有共同祖
一支ꎬ但自展支持率都很低ꎮ ＮＪ 树表明狭叶润楠、 先 Ｖ１ꎬ黄绒润楠和基脉润楠具有共同祖先 Ｖ２ꎬ滇
建润楠和红楠单独成为一支( 支持率为 １.００)ꎬ其 润楠、黔桂润楠、薄叶润楠、粗壮润楠和琼桂润楠
他的 ９ 个物种聚为一支且无明显的进化关系ꎬ支 具有共同祖先 Ｖ３ꎬ刨花润楠和华润楠具有共同祖
持率很低ꎮ 这两种基于距离矩阵方法给出的结果 先 Ｖ４ꎬ建润楠、红楠和狭叶润楠具有共同祖先 Ｖ５ꎮ
基本一致ꎬ红楠、建润楠和狭叶润楠聚为一支的支 从 Ｖ１ 到 Ｖ２ 经过 ２ 步突变ꎬ从 Ｖ２ 到 Ｖ３ 以及从 Ｖ３
持率都很高ꎬ其他物种没有明显进化关系( 图 ２)ꎮ 到 Ｖ４ 都只经历了 １ 步突变ꎬ而从 Ｖ２ 到 Ｖ５ 经历了
ＭＰ 树显示建润楠、红楠和狭叶润楠聚为一支ꎬ剩 １４ 步突变ꎮ
下的其他物种聚为一支ꎮ 贝叶斯系统树显示短序润楠单独成为一支ꎬ
２.２ 微卫星位点分析 红楠、建润楠和狭叶润楠聚为一支ꎬ黄绒润楠和基
２.２.１ 中性位点及哑等位基因分析 使用 １０ 个微卫 脉润楠聚为一支ꎬ华润楠和刨花润楠聚为一支ꎬ剩
星位点对 ２３４ 个个体进行 ＰＣＲ 扩增ꎬ均获得清晰的 下的 ５ 个物种聚为一支ꎬ但这些系统树节点的后
产物ꎮ 对异常位点的检测发现ꎬ有 ４ 个位点分布于 验概率值都不高ꎮ 基于樟科系统进化树确定的润
中性区间(９５％置信区间) 之外ꎮ 这些位点的遗传 楠属和楠属分化时间为 ２ ３００ 万年前(２３ Ｍａ)ꎬ润
分化系数 Ｆ 大于中性标记模拟的 Ｆ 值的 ０.９７５ 分 楠属的树冠时间( ｃｒｏｗｎ ａｇｅ) 为 １０.７５ Ｍａ( 图 ３:
ＳＴ ＳＴ
位数ꎬ表明这些位点可能受到歧化选择的作用ꎮ Ｂ)ꎮ 根据这 ２ 个估计时间对润楠属系统树定年ꎬ
１０ 个位点两两之间均不存在显著的连锁关 结果表明 Ｖ１ 形成时间大约为 ７.６９ ＭａꎬＶ２ 形成的
系ꎬ表明它们为独立遗传的ꎮ 对哑等位基因的检 时间大约为 ４.６５ Ｍａꎬ Ｖ３ 形成的时间大约为 ３.９９
测表明ꎬ１０ 个微卫星位点中有 ３ 个可能存在哑等 Ｍａꎬ Ｖ４ 形成的时间大约为 １.６１ ＭａꎬＶ５ 的形成时
位基因ꎬ对这 ３ 个可能存在哑等位基因的位点的 间大约为 ２.５８ Ｍａꎮ
等位 基 因 频 率 校 正 后 用 于 后 续 遗 传 多 样 性 的 使用 １０ 个中性微卫星位点对 １２ 个润楠属物
计算ꎮ 种的个体进行聚类ꎬ 最佳 Ｋ 值为 ２、５ 和 ８ (ΔＫ 值
