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

1 6 9 0                               广  西  植  物                                          42 卷
       如果物种 A、B 存在共享的多倍化事件ꎬ那么
   这次多倍化事件在不同物种中发生的时间应该是
   相同的ꎬKs 峰值也应该是相等的( 图 3:A)ꎮ 黄色
   方块代表两个物种共享的多倍化ꎬ即 Ks                     = Ks ꎬ
                                        AA    BB
   对应的时间范围为物种 A、B 从多倍化事件到当前
   的时间点(绿色的大括号)ꎮ 由于物种不同的进化
   速率ꎬ因此现实情况下的 Ks 和 Ks 并不相等ꎮ 假
                            AA    BB
   设多倍化事件之后物种 A 和 B 有各自的进化速率
   分别为 v 和 v ꎬO 是物种 A、B 的分化节点ꎬ从多倍                      A. 共享多倍化事件ꎻ B. 共享早期分化ꎮ
           A    B
   化事件到分歧点 Oꎬ物种 A、B 的祖先拥有的进化速                          A. Shared polyploidy eventsꎻ B. Shared early divergence.
   率为 vꎮ 那么ꎬ物种 A 的进化速率 v 要想恢复到 v                               图 3  Ks 分布矫正方法的原理
                                   A
                               v                      Fig. 3  Principle of the Ks distribution correction method
   就要乘以它的矫正系数为 λ =               ꎮ 同理ꎬ物种 B 的
                            A
                               v
                                A
                  v                                  心被子植物祖先可能发生了快速辐射分化( Yang
   矫正系数为 λ =        ꎮ 因而ꎬ物种 A、B 间分化的 Ks
                  v                                  et al.ꎬ 2020)ꎮ 因此ꎬ在矫正过程中ꎬ以无油樟目
               B                               AB
                   B
   矫正后为 Ks           =Ks λ λ (Yang et al.ꎬ 2020)ꎮ    为作为参考ꎬ不讨论它和睡莲目的关系ꎬ认为五大
              AB-correction  AB  A  B
       如果两个物种 A、B 虽不存在共享的多倍化事                        分支的分化时间尺度在同一个时间范围内ꎮ 基于
   件但存在共享的早期分化事件ꎬ就通过寻找外类                             核心真 双 子 叶 植 物 共 享 的 γ 事 件ꎬ 时 间 范 围 为
   群来辅助矫正(图 3:B)ꎮ 物种 C、D、E 是外类群ꎬ                     115 ~ 130 百万年( Million years agoꎬMya)ꎬ对 44 个
   物种 C 和 D 的祖先在 P 点与物种 A、B 的祖先分                     被子植物基因组( 表 1) 进行了时间尺度矫正( 图
   化ꎬ所以物种 C 与 A、B 间的 Ks 峰值应该相等ꎬ物                     4)ꎮ 从矫正后的时间尺度来看ꎬ被子植物在 130
   种 D 与 A、B 间的 Ks 峰值也应该相等ꎬ即 Ks                 =    百万年前附近ꎬ单子叶植物、真双子叶植物、木兰
                                             CA
   Ks ꎬKs   = Ks ꎮ 同样ꎬ由于物种间不同的进化速                    类植物祖先都发生了快速辐射进化ꎬ与 Zhang 等
     CB   DA    DB
   率ꎬ因此现实情况下它们大多不相等ꎮ 按照前面                            (2020b) 的结论一致ꎮ 此外ꎬ在早白垩世(130 百
   的假设ꎬ                                              万年)时期ꎬ白垩纪-古新世( K-Pg) 边界时期(66
       Ks        Ks λ λ    Ks λ       λ   Ks         百万年) 和中新世(20 百万年ꎬ靠近冰川期) 很多
         CA-correction  CA  C  A  CA  A  A  CB
                =        =       = 1ꎬ即   =    ꎮ
       Ks        Ks λ λ    Ks λ        λ  Ks         被子植物发生的多倍化事件ꎬ研究发现 WGD 的时
         CB-correction  CB  C  B  CB  B  B  CA
   同理ꎬ                                               间在被子植物的系统发育中并不是随机分布与
       Ks        Ks λ λ    Ks λ        λ   Ks
         DA-correction  DA  D  A  DA  A  A   DB      Wu 等(2020)的结论一致ꎮ
                =         =      = 1ꎬ即   =    ꎮ
       Ks        Ks λ λ    Ks λ        λ   Ks            尽管不同物种的进化速率数值显著不同ꎬ但
         DB-correction  DB  D  B  DB  B  B   DA
       当选取的外类群越多ꎬ获取的 λ 和 λ 的关
                                     A     B         是同一类群中的进化速率往往具有部分一致性ꎮ
                                             λ       由矫正方法可知ꎬ矫正后的 Ks 峰值应该相等ꎮ 因
   系越准确ꎮ 取平均值表示它们之间的关系                        A  =
                                             λ
                                              B      此ꎬKs 峰值越大ꎬ表明进化速率越快ꎮ 对木兰类植
         Ks   Ks                                     物、真双子叶植物和单子叶植物与无油樟的 Ks 峰
   mean(   CB ꎬ  DB ꎬ...)ꎮ
         Ks   Ks                                     值的比较发现ꎬ木兰类植物( 大多数为木本) 进化
           CA   DA
   2.3 被子植物系统发育树时间矫正                                 速率最慢ꎬ真双子叶植物(大多数为灌木) 次之ꎬ单
       目前ꎬ很多用系统发育树的方法推测被子植                           子叶植物(大多数为草本) 进化速率最快( 表 3)ꎬ
   物的演化时间ꎬ认为被子植物的起源为三叠纪 225                          这与多年生木本植物比草本植物的分子进化速率
   百万年至 240 百万年前( Magallónꎬ 2010)ꎬ这与起                慢的结论相符( Lanfear et al.ꎬ 2013)ꎮ 此外ꎬ对多
   传粉作用的核心植食性鳞翅目昆虫的起源时间                              倍化事件发生的时间与矫正前后的 Ks 峰值比较
   (约 230 百万年前)一致(Li et al.ꎬ 2019)ꎮ 由于无              (图 5)发现ꎬ矫正前的 Ks 峰值与时间并不是线性
   油樟目和睡莲目、核心被子植物五大分支之间的                             关系ꎮ 随着 Ks 峰值的增大ꎬ多倍化事件发生的时
   关系仍然没有完全解析ꎬ且已有多个证据暗示核                             间并没有更古老ꎻ 由于矫正后的 Ks 峰值与时间成
   67   68   69   70   71   72   73   74   75   76   77