Page 55 - 《广西植物》2020年第1期
P. 55
1 期 金鑫等: 基于 5 993 个核基因的被子植物系统发育关系研究 5 1
续表 1
物种 数据来源类型 缩写 目 数据来源
Species Data source type Abbreviation Order Data source
锦葵叶刺核藤 Pyrenacantha malvifolia Transcriptome QZZU 茶茱萸目 Icacinales onekp
钩药茶 Oncotheca balansae Transcriptome PVGM 茶茱萸目 Icacinales onekp
美国蜡梅 Calycanthus floridus Transcriptome FALI 樟目 Laurales onekp
美檫树 Sassafras albidum Transcriptome ABSS 樟目 Laurales onekp
阿福花状旱叶草Xerophyllum asphodeloides Transcriptome AFLV 百合目 Liliales onekp
锯齿绿荆棘 Smilax bona-nox Transcriptome MWYQ 百合目 Liliales onekp
刺果番荔枝 Annona muricata Transcriptome YZRI 木兰目 Magnoliales onekp
肉豆蔻 Myristica fragrans Transcriptome OBPL 木兰目 Magnoliales onekp
黄瑞香 Daphne giraldii Transcriptome PUDI 锦葵目 Malvales onekp
圆叶肋果萍蓬草 Nuphar advena Transcriptome WTKZ 睡莲目 Nymphaeales onekp
露兜树 Xerophyta villosa Transcriptome QOXT 露兜树目 Pandanales onekp
马蹄香 Saruma henryi Transcriptome QDVW 胡椒目 Piperales onekp
墨西哥胡椒 Piper auritum Transcriptome MUNP 胡椒目 Piperales onekp
银桦 Grevillea robusta Transcriptome GRRW 山龙眼目 Proteales onekp
框东坚果 Santalum acuminatum Transcriptome RSPO 檀香目 Santalales onekp
昆栏树 Trochodendron aralioides Transcriptome SWOH 昆栏树目 Trochodendrales onekp
蔓茎刺球果 Krameria lanceolata Transcriptome ZHMB 蒺藜目 Zygophyllales onekp
蒺藜 Tribulus eichlerianus Transcriptome KVAY 蒺藜目 Zygophyllales onekp
露水草 Cyanotis arachnoidea RNA-SEQ TCYTR 鸭跖草目 Commelinales https:/ / www.ncbi.nlm.nih.gov/
sra/ SRP144398
无叶莲 Petrosavia sakurai RNA-SEQ TPETR 无叶莲科 Petrosaviaceae This study
1.4 系统发育进化树构建 (Lanfear et al.ꎬ 2009)对串联序列进行分区和进化
我们采用两种方法串联法 ( concatenation) 和 模型检测ꎬ从而设置较合理的分区和为每个分区
溯祖法( coalescence)ꎬ并分别使用 CDS 序列和氨 选择合理的进化模型ꎮ 对 CDS 序列检测下列的四
基酸序列构建进化树ꎮ 无论是 CDS 序列还是蛋白 个分区策略( 表 2):no partitioningꎬ partitioning by
质序 列ꎬ 都 使 用 PRANK v. 170427 软 件 ( http: / / each codon position (three partitions)ꎬ partitioning by
wasabiapp.org / software / prank / ) 进 行 多 序 列 比 对ꎬ gene 和 partitioning by each codon position within
使用 PHYUTILITY v2. 2. 6 软 件 ( Smith & Dunnꎬ each geneꎮ 对蛋白质序列检测下列两个分区策略:
2008)修剪缺失率大于 70% 的位点ꎬ其中 CDS 序 no partitioning 和 partitioning by geneꎮ 参数设置如
列需去除长度小于 300 个碱基的序列ꎬ蛋白质序 下:branch lengths = linkedꎻmodel_selection = aiccꎻ
列需去除长度小于 100 个氨基酸的序列ꎮ search = userꎻmodels = GTRꎬ GTR+Gꎬ GTR+I+G
溯祖法ꎬ先对每个基因使用 RAxML v8.2.11 软 (CDS 序列)或者 models = LG+Gꎬ LG+I+Gꎬ WAG+
件( 默 认 参 数) ( Stamatakisꎬ 2014) 画 树ꎬ 再 使 用 Gꎬ WAG+I+G(蛋白质序列)ꎮ 再使用 iqtree v1.5.5
ASTRAL v5.5.9 软件(Zhang et al.ꎬ 2017)处理所有 软 件 画 树 [ 1000 ultrafast bootstrap replicates ( Von
基因树ꎬ得到共有树ꎬ参数设置“ -t 1--gene-only” Haeseler et al.ꎬ 2013)ꎬ-spp 设置最优分区策略]ꎬ
以获得 bootstrap 值和基因支持率ꎬ枝长使用 iqtree 基因支持率使用 ASTRAL v5.5.9 软件( -t 1) 获得ꎮ
v1.5.5 软件(Nguyen et al.ꎬ 2015)获得ꎮ 最后使用软件 Evolview v2(He et al.ꎬ 2016) 对获得
串联 法ꎬ 先 使 用 PartitionFinder v2. 1. 1 软 件 的所有进化树进行美化ꎮ