Page 50 - 广西植物2024年1期
P. 50
4 6 广 西 植 物 44 卷
兰和景洪白蝶兰的简单重复序列( SSRs)ꎬ单核苷 景洪白蝶兰 3.37 Gb 和 1.07 Gb 的叶绿体基因组
酸的最小重复单位设置为 10 个ꎬ二核苷酸的最小 序列数据ꎮ 在组装之后通过 Bandage v0.8.1 软件
重复单位设置为 5 个ꎬ三核苷酸的最小重复单位 检查ꎬ判断组装后的数据可用ꎮ 基因组结构为闭
设置为 4 个ꎬ四核苷酸、五核苷酸和六核苷酸的最 合的环状四分体结构ꎬ与典型的被子植物叶绿体
小重复单位设置为 3 个ꎮ 基因组结构相同ꎮ
1.5 序列差异分析 2 个白蝶兰属植物的叶绿体基因组的基因图谱
使用 Geneious prime 2022.0.2 的插件 MAFFT 如图 2 所示ꎮ 龙头兰叶绿体基因组全长 154 407
Alignment 将 2 种白蝶兰属的叶绿体基因组序列进 bpꎬ包含一对 26 550 bp 的 IR 区ꎬ由 84 204 bp 的
行比对ꎬ通过 BioEdit 软件编辑整理ꎬ使用 DnaSP LSC 区和 17 103 bp 的 SSC 区分开ꎮ 景洪白蝶兰叶
v6.12.03 计算单核苷酸变异( SNPs) 和插入缺失 绿体基因组全长 153 891 bpꎬ包含一对 26 523 bp 的
(InDels)的数量ꎬ每 100 bp 计算突变率ꎬ并评估 2 IR 区ꎬ由 83 756 bp 的 LSC 区和 17 089 bp 的 SSC 区
种白蝶兰属植物的叶绿体基因组的核苷酸多样性 分开( 表 1)ꎮ 龙头兰叶绿体基因组的 GC 含量为
(Pi)值ꎬ当 Pi 值高于 0. 030 0 时ꎬ将其定义为高 36.5%ꎬ景洪白蝶兰的 GC 含量为 36.6%ꎬ二者都是
变区ꎮ 在 IR 区的 GC 含量最高(43.0%)ꎬ并且高于基因组
1.6 比较基因组分析 的 GC 含量ꎮ 龙头兰和景洪白蝶兰的叶绿体基因组
从 NCBI 上下载已公布的 1 种白蝶兰属和 3 均注释了 111 个唯一基因( unique genes)ꎬ包含 77
种玉凤花属的叶绿体基因组ꎬ包括狭叶白蝶兰、鹅 个蛋白质编码基因、30 个转运 RNA(tRNA) 基因和
毛玉凤 花 ( H. dentateꎬ OK012095)、 线 叶 十 字 兰、 4 个核糖体 RNA(rRNA)基因(表 2)ꎮ 其中ꎬ18 个基
H. cruciformisꎬ与本研究组装的龙头兰和景洪白蝶 因含有 1 个内含子ꎬ2 个基因(pafI 和 clpP 1)含有 2
兰ꎬ利用 mVISTA 程序(Frazer et al.ꎬ2004))对 3 种 个内含子(表 2)ꎮ 在所有基因中ꎬ有 19 个基因位于
白蝶兰属和 3 种玉凤花属的叶绿体基因组序列进 反向 重 复 区 且 具 有 2 个 拷 贝 ( ndhB、 rps7、 rps12、
行比 较ꎮ 通 过 R v4. 1. 3 软 件 运 行 IRscope 脚 本 rps19、rpl2、rpl23、rrn4.5、rrn5、rrn16、rrn23、trnA ̄UGC、
(Amiryousefi et al.ꎬ2018)ꎬ 比 较 6 个 物 种 中 LSC trnH ̄GUG、 trnI ̄CAU、 trnI ̄GAU、 trnL ̄CAA、 ycf2、 trnN ̄
区、SSC 区和 IR 区的边界基因ꎮ GUU、trnR ̄ACG 和 trnV ̄GACꎬ表 2)ꎮ 注释的龙头兰和
1.7 系统发育分析 景洪白蝶兰的叶绿体基因组信息已提交至 GenBank
除了已下载的上述 4 种白蝶兰属和玉凤花属 上ꎬ序列号分别为 OP435916 和 OP435917ꎮ
植物以外ꎬ再从 NCBI 上下载另外 4 种已公布的玉 2.2 叶绿体基因组中的简单重复序列特征
凤花属植物的叶绿体基因组序列ꎬ包括丝瓣玉凤 在龙头兰的叶绿体基因组中共检测到 94 个
花 ( H. fordiiꎬ NC _ 026775 )、 毛 葶 玉 凤 花 ( H. SSRsꎬ包括 4 种类型的 SSRsꎬ单核苷酸、二核苷酸、
ciliolarisꎬMN495594)、H. chejuensis( NC_046821)、 三核苷酸和四核苷酸重复序列分别为 74、16、1、3
落地金钱( H. aitchisoniiꎬMW316693)ꎬ以红门兰亚 个ꎬ其 中 54 个 SSRs 位 于 基 因 间 隔 区 ( intergenic
族的缘毛鸟足兰( Satyrium nepalense var. ciliatumꎬ spacerꎬIGS)、17 个 SSRs 位于蛋白质编码区(coding
MN497244)为外类群( Jin et al.ꎬ2017)ꎬ将龙头兰 sequenceꎬCDS)、23 个 SSRs 位于内含子上(图 3)ꎮ
和景洪白蝶兰的叶绿体基因组序列与这些近缘物 在景洪白蝶兰的叶绿体基因组中共检测到 92 个
种进行系统发育关系分析ꎮ 通过 PhyloSuite v1.2.2 SSRsꎬ包括 5 种类型的 SSRsꎬ单核苷酸、二核苷酸、
平台(Zhang et al.ꎬ2020) 查找最适模型后使用 IQ ̄ 三核苷酸、四核苷酸和五核苷酸重复序列分别为
TREE 构建 Maximum Likelihood( ML) 系统发育树ꎮ 72、15、1、3、1 个ꎬ其中ꎬ58 个 SSRs 位于 IGS 中ꎬ15
通过 FigTree v1.4.0 软件编辑系统发育树ꎮ 个 SSRs 位于 CDS 中ꎬ19 个 SSRs 位于内含子中(图
3)ꎮ 在单核苷酸 SSRs 中ꎬ龙头兰和景洪白蝶兰分
2 结果与分析 别有 3 个和 4 个 C / G 型ꎬ其余均为 A/ T 型ꎮ 大多数
二核苷酸 SSRs 为 AT/ TA 型ꎮ 这一结果与叶绿体基
2.1 龙头兰和景洪白蝶兰叶绿体基因组的基本特征 因组的 SSRs 通常由 A/ T 组成而很少包含 C / G 的观
在去除低质量 reads 后ꎬ分别得到了龙头兰和 点一致(Kuang et al.ꎬ2011ꎻ陈模舜和杨仲毅ꎬ2022)ꎮ