Page 81 - 《广西植物》2020年第1期
P. 81
1 期 李依容等: 基于叶绿体基因证据的民族药滇白珠复合群系统发育关系 7 7
续表 1
居群编号 采集号 个体数 海拔
变种名称 采集地 纬度 经度
Population Collection Number of Altitude
Name of variety Sampling location Latitude Longitude
code number samples (m)
YHA 中国贵州印江 LL ̄2017-07 3 27°57′ N 108°36′ E 815
Yinjiangꎬ Guizhouꎬ China
秃果白珠 HXQ 中国台湾玉山 HXQ ̄1001 3 23°27′ N 120°54′ E 1 774
G. leucocarpa Yushanꎬ Taiwanꎬ China
var. psilocarpa
F 菲律宾南达沃棉兰老岛 Penneys 2355 3 6°59′ N 125°16′ E 2 825
Mindanaoꎬ Davao del Surꎬ Philippines
白果白珠原变种 WH 马来西亚鼓亨云顶高原 WH ̄2011 3 3°25′ N 101°47′ E 1 748
G. leucocarpa Genting Highlandsꎬ Pahangꎬ Malaysia
var. leucocarpa
芳香白珠(外类群) FX 中国云南大理苍山 LL ̄2019-04 3 25°48′ N 100°06′ E 2 179
G. fragrantissima Cangshanꎬ Daliꎬ Yunnanꎬ China
(outgroup)
计 ( NanoDrop ® ND ̄1000ꎬ Thermo Fisher Scientificꎬ al.ꎬ 2008)ꎮ 利 用 软 件 jmodeltest 2. 110 ( Darribaꎬ
DelawareꎬUSA)测定 DNA 的浓度和纯度ꎬ用 ddH O 2012)ꎬ分别检测到两个叶绿体基因片段的最佳碱
2
 ̄1 基替代模型都为 GTRꎮ 最大似然法依据 DNA 序列
稀释至 30 ~ 50 ngμL 以备后续 PCR 扩增使用ꎮ
扩增体 系: 10 μL PCR Master Mix ( 2 ×) ( Thermo 突变是随机的这种假设ꎬ将构建系统树所用到的分
 ̄1 枝长度与推断肯定发生的额外突变结合起来ꎬ利用
Scientific)ꎬ9 μL ddH Oꎬ正反向引物(10 ngμL )
2
各 0.4 μLꎬ1 μL 模板 DNAꎮ PCR 扩增反应程序: 特定突变概率的所有假设共同组成一个演化模型
94 ℃ 预 变 性 4 minꎻ 94 ℃ 变 性 45 sꎬ 52 ℃ (model of evolution)ꎬ选取具有最大可能性的树作为
(trnL ̄rpl32) / 60 ℃(rpl33 ̄psaJ)退火 1 minꎬ72 ℃ 延 引证ꎬ 分 析 其 系 统 进 化 关 系ꎬ 并 计 算 支 持 率
伸 45 sꎬ循环 35 次ꎻ72 ℃ 延伸 10 minꎮ 取 2 μL 扩 (BootstrapꎬBP)ꎮ 用贝叶斯法构建系统进化关系树
增产物进行 1.5% 琼脂糖凝胶电泳ꎬ检测扩增产物 时ꎬ需对 4 个平行进化链进行1 000 万代运算ꎬ每
的产量和质量ꎮ 剩余的 PCR 扩增产物纯化后ꎬ采用 10 000代保存一棵进化树ꎬ分别去掉最初和最后的
5.5 μL 的反应体积ꎬ包括 1. 5 μL 的测序混合物 25%系统进化树作为“burn ̄in”ꎮ 每个分支的支持率
(BigDye)、1 μL的测序引物、1 μL 的纯化 PCR 产物 由保存的进化树中出现的后验概率(posterior proba ̄
和 2 μL 的双蒸水ꎮ 反应程序:95 ℃ 预变性 40 sꎻ bilityꎬPP)表示ꎬ后验概率值 PP<0.95 视为低或较弱
95 ℃ 20 s、50 ℃ 10 s、60 ℃ 1 min 30 sꎬ共 35 个循 的支持ꎬ0.95 ~ 1.00 视为高或好的支持ꎮ 两种建树
环进行测序反应ꎮ 反应产物经过 75%乙醇和无水 方法 都 通 过 CIPRES 网 上 建 树 平 台 ( https:/ /
乙醇纯化后ꎬ加入 10 μL 的双蒸水充分溶解ꎬ用 ABI www.phylo.org / )实现ꎬ最后得到的基于 ML 和 BI 法
3730 xl 自 动 测 序 仪 ( Applied Biosystemsꎬ Foster 的严 格 一 致 树 用 Figtree 1. 43 软 件 ( http:/ / tree.
CityꎬCaliforniaꎬUSA)进行测序ꎮ bio.ed.ac.uk / software / figtree / )查看ꎮ 为了直观地判
1.4 数据分析 断滇白珠复合群的系统演化关系ꎬ用 Splitstree4 软
所得序列用 Sequencher 软件进行拼接ꎮ 拼接好 件运用邻接网法(Neighbor ̄Net) 构建系统发育网络
的序列首先逐一对每一个变异位点的原始峰图进行 (Bryantꎬ2003)ꎬBootstrapping 通过 1 000 次重复ꎬ支
校对ꎬ确认无误后用 Se ̄al 软件进行比对构建矩阵ꎬ 持率大于 50%的显示在支上ꎮ
根据头尾的缺失情况对序列做适当调整ꎬ矩阵中的
插入/ 缺失(indel)处理为单次进化事件ꎮ 为解析目 2 结果与分析
前所获各居群间的系统进化关系ꎬ我们对两个叶绿
体片段 rpl33 ̄psaJ 和 trnL ̄rpl32 进行了联合分析ꎮ 本 2.1 从滇白珠复合群 241 个样本和外类群芳香白
研究用最大似然法(Maximum LikelihoodꎬML) 和贝 珠 3 个样本获得 rpl33 ̄psaJ 和 trnL ̄rpl32 片段
叶斯法(Bayesian InferenceꎬBI) 构建系统树ꎬ以比较 rpl33 ̄psaJ 片段比对矩阵长 678 bpꎬ共有 241
不同拓扑结构的异同ꎬ并计算支持率(Stamatakis et 条序列ꎬ一致位点( conserved site) 有 665 个ꎬ变异