Page 58 - 《广西植物》2020年第1期
P. 58

5 4                                   广  西  植  物                                         40 卷
   树ꎬ并比较它们之间的不同( 图 4)ꎬ以评估树的稳                         类植物的姊妹群”ꎬ与 APG IV 一致(图 4)ꎮ
   定性ꎮ CDS 序列和蛋白质序列ꎬ都分别使用五个                          2.3 分化时间估计
   数据集ꎬ总共构建 20 棵树(5 棵 CDS 串联法树ꎬ5                         基于 742 个基因 CDS 序列串联方法建立的进
   棵 CDS 溯祖法树ꎬ5 棵 AA 串联法树和 5 棵 AA 溯                  化树ꎬ我们估计了被子植物的分化时间( 图 6)ꎮ
   祖法 树 )ꎮ 这 5 个 数 据 集 分 别 包 含 5 928 个               我们认为被子植物的起源时间为 237.78 百万年前
   orthologs( ≥50 samples)、3 384 个 orthologs( ≥70    (95%置信区间为 202.6 ~ 278.08)ꎬ与主流观点认
   samples)、1 791 个 orthologs(≥80 samples)、742 个     为的 225 百万年至 240 百万年前一致 ( Magallonꎬ
   orthologs( ≥ 85 samples) 及 42 个 orthologs ( ≥ 89  2010ꎻ Smith et al.ꎬ 2010ꎻ Zeng et al.ꎬ 2014)ꎮ 木
   samples)ꎮ                                         兰类植物与单子叶植物和真双子叶植物的分化时
       这 20 棵进化树主要是为了进一步确定图 1 中                      间约为 166.11 百万年前ꎻ五桠果科与蔷薇类和菊
   五类被子植物间演化关系和真双子叶植物内部各                             类植物的分化时间约为 124.23 百万年前ꎻ蔷薇类
   目间系统发育关系ꎮ 这些进化树中的大多数ꎬ是                            植物与菊类植物的分化时间约为 116.98 百万年

   与使用 742 个基因 CDS 序列( 共 4 069 848 位点)               前ꎻ唇形类植物( Lmiids) 与桔梗类植物( Campanu ̄
   串联方法建立的进化树高度一致的( 图 5) ( 使用                        lids)的分化时间约为 102.37 百万年前ꎮ

   3 384个基因 AA 序列建立的进化树ꎬ和使用 1 791
   个基因 AA 序列建立的进化树ꎬ也是相同的最佳拓                          3  讨论与结论

   扑结构)ꎮ
   2.2.1 木兰类植物、单子叶植物及双子叶植物间演化                            长期以来ꎬ被子植物的系统发育关系重建ꎬ都
   关系  无论核酸序列还是蛋白质序列ꎬ使用串联法                           是使用质体基因、线粒体基因或少数保守的单拷
   和溯祖法建立的进化树基本都支持拓扑结构[(真                            贝核基因ꎮ Yang & Smith(2014) 报道了一种基于

   双子叶植物ꎬ单子叶植物)ꎬ木兰类植物](图 4)ꎮ                         系统进化树的同源基因聚类及去旁系同源基因的
   2.2.2 金栗兰科与金鱼藻科  我们的研究表示ꎬ金                        方法ꎬ我们使用此种方法对收集的 88 种植物核基
   鱼藻科是真双子叶植物的姊妹群ꎬ这与前人的研                             因集进行聚类ꎬ共获得了多达 5 993 个 one-to-one
   究结果一致(图 4)ꎮ 但金栗兰科是所有被子植物                          基因家族ꎬ并从这个数据集里面截取各种大小的
   (除 ANITA 外) 的基底旁系群ꎬ这与 APG IV 认为                   数据进行进化树重建ꎬ以测定进化树的稳定性ꎮ
   的“金栗兰科是木兰类植物的姊妹群”是不同的ꎮ                                获得比以前更多的核基因家族后ꎬ制约系统
   2.2.3 双子叶植物内部各目的系统发育关系  我                         演化关系构建的另一个因素就是大量的计算资源
   们的研究认为ꎬ五桠果科是蔷薇类植物和菊类植                             和计算时间ꎮ 构建系统进化树时ꎬ一般需要设置
   物共同的姊妹群ꎬ虎耳草目是蔷薇类植物的姊妹                             bootstrap 值(100 ~ 1 000) 迭代ꎬ此步骤非常耗费计
   群ꎬ这都与 APG IV 一致(图 4)ꎮ                             算时间ꎮ Nguyen et al.(2015)发表的软件 iqtreeꎬ采
       APG IV 认为“檀香目和石竹目是菊类植物的                       用 ultrafast bootstrap approximation( UFBoot) 方法获
   姊妹群”ꎬ而我们的研究否定了这一结论:20 棵进                          得 bootstrap 值(Von Haeseler et al.ꎬ 2013)ꎬ比 RAx ̄
   化树中ꎬ所有结果都支持“ 石竹目是蔷薇类植物的                           ML 软件的传统方法ꎬ计算速度快 10 ~ 40 倍ꎬ并且
   姊妹群”ꎻ大部分支持“ 檀香目是蔷薇类植物的姊                           获得的 bootstrap 值更精确ꎮ
   妹群”ꎬ这与 Zeng et al.( 2017) 的研究结果一致ꎻ                    我们使用多达 5 993 个 one-to-one 基因家族
   少部分支持“ 檀香目是蔷薇类植物和菊类植物共                            构建的进化树ꎬ与 APG IV 报道的主要差异为檀香
   同的姊妹群”(图 4)ꎮ                                      目和石竹目在系统发育树中的位置ꎬ本研究认为
       APG IV 认为“ 智利藤目 是 菊 类 植 物 的 姊 妹               “檀香目 和 石 竹 目 是 蔷 薇 类 植 物 的 姊 妹 群”ꎬ 而
   群”ꎬ而我们的研究只有少部分支持这一结论ꎮ 使                           APG IV 认为“檀香目和石竹目是菊类植物的姊妹
   用蛋白质序列建立的进化树ꎬ无论串联还是溯祖                             群”ꎮ 可能原因有以下两个:一是基因数目的增
   法ꎬ都支持“智利藤目是蔷薇类植物和菊类植物共                            多ꎻ二是本研究所选 88 个植物只有一半使用的基
   同的姊妹群”ꎮ 使用核酸序列建立的进化树ꎬ随着                           因组序列ꎬ另一半为转录组序列ꎬ而转录组序列一
   基因数目的增多ꎬ逐渐转变为支持“ 智利藤目是菊                           般存在大量的基因缺失(即未表达基因较多)ꎮ
   53   54   55   56   57   58   59   60   61   62   63