Page 38 - 《广西植物》2023年第1期
P. 38

3 4                                    广  西  植  物                                         43 卷
            方法具有一定的局限性ꎬ尤其是近缘种、易混种在                             基因的注释依据参考序列和开放阅读框确定起始
            没有花和果等关键特征的情况下分类鉴定难度比                              密码子和终止子位置ꎮ
            较大ꎬ因此容易出现错误鉴定ꎮ 随着 DNA 测序技                          1.3 系统发育分析
            术的快速发展ꎬ分子鉴定方法因具有方便、快捷和                                 在 Geneious 软件中提取新测 29 份样品的叶绿

            准确等特点而被广泛应用于物种和药材鉴定( Yu                            体基因片段 psbD ̄trnT 和 trnL ̄trnFꎬ以及核糖体基因
            et al.ꎬ 2021)ꎮ 本研究以澄清傣药“ 傣百解” 基原                   序列 ITS ( 包 括 ITS1、 5. 8S、 ITS2 三 个 区)ꎬ 并 从

            植物为目的ꎬ利用栽培的“傣百解” 基原植物材料、                           GenBank 上下载已有的牛奶菜属 42 种、南山藤属 1
            通光散以及牛奶菜族的其他植物样品ꎬ采用形态                              种以及匙羹藤属 2 种的 psbD ̄trnT、trnL ̄trnT 以及 ITS
            结合分子的方法ꎬ通过整合形态性状比较和系统                              序列ꎬ 以 吊 灯 花 属 ( Ceropegia) 植 物 锥 顶 吊 灯 花
            发育树分析ꎬ探讨以下问题:(1) “ 傣百解” 基原植                        (Ceropegia nilotica)作为外类群(表 2)ꎮ 对 3 个基因

            物是否为通光散ꎻ(2) “傣百解” 的近缘种有哪些ꎻ                         片段的序列分别利用 MAFFT 7.450 软件进行多序
            (3) 如何鉴别区分“傣百解”及其近缘种ꎮ                              列自动比对 (Katoh et al.ꎬ 2013)ꎬ输入到 Geneious
                                                               软件将模糊的区域或空位进行删除或调整ꎮ 使用
            1  材料与方法                                           贝叶斯推断法(Bayesian inferenceꎬ BI)、最大似然法
                                                               (maximum likelihoodꎬ ML) 和最大简约法(maximum
            1.1 材料                                             parsimonyꎬ MP)对 3 个片段分别进行单独和联合重
                 新收集了夹竹桃科牛奶菜属 5 种共 17 份样                       建ꎮ 对 3 个基因 片 段 矩 阵 分 别 利 用 jModeltest 2
            品ꎬ包括栽培于中国医学科学院药用植物所云南                              2.1.6ꎬ选择贝叶斯信息量准则(Bayesian information
            分所南药 园 “ 傣 百 解” 样 品 4 份ꎬ 以 及 南 山 藤 属               criterionꎬ BIC)(Darriba et al.ꎬ 2012) 筛选出序列矩
            (Dregea) 2 种 7 份样品和匙羹藤属( Gymnema) 1                阵最 佳 核 苷 酸 替 代 模 型ꎬ 并 采 用 Mrbayes 3. 2. 7
            种 1 份样品ꎮ 具体样品来源和凭证信息如表 1 所                         (Huelsenbeck & Ronquistꎬ 2001)进行贝叶斯系统发
            示ꎮ 采集植物样品的新鲜叶片经变色硅胶进行快                             生分析ꎮ 采用最大似然法时选用 RAxML ̄HPC2 on
            速干燥ꎬ并低温保存ꎬ以备实验所需ꎮ                                  XSEDE 8.2.12 工具(Alexandrosꎬ 2014) ꎬ选择 GTR+
            1.2 DNA 提取、测序和数据处理                                 CAT 模型ꎬ并设置 1 000 次靴带值估算分支的支持
                 对用硅胶干燥后的样品ꎬ先利用改良的十六                           率(bootstrap support valueꎬ BS)ꎬ使用 CIPRES 网站
            烷基三甲基溴化铵(CTAB)法提取总 DNA (Doyle &                    在线构建 ML 树ꎮ 采用最大简约法时使用 MEGA
            Doyleꎬ 1987)ꎬ再使用琼脂糖凝胶电泳和 Nanodrop                  10.2.6 软件(Kumar et al.ꎬ 2018)ꎬ使用子树修剪和
            分光光度计评估 DNA 质量ꎮ 对检测合格的 DNA                         嫁接(Subtree ̄Pruning ̄Regrafitingꎬ SPR)算法进行分
            样品ꎬ使用 Covaris 超声波破碎仪随机打断ꎬ随后通                       析获得 MP 树(Nei & Kumarꎬ 2000)ꎮ 构建系统进化
            过末端修复、加 A 尾、加测序接头、纯化、PCR 扩                         树时使用 FigTree 1.4.1 (Price et al.ꎬ 2009) 进行美化ꎮ
            增、片段筛选等步骤ꎬ完成片段长度为 350 bp 的测                        1.4 形态特征比较
            序文库制备ꎮ 对构建好的测序文库使用 Qubit 2.0                           将采集到的“傣百解”样品与通光散的模式标
            对浓度进行检测ꎬ定量范围为 2 ~ 1 000 ngꎬ依据溶                     本进行形态观察ꎬ并进一步比较分析在系统树上

            度进行多样品混合ꎬ并利用 Illumina Novaseq 6000                 与“傣百解”亲缘关系最近的物种ꎬ观察和比较“ 傣
            进行双向 150 bp 读长测序ꎮ 对测序原始数据( raw                     百解”基原植物与其姐妹种之间的形态特征差异ꎮ
            data) 直 接 使 用 GetOrganelle 软 件 包 ( Jin et al.ꎬ
            2020) 进行叶绿体基因组和核糖体 DNA 从头( de                      2  结果与分析
            novo)组装ꎮ 选取长春花( Catharanthus roseus) 的
            叶 绿 体 基 因 组 ( KC561139 ) 和 核 糖 体 基 因               2.1 基因序列特征

            (HQ130657) 序 列 作 为 参 考 序 列ꎬ 借 助 Geneious               核糖体 ITS 矩阵的长度为 768 bpꎬ包含 334 个
            (Kearse et al.ꎬ 2012) 对新测物种进行初步注释ꎬ                 变异位点(variable sites)ꎬ其中 202 个为位系统发
            相似 性 参 数 设 置 为 70%ꎬ 并 结 合 ORF ( open               育信 息 位 点 ( parsimony ̄informative site)ꎮ 叶 绿 体
            reading frame) 手工调整注释结果ꎬ其中蛋白编码                     片段 psbD ̄trnT 矩阵的长度为 1 676 bpꎬ 包含 247
   33   34   35   36   37   38   39   40   41   42   43