Page 129 - 《广西植物》2022年第10期
P. 129

10 期                   陈丽琼等: 质体系统发育基因组解析旋花科系统发育关系                                         1 7 4 7

           表 3  AU 检验和 SH 检验的结果                      旋花科的界定混乱ꎬ特别是菟丝子属和马蹄金族
         Table 3  Results of AU test and SH test     (Austinꎬ 1973)ꎮ 多数学者主张将寄生习性的菟
                                                     丝子 属 处 理 为 旋 花 科 内 单 属 族 ( Choisyꎬ 1834ꎻ
    拓扑结构约束         最大似然             近似无   SH 检验
    Topology        估计值       deltaL  偏检验            Bentham & Hookerꎬ 1873ꎻ Peterꎬ 1897ꎻ Austinꎬ
                                           SH test
    constraint       -LnL           AU test
                                                     1998)或亚科(Peterꎬ 1891)ꎬ少数学者主张将菟丝
    WCG 拓扑结构       721 494.39  —      —     —        子属处理为单属科( Lindleyꎬ 1853ꎻRobertyꎬ 1952ꎬ
    WCG topology
                                                     1964ꎻAustinꎬ 1973)ꎮ 另外ꎬ还有学者将马蹄金族
    LSC 拓扑结构       721 516.48  37.24  0.08  0.45
    LSC topology                                     作为 独 立 的 科 ( Dumortierꎬ 1829ꎻ Austinꎬ 1973ꎻ
    SSC 拓扑结构       723 779.52  339.44  0.00  0.00    Stefanovi et al.ꎬ 2003)ꎮ 本研究的质体基因组系统
    SSC topology
                                                     发育分析结果表明ꎬ现在划定的旋花科是一个“ 自
    IR 拓扑结构        726 774.28  512.99  0.00  0.00
    IR topology                                      然”类群ꎬ即菟丝子属与马蹄金族均应包括在旋花
                                                     科内ꎬ不应独立成科ꎬ这与 Stefanovi 等(2002) 基于
                                                     4 个质体基因组片段系统发育结果一致ꎮ 同时ꎬ该
   WCG ̄/ CDS ̄PP = 1.00 / 1.00)ꎬ剩下的 5 个族分为 2
                                                     结果支持 Stefanovi 等(2003) 对旋花科的科下等
   个分支在所有的分析中都得到了强支持ꎬ即旋花
                                                     级划分ꎬ即菟丝子科并入旋花科内为菟丝子族ꎬ与
   族、鱼黄草族和番薯族组成 Clade Iꎬ盐帚花族和马
                                                     马蹄金族一起被置于旋花亚科内ꎮ
   蹄金族组成 Clade Ⅱꎮ 在 Clade I 中ꎬ旋花族和番
                                                     3.2 旋花亚科的系统发育关系
   薯族各自形成单系群ꎬ嵌套在鱼黄草族中ꎬ且 LSC
                                                         虽然现有的分子系统学研究已构建了旋花科
   矩阵分析时和菟丝子族聚在一起( BS = 60)ꎮ 但
                                                     2 亚科 12 族的系统框架ꎬ但旋花亚科内主要分支
   是ꎬIR 矩阵(BS = 93)和 SSC 矩阵( BS = 61) 支持菟
                                                     的系统发育关系尚未厘清ꎬ特别是心被藤族、丁公
   丝子族与盐帚花族+马蹄金族为姊妹类群ꎮ 值得
                                                     藤 族 和 菟 丝 子 族 的 系 统 位 置 ( 李 攀ꎬ 2020 )ꎮ
   关注的是ꎬ鱼黄草属 4 个种分散在 3 个不同的分支
                                                     Stefanovi 等(2002)基于 4 个质体基因数据推断丁
   上ꎬ即海南山猪菜( Merremia hainanensis) 与旋花族
                                                     公藤族和心被藤族各为一个分支ꎬ依次位于其他
   聚在 一 起 ( BS ≥ 96ꎬ PP = 1. 00 )ꎬ 五 叶 鱼 黄 藤
                                                     旋花亚科的基部ꎮ 随后增加线粒体基因 atpA、核
   (Merremia quinquefolia) 与番薯族聚在一起( BS≥
                                                     基因 RPB2 和质体基因 rpl2 后ꎬ心被藤族、丁公藤
   98ꎬPP = 1.00)ꎬ剩下 2 种和其他 3 种鱼黄草族的
                                                     族和旋花亚科其他族形成“梳齿”关系(Stefanovi &
   属种一起(BS = 100ꎬPP = 1.00)ꎮ 番薯属是一个并
                                                     Olmsteadꎬ 2004)ꎬ 无 法 确 定 系 统 位 置ꎮ Refulio ̄
   系类群ꎬ所有分析都形成了 2 个分支( BS = 100ꎬ
                                                     Rodriguez 和 Olmstead(2014)联合 9 个质体基因和
   PP = 1.00)ꎬ其中 1 支包括了另 4 属 5 种植物ꎮ                  1 个线粒体基因序列的分析结果较弱地支持心被
   2.4 系统发生冲突检测
                                                     藤族和丁公藤族聚为一支ꎬ位于旋花亚科基部ꎮ
       由于 WCG 矩阵和 CDS 矩阵与其他 3 个矩阵
                                                     由于该研究旋花科取样较少ꎬ这 2 个族各仅有 1
   构建的系统树存在明显的拓扑结构差异ꎬ因此利
                                                     种ꎬ支持率较低ꎬ且缺少菟丝子族、盐帚花族等 6
   用 AU 检验和 SH 检验进行冲突显著性评估( 表                        族的代表ꎬ因此二者的亲缘关系和系统位置没有
   3)ꎮ 检测结果表明ꎬWCG 矩阵与 LSC 矩阵的拓扑                      完全澄清ꎮ 本研究中的 4 组数据矩阵( 除 IR 矩阵
   结构不存在明显的拓扑结构冲突ꎬ即菟丝子族与                             外)分析结果均强支持心被藤族和丁公藤族为姊
   Clade I 构 成 姊 妹 关 系 存 在 可 能 ( P     = 0.083 7ꎬ
                                     AU              妹关系ꎬ且可能是最早从旋花亚科内分化出来ꎮ
   P = 0.451)ꎬ但拒绝 IR 区和 SSC 区的拓扑结构关                  由于 IR 区和 SSC 区矩阵的 ML 树在这些分支的支
    SH
   系(P 值均小于 0.001)(表 3)ꎮ                             持率较弱ꎬ因此 AU 检验和 SH 检验均显示两种拓
                                                     扑结构与 WCG 矩阵的拓扑结构存在明显冲突ꎬ这
   3  讨论与结论                                          可能是 IR 区和 SSC 区基因的核苷酸替换速率不

                                                     一致所造成的(Zhu et al.ꎬ 2016)ꎮ
   3.1 广义旋花科的界定                                          针对菟丝子族在旋花科内的系统位置ꎬ本研
       传统分类上ꎬ研究者们对性状加权不同ꎬ使得                          究结果显示ꎬ首先 WCG 矩阵和 CDS 矩阵强支持
   124   125   126   127   128   129   130   131   132   133   134